[ncl-talk] Incorrect reading of dimensions from he5

Rebecca Buchholz buchholz at ucar.edu
Thu Nov 20 11:28:14 MST 2014


Dear NCL team,

I'm having difficulty reading some MOPITT V6 profile data. I am finding
zeros where there should not be any. The issue seems to be coming from NCL
incorrectly reading the dimensions of the hdf 5 file (he5). It happens for
every file but here is an example:

The file:
ftp://l5eil01.larc.nasa.gov/MOPITT/MOP02T.006/2010.09.07/MOP02T-20100907-L2V16.2.1.he5

I have attached an example ncl file script (test3.ncl) that shows how I
read in the file.

Differences:
---------------------------------------------------------------------------
*Using hdfview these are the correct values for the profile dimension:*
OBJECT = Dimension_4
                     DimensionName="nPrs"
                     Size =9
END_OBJECT=Dimension_4


*Dimensions from an ncl_filedump gives:*
nPrs_MOP02 = 10
-----------------------------------------------------------------------------

All the other dimensions are identical in size, except for nTime, which in
hdfview has a size of 2 and ncl_filedump gives a size of 208903 (the number
of time values).

Also, the named dimensions of some variables are incorrect. Specifically
for "RetrievalAveragingKernelMatrix" and "RetrievalErrorCovarianceMatrix"
the named dimensions are not the same:
*hdfview:* (ntime, nPrs2, nPrs2)
*ncl_filedump:* (ntime_MOP02, nPrs2_MOP02, nPrs_MOP02)
The averaging kernel and co-variance matrices should be 10x10 i.e. nPrs2,
nPrs2.

Does anyone have any idea why NCL is incorrectly reading in the dimensions
of nPrs and how I can fix it?

Specs: I'm running NCL version 6.1.2
uname-a gives
Linux mopfl 2.6.32-431.29.2.el6.x86_64 #1 SMP Tue Sep 9 13:45:55 CDT 2014
x86_64 x86_64 x86_64 GNU/Linux

Kind Regards,
Rebecca
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141120/a995483f/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test3.ncl
Type: text/x-ncl
Size: 2519 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141120/a995483f/attachment.bin 


More information about the ncl-talk mailing list