[ncl-talk] Fw: Can't got value of variables other than Lat and Lon

Mary Haley haley at ucar.edu
Wed Aug 30 18:12:03 MDT 2017


Xiaoli,

Thanks for providing the file offline.

I was able to read this data using NCL V6.4.0. There is indeed a problem if
you read it with NCL V6.3.0.

I've attached the simple script I used, to try two reading two different
precip variables, using two different methods for reading the group data.
Both methods resulted in the same min/max being printed:

min=0   max=265.672
min=0   max=294.619
min=0   max=265.672
min=0   max=294.619
difference between two methods for HS group : min=0   max=0
difference between two methods for NS group : min=0   max=0

If you are running NCL V6.3.0, then I suggest that you upgrade to 6.4.0.
The easy way to upgrade may be with conda:

http://www.ncl.ucar.edu/Download/conda.shtml

--Mary




On Tue, Aug 29, 2017 at 2:25 PM, Liu, Xiaoli <xiaoli.liu at ou.edu> wrote:

>
>
>
> ------------------------------
> *From:* Liu, Xiaoli <xiaoli.liu at ou.edu>
> *Sent:* Wednesday, August 30, 2017 6:18 AM
> *To:* ncl-talk at ucar.edu
> *Subject:* Can't got value of variables other than Lat and Lon
>
>
> Hello,
>
> I want use the data of GPM 2A.GPM.DPR.V7 with NCL. I found that I can get
> the value of Lat and lon correctlly. But as to the other variables, the
> value got is very unusal.
>
> Here is my script:
>
>  dir="/scratch/liu4119/gpm/data/"
>   FILES=systemfunc(" ls -1 "+dir+"*.HDF5")
>   numFILES = dimsizes(FILES)
> ;  print("numFILES = " + numFILES)
> ;  print(FILES)
> ;  print (" ")
>   do i=32,32
>     h5_file = addfile(FILES(i),"r")
> ;  g = h5_file=>/NS
> ;   print(g)
> ;  g1=g=>/SLV
> ;   print(g1)
> ;  list_filevars(g)
>   data_raw = h5_file->/NS/SLV/precipRateESurface
>   print(data_raw(5396,:))
>   data = data_raw
>   data at long_name = "precipRate"
>   data at _FillValue = -9999.9
>   data at units = data_raw at Units
>   longitude=h5_file->/NS/Longitude
>   printVarSummary(longitude)
>   latitude=h5_file->/NS/Latitude
>   print(latitude(5396,:))
>   print(longitude(5396,:))
>   printVarSummary(latitude)
>
>
>
>
> Thanks!
>
> Xiaoli
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170830/92757b79/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: data_raw.png
Type: image/png
Size: 25095 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170830/92757b79/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: longitude.png
Type: image/png
Size: 32193 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170830/92757b79/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: latitude.png
Type: image/png
Size: 33981 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170830/92757b79/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gpm_simple_read.ncl
Type: application/octet-stream
Size: 1033 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170830/92757b79/attachment.obj>


More information about the ncl-talk mailing list