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

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Wed Aug 31 2011 - 13:43:48 MDT

Maybe this is a question better asked of the CALIPSO provider.

Unfortunately, CALIPSO created hdf files are lacking in information.
Really, the whole idea of HDF and NETCDF is that they are self
describing. No external information should be necessary.
After all the $$$ spent on CALIPSO somebody should have been
looking at the output to see if it was informative to
non-experts. :-(

You can send a sample file to:
ftp ftp.cgd.ucar.edu
anonymous
Use your email address for the password
cd incoming
put CAL_LID_L1-ValStage1-V3-01.2011-05-01T05-00-22ZN.hdf
quit

Send email to ncl-talk (you must subscribe first) with the following
information:

     * The version of NCL you are running (ncl -V).
     * The type of system you are on (uname -a).

On 08/31/2011 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

-- 
======================================================
Dennis J. Shea                  tel: 303-497-1361    |
P.O. Box 3000                   fax: 303-497-1333    |
Climate Analysis Section                             |
Climate & Global Dynamics Div.                       |
National Center for Atmospheric Research             |
Boulder, CO  80307                                   |
USA                        email: shea 'at' ucar.edu |
======================================================
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Aug 31 13:43:54 2011

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