Hi,
I haven't used Calipso, but I guess it is pretty much the same as Cloudsat:
There are two modes to open these files hdf and hdfeos. Just modify the extension when using addfile, e.g.
eos = addfile(satdir+filename+".hdfeos","r")
hdf = addfile(satdir+filename+".hdf","r")
note that both are in fact the same file. This is described is several posts to ncl-talk.
Some variable can only be read from hdf
dbz_short = hdf->Radar_Reflectivity
dbz = short2flt(dbz_short)
others from eos
utc = eos->UTC_start_2B_GEOPROF
tai = eos->TAI_start_2B_GEOPROF
time = eos->Profile_time_2B_GEOPROF
lat = eos->Latitude_2B_GEOPROF
lon = eos->Longitude_2B_GEOPROF
you may use
eNames = getfilevarnames(eos) ; get names of variables on file
vNames = getfilevarnames(hdf) ; get names of variables on file
print(eNames)
print(vNames)
or ncl_filedump.
To plot the fields, I recommend "TriangularMesh" and "Rasterfill"
res = True
res@cnFillOn = True ; color fill
res@cnFillMode = "RasterFill" ; raster mode
res@cnLinesOn = False ; Turn off contour lines
res@trGridType = "TriangularMesh"
Use the time and height coordinates to get the "quicklooks"
res@sfXArray = time2d
res@sfYArray = height2d
and you are done
plot = gsn_csm_contour(wks,y, res) ; create plot
Best, Axel
________________________________________
Von: ncl-talk-bounces@ucar.edu [ncl-talk-bounces@ucar.edu]" im Auftrag von "Glauber Mariano [glaubermariano@gmail.com]
Gesendet: Mittwoch, 31. August 2011 20:23
Bis: ncl-talk@ucar.edu
Betreff: [ncl-talk] Calipso Data - how can I get the height?
Hi NCL users,
I got the calipso data and I'm trying to plot of Total_Attenuated_Backscatter_532 in a xy-plot. I intend to do just an average of the 50 lat/lon.
How can I get the height of the data? Because I don't have a var named "height" as the "latitude" and "longitude".
My script is below:
--------------------------------------------------------------------------------------------------------------
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
calipsofile = addfile("CAL_LID_L1-ValStage1-V3-01.2011-05-01T05-00-22ZN.hdf","r")
tab = calipsofile->$"Total_Attenuated_Backscatter_532"$
printVarSummary(tab)
avgtab = dim_avg(tab(fakeDim41|1:582, fakeDim40|1:50))
wks = gsn_open_wks( "pdf", "graficos_calipso_mean" )
res = True
res@trYReverse = True
res@trXMinF = 0.
grafico1 = gsn_xy( wks, avgtab(:), ispan(1,582,1), res )
--------------------------------------------------------------------------------------------------------------
Someone could help me?
Best regards,
Glauber Lopes Mariano
Email: <mailto:glauber.mariano@cptec.inpe.br> glaubermariano@gmail.com<mailto:glaubermariano@gmail.com>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Aug 31 12:41:30 2011
This archive was generated by hypermail 2.1.8 : Wed Sep 07 2011 - 10:58:58 MDT