Re: Calipso Data - how can I get the height?

From: Seifert Axel <Axel.Seifert_at_nyahnyahspammersnyahnyah>
Date: Wed Aug 31 2011 - 12:41:21 MDT


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

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: []&quot; im Auftrag von &quot;Glauber Mariano []
Gesendet: Mittwoch, 31. August 2011 20:23
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"$

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: <><>
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
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