¿Cómo obtener datos de aerosoles CALIPSO de forma cómoda y eficaz? Además de descargar los datos originales del sitio web oficial.
El código es el siguiente:
cargar "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
cargar
" $NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
cargar
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
cargar
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
cargar
"$NCARG_ROOT/lib /ncarg/ nclscripts/csm/skewt_func.ncl"
cargar
"$NCARG_ROOT/lib/ncarg/nclscripts/csm/wind_rose.ncl"
comenzar
dir=("./DATA/")
fi=addfile(dir "./CAL_LID_L2_VFM-ValStage1-V3-01.2007-02-06T16-29-28ZN.hdf" ,"r ")
lat
= fi-gt;Latitud(:,:)
lon = fi-gt;Longitud(:,:)
tiempo =
fi-gt;Profile_UTC_Time(:,:)
aerosol =
fi-gt;Feature_Classification_Flags(:, :)
aerosol@units
="categroy"
wks = gsn_open_wks("png", "gsn_contour"); abre un ps
archivo
gsn_define_colormap(wks, "rainbow")
res =
True
res@gsnMaximize = True
res@gsnLeftString =
"Máscara de característica vertical"
res@tiYAxisString = "Altitud
(km)"
res @tiYAxisFontHeightF = 0.02
res@cnLinesOn =
False
res@cnFillOn = True
res@cnFillMode =
"RasterFill" ; ¡Mucho más rápido!
res@gsnAddCyclic =
False
res
@cnLevelSelectionMode = "Niveles Explícitos"
res@cnLevels
= (/0, 1, 2, 3, 4, 5, 6/)
res@ gsnSpreadColors =
True
res@gsnSpreadColorStart = 1
res@gsnSpreadColorEnd = -3; no
incluye gris o blanco; al final
;---Marcar recursos
latstep
= 10
res@tmXBMode = "Explicit"
res@tmXBValues
= lat(::latstep)
res@tmXBLabels =
sprintf("5.2f", lat(: :latstep )) "~C~" \
sprintf("5.2f",lon(::latstep))
;---Tickmark
recursos
res@lbOrientation = "Vertical"
res@lbLabelAlignment =
"BoxCenters" ； El valor predeterminado es "InteriorEdges"
res @lbLabelAutoStride =
True
res@lbLabelFontHeightF = 0.015
res@lbLabelStrings = (/"No
determinado", "Limpiar Marino", "Polvo", "Continental contaminado", \
"Continental limpio", "Polvo contaminado", "Humo", "Otros"/)
plot =
gsn_csm_contour(wks, aerosol(:,:),res Dibujar un trazado de contorno.
end