[ncl-talk] set_dimension_query
Kunal Bali
kunal.bali9 at gmail.com
Wed Oct 4 03:05:55 MDT 2017
Hello Rick sir, I tried your suggestion, but it did not work.
because lons and lats are not in the
*a*
*That's why it shows the error variable (lons) is not in file (a)*
So, I removed the part from the given script (below) and then tried. But
still not succeed in generating the netcdf file.
please let me know if you have any more suggestions.
;----------------------------------------------------------------------
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
;---Read data
a = addfile("/Users/Pushp/Desktop/test/
*MAIACTAOT.h01v03.20003660525.hdf*","r")
var = short2flt(a->Optical_Depth_055(:,:))
x = fspan(-3450000.+500., -2250000.-500., 1200)
y = fspan(400000.-500., -800000.+500., 1200)
x := conform_dims((/1200, 1200/), x, 1)
y := conform_dims((/1200, 1200/), y, 0)
src = "+proj=aea +lat_0=23.0 +lat_1=21d03m +lat_2=41 +lon_0=102
+ellps=clrk66"
dst = "+proj=latlon +ellps=sphere"
x := flt2dble(ndtooned(x))
y := flt2dble(ndtooned(y))
z = x ; A quick-n-dirty way to get an empty Z array of same
dimensions/type
z = 0. ;
transform_coordinate(src, dst, x, y, z)
; convert back to curvilinear coord arrays
lon = onedtond(x, (/1200, 1200/))
lat = onedtond(y, (/1200, 1200/))
; print(lat) ; check that these are the correct coordinates
; print(lon) ; check that these are the correct coordinates
lat at units = "degrees_north"
lon at units = "degrees_east"
lat!0 = "lat"
; lat&lat = lat
lon!0 = "lon"
; lon&lon = lon
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&lon = lon
var at lat2d = lat
var at lon2d = lon
system("/bin/rm -f simple_new.nc")
ncdf = addfile("simple_new.nc" ,"c") ; open output netCDF file
;===================================================================
; create global attributes of the file (optional)
;===================================================================
fAtt = True ; assign file attributes
fAtt at title = "NCL Simple Approach to netCDF Creation"
fAtt at source_file = "original-file.nc"
fAtt at Conventions = "None"
fAtt at creation_date = systemfunc ("date")
fileattdef( ncdf, fAtt ) ; copy file attributes
ncdf->var = var
end
Kunal Bali
Research Scholar
Radio & Atmospheric Science Division
CSIR - National Physical Laboratory
New Delhi - 110012
India
On Tue, Oct 3, 2017 at 9:24 PM, Rick Brownrigg <brownrig at ucar.edu> wrote:
> Hi Kunal,
>
> Sorry for the delay in response -- I have been travelling. I am not
> expert at the issue of 2D coordinates, and have not had a chance to
> experiment/explore the issue. My very limited understanding is that with
> 2D, you don't need to do this sort of stuff:
>
> 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&lon = lon
>
> Rather, you would write to NetCDF your variable and the 2D coordinate
> arrays computed via transform_coordinate() as just regular variables. When
> you want to read and plot that NetCDF, you need to read the variable along
> the 2D arrays, and assign them as attributes to the variable, named very
> specifically lat2d and lon2d:
>
> myVar = f->var
> lons = f->lons ; doesn't matter what the actual names are
> lats = f->lats ; ditto
>
> myVar at lat2d = lats
> myVar at lon2d = lons
>
> plot = gsn_xxxxxx(wks, myVar,....)
>
> Again, this is my limited understand and untested. There's more
> information at thes page, and in particular example #3:
>
> http://www.ncl.ucar.edu/Applications/plot_data_on_map.shtml
>
> http://www.ncl.ucar.edu/Applications/plot_data_on_map.shtml#ex3
>
>
> It looks like your 2D coordinate variables may need a "degrees_north"
> ,etc. attribute (see the example).
>
> One of my colleagues further offers this guidance:
>
> "Also to be CF compliant and for ease of use as well, all variables
> that use these coordinates should have a "coordinates" attribute that
> lists each of these coordinates by name separated by spaces. See the
> CF document for details."
>
> I hope that helps...
> Rick
>
> On Sat, Sep 30, 2017 at 1:25 PM, Kunal Bali <kunal.bali9 at gmail.com> wrote:
>
>> Thank You so much, Rick Sir for providing all these valuable information.
>> It's little complex to me but I am trying to understand.
>>
>> However, I am trying to convert this hdf file into netcdf by
>> incorporating your code. But file did not convert with Geo2D. please check
>> the code below if you have time.
>>
>> begin
>> ;---Read data
>> a = addfile("/home/kunal/mishra_sir/MAIACTAOT.h00v02.20000570505
>> .hdf","r")
>> var = short2flt(a->Optical_Depth_055(:,:))
>>
>> x = fspan(-3450000.+500., -2250000.-500., 1200)
>> y = fspan(1600000.-500., 400000.+500., 1200)
>>
>> x := conform_dims((/1200, 1200/), x, 1)
>> y := conform_dims((/1200, 1200/), y, 0)
>>
>>
>> src = "+proj=aea +lat_0=23.0 +lat_1=21d03m +lat_2=41 +lon_0=102
>> +ellps=clrk66"
>> dst = "+proj=latlon +ellps=sphere"
>> x := flt2dble(ndtooned(x))
>> y := flt2dble(ndtooned(y))
>> z = x ; A quick-n-dirty way to get an empty Z array of same
>> dimensions/type
>> z = 0. ;
>> transform_coordinate(src, dst, x, y, z)
>>
>> ; convert back to curvilinear coord arrays
>> lon = onedtond(x, (/1200, 1200/))
>> lat = onedtond(y, (/1200, 1200/))
>>
>> lat at units = "degrees_north"
>> lat!0 = "lat"
>> ; lat&lat = lat
>> lon at units = "degrees_east"
>> lon!0 = "lon"
>> ; lon&lon = lon
>>
>> 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&lon = lon
>>
>>
>>
>> ; system("/bin/rm -f simple1.nc")
>> ncdf = addfile("simple3.nc" ,"c") ; open output netCDF file
>>
>> ;===================================================================
>> ; create global attributes of the file (optional)
>> ;===================================================================
>> fAtt = True ; assign file attributes
>> fAtt at title = "NCL Simple Approach to netCDF Creation"
>> fAtt at source_file = "original-file.nc"
>> fAtt at Conventions = "None"
>> fAtt at creation_date = systemfunc ("date")
>> fileattdef( ncdf, fAtt ) ; copy file attributes
>>
>> ;===================================================================
>> ; make time an UNLIMITED dimension; recommended for most applications
>> ;===================================================================
>>
>> ncdf->var = var
>>
>>
>> end
>>
>>
>>
>>
>>
>> Kunal Bali
>>
>>
>>
>>
>> On Sat, Sep 30, 2017 at 10:15 PM, Rick Brownrigg <brownrig at ucar.edu>
>> wrote:
>>
>>> Hi,
>>>
>>> This file is an interesting and challenging one to geo-reference.
>>> Without proper documentation on how the file is constructed/encoded, all of
>>> what follows is speculation. But I'm fairly confident it is sound.
>>>
>>> In the group GRID_1, there are these parameters:
>>>
>>> GridName="grid1km"
>>> XDim=1200
>>> YDim=1200
>>> UpperLeftPointMtrs=(-3450000.000000,1600000.000000)
>>> LowerRightMtrs=(-2250000.000000,400000.000000)
>>> Projection=GCTP_ALBERS
>>> ProjParams=(0,0,21030000,41000000,102000000,23000000,0,0,0,0,0,0,0)
>>> SphereCode=12
>>> GridOrigin=HDFE_GD_UL
>>> GROUP=Dimension
>>> END_GROUP=Dimension
>>>
>>> The name suggests a 1KM grid, and as such, by its resolution, it should
>>> span 1200KM. Indeed, from the corner points:
>>>
>>> In X: (-3450000 - -2250000) = a span of 1200000meters; divided by
>>> xres=1200, yields 1000meter grid cells
>>> In Y: (1600000 - 400000) = 1200000meters; div by yres=1200 yields
>>> 1000meter grid cells
>>>
>>> So everything is consistent so far. The GCTP_ALBERS suggest the grid was
>>> constructed based upon an ALBERS projection, using the GCTP library. The
>>> ProjParams give further details:
>>>
>>> ProjParams=(0,0,21030000,41000000,102000000,23000000,0,0,0,0,0,0,0)
>>>
>>> Again, without documentation, this is speculation, but looking at
>>> pg13/appendexD in:
>>>
>>> https://www.ngs.noaa.gov/PUBS_LIB/GeneralCartographicTransfo
>>> rmationPackage_v2_TR_NOS124_CGS9.pdf
>>>
>>> for the Albers projection, the first 2 values appear to be semimajor
>>> axis of ellipsoid, and eccentricity. Both are zero, implying defaults of
>>> Clarke's 1866 ellipsoid and spherical earth (which seems to be an oymoron?
>>> I'm confused about that)
>>>
>>> The next 6 values might be:
>>>
>>> lat of 1st standard parallel
>>> lat of 2nd standard parallel
>>> lon of central meridian
>>> lat of projection's origin
>>> false easting (zero in this case)
>>> false northing (zero)
>>>
>>> the remaining zeros don't seem to apply.
>>>
>>> NCL does not have an public interface to GCTP, but there is an
>>> undocumented interface to the well-known proj4 library. The trick here is
>>> to use it to convert from meters-on-the ground units into lat/lon. The
>>> first step is to construct arrays representing the grid's position in
>>> meters. In the X direction, if we did
>>>
>>> x = fspan(-3450000., -2250000, 1200)
>>>
>>> we get an array of 1201, which is wrong. I speculate that for a grid of
>>> 1200, those values represent the total *spatial extent* of the grid, but
>>> that the data are based upon cell centers, not corners. Thus something like:
>>>
>>> x = fspan(-3450000.+500., -2250000.-500., 1200)
>>> y = fspan(1600000.-500., 400000.+500., 1200)
>>>
>>> A 1km grid cannot obviously be constant spacing in longitude, so
>>> ultimately, your coordinate arrays are going to be curvilinear:
>>>
>>> x := conform_dims((/1200, 1200/), x, 1)
>>> y := conform_dims((/1200, 1200/), y, 0)
>>>
>>> NOTE: THIS ASSUMES THE GRID IS ROW MAJOR, WEST-->EAST, NORTH-->SOUTH
>>>
>>> The undocumented NCL function "transform_coordinates" takes a proj4
>>> source-porjection string, dest-proj string, x, y, z 1D arrays of type
>>> double, and returns the results in the x,y,z arrays. I won't go into the
>>> details of the proj-strings. We want something like:
>>>
>>> src = "+proj=aea +lat_0=23.0 +lat_1=21d03m +lat_2=41 +lon_0=102
>>> +ellps=clrk66"
>>> dst = "+proj=latlon +ellps=sphere"
>>> x := flt2dble(ndtooned(x))
>>> y := flt2dble(ndtooned(y))
>>> z = x ; A quick-n-dirty way to get an empty Z array of same
>>> dimensions/type
>>> z = 0. ;
>>> transform_coordinate(src, dst, x, y, z)
>>>
>>> ; convert back to curvilinear coord arrays
>>> lon = onedtond(x, (/1200, 1200/))
>>> lat = onedtond(y, (/1200, 1200/))
>>>
>>>
>>> From here, you have to associate the lat/lon coordinate arrays with the
>>> variables you want to plot, as Adam and other previously described.
>>>
>>> Attached is a script and example plot of the resultant grid. Please let
>>> the group know if this appears to be the correct area-of-interest.
>>>
>>> Rick
>>>
>>>
>>>
>>> On Fri, Sep 29, 2017 at 11:45 AM, Rick Brownrigg <brownrig at ucar.edu>
>>> wrote:
>>>
>>>> Hi,
>>>>
>>>> I do not know for certain without examining the file, but it looks like
>>>> those values are map-projection parameters, from which curvilinear grid was
>>>> generated as equal-area grid cells in lat/lon spac.
>>>>
>>>> Rick
>>>>
>>>> On Fri, Sep 29, 2017 at 10:16 AM, Kunal Bali <kunal.bali9 at gmail.com>
>>>> wrote:
>>>>
>>>>> Thanks for all the suggestions and help Adam Sir. Now the data is
>>>>> showing over the Indian region.
>>>>> but one thing I need to confirm here.
>>>>>
>>>>> Some values in the form of UpperLeftPointMtrs and LowerRightMtrs are
>>>>> given in the file ((in the original hdf file) ) description along with X
>>>>> YDim.
>>>>>
>>>>> So, does these value interfer while creating the lat lon dimensions
>>>>> during conversion netcdf from hdf?
>>>>>
>>>>> Sould these values be consider or we can ignore these values?
>>>>>
>>>>>
>>>>> group: latlon {
>>>>> variables:
>>>>> short _HDFEOS_CRS;
>>>>> :Projection = "GCTP_ALBERS";
>>>>> :*UpperLeftPointMtrs* = -3450000.0, 1600000.0; // double
>>>>> *:LowerRightMtrs* = -2250000.0, 400000.0; // double
>>>>> :ProjParams = 0.0, 0.0, 2.103E7, 4.1E7, 1.02E8, 2.3E7, 0.0,
>>>>> 0.0, 0.0, 0.0, 0.0, 0.0, 0.0; // double
>>>>> :SphereCode = "12";
>>>>>
>>>>> group: Data_Fields {
>>>>> dimensions:
>>>>> YDim = 1200;
>>>>> XDim = 1200;
>>>>> variables:
>>>>> float lat(YDim=1200, XDim=1200);
>>>>>
>>>>> float lon(YDim=1200, XDim=1200);
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> Regards
>>>>> Kunal Bali
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> On Thu, Sep 28, 2017 at 9:27 PM, Adam Phillips <asphilli at ucar.edu>
>>>>> wrote:
>>>>>
>>>>>> Hi Kunal,
>>>>>> As the documentation states, latGlobeFo (lonGlobeFo) generates
>>>>>> latitudes (longitudes) and associated metadata for a *global* fixed
>>>>>> offset grid.
>>>>>> https://www.ncl.ucar.edu/Document/Functions/Contributed/latG
>>>>>> lobeFo.shtml
>>>>>> https://www.ncl.ucar.edu/Document/Functions/Contributed/lonG
>>>>>> lobeFo.shtml
>>>>>> If you need to create coordinates for a regional grid, then you will
>>>>>> have to do it by hand. For example, if you know your data are equally
>>>>>> spaced in the both the latitude and longitude direction from 0:40N, 60:100E
>>>>>> (for example), you can set up your coordinates like this:
>>>>>>
>>>>>> lat = fspan(0,40,1200)
>>>>>> lon = fspan(60,100,1200)
>>>>>> print(lat) ; check that these are the correct coordinates
>>>>>> print(lon) ; check that these are the correct coordinates
>>>>>> lat at units = "degrees_north"
>>>>>> lat!0 = "lat"
>>>>>> lat&lat = lat
>>>>>> lon at units = "degrees_east"
>>>>>> lon!0 = "lon"
>>>>>> lon&lon = lon
>>>>>> 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&lon = lon
>>>>>>
>>>>>> The above is similar to what latGlobeFo/lonGlobeFo do for global
>>>>>> grids. Hope that helps. If you have further questions please respond to
>>>>>> ncl-talk.
>>>>>> Adam
>>>>>>
>>>>>> On Tue, Sep 26, 2017 at 10:26 AM, Kunal Bali <kunal.bali9 at gmail.com>
>>>>>> wrote:
>>>>>>
>>>>>>> Thanks for all the suggestions Dennis sir.
>>>>>>> But now I have come to know that the dimensions
>>>>>>>
>>>>>>> nlat=1200
>>>>>>> nlon=1200
>>>>>>>
>>>>>>> are belongs to only of Indian region.
>>>>>>>
>>>>>>> And the code I am using below is not displaying the data over Indian
>>>>>>> region.
>>>>>>>
>>>>>>> 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
>>>>>>>
>>>>>>> So please help me regarding this issue.
>>>>>>>
>>>>>>> regards
>>>>>>> Kunal Bali
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> On Tue, Sep 19, 2017 at 10:36 PM, Dennis Shea <shea at ucar.edu> wrote:
>>>>>>>
>>>>>>>> Ditto what Marston posted.
>>>>>>>>
>>>>>>>> ===
>>>>>>>> Please look at NCL's FAQ.
>>>>>>>>
>>>>>>>> [1] Click "Support" at top of page
>>>>>>>> [2] Click FAQ
>>>>>>>>
>>>>>>>> I think it would be instructive to look at all of the topics.
>>>>>>>>
>>>>>>>> For example:
>>>>>>>>
>>>>>>>> *Error messages and other issues*
>>>>>>>> Look at #7
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> On Tue, Sep 19, 2017 at 3:03 AM, Marston Johnston <
>>>>>>>> shejo284 at gmail.com> wrote:
>>>>>>>>
>>>>>>>>> It tells you in the warning message, exactly what to do to remove
>>>>>>>>> the warning.
>>>>>>>>>
>>>>>>>>> If you want to improve you coding skills, it would better to
>>>>>>>>> understand why the warning is occurring in the first place: dimension 0 and
>>>>>>>>> 1 of the variable you are passing has no names or the names do not match
>>>>>>>>> the variable on the lhs.
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> /M
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
>>>>>>>>>
>>>>>>>>> Marston S. Ward, PhD
>>>>>>>>>
>>>>>>>>> Department of Earth Sciences
>>>>>>>>>
>>>>>>>>> University of Gothenburg, Sweden
>>>>>>>>>
>>>>>>>>> Email: marston.johnston at gu.se
>>>>>>>>>
>>>>>>>>> SkypeID: marston.johnston
>>>>>>>>>
>>>>>>>>> Phone: +46-31-7864901 <+46%2031%20786%2049%2001>
>>>>>>>>>
>>>>>>>>> Only the fruitful thing is true!
>>>>>>>>>
>>>>>>>>> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> *From: *ncl-talk <ncl-talk-bounces at ucar.edu> on behalf of Kunal
>>>>>>>>> Bali <kunal.bali9 at gmail.com>
>>>>>>>>> *Date: *Tuesday, 19 September 2017 at 10:50
>>>>>>>>> *Cc: *"ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
>>>>>>>>> *Subject: *Re: [ncl-talk] set_dimension_query
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Thank you for providing this information.
>>>>>>>>>
>>>>>>>>> It's working. The script given below is now able to convert
>>>>>>>>> multiple files into .netcdf with the desired file name.
>>>>>>>>> However, on the same time, it also shows some warnings. But noted
>>>>>>>>> that, I am successfully creating the multiple files into the netcdf. I just
>>>>>>>>> need to know that how can I remove these warnings?
>>>>>>>>>
>>>>>>>>> warning:VarVarWrite: Dimension names for dimension number (0)
>>>>>>>>> don't match, assigning name of rhs dimension to lhs and overwriting
>>>>>>>>> coordinate variable, use "(/../)" if this change is not desired
>>>>>>>>> warning:VarVarWrite: Dimension names for dimension number (1)
>>>>>>>>> don't match, assigning name of rhs dimension to lhs and overwriting
>>>>>>>>> coordinate variable, use "(/../)" if this change is not desired
>>>>>>>>> warning:["Execute.c":8640]:Execute: Error occurred at or near
>>>>>>>>> line 22 in file netcdf2.ncl
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> ;-----------------------------------------------------------
>>>>>>>>> -----------
>>>>>>>>>
>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>>>>>>>>
>>>>>>>>> begin
>>>>>>>>> ;---Read data
>>>>>>>>>
>>>>>>>>> diri = "/Users/Pushp/Desktop/test/" ; input directory
>>>>>>>>> fili = systemfunc("cd "+diri+" ; ls MAIACTAOT.h00v02*hdf")
>>>>>>>>> nfili = dimsizes(fili)
>>>>>>>>> print("nfili="+nfili)
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> dirnc = "/Users/Pushp/Desktop/test/" ; output (netCDF)
>>>>>>>>> directory
>>>>>>>>>
>>>>>>>>> do nf=0,nfili-1
>>>>>>>>> pthi = diri+fili(nf)
>>>>>>>>> f = addfile(pthi,"r")
>>>>>>>>>
>>>>>>>>> var = short2flt(f->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&lon = lon
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> filroot = str_get_cols(fili(nf), 0, 27) ; eg:
>>>>>>>>> "MAIACTAOT.h00v02.20000570505"
>>>>>>>>> filnc = filroot+".nc"
>>>>>>>>> pthnc = dirnc + filnc
>>>>>>>>>
>>>>>>>>> system("/bin/rm -f "+pthnc)
>>>>>>>>> ncdf = addfile(pthnc,"c") ; open new netCDF file
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> end do ; end 'nf' loop
>>>>>>>>>
>>>>>>>>> ; system("/bin/rm -f simple2.nc")
>>>>>>>>> ; ncdf = addfile("simple2.nc" ,"c") ; open output netCDF
>>>>>>>>> file
>>>>>>>>>
>>>>>>>>> ;===========================================================
>>>>>>>>> ========
>>>>>>>>> ; create global attributes of the file (optional)
>>>>>>>>> ;===========================================================
>>>>>>>>> ========
>>>>>>>>> fAtt = True ; assign file
>>>>>>>>> attributes
>>>>>>>>> fAtt at title = "NCL Simple Approach to netCDF
>>>>>>>>> Creation"
>>>>>>>>> fAtt at source_file = "original-file.nc"
>>>>>>>>> fAtt at Conventions = "None"
>>>>>>>>> fAtt at creation_date = systemfunc ("date")
>>>>>>>>> fileattdef( ncdf, fAtt ) ; copy file attributes
>>>>>>>>>
>>>>>>>>> ;===========================================================
>>>>>>>>> ========
>>>>>>>>> ; make time an UNLIMITED dimension; recommended for most
>>>>>>>>> applications
>>>>>>>>> ;===========================================================
>>>>>>>>> ========
>>>>>>>>> filedimdef(ncdf,"time",-1,True)
>>>>>>>>>
>>>>>>>>> ; ncdf->var&lat = lat
>>>>>>>>> ; ncdf->var&lon = lon
>>>>>>>>> ncdf->var = var
>>>>>>>>>
>>>>>>>>> end
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Kunal Bali
>>>>>>>>>
>>>>>>>>> Research Scholar
>>>>>>>>>
>>>>>>>>> Radio & Atmospheric Science Division
>>>>>>>>>
>>>>>>>>> CSIR - National Physical Laboratory
>>>>>>>>>
>>>>>>>>> New Delhi - 110012
>>>>>>>>>
>>>>>>>>> India
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Mon, Sep 18, 2017 at 9:04 PM, Dennis Shea <shea at ucar.edu>
>>>>>>>>> wrote:
>>>>>>>>>
>>>>>>>>> Loop over the files ...
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> diri = "/media/Local Disk/NPL/MODIS_FPC/" ; input directory
>>>>>>>>>
>>>>>>>>> fili = systemfunc("cd "+diri +; ls MAIACTAOT.h00v02*hdf")
>>>>>>>>>
>>>>>>>>> nfili = dimsizes(fili)
>>>>>>>>>
>>>>>>>>> print("nfili="+nfili)
>>>>>>>>>
>>>>>>>>> dirnc = "/media/Local Disk/NPL/MODIS_FPC/" ; output (netCDF)
>>>>>>>>> directory
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> do nf=0,nfili-1
>>>>>>>>>
>>>>>>>>> pthi = diri+fili(nf)
>>>>>>>>>
>>>>>>>>> f = addfile(pthi,"r")
>>>>>>>>>
>>>>>>>>> ....
>>>>>>>>>
>>>>>>>>> filroot = *str_get_cols*(fili(nf), 0, 27) ; eg:
>>>>>>>>> "MAIACTAOT.h00v02.20000570505"
>>>>>>>>>
>>>>>>>>> filnc = filroot+".nc"
>>>>>>>>>
>>>>>>>>> pthnc = dirnc + filnc
>>>>>>>>>
>>>>>>>>> system("/bin/rm -f "+pthnc)
>>>>>>>>> ncdf = addfile(pthnc,"c") ; open new netCDF file
>>>>>>>>>
>>>>>>>>> ....
>>>>>>>>>
>>>>>>>>> end do ; end 'nf' loop
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> =====
>>>>>>>>>
>>>>>>>>> If there is something you do not understand use some print
>>>>>>>>> statements.
>>>>>>>>>
>>>>>>>>> Become familiar with the 'string' category. Look at all the
>>>>>>>>> functions.
>>>>>>>>>
>>>>>>>>> https://www.ncl.ucar.edu/Document/Functions/string.shtml
>>>>>>>>>
>>>>>>>>> Please read the documentation for the function used above
>>>>>>>>>
>>>>>>>>> https://www.ncl.ucar.edu/Document/Functions/Built-in/str_get
>>>>>>>>> _cols.shtml
>>>>>>>>>
>>>>>>>>> Good Luck
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Mon, Sep 18, 2017 at 6:11 AM, Kunal Bali <kunal.bali9 at gmail.com>
>>>>>>>>> wrote:
>>>>>>>>>
>>>>>>>>> Thnks for the reply.
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> that will change the name of only one file. isn't it. but what
>>>>>>>>> about the other files
>>>>>>>>>
>>>>>>>>> I have 365 files in one directory and I need to change all the
>>>>>>>>> file name one by one.
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Kunal Bali
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Mon, Sep 18, 2017 at 5:14 PM, Dennis Shea <shea at ucar.edu>
>>>>>>>>> wrote:
>>>>>>>>>
>>>>>>>>> Really?
>>>>>>>>>
>>>>>>>>> Change:
>>>>>>>>> system("/bin/rm -f simple.nc")
>>>>>>>>> ncdf = addfile("simple.nc" ,"c")
>>>>>>>>>
>>>>>>>>> To:
>>>>>>>>> system("/bin/rm -f MAIACTAOT.h00v02.20000570505.nc
>>>>>>>>> <http://simple.nc>")
>>>>>>>>> ncdf = addfile("MAIACTAOT.h00v02.20000570505.nc" ,"c")
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Mon, Sep 18, 2017 at 12:37 AM, Kunal Bali <
>>>>>>>>> kunal.bali9 at gmail.com> wrote:
>>>>>>>>>
>>>>>>>>> Thanks, it worked.
>>>>>>>>>
>>>>>>>>> one more question is that.
>>>>>>>>>
>>>>>>>>> I used the code (given below). So it gives simple.nc name as an
>>>>>>>>> output file.
>>>>>>>>>
>>>>>>>>> ncdf = addfile("simple.nc" ,"c")
>>>>>>>>>
>>>>>>>>> If I want to keep the original file name with the output file name
>>>>>>>>> then what should I do?
>>>>>>>>>
>>>>>>>>> I mean the original file name is MAIACTAOT.h00v02.20003660700.hdf
>>>>>>>>>
>>>>>>>>> Now I want to create my netcdf file name as
>>>>>>>>> MAIACTAOT.h00v02.20003660700.nc NOT simple.nc.
>>>>>>>>>
>>>>>>>>> I can not simply write simple.nc file name because I have many
>>>>>>>>> files.
>>>>>>>>>
>>>>>>>>> e.g I have one directory having 365 .hdf file. So I need to
>>>>>>>>> convert all the hdf file to netcdf file at once with the same name of
>>>>>>>>> original file names.
>>>>>>>>>
>>>>>>>>> I hope you have got my query.
>>>>>>>>>
>>>>>>>>> please let me know that too.
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Thank You
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Kunal Bali
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Thu, Sep 14, 2017 at 9:09 PM, Adam Phillips <asphilli at ucar.edu>
>>>>>>>>> wrote:
>>>>>>>>>
>>>>>>>>> Hi Kunal,
>>>>>>>>>
>>>>>>>>> It looks to me like your latitudes are flipped. You are setting up
>>>>>>>>> and assigning your coordinate lat/lon variables in this coding:
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> 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&lon = lon
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Assuming I'm right on your latitudes being the issue, one of the
>>>>>>>>> following two modifications should work:
>>>>>>>>>
>>>>>>>>> Delete this line:
>>>>>>>>>
>>>>>>>>> lat = lat(::-1)
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> If the above doesn't fix it, try this:
>>>>>>>>>
>>>>>>>>> Change this:
>>>>>>>>>
>>>>>>>>> lat = lat(::-1)
>>>>>>>>>
>>>>>>>>> to this:
>>>>>>>>>
>>>>>>>>> lat = lat(::-1)
>>>>>>>>>
>>>>>>>>> lat&lat = lat
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> If you continue to have issues please respond to ncl-talk.
>>>>>>>>>
>>>>>>>>> Adam
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Wed, Sep 13, 2017 at 12:35 PM, Kunal Bali <
>>>>>>>>> kunal.bali9 at gmail.com> wrote:
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> For direct output as a netcdf format, I used the script given
>>>>>>>>> below. It produced the netcdf file easily. But the netcdf file and original
>>>>>>>>> hdf file both showing the different results. The data pattern is shifted. I
>>>>>>>>> mean it may be related to the lat lon position. please see the attached
>>>>>>>>> file, you will understand. And please let me know which one is correct.
>>>>>>>>>
>>>>>>>>> ;-----------------------------------------------------------
>>>>>>>>> -----------
>>>>>>>>>
>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>>>>>>>>>
>>>>>>>>> begin
>>>>>>>>> ;---Read data
>>>>>>>>> a = addfile("/media/Local Disk/NPL/MODIS_FPC/MAIACTAOT.h
>>>>>>>>> 00v02.20000570505.hdf","r")
>>>>>>>>> var = short2flt(a->Optical_Depth_055(:,:))
>>>>>>>>> ; short2flt(a[:]->noxfire)
>>>>>>>>> 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&lon = lon
>>>>>>>>>
>>>>>>>>> system("/bin/rm -f simple.nc")
>>>>>>>>> ncdf = addfile("simple.nc" ,"c") ; open output
>>>>>>>>> netCDF file
>>>>>>>>>
>>>>>>>>> ;===========================================================
>>>>>>>>> ========
>>>>>>>>> ; create global attributes of the file (optional)
>>>>>>>>> ;===========================================================
>>>>>>>>> ========
>>>>>>>>> fAtt = True ; assign file
>>>>>>>>> attributes
>>>>>>>>> fAtt at title = "NCL Simple Approach to netCDF
>>>>>>>>> Creation"
>>>>>>>>> fAtt at source_file = "original-file.nc"
>>>>>>>>> fAtt at Conventions = "None"
>>>>>>>>> fAtt at creation_date = systemfunc ("date")
>>>>>>>>> fileattdef( ncdf, fAtt ) ; copy file attributes
>>>>>>>>>
>>>>>>>>> ;===========================================================
>>>>>>>>> ========
>>>>>>>>> ; make time an UNLIMITED dimension; recommended for most
>>>>>>>>> applications
>>>>>>>>> ;===========================================================
>>>>>>>>> ========
>>>>>>>>> filedimdef(ncdf,"time",-1,True)
>>>>>>>>> ncdf->var = var
>>>>>>>>>
>>>>>>>>> end
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> Kunal Bali
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> On Wed, Sep 13, 2017 at 11:27 PM, Kunal Bali <
>>>>>>>>> kunal.bali9 at gmail.com> wrote:
>>>>>>>>>
>>>>>>>>> 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
>>>>>>>>> <%28303%29%20497-1726>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> _______________________________________________
>>>>>>>>> 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
>>>>>>>>> <%28303%29%20497-1726>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> _______________________________________________
>>>>>>>>> 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
>>>>>>>>> <%28303%29%20497-1726>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> _______________________________________________
>>>>>>>>> ncl-talk mailing list
>>>>>>>>> ncl-talk at ucar.edu
>>>>>>>>> List instructions, subscriber options, unsubscribe:
>>>>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> _______________________________________________
>>>>>>>>> ncl-talk mailing list
>>>>>>>>> ncl-talk at ucar.edu
>>>>>>>>> List instructions, subscriber options, unsubscribe:
>>>>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> _______________________________________________ ncl-talk mailing
>>>>>>>>> list ncl-talk at ucar.edu List instructions, subscriber options,
>>>>>>>>> unsubscribe: http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>>>>>
>>>>>>>>> _______________________________________________
>>>>>>>>> ncl-talk mailing list
>>>>>>>>> ncl-talk at ucar.edu
>>>>>>>>> List instructions, subscriber options, unsubscribe:
>>>>>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>>>>>
>>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>> _______________________________________________
>>>>>>> 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
>>>>>> <%28303%29%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
>>>>>
>>>>>
>>>>
>>>
>>
>> _______________________________________________
>> 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/20171004/ae47c800/attachment.html>
More information about the ncl-talk
mailing list