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

Mary Haley haley at ucar.edu
Tue Aug 29 16:02:50 MDT 2017


Xiaoli,

Just because you have a bunch of values equal to 0.0, it doesn't mean your
data is incorrect.
​ You should ​
also look at the min/max of your
​whole ​
data
​variable ​
with:

printMinMax(data_raw,0)

Also, to be on the safe side, rather than hard-coding data at _FillValue to
-9999.9, you should use the attribute that is attached to data:

data at _FillValue = data at CodeMissingValue
​data_raw = data​

Please do this before the "printMinMax" call.

​If you still think something is wron​g with your data after looking at the
printMinMax output, then respond back to this message and cc
ncl-talk at ucar.edu.  It helps if you can provide more information about why
you think the data is incorrect.

Also, if you continue to have problems, then it also helps if you can
provide your data file. You can put it on our anonymous ftp:

http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP

​Thanks,
​

--Mary

On Tue, Aug 29, 2017 at 3:38 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. When I use orbit view, I found there are
> normal values for these variables.
>
> 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
>
> _______________________________________________
> 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/20170829/d43fa209/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/20170829/d43fa209/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/20170829/d43fa209/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/20170829/d43fa209/attachment-0002.png>


More information about the ncl-talk mailing list