[ncl-talk] set_dimension_query

Kunal Bali kunal.bali9 at gmail.com
Wed Sep 13 11:57:45 MDT 2017 

No problem, I sorted out.

Just changed short to float.

Kunal Bali
Research Scholar
Radio & Atmospheric Science Division
CSIR - National Physical Laboratory
New Delhi - 110012
India





On Wed, Sep 13, 2017 at 11:24 PM, Kunal Bali <kunal.bali9 at gmail.com> wrote:

> Thanks for providing this information.
>
> Also, I would like to mention that when I plot this data. The values are
> not in the domain. I mean the real values lie in-between 0 to 1 but here it
> is reaching to 400.
>
> So, how to correct the values?
>
> Description of the variable is
>
> short Optical_Depth_055(YDim=1200, XDim=1200);
>   :long_name = "AOT at 0.55 micron";
>   :scale_factor = 0.001; // double
>   :add_offset = 0.0; // double
>   :unit = "None";
>   :_FillValue = -28672S; // short
>   :valid_range = -100S, 5000S; // short
>
>
>
>
> Kunal Bali
>
>
>
>
>
>
> On Wed, Sep 13, 2017 at 10:28 PM, Adam Phillips <asphilli at ucar.edu> wrote:
>
>> Hi Kunal,
>> Yes, clicking on the output netCDF link from the Applications page:
>> http://www.ncl.ucar.edu/Applications/o-netcdf.shtml
>> I would recommend following the inefficient method #1. Unless you are
>> writing a file with many large variables, the inefficient method works just
>> fine.
>> http://www.ncl.ucar.edu/Applications/method_1.shtml
>> Good luck,
>> Adam
>>
>> On Wed, Sep 13, 2017 at 10:25 AM, Kunal Bali <kunal.bali9 at gmail.com>
>> wrote:
>>
>>> Thank you so much, it worked
>>>
>>> I want to know one more thing.
>>> After rearranging the dimensions can we now convert (or write) this
>>> arranged file into netcdf format?
>>>
>>>
>>>
>>> Kunal Bali
>>>
>>>
>>>
>>>
>>>
>>>
>>> On Wed, Sep 13, 2017 at 8:55 PM, Adam Phillips <asphilli at ucar.edu>
>>> wrote:
>>>
>>>> Hi Kunal,
>>>> I think you just need to rearrange the order of your lines and tweak a
>>>> couple of lines. As the error message states, lat is not defined in your
>>>> 3rd line and you are referring to it as if it is.Try this:
>>>> a = addfile("/media/Local Disk/NPL/MODIS_FPC/MAIACTAOT.h
>>>> 00v02.20000570505.hdf","r")
>>>> var  = a->Optical_Depth_055(:,:)
>>>> nlat = 1200
>>>> nlon = 1200
>>>> lat = latGlobeFo(nlat, "lat", "latitude", "degrees_north")
>>>> lon = lonGlobeFo(nlon, "lon", "longitude", "degrees_east")
>>>> lat = lat(::-1)
>>>> lon = (/ lon - 180. /)  ; subtract 180 from all values
>>>> lon&lon = lon           ; update coordinates
>>>>
>>>> var!0 = "lat"    ; you can name these dimensions whatever you want..
>>>> var!1 = "lon"
>>>> var&lat = lat     ; but make sure you refer to the correct named
>>>> dimensions
>>>> var&lat = lon
>>>> ;var&XDim_grid1km = lat
>>>> ;var&YDim_grid1km = lon
>>>>
>>>> If you have any further questions please respond to the ncl-talk email
>>>> list.
>>>> Adam
>>>>
>>>>
>>>>
>>>> On Wed, Sep 13, 2017 at 8:38 AM, Kunal Bali <kunal.bali9 at gmail.com>
>>>> wrote:
>>>>
>>>>> Dear NCL
>>>>>
>>>>> I have a file with the variable summary. The dimensions are in 2D.
>>>>>
>>>>> ncl 2>  printVarSummary(var)
>>>>>
>>>>> Variable: var
>>>>> Type: short
>>>>> Total Size: 2880000 bytes
>>>>>             1440000 values
>>>>> Number of Dimensions: 2
>>>>> Dimensions and sizes:    [*YDim_grid1km | 1200] x [XDim_grid1km |
>>>>> 1200]*
>>>>> Coordinates:
>>>>> Number Of Attributes: 7
>>>>>   long_name :    AOT at 0.55 micron
>>>>>   scale_factor :    0.001
>>>>>   add_offset :       0
>>>>>   unit :    None
>>>>>   _FillValue :    -28672
>>>>>   valid_range :    ( -100, 5000 )
>>>>>   hdf_name :    Optical_Depth_055
>>>>>
>>>>>
>>>>> I am trying to read the dimensions of this file as
>>>>>
>>>>> begin
>>>>> ;---Read data
>>>>>          a = addfile("/media/Local Disk/NPL/MODIS_FPC/MAIACTAOT.h
>>>>> 00v02.20000570505.hdf","r")
>>>>>
>>>>>             var  = a->Optical_Depth_055(:,:)
>>>>>
>>>>>          var&XDim_grid1km = lat
>>>>>          var&YDim_grid1km = lon
>>>>>
>>>>>         nlat = 1200
>>>>>          nlon = 1200
>>>>>             lat = latGlobeFo(nlat, "lat", "latitude", "degrees_north")
>>>>>         lon = lonGlobeFo(nlon, "lon", "longitude", "degrees_east")
>>>>>             lat = lat(::-1)
>>>>>             lon = (/ lon - 180. /)  ; subtract 180 from all values
>>>>>             lon&lon = lon           ; update coordinates
>>>>>
>>>>>
>>>>>         var!0 = "lat"
>>>>>         var!1 = "lon"
>>>>>
>>>>>
>>>>> But the error appeared as
>>>>> *fatal:Variable (lat) is undefined*
>>>>>
>>>>>
>>>>> So, could anyone please let me know that how to read dimension of this
>>>>> file.
>>>>>
>>>>> Thank You
>>>>>
>>>>>
>>>>> Regards
>>>>> Kunal Bali
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> ncl-talk at ucar.edu
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>
>>>>>
>>>>
>>>>
>>>> --
>>>> Adam Phillips
>>>> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
>>>> www.cgd.ucar.edu/staff/asphilli/   303-497-1726 <(303)%20497-1726>
>>>>
>>>> <http://www.cgd.ucar.edu/staff/asphilli>
>>>>
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>>
>> --
>> Adam Phillips
>> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
>> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>>
>> <http://www.cgd.ucar.edu/staff/asphilli>
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170913/d74213ff/attachment.html>

More information about the ncl-talk mailing list