[ncl-talk] TRMM 2A12 data

Dennis Shea shea at ucar.edu
Tue May 31 15:17:15 MDT 2016


I don't think people associated with ncl-talk have time to read a large
document. We have our own jobs to do.
---

Look at your data. You must look at the file's contents.

Presumably you have something like:

%> ncl_filedump 2A12.20100101.69095.7.HDF -v heightLayerTop | less

Many part snipped out.

  dimensions:
      nlayer = 28
      ncluster = 100
      nfindex = 13
      nspecies = 6
      nscan = 2921  // unlimited
      npixel = 208

    float heightLayerTop ( nlayer )
         units :        km
         hdf_name :     heightLayerTop

     float Latitude ( nscan, npixel )
         units :        degrees
         hdf_name :     Latitude

      float Longitude ( nscan, npixel )
         units :        degrees
         hdf_name :     Longitude

       float surfaceRain ( nscan, npixel )
         units :        mm/hr
         hdf_name :     surfaceRain

      float cloudWaterPath ( nscan, npixel )    <==== *no layers*
         units :        kg/m^2
         hdf_name :     cloudWaterPath

     float cluster ( ncluster, nlayer, nfindex, nspecies )
         hdf_name :     cluster

     byte clusterNumber ( nscan, npixel, nspecies )
         hdf_name :     clusterNumber

      float clusterScale ( nscan, npixel, nspecies )
         hdf_name :     clusterScale


Variables:

Variable No. 0
     heightLayerTop:

  0.5000       1.000       1.500       2.000       2.500       3.000
   3.500       4.000       4.500       5.000       5.500       6.000
   6.500       7.000       7.500       8.000       8.500       9.000
   9.500       10.00       11.00       12.00       13.00       14.00
   15.00       16.00       17.00       18.00

=====================================
So ...
 layer index 0 ==> 0.5000
 layer index 1 ==> 1.000
etc

====================================

   LAYER = 1                 ; you specify

   fa12 = addfile("....hdf","r")
   x_1km     = fa12->cluster(?, LAYER, ?,?)
   printVarSummary(x_1km)

 ======================================

The Latitude and Longitude units must be changes to degrees_north and
degrees_east





On Sat, May 28, 2016 at 12:11 PM, Geeta Geeta <geetag54 at yahoo.com> wrote:

> Hi .
> i have to make spatial plots of each hydrometeor available in 2A12 TRMM
> data.
> The document says that  (Pls refer to Pg 119 of the document attached).
>
> Profile Value = clusterScale * cluster(S,F,L,C)
> These arrays are of different sizes. How can I get the spatial plot of say
> CLOUDWater at a height of 1.0km above msl.
> I am able to plot surface_rain variable but not the variable which is
> extending in the vertical
>
> Pls suggest.
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160531/1920877f/attachment.html 


More information about the ncl-talk mailing list