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

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Fri Sep 02 2011 - 08:30:03 MDT

As noted in previous postings to ncl-talk,
NCL can't read vgroup data on HDF files, so we got these
values via hdfview and wrote them to an ascii file. (attached)

The variable has a "fillvalue" attribute but this is not recognized
by many software tools including NCL which understands "_FillValue".
This was added to the imported variable.

Good Luck

On 8/31/11 12:23 PM, Glauber Mariano wrote:
> 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

Received on Fri Sep 2 08:30:09 2011

This archive was generated by hypermail 2.1.8 : Wed Sep 07 2011 - 10:58:58 MDT