<div dir="ltr"><div><div>Dear David Sir, Thank you so much for this information. It worked nicely. Now I am able to create and plot <b>(PLOT_1)</b> netcdf file easily. <br><br></div>But, I also got some suggestions by Rick sir and the script which he attached from that script my plot <b>(PLOT_2)</b> just shifted. So I am little confused that which one is correct. <br><br></div>And the grid_1km structure for the given file also attached. (( the hdf file was attached in previous mail)). MAIACTAOT.<b>h00v02</b>.20000570505.hdf<br><div><br></div><div>the script I used for creating netcdf is also attached.<br></div><div><br></div><div>please provide some information on that as well. <br></div><div><br> </div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Kunal Bali<br></div><br><div><br></div><div><br></div><div><br></div><div><p style="margin:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></p></div></div></div></div></div></div>
<br><div class="gmail_quote">On Thu, Oct 5, 2017 at 1:57 AM, Rick Brownrigg <span dir="ltr"><<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">HI,<div><br></div><div>Just add to what Dave said:</div><div><br></div><div>Attached is a script that shows how to put these things together when plotting. The script generates a contour plot of the variable "var" over a map. Note that most of the data are missing, and only a triangular shaped piece shows up in the southwest corner of the overall grid. You can see this by commenting out the lines noted in the script (BTW: this is true for the file you originally sent me "MAIACTAOT.h00<font color="#ff0000">v02</font>.20000570505.<wbr>hdf, I don't know about any others).</div><div><br></div><div>You also may still have a scaling issue to deal with in the conversion script. Note that the original variable that "var" is derived from contains scaling and offset attributes, The offset is zero, but I **think** (not certain) you want to multiply the data by that value. </div><span class="gmail-HOEnZb"><font color="#888888"><div><br></div><div>Rick</div></font></span></div><div class="gmail-HOEnZb"><div class="gmail-h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Oct 4, 2017 at 11:08 AM, David Brown <span dir="ltr"><<a href="mailto:dbrown@ucar.edu" target="_blank">dbrown@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Kunai,<br>
The lat2d and lon2d attributes are only recognized by the NCL graphics<br>
routines. NetCDF only supports attributes with no more than 1<br>
dimension. When you have 2d coordinates you should just write them<br>
directly to the NetCDF file as separate variables. Furthermore your 2d<br>
coordinates should have the same dimension names as the data variable<br>
and to be CF compliant the data variable should have a coordinates<br>
attribute that lists the names of the coordinate variables.<br>
<br>
So if you want the coordinates to be named "lat" and "lon", then you<br>
should introduce two new dimension names, for example, "x" and "y",<br>
that will be shared by the data variable and the coordinate variables:<br>
<br>
So replace this code:<br>
<span><br>
 lat!0 = "lat"<br>
 ; lat&lat = lat<br>
 lon!0 = "lon"<br>
; lon&lon = lon<br>
<br>
 var!0 = "lat"    ; you can name these dimensions whatever you want..<br>
 var!1 = "lon"<br>
; var&lat = lat     ; but make sure you refer to the correct named dimensions<br>
; var&lon = lon<br>
<br>
 var@lat2d = lat<br>
 var@lon2d = lon<br>
<br>
</span>with<br>
<br>
lat!0 = "y"<br>
lat!1 = "x"<br>
lon!0 = "y"<br>
lon!1 = "x"<br>
var!0 = "x'<br>
var!1 = "y"<br>
<br>
var@coordinates = "lat lon"<br>
<br>
then add after (or before - it shouldn't matter) this line:<br>
<br>
ncdf->var = var<br>
<br>
ncdf->lat = lat<br>
ncfd->lon = lon<br>
<br>
<br>
Hope this helps.<br>
 -dave<br>
<div class="gmail-m_-2667434472774276178HOEnZb"><div class="gmail-m_-2667434472774276178h5"><br>
<br>
<br>
<br>
On Wed, Oct 4, 2017 at 3:05 AM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<br>
> Hello Rick sir, I tried your suggestion, but it did not work.<br>
><br>
> because lons and lats are not in the a<br>
> That's why it shows the error variable (lons) is not in file (a)<br>
> So, I removed the part from the given script (below) and then tried. But<br>
> still not succeed in generating the netcdf file.<br>
><br>
> please let me know if you have any more suggestions.<br>
><br>
> ;-----------------------------<wbr>------------------------------<wbr>-----------<br>
><br>
> load "$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/gsn_code.ncl"<br>
> load "$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/gsn_csm.ncl"<br>
> load "$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/contributed.ncl"<br>
><br>
> begin<br>
> ;---Read data<br>
>       a =<br>
> addfile("/Users/Pushp/Desktop/<wbr>test/MAIACTAOT.h01v03.20003660<wbr>525.hdf","r")<br>
>       var  = short2flt(a->Optical_Depth_055<wbr>(:,:))<br>
><br>
>       x = fspan(-3450000.+500., -2250000.-500., 1200)<br>
>       y = fspan(400000.-500., -800000.+500., 1200)<br>
><br>
>       x := conform_dims((/1200, 1200/), x, 1)<br>
>       y := conform_dims((/1200, 1200/), y, 0)<br>
><br>
><br>
>       src = "+proj=aea +lat_0=23.0 +lat_1=21d03m +lat_2=41 +lon_0=102<br>
> +ellps=clrk66"<br>
>       dst = "+proj=latlon +ellps=sphere"<br>
>       x := flt2dble(ndtooned(x))<br>
>       y := flt2dble(ndtooned(y))<br>
>       z = x   ; A quick-n-dirty way to get an empty Z array of same<br>
> dimensions/type<br>
>       z = 0.   ;<br>
>       transform_coordinate(src, dst, x, y, z)<br>
><br>
>     ; convert back to curvilinear coord arrays<br>
>       lon = onedtond(x, (/1200, 1200/))<br>
>       lat = onedtond(y, (/1200, 1200/))<br>
><br>
>     ; print(lat)   ; check that these are the correct coordinates<br>
>     ; print(lon)   ; check that these are the correct coordinates<br>
><br>
>       lat@units = "degrees_north"<br>
>       lon@units = "degrees_east"<br>
><br>
>       lat!0 = "lat"<br>
>     ; lat&lat = lat<br>
>       lon!0 = "lon"<br>
>     ; lon&lon = lon<br>
><br>
>       var!0 = "lat"    ; you can name these dimensions whatever you want..<br>
>       var!1 = "lon"<br>
>     ; var&lat = lat     ; but make sure you refer to the correct named<br>
> dimensions<br>
>     ; var&lon = lon<br>
><br>
>       var@lat2d = lat<br>
>       var@lon2d = lon<br>
><br>
>       system("/bin/rm -f <a href="http://simple_new.nc" rel="noreferrer" target="_blank">simple_new.nc</a>")<br>
>       ncdf = addfile("<a href="http://simple_new.nc" rel="noreferrer" target="_blank">simple_new.nc</a>" ,"c")  ; open output netCDF file<br>
><br>
>     ;============================<wbr>==============================<wbr>=========<br>
>     ; create global attributes of the file (optional)<br>
>     ;============================<wbr>==============================<wbr>=========<br>
>        fAtt               = True            ; assign file attributes<br>
>        fAtt@title         = "NCL Simple Approach to netCDF Creation"<br>
>        fAtt@source_file   =  "<a href="http://original-file.nc" rel="noreferrer" target="_blank">original-file.nc</a>"<br>
>        fAtt@Conventions   = "None"<br>
>        fAtt@creation_date = systemfunc ("date")<br>
>        fileattdef( ncdf, fAtt )            ; copy file attributes<br>
><br>
>        ncdf->var = var<br>
><br>
> end<br>
><br>
><br>
> Kunal Bali<br>
> Research Scholar<br>
> Radio & Atmospheric Science Division<br>
> CSIR - National Physical Laboratory<br>
> New Delhi - 110012<br>
> India<br>
><br>
><br>
><br>
><br>
><br>
> On Tue, Oct 3, 2017 at 9:24 PM, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>> wrote:<br>
>><br>
>> Hi Kunal,<br>
>><br>
>> Sorry for the delay in response -- I have been travelling.  I am not<br>
>> expert at the issue of 2D coordinates, and have not had a chance to<br>
>> experiment/explore the issue. My very limited understanding is that with 2D,<br>
>> you don't need to do this sort of stuff:<br>
>><br>
>>      var!0 = "lat"    ; you can name these dimensions whatever you want..<br>
>>         var!1 = "lon"<br>
>> ;        var&lat = lat     ; but make sure you refer to the correct named<br>
>> dimensions<br>
>> ;        var&lon = lon<br>
>><br>
>> Rather, you would write to NetCDF your variable and the 2D coordinate<br>
>> arrays computed via transform_coordinate() as just regular variables.  When<br>
>> you want to read and plot that NetCDF, you need to read the variable along<br>
>> the 2D arrays, and assign them as attributes to the variable, named very<br>
>> specifically lat2d and lon2d:<br>
>><br>
>>      myVar = f->var<br>
>>      lons = f->lons   ; doesn't matter what the actual names are<br>
>>     lats = f->lats      ;   ditto<br>
>><br>
>>     myVar@lat2d = lats<br>
>>     myVar@lon2d = lons<br>
>><br>
>>     plot = gsn_xxxxxx(wks, myVar,....)<br>
>><br>
>> Again, this is my limited understand and untested.  There's more<br>
>> information at thes page, and in particular example #3:<br>
>><br>
>> <a href="http://www.ncl.ucar.edu/Applications/plot_data_on_map.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/plot_data_on_map.shtml</a><br>
>><br>
>> <a href="http://www.ncl.ucar.edu/Applications/plot_data_on_map.shtml#ex3" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/plot_data_on_map.shtml#<wbr>ex3</a><br>
>><br>
>><br>
>> It looks like your 2D coordinate variables may need a "degrees_north"<br>
>> ,etc. attribute (see the example).<br>
>><br>
>> One of my colleagues further offers this guidance:<br>
>><br>
>> "Also to be CF compliant and for ease of use as well, all variables<br>
>> that use these coordinates should have a "coordinates" attribute that<br>
>> lists each of these coordinates by name separated by spaces. See the<br>
>> CF document for details."<br>
>><br>
>> I hope that helps...<br>
>> Rick<br>
>><br>
>> On Sat, Sep 30, 2017 at 1:25 PM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<br>
>>><br>
>>> Thank You so much, Rick Sir for providing all these valuable information.<br>
>>> It's little complex to me but I am trying to understand.<br>
>>><br>
>>> However, I am trying to convert this hdf file into netcdf by<br>
>>> incorporating your code. But file did not convert with Geo2D. please check<br>
>>> the code below if you have time.<br>
>>><br>
>>> begin<br>
>>> ;---Read data<br>
>>>         a =<br>
>>> addfile("/home/kunal/mishra_si<wbr>r/MAIACTAOT.h00v02.20000570505<wbr>.hdf","r")<br>
>>>         var  = short2flt(a->Optical_Depth_055<wbr>(:,:))<br>
>>><br>
>>>       x = fspan(-3450000.+500., -2250000.-500., 1200)<br>
>>>       y = fspan(1600000.-500., 400000.+500., 1200)<br>
>>><br>
>>>       x := conform_dims((/1200, 1200/), x, 1)<br>
>>>       y := conform_dims((/1200, 1200/), y, 0)<br>
>>><br>
>>><br>
>>>       src = "+proj=aea +lat_0=23.0 +lat_1=21d03m +lat_2=41 +lon_0=102<br>
>>> +ellps=clrk66"<br>
>>>       dst = "+proj=latlon +ellps=sphere"<br>
>>>       x := flt2dble(ndtooned(x))<br>
>>>       y := flt2dble(ndtooned(y))<br>
>>>       z = x   ; A quick-n-dirty way to get an empty Z array of same<br>
>>> dimensions/type<br>
>>>       z = 0.   ;<br>
>>>       transform_coordinate(src, dst, x, y, z)<br>
>>><br>
>>>       ; convert back to curvilinear coord arrays<br>
>>>       lon = onedtond(x, (/1200, 1200/))<br>
>>>       lat = onedtond(y, (/1200, 1200/))<br>
>>><br>
>>>      lat@units = "degrees_north"<br>
>>>     lat!0 = "lat"<br>
>>> ;    lat&lat = lat<br>
>>>     lon@units = "degrees_east"<br>
>>>     lon!0 = "lon"<br>
>>> ;    lon&lon = lon<br>
>>><br>
>>>         var!0 = "lat"    ; you can name these dimensions whatever you<br>
>>> want..<br>
>>>         var!1 = "lon"<br>
>>> ;        var&lat = lat     ; but make sure you refer to the correct named<br>
>>> dimensions<br>
>>> ;        var&lon = lon<br>
>>><br>
>>><br>
>>><br>
>>>             ;     system("/bin/rm -f <a href="http://simple1.nc" rel="noreferrer" target="_blank">simple1.nc</a>")<br>
>>>              ncdf = addfile("<a href="http://simple3.nc" rel="noreferrer" target="_blank">simple3.nc</a>" ,"c")  ; open output netCDF file<br>
>>><br>
>>>     ;============================<wbr>==============================<wbr>=========<br>
>>>     ; create global attributes of the file (optional)<br>
>>>     ;============================<wbr>==============================<wbr>=========<br>
>>>        fAtt               = True            ; assign file attributes<br>
>>>        fAtt@title         = "NCL Simple Approach to netCDF Creation"<br>
>>>        fAtt@source_file   =  "<a href="http://original-file.nc" rel="noreferrer" target="_blank">original-file.nc</a>"<br>
>>>        fAtt@Conventions   = "None"<br>
>>>        fAtt@creation_date = systemfunc ("date")<br>
>>>        fileattdef( ncdf, fAtt )            ; copy file attributes<br>
>>><br>
>>>     ;============================<wbr>==============================<wbr>=========<br>
>>>     ; make time an UNLIMITED dimension; recommended  for most<br>
>>> applications<br>
>>>     ;============================<wbr>==============================<wbr>=========<br>
>>><br>
>>>        ncdf->var = var<br>
>>><br>
>>><br>
>>> end<br>
>>><br>
>>><br>
>>><br>
>>><br>
>>><br>
>>> Kunal Bali<br>
>>><br>
>>><br>
>>><br>
>>><br>
>>> On Sat, Sep 30, 2017 at 10:15 PM, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>><br>
>>> wrote:<br>
>>>><br>
>>>> Hi,<br>
>>>><br>
>>>> This file is an interesting and challenging one to geo-reference.<br>
>>>> Without proper documentation on how the file is constructed/encoded, all of<br>
>>>> what follows is speculation. But I'm fairly confident it is sound.<br>
>>>><br>
>>>> In the group GRID_1, there are these parameters:<br>
>>>><br>
>>>>   GridName="grid1km"<br>
>>>>   XDim=1200<br>
>>>>   YDim=1200<br>
>>>>   UpperLeftPointMtrs=(-3450000.<wbr>000000,1600000.000000)<br>
>>>>   LowerRightMtrs=(-2250000.0000<wbr>00,400000.000000)<br>
>>>>   Projection=GCTP_ALBERS<br>
>>>>   ProjParams=(0,0,21030000,4100<wbr>0000,102000000,23000000,0,0,0,<wbr>0,0,0,0)<br>
>>>>   SphereCode=12<br>
>>>>   GridOrigin=HDFE_GD_UL<br>
>>>>   GROUP=Dimension<br>
>>>>   END_GROUP=Dimension<br>
>>>><br>
>>>> The name suggests a 1KM grid, and as such, by its resolution, it should<br>
>>>> span 1200KM. Indeed, from the corner points:<br>
>>>><br>
>>>> In X:  (-3450000 - -2250000) = a span of 1200000meters; divided by<br>
>>>> xres=1200, yields 1000meter grid cells<br>
>>>> In Y:  (1600000 - 400000) = 1200000meters; div by yres=1200 yields<br>
>>>> 1000meter grid cells<br>
>>>><br>
>>>> So everything is consistent so far. The GCTP_ALBERS suggest the grid was<br>
>>>> constructed based upon an ALBERS projection, using the GCTP library. The<br>
>>>> ProjParams give further details:<br>
>>>><br>
>>>>   ProjParams=(0,0,21030000,4100<wbr>0000,102000000,23000000,0,0,0,<wbr>0,0,0,0)<br>
>>>><br>
>>>> Again, without documentation, this is speculation, but looking at<br>
>>>> pg13/appendexD in:<br>
>>>><br>
>>>><br>
>>>> <a href="https://www.ngs.noaa.gov/PUBS_LIB/GeneralCartographicTransformationPackage_v2_TR_NOS124_CGS9.pdf" rel="noreferrer" target="_blank">https://www.ngs.noaa.gov/PUBS_<wbr>LIB/GeneralCartographicTransfo<wbr>rmationPackage_v2_TR_NOS124_<wbr>CGS9.pdf</a><br>
>>>><br>
>>>> for the Albers projection, the first 2 values appear to be semimajor<br>
>>>> axis of ellipsoid, and eccentricity. Both are zero, implying defaults of<br>
>>>> Clarke's 1866 ellipsoid and spherical earth (which seems to be an oymoron?<br>
>>>> I'm confused about that)<br>
>>>><br>
>>>> The next 6 values might be:<br>
>>>><br>
>>>>   lat of 1st standard parallel<br>
>>>>   lat of 2nd standard parallel<br>
>>>>   lon of central meridian<br>
>>>>   lat of projection's origin<br>
>>>>   false easting (zero in this case)<br>
>>>>   false northing (zero)<br>
>>>><br>
>>>> the remaining zeros don't seem to apply.<br>
>>>><br>
>>>> NCL does not have an public interface to GCTP, but there is an<br>
>>>> undocumented interface to the well-known proj4 library. The trick here is to<br>
>>>> use it to convert from meters-on-the ground units into lat/lon.  The first<br>
>>>> step is to construct arrays representing the grid's position in meters. In<br>
>>>> the X direction, if we did<br>
>>>><br>
>>>>   x = fspan(-3450000., -2250000, 1200)<br>
>>>><br>
>>>> we get an array of 1201, which is wrong. I speculate that for a grid of<br>
>>>> 1200, those values represent the total *spatial extent* of the grid, but<br>
>>>> that the data are based upon cell centers, not corners. Thus something like:<br>
>>>><br>
>>>>   x = fspan(-3450000.+500., -2250000.-500., 1200)<br>
>>>>   y = fspan(1600000.-500., 400000.+500., 1200)<br>
>>>><br>
>>>> A 1km grid cannot obviously be constant spacing in longitude, so<br>
>>>> ultimately, your coordinate arrays are going to be curvilinear:<br>
>>>><br>
>>>>   x := conform_dims((/1200, 1200/), x, 1)<br>
>>>>   y := conform_dims((/1200, 1200/), y, 0)<br>
>>>><br>
>>>> NOTE: THIS ASSUMES THE GRID IS ROW MAJOR, WEST-->EAST, NORTH-->SOUTH<br>
>>>><br>
>>>> The undocumented NCL function "transform_coordinates" takes a proj4<br>
>>>> source-porjection string, dest-proj string, x, y, z 1D arrays of type<br>
>>>> double, and returns the results in the x,y,z arrays. I won't go into the<br>
>>>> details of the proj-strings. We want something like:<br>
>>>><br>
>>>> src = "+proj=aea +lat_0=23.0 +lat_1=21d03m +lat_2=41 +lon_0=102<br>
>>>> +ellps=clrk66"<br>
>>>> dst = "+proj=latlon +ellps=sphere"<br>
>>>> x := flt2dble(ndtooned(x))<br>
>>>> y := flt2dble(ndtooned(y))<br>
>>>> z = x   ; A quick-n-dirty way to get an empty Z array of same<br>
>>>> dimensions/type<br>
>>>> z = 0.  ;<br>
>>>> transform_coordinate(src, dst, x, y, z)<br>
>>>><br>
>>>> ; convert back to curvilinear coord arrays<br>
>>>> lon = onedtond(x, (/1200, 1200/))<br>
>>>> lat = onedtond(y, (/1200, 1200/))<br>
>>>><br>
>>>><br>
>>>> From here, you have to associate the lat/lon coordinate arrays with the<br>
>>>> variables you want to plot, as Adam and other previously described.<br>
>>>><br>
>>>> Attached is a script and example plot of the resultant grid. Please let<br>
>>>> the group know if this appears to be the correct area-of-interest.<br>
>>>><br>
>>>> Rick<br>
>>>><br>
>>>><br>
>>>><br>
>>>> On Fri, Sep 29, 2017 at 11:45 AM, Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>><br>
>>>> wrote:<br>
>>>>><br>
>>>>> Hi,<br>
>>>>><br>
>>>>> I do not know for certain without examining the file, but it looks like<br>
>>>>> those values are map-projection parameters, from which curvilinear grid was<br>
>>>>> generated as equal-area grid cells in lat/lon spac.<br>
>>>>><br>
>>>>> Rick<br>
>>>>><br>
>>>>> On Fri, Sep 29, 2017 at 10:16 AM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>><br>
>>>>> wrote:<br>
>>>>>><br>
>>>>>> Thanks for all the suggestions and help Adam Sir. Now the data is<br>
>>>>>> showing over the Indian region.<br>
>>>>>> but one thing I need to confirm here.<br>
>>>>>><br>
>>>>>> Some values in the form of UpperLeftPointMtrs and LowerRightMtrs are<br>
>>>>>> given in the file ((in the original hdf file) ) description along with X<br>
>>>>>> YDim.<br>
>>>>>><br>
>>>>>> So, does these value interfer while creating the lat lon dimensions<br>
>>>>>> during conversion netcdf from hdf?<br>
>>>>>><br>
>>>>>> Sould these values be consider or we can ignore these values?<br>
>>>>>><br>
>>>>>><br>
>>>>>> group: latlon {<br>
>>>>>>     variables:<br>
>>>>>>       short _HDFEOS_CRS;<br>
>>>>>>         :Projection = "GCTP_ALBERS";<br>
>>>>>>         :UpperLeftPointMtrs = -3450000.0, 1600000.0; // double<br>
>>>>>>         :LowerRightMtrs = -2250000.0, 400000.0; // double<br>
>>>>>>         :ProjParams = 0.0, 0.0, 2.103E7, 4.1E7, 1.02E8, 2.3E7, 0.0,<br>
>>>>>> 0.0, 0.0, 0.0, 0.0, 0.0, 0.0; // double<br>
>>>>>>         :SphereCode = "12";<br>
>>>>>><br>
>>>>>>     group: Data_Fields {<br>
>>>>>>       dimensions:<br>
>>>>>>         YDim = 1200;<br>
>>>>>>         XDim = 1200;<br>
>>>>>>       variables:<br>
>>>>>>         float lat(YDim=1200, XDim=1200);<br>
>>>>>><br>
>>>>>>         float lon(YDim=1200, XDim=1200);<br>
>>>>>><br>
>>>>>><br>
>>>>>><br>
>>>>>><br>
>>>>>><br>
>>>>>><br>
>>>>>> Regards<br>
>>>>>> Kunal Bali<br>
>>>>>><br>
>>>>>><br>
>>>>>><br>
>>>>>><br>
>>>>>> On Thu, Sep 28, 2017 at 9:27 PM, Adam Phillips <<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>><br>
>>>>>> wrote:<br>
>>>>>>><br>
>>>>>>> Hi Kunal,<br>
>>>>>>> As the documentation states, latGlobeFo (lonGlobeFo) generates<br>
>>>>>>> latitudes (longitudes) and associated metadata for a global fixed offset<br>
>>>>>>> grid.<br>
>>>>>>><br>
>>>>>>> <a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/latGlobeFo.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Docum<wbr>ent/Functions/Contributed/<wbr>latGlobeFo.shtml</a><br>
>>>>>>><br>
>>>>>>> <a href="https://www.ncl.ucar.edu/Document/Functions/Contributed/lonGlobeFo.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Docum<wbr>ent/Functions/Contributed/<wbr>lonGlobeFo.shtml</a><br>
>>>>>>> If you need to create coordinates for a regional grid, then you will<br>
>>>>>>> have to do it by hand. For example, if you know your data are equally spaced<br>
>>>>>>> in the both the latitude and longitude direction from 0:40N, 60:100E (for<br>
>>>>>>> example), you can set up your coordinates like this:<br>
>>>>>>><br>
>>>>>>> lat = fspan(0,40,1200)<br>
>>>>>>> lon = fspan(60,100,1200)<br>
>>>>>>> print(lat)   ; check that these are the correct coordinates<br>
>>>>>>> print(lon)   ; check that these are the correct coordinates<br>
>>>>>>> lat@units = "degrees_north"<br>
>>>>>>> lat!0 = "lat"<br>
>>>>>>> lat&lat = lat<br>
>>>>>>> lon@units = "degrees_east"<br>
>>>>>>> lon!0 = "lon"<br>
>>>>>>> lon&lon = lon<br>
>>>>>>> var!0 = "lat"    ; you can name these dimensions whatever you want..<br>
>>>>>>> var!1 = "lon"<br>
>>>>>>> var&lat = lat     ; but make sure you refer to the correct named<br>
>>>>>>> dimensions<br>
>>>>>>> var&lon = lon<br>
>>>>>>><br>
>>>>>>> The above is similar to what latGlobeFo/lonGlobeFo do for global<br>
>>>>>>> grids. Hope that helps. If you have further questions please respond to<br>
>>>>>>> ncl-talk.<br>
>>>>>>> Adam<br>
>>>>>>><br>
>>>>>>> On Tue, Sep 26, 2017 at 10:26 AM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>><br>
>>>>>>> wrote:<br>
>>>>>>>><br>
>>>>>>>> Thanks for all the suggestions Dennis sir.<br>
>>>>>>>> But now I have come to know that the dimensions<br>
>>>>>>>><br>
>>>>>>>> nlat=1200<br>
>>>>>>>> nlon=1200<br>
>>>>>>>><br>
>>>>>>>> are belongs to only of Indian region.<br>
>>>>>>>><br>
>>>>>>>> And the code I am using below is not displaying the data over Indian<br>
>>>>>>>> region.<br>
>>>>>>>><br>
>>>>>>>>  lat = latGlobeFo(nlat, "lat", "latitude", "degrees_north")<br>
>>>>>>>>  lon = lonGlobeFo(nlon, "lon", "longitude", "degrees_east")<br>
>>>>>>>> lat = lat(::-1)<br>
>>>>>>>> lon = (/ lon - 180. /)  ; subtract 180 from all values<br>
>>>>>>>> lon&lon = lon<br>
>>>>>>>><br>
>>>>>>>> So please help me regarding this issue.<br>
>>>>>>>><br>
>>>>>>>> regards<br>
>>>>>>>> Kunal Bali<br>
>>>>>>>><br>
>>>>>>>><br>
>>>>>>>><br>
>>>>>>>><br>
>>>>>>>><br>
>>>>>>>><br>
>>>>>>>> On Tue, Sep 19, 2017 at 10:36 PM, Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<br>
>>>>>>>>><br>
>>>>>>>>> Ditto what Marston posted.<br>
>>>>>>>>><br>
>>>>>>>>> ===<br>
>>>>>>>>> Please look at NCL's FAQ.<br>
>>>>>>>>><br>
>>>>>>>>> [1] Click "Support" at top of page<br>
>>>>>>>>> [2] Click FAQ<br>
>>>>>>>>><br>
>>>>>>>>> I think it would be instructive to look at all of the topics.<br>
>>>>>>>>><br>
>>>>>>>>> For example:  Error messages and other issues<br>
>>>>>>>>><br>
>>>>>>>>> Look at #7<br>
>>>>>>>>><br>
>>>>>>>>><br>
>>>>>>>>><br>
>>>>>>>>> On Tue, Sep 19, 2017 at 3:03 AM, Marston Johnston<br>
>>>>>>>>> <<a href="mailto:shejo284@gmail.com" target="_blank">shejo284@gmail.com</a>> wrote:<br>
>>>>>>>>>><br>
>>>>>>>>>> It tells you in the warning message, exactly what to do to remove<br>
>>>>>>>>>> the warning.<br>
>>>>>>>>>><br>
>>>>>>>>>> If you want to improve you coding skills, it would better to<br>
>>>>>>>>>> understand why the warning is occurring in the first place: dimension 0 and<br>
>>>>>>>>>> 1 of the variable you are passing has no names or the names do not match the<br>
>>>>>>>>>> variable on the lhs.<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> /M<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<br>
>>>>>>>>>><br>
>>>>>>>>>> Marston S. Ward, PhD<br>
>>>>>>>>>><br>
>>>>>>>>>> Department of Earth Sciences<br>
>>>>>>>>>><br>
>>>>>>>>>> University of Gothenburg, Sweden<br>
>>>>>>>>>><br>
>>>>>>>>>> Email: <a href="mailto:marston.johnston@gu.se" target="_blank">marston.johnston@gu.se</a><br>
>>>>>>>>>><br>
>>>>>>>>>> SkypeID: marston.johnston<br>
>>>>>>>>>><br>
>>>>>>>>>> Phone: <a href="tel:%2B46-31-7864901" value="+46317864901" target="_blank">+46-31-7864901</a><br>
>>>>>>>>>><br>
>>>>>>>>>> Only the fruitful thing is true!<br>
>>>>>>>>>><br>
>>>>>>>>>> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> From: ncl-talk <<a href="mailto:ncl-talk-bounces@ucar.edu" target="_blank">ncl-talk-bounces@ucar.edu</a>> on behalf of Kunal Bali<br>
>>>>>>>>>> <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>><br>
>>>>>>>>>> Date: Tuesday, 19 September 2017 at 10:50<br>
>>>>>>>>>> Cc: "<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
>>>>>>>>>> Subject: Re: [ncl-talk] set_dimension_query<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> Thank you for providing this information.<br>
>>>>>>>>>><br>
>>>>>>>>>> It's working. The script given below is now able to convert<br>
>>>>>>>>>> multiple files into .netcdf with the desired file name.<br>
>>>>>>>>>> However, on the same time, it also shows some warnings. But noted<br>
>>>>>>>>>> that, I am successfully creating the multiple files into the netcdf. I just<br>
>>>>>>>>>> need to know that how can I remove these warnings?<br>
>>>>>>>>>><br>
>>>>>>>>>> warning:VarVarWrite: Dimension names for dimension number (0)<br>
>>>>>>>>>> don't match, assigning name of rhs dimension to lhs and overwriting<br>
>>>>>>>>>> coordinate variable, use "(/../)" if this change is not desired<br>
>>>>>>>>>> warning:VarVarWrite: Dimension names for dimension number (1)<br>
>>>>>>>>>> don't match, assigning name of rhs dimension to lhs and overwriting<br>
>>>>>>>>>> coordinate variable, use "(/../)" if this change is not desired<br>
>>>>>>>>>> warning:["Execute.c":8640]:Exe<wbr>cute: Error occurred at or near line<br>
>>>>>>>>>> 22 in file netcdf2.ncl<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> ;-----------------------------<wbr>------------------------------<wbr>-----------<br>
>>>>>>>>>><br>
>>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/gsn_code.ncl"<br>
>>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/gsn_csm.ncl"<br>
>>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/contributed.ncl"<br>
>>>>>>>>>><br>
>>>>>>>>>> begin<br>
>>>>>>>>>> ;---Read data<br>
>>>>>>>>>><br>
>>>>>>>>>>         diri = "/Users/Pushp/Desktop/test/"   ; input directory<br>
>>>>>>>>>>         fili  = systemfunc("cd "+diri+" ; ls<br>
>>>>>>>>>> MAIACTAOT.h00v02*hdf")<br>
>>>>>>>>>>         nfili = dimsizes(fili)<br>
>>>>>>>>>>         print("nfili="+nfili)<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>>         dirnc = "/Users/Pushp/Desktop/test/"   ; output (netCDF)<br>
>>>>>>>>>> directory<br>
>>>>>>>>>><br>
>>>>>>>>>>         do nf=0,nfili-1<br>
>>>>>>>>>>         pthi = diri+fili(nf)<br>
>>>>>>>>>>         f    = addfile(pthi,"r")<br>
>>>>>>>>>><br>
>>>>>>>>>>         var  = short2flt(f->Optical_Depth_055<wbr>(:,:))<br>
>>>>>>>>>><br>
>>>>>>>>>>         nlat = 1200<br>
>>>>>>>>>>         nlon = 1200<br>
>>>>>>>>>>         lat  = latGlobeFo(nlat, "lat", "latitude",<br>
>>>>>>>>>> "degrees_north")<br>
>>>>>>>>>>         lon  = lonGlobeFo(nlon, "lon", "longitude",<br>
>>>>>>>>>> "degrees_east")<br>
>>>>>>>>>>     ;    lat  = lat(::-1)<br>
>>>>>>>>>>         lon  = (/ lon - 180. /)  ; subtract 180 from all values<br>
>>>>>>>>>>         lon&lon = lon           ; update coordinates<br>
>>>>>>>>>><br>
>>>>>>>>>>         var!0 = "lat"    ; you can name these dimensions whatever<br>
>>>>>>>>>> you want..<br>
>>>>>>>>>>         var!1 = "lon"<br>
>>>>>>>>>>         var&lat = lat     ; but make sure you refer to the correct<br>
>>>>>>>>>> named dimensions<br>
>>>>>>>>>>         var&lon = lon<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>>      filroot = str_get_cols(fili(nf), 0, 27)  ; eg:<br>
>>>>>>>>>> "MAIACTAOT.h00v02.20000570505"<br>
>>>>>>>>>>      filnc    = filroot+".nc"<br>
>>>>>>>>>>      pthnc  = dirnc + filnc<br>
>>>>>>>>>><br>
>>>>>>>>>>      system("/bin/rm -f "+pthnc)<br>
>>>>>>>>>>      ncdf = addfile(pthnc,"c")  ; open new netCDF file<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>>      end do    ; end 'nf' loop<br>
>>>>>>>>>><br>
>>>>>>>>>>    ;     system("/bin/rm -f <a href="http://simple2.nc" rel="noreferrer" target="_blank">simple2.nc</a>")<br>
>>>>>>>>>>    ;     ncdf = addfile("<a href="http://simple2.nc" rel="noreferrer" target="_blank">simple2.nc</a>" ,"c")  ; open output netCDF<br>
>>>>>>>>>> file<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> ;=============================<wbr>==============================<wbr>========<br>
>>>>>>>>>>     ; create global attributes of the file (optional)<br>
>>>>>>>>>><br>
>>>>>>>>>> ;=============================<wbr>==============================<wbr>========<br>
>>>>>>>>>>        fAtt               = True            ; assign file<br>
>>>>>>>>>> attributes<br>
>>>>>>>>>>        fAtt@title         = "NCL Simple Approach to netCDF<br>
>>>>>>>>>> Creation"<br>
>>>>>>>>>>        fAtt@source_file   =  "<a href="http://original-file.nc" rel="noreferrer" target="_blank">original-file.nc</a>"<br>
>>>>>>>>>>        fAtt@Conventions   = "None"<br>
>>>>>>>>>>        fAtt@creation_date = systemfunc ("date")<br>
>>>>>>>>>>        fileattdef( ncdf, fAtt )            ; copy file attributes<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> ;=============================<wbr>==============================<wbr>========<br>
>>>>>>>>>>     ; make time an UNLIMITED dimension; recommended  for most<br>
>>>>>>>>>> applications<br>
>>>>>>>>>><br>
>>>>>>>>>> ;=============================<wbr>==============================<wbr>========<br>
>>>>>>>>>>        filedimdef(ncdf,"time",-1,True<wbr>)<br>
>>>>>>>>>><br>
>>>>>>>>>> ;       ncdf->var&lat = lat<br>
>>>>>>>>>> ;       ncdf->var&lon = lon<br>
>>>>>>>>>>        ncdf->var = var<br>
>>>>>>>>>><br>
>>>>>>>>>> end<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> Kunal Bali<br>
>>>>>>>>>><br>
>>>>>>>>>> Research Scholar<br>
>>>>>>>>>><br>
>>>>>>>>>> Radio & Atmospheric Science Division<br>
>>>>>>>>>><br>
>>>>>>>>>> CSIR - National Physical Laboratory<br>
>>>>>>>>>><br>
>>>>>>>>>> New Delhi - 110012<br>
>>>>>>>>>><br>
>>>>>>>>>> India<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> On Mon, Sep 18, 2017 at 9:04 PM, Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
>>>>>>>>>> wrote:<br>
>>>>>>>>>><br>
>>>>>>>>>> Loop over the files ...<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>>  diri = "/media/Local Disk/NPL/MODIS_FPC/"   ; input directory<br>
>>>>>>>>>><br>
>>>>>>>>>>  fili  = systemfunc("cd "+diri +; ls MAIACTAOT.h00v02*hdf")<br>
>>>>>>>>>><br>
>>>>>>>>>>  nfili = dimsizes(fili)<br>
>>>>>>>>>><br>
>>>>>>>>>>  print("nfili="+nfili)<br>
>>>>>>>>>><br>
>>>>>>>>>>  dirnc = "/media/Local Disk/NPL/MODIS_FPC/"   ; output (netCDF)<br>
>>>>>>>>>> directory<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>>  do nf=0,nfili-1<br>
>>>>>>>>>><br>
>>>>>>>>>>       pthi = diri+fili(nf)<br>
>>>>>>>>>><br>
>>>>>>>>>>       f      = addfile(pthi,"r")<br>
>>>>>>>>>><br>
>>>>>>>>>>       ....<br>
>>>>>>>>>><br>
>>>>>>>>>>      filroot = str_get_cols(fili(nf), 0, 27)  ; eg:<br>
>>>>>>>>>> "MAIACTAOT.h00v02.20000570505"<br>
>>>>>>>>>><br>
>>>>>>>>>>      filnc    = filroot+".nc"<br>
>>>>>>>>>><br>
>>>>>>>>>>      pthnc  = dirnc + filnc<br>
>>>>>>>>>><br>
>>>>>>>>>>      system("/bin/rm -f "+pthnc)<br>
>>>>>>>>>>      ncdf = addfile(pthnc,"c")  ; open new netCDF file<br>
>>>>>>>>>><br>
>>>>>>>>>>      ....<br>
>>>>>>>>>><br>
>>>>>>>>>> end do    ; end 'nf' loop<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> =====<br>
>>>>>>>>>><br>
>>>>>>>>>> If there is something you do not understand use some print<br>
>>>>>>>>>> statements.<br>
>>>>>>>>>><br>
>>>>>>>>>> Become familiar with the 'string' category. Look at all the<br>
>>>>>>>>>> functions.<br>
>>>>>>>>>><br>
>>>>>>>>>>    <a href="https://www.ncl.ucar.edu/Document/Functions/string.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Docum<wbr>ent/Functions/string.shtml</a><br>
>>>>>>>>>><br>
>>>>>>>>>> Please read the documentation for the function used above<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/str_get_cols.shtml" rel="noreferrer" target="_blank">https://www.ncl.ucar.edu/Docum<wbr>ent/Functions/Built-in/str_<wbr>get_cols.shtml</a><br>
>>>>>>>>>><br>
>>>>>>>>>> Good Luck<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> On Mon, Sep 18, 2017 at 6:11 AM, Kunal Bali<br>
>>>>>>>>>> <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<br>
>>>>>>>>>><br>
>>>>>>>>>> Thnks for the reply.<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> that will change the name of only one file. isn't it. but what<br>
>>>>>>>>>> about the other files<br>
>>>>>>>>>><br>
>>>>>>>>>> I have 365 files in one directory and I need to change all the<br>
>>>>>>>>>> file name one by one.<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> Kunal Bali<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> On Mon, Sep 18, 2017 at 5:14 PM, Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>><br>
>>>>>>>>>> wrote:<br>
>>>>>>>>>><br>
>>>>>>>>>> Really?<br>
>>>>>>>>>><br>
>>>>>>>>>> Change:<br>
>>>>>>>>>>       system("/bin/rm -f <a href="http://simple.nc" rel="noreferrer" target="_blank">simple.nc</a>")<br>
>>>>>>>>>>       ncdf = addfile("<a href="http://simple.nc" rel="noreferrer" target="_blank">simple.nc</a>" ,"c")<br>
>>>>>>>>>><br>
>>>>>>>>>> To:<br>
>>>>>>>>>>       system("/bin/rm -f <a href="http://MAIACTAOT.h00v02.20000570505.nc" rel="noreferrer" target="_blank">MAIACTAOT.h00v02.20000570505.n<wbr>c</a>")<br>
>>>>>>>>>>       ncdf = addfile("<a href="http://MAIACTAOT.h00v02.20000570505.nc" rel="noreferrer" target="_blank">MAIACTAOT.h00v02.2000<wbr>0570505.nc</a>" ,"c")<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> On Mon, Sep 18, 2017 at 12:37 AM, Kunal Bali<br>
>>>>>>>>>> <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<br>
>>>>>>>>>><br>
>>>>>>>>>> Thanks, it worked.<br>
>>>>>>>>>><br>
>>>>>>>>>> one more question is that.<br>
>>>>>>>>>><br>
>>>>>>>>>> I used the code (given below). So it gives <a href="http://simple.nc" rel="noreferrer" target="_blank">simple.nc</a> name as an<br>
>>>>>>>>>> output file.<br>
>>>>>>>>>><br>
>>>>>>>>>> ncdf = addfile("<a href="http://simple.nc" rel="noreferrer" target="_blank">simple.nc</a>" ,"c")<br>
>>>>>>>>>><br>
>>>>>>>>>> If I want to keep the original file name with the output file name<br>
>>>>>>>>>> then what should I do?<br>
>>>>>>>>>><br>
>>>>>>>>>> I mean the original file name is MAIACTAOT.h00v02.20003660700.h<wbr>df<br>
>>>>>>>>>><br>
>>>>>>>>>> Now I want to create my netcdf file name  as<br>
>>>>>>>>>> <a href="http://MAIACTAOT.h00v02.20003660700.nc" rel="noreferrer" target="_blank">MAIACTAOT.h00v02.20003660700.n<wbr>c</a> NOT <a href="http://simple.nc" rel="noreferrer" target="_blank">simple.nc</a>.<br>
>>>>>>>>>><br>
>>>>>>>>>> I can not simply write <a href="http://simple.nc" rel="noreferrer" target="_blank">simple.nc</a> file name because I have many<br>
>>>>>>>>>> files.<br>
>>>>>>>>>><br>
>>>>>>>>>> e.g I have one directory having 365 .hdf file. So I need to<br>
>>>>>>>>>> convert all the hdf file to netcdf file at once with the same name of<br>
>>>>>>>>>> original file names.<br>
>>>>>>>>>><br>
>>>>>>>>>> I hope you have got my query.<br>
>>>>>>>>>><br>
>>>>>>>>>> please let me know that too.<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> Thank You<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> Kunal Bali<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> On Thu, Sep 14, 2017 at 9:09 PM, Adam Phillips <<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>><br>
>>>>>>>>>> wrote:<br>
>>>>>>>>>><br>
>>>>>>>>>> Hi Kunal,<br>
>>>>>>>>>><br>
>>>>>>>>>> It looks to me like your latitudes are flipped. You are setting up<br>
>>>>>>>>>> and assigning your coordinate lat/lon variables in this coding:<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> nlat = 1200<br>
>>>>>>>>>> nlon = 1200<br>
>>>>>>>>>> lat = latGlobeFo(nlat, "lat", "latitude", "degrees_north")<br>
>>>>>>>>>><br>
>>>>>>>>>> lon = lonGlobeFo(nlon, "lon", "longitude", "degrees_east")<br>
>>>>>>>>>> lat = lat(::-1)<br>
>>>>>>>>>> lon = (/ lon - 180. /)  ; subtract 180 from all values<br>
>>>>>>>>>><br>
>>>>>>>>>> lon&lon = lon           ; update coordinates<br>
>>>>>>>>>><br>
>>>>>>>>>> var!0 = "lat"    ; you can name these dimensions whatever you<br>
>>>>>>>>>> want..<br>
>>>>>>>>>> var!1 = "lon"<br>
>>>>>>>>>> var&lat = lat     ; but make sure you refer to the correct named<br>
>>>>>>>>>> dimensions<br>
>>>>>>>>>> var&lon = lon<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> Assuming I'm right on your latitudes being the issue, one of the<br>
>>>>>>>>>> following two modifications should work:<br>
>>>>>>>>>><br>
>>>>>>>>>> Delete this line:<br>
>>>>>>>>>><br>
>>>>>>>>>> lat = lat(::-1)<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> If the above doesn't fix it, try this:<br>
>>>>>>>>>><br>
>>>>>>>>>> Change this:<br>
>>>>>>>>>><br>
>>>>>>>>>> lat = lat(::-1)<br>
>>>>>>>>>><br>
>>>>>>>>>> to this:<br>
>>>>>>>>>><br>
>>>>>>>>>> lat = lat(::-1)<br>
>>>>>>>>>><br>
>>>>>>>>>> lat&lat = lat<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> If you continue to have issues please respond to ncl-talk.<br>
>>>>>>>>>><br>
>>>>>>>>>> Adam<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> On Wed, Sep 13, 2017 at 12:35 PM, Kunal Bali<br>
>>>>>>>>>> <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> For direct output as a netcdf format, I used the script given<br>
>>>>>>>>>> below. It produced the netcdf file easily. But the netcdf file and original<br>
>>>>>>>>>> hdf file both showing the different results. The data pattern is shifted. I<br>
>>>>>>>>>> mean it may be related to the lat lon position. please see the attached<br>
>>>>>>>>>> file, you will understand.  And please let me know which one is correct.<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> ;-----------------------------<wbr>------------------------------<wbr>-----------<br>
>>>>>>>>>><br>
>>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/gsn_code.ncl"<br>
>>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/gsn_csm.ncl"<br>
>>>>>>>>>> load "$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/contributed.ncl"<br>
>>>>>>>>>><br>
>>>>>>>>>> begin<br>
>>>>>>>>>> ;---Read data<br>
>>>>>>>>>>         a = addfile("/media/Local<br>
>>>>>>>>>> Disk/NPL/MODIS_FPC/MAIACTAOT.h<wbr>00v02.20000570505.hdf","r")<br>
>>>>>>>>>>         var  = short2flt(a->Optical_Depth_055<wbr>(:,:))<br>
>>>>>>>>>>             ;          short2flt(a[:]->noxfire)<br>
>>>>>>>>>>         nlat = 1200<br>
>>>>>>>>>>         nlon = 1200<br>
>>>>>>>>>>         lat = latGlobeFo(nlat, "lat", "latitude", "degrees_north")<br>
>>>>>>>>>>         lon = lonGlobeFo(nlon, "lon", "longitude", "degrees_east")<br>
>>>>>>>>>>         lat = lat(::-1)<br>
>>>>>>>>>>         lon = (/ lon - 180. /)  ; subtract 180 from all values<br>
>>>>>>>>>>         lon&lon = lon           ; update coordinates<br>
>>>>>>>>>><br>
>>>>>>>>>>         var!0 = "lat"    ; you can name these dimensions whatever<br>
>>>>>>>>>> you want..<br>
>>>>>>>>>>         var!1 = "lon"<br>
>>>>>>>>>>         var&lat = lat     ; but make sure you refer to the correct<br>
>>>>>>>>>> named dimensions<br>
>>>>>>>>>>         var&lon = lon<br>
>>>>>>>>>><br>
>>>>>>>>>>                  system("/bin/rm -f <a href="http://simple.nc" rel="noreferrer" target="_blank">simple.nc</a>")<br>
>>>>>>>>>>              ncdf = addfile("<a href="http://simple.nc" rel="noreferrer" target="_blank">simple.nc</a>" ,"c")  ; open output<br>
>>>>>>>>>> netCDF file<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> ;=============================<wbr>==============================<wbr>========<br>
>>>>>>>>>>     ; create global attributes of the file (optional)<br>
>>>>>>>>>><br>
>>>>>>>>>> ;=============================<wbr>==============================<wbr>========<br>
>>>>>>>>>>        fAtt               = True            ; assign file<br>
>>>>>>>>>> attributes<br>
>>>>>>>>>>        fAtt@title         = "NCL Simple Approach to netCDF<br>
>>>>>>>>>> Creation"<br>
>>>>>>>>>>        fAtt@source_file   =  "<a href="http://original-file.nc" rel="noreferrer" target="_blank">original-file.nc</a>"<br>
>>>>>>>>>>        fAtt@Conventions   = "None"<br>
>>>>>>>>>>        fAtt@creation_date = systemfunc ("date")<br>
>>>>>>>>>>        fileattdef( ncdf, fAtt )            ; copy file attributes<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> ;=============================<wbr>==============================<wbr>========<br>
>>>>>>>>>>     ; make time an UNLIMITED dimension; recommended  for most<br>
>>>>>>>>>> applications<br>
>>>>>>>>>><br>
>>>>>>>>>> ;=============================<wbr>==============================<wbr>========<br>
>>>>>>>>>>        filedimdef(ncdf,"time",-1,True<wbr>)<br>
>>>>>>>>>>        ncdf->var = var<br>
>>>>>>>>>><br>
>>>>>>>>>> end<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> Kunal Bali<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> On Wed, Sep 13, 2017 at 11:27 PM, Kunal Bali<br>
>>>>>>>>>> <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<br>
>>>>>>>>>><br>
>>>>>>>>>> No problem, I sorted out.<br>
>>>>>>>>>><br>
>>>>>>>>>> Just changed short to float.<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> Kunal Bali<br>
>>>>>>>>>><br>
>>>>>>>>>> Research Scholar<br>
>>>>>>>>>><br>
>>>>>>>>>> Radio & Atmospheric Science Division<br>
>>>>>>>>>><br>
>>>>>>>>>> CSIR - National Physical Laboratory<br>
>>>>>>>>>><br>
>>>>>>>>>> New Delhi - 110012<br>
>>>>>>>>>><br>
>>>>>>>>>> India<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> On Wed, Sep 13, 2017 at 11:24 PM, Kunal Bali<br>
>>>>>>>>>> <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<br>
>>>>>>>>>><br>
>>>>>>>>>> Thanks for providing this information.<br>
>>>>>>>>>><br>
>>>>>>>>>> Also, I would like to mention that when I plot this data. The<br>
>>>>>>>>>> values are not in the domain. I mean the real values lie in-between 0 to 1<br>
>>>>>>>>>> but here it is reaching to 400.<br>
>>>>>>>>>><br>
>>>>>>>>>> So, how to correct the values?<br>
>>>>>>>>>><br>
>>>>>>>>>> Description of the variable is<br>
>>>>>>>>>><br>
>>>>>>>>>> short Optical_Depth_055(YDim=1200, XDim=1200);<br>
>>>>>>>>>>   :long_name = "AOT at 0.55 micron";<br>
>>>>>>>>>>   :scale_factor = 0.001; // double<br>
>>>>>>>>>>   :add_offset = 0.0; // double<br>
>>>>>>>>>>   :unit = "None";<br>
>>>>>>>>>>   :_FillValue = -28672S; // short<br>
>>>>>>>>>>   :valid_range = -100S, 5000S; // short<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> Kunal Bali<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> On Wed, Sep 13, 2017 at 10:28 PM, Adam Phillips<br>
>>>>>>>>>> <<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>> wrote:<br>
>>>>>>>>>><br>
>>>>>>>>>> Hi Kunal,<br>
>>>>>>>>>><br>
>>>>>>>>>> Yes, clicking on the output netCDF link from the Applications<br>
>>>>>>>>>> page:<br>
>>>>>>>>>><br>
>>>>>>>>>> <a href="http://www.ncl.ucar.edu/Applications/o-netcdf.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/o-netcdf.shtml</a><br>
>>>>>>>>>><br>
>>>>>>>>>> I would recommend following the inefficient method #1. Unless you<br>
>>>>>>>>>> are writing a file with many large variables, the inefficient method works<br>
>>>>>>>>>> just fine.<br>
>>>>>>>>>><br>
>>>>>>>>>> <a href="http://www.ncl.ucar.edu/Applications/method_1.shtml" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/method_1.shtml</a><br>
>>>>>>>>>><br>
>>>>>>>>>> Good luck,<br>
>>>>>>>>>><br>
>>>>>>>>>> Adam<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> On Wed, Sep 13, 2017 at 10:25 AM, Kunal Bali<br>
>>>>>>>>>> <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<br>
>>>>>>>>>><br>
>>>>>>>>>> Thank you so much, it worked<br>
>>>>>>>>>><br>
>>>>>>>>>> I want to know one more thing.<br>
>>>>>>>>>> After rearranging the dimensions can we now convert (or write)<br>
>>>>>>>>>> this arranged file into netcdf format?<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> Kunal Bali<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> On Wed, Sep 13, 2017 at 8:55 PM, Adam Phillips <<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>><br>
>>>>>>>>>> wrote:<br>
>>>>>>>>>><br>
>>>>>>>>>> Hi Kunal,<br>
>>>>>>>>>><br>
>>>>>>>>>> I think you just need to rearrange the order of your lines and<br>
>>>>>>>>>> tweak a couple of lines. As the error message states, lat is not defined in<br>
>>>>>>>>>> your 3rd line and you are referring to it as if it is.Try this:<br>
>>>>>>>>>><br>
>>>>>>>>>> a = addfile("/media/Local<br>
>>>>>>>>>> Disk/NPL/MODIS_FPC/MAIACTAOT.h<wbr>00v02.20000570505.hdf","r")<br>
>>>>>>>>>> var  = a->Optical_Depth_055(:,:)<br>
>>>>>>>>>> nlat = 1200<br>
>>>>>>>>>> nlon = 1200<br>
>>>>>>>>>> lat = latGlobeFo(nlat, "lat", "latitude", "degrees_north")<br>
>>>>>>>>>><br>
>>>>>>>>>> lon = lonGlobeFo(nlon, "lon", "longitude", "degrees_east")<br>
>>>>>>>>>> lat = lat(::-1)<br>
>>>>>>>>>> lon = (/ lon - 180. /)  ; subtract 180 from all values<br>
>>>>>>>>>> lon&lon = lon           ; update coordinates<br>
>>>>>>>>>><br>
>>>>>>>>>> var!0 = "lat"    ; you can name these dimensions whatever you<br>
>>>>>>>>>> want..<br>
>>>>>>>>>> var!1 = "lon"<br>
>>>>>>>>>><br>
>>>>>>>>>> var&lat = lat     ; but make sure you refer to the correct named<br>
>>>>>>>>>> dimensions<br>
>>>>>>>>>><br>
>>>>>>>>>> var&lat = lon<br>
>>>>>>>>>><br>
>>>>>>>>>> ;var&XDim_grid1km = lat<br>
>>>>>>>>>> ;var&YDim_grid1km = lon<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> If you have any further questions please respond to the ncl-talk<br>
>>>>>>>>>> email list.<br>
>>>>>>>>>><br>
>>>>>>>>>> Adam<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> On Wed, Sep 13, 2017 at 8:38 AM, Kunal Bali<br>
>>>>>>>>>> <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<br>
>>>>>>>>>><br>
>>>>>>>>>> Dear NCL<br>
>>>>>>>>>><br>
>>>>>>>>>> I have a file with the variable summary. The dimensions are in 2D.<br>
>>>>>>>>>><br>
>>>>>>>>>> ncl 2>  printVarSummary(var)<br>
>>>>>>>>>><br>
>>>>>>>>>> Variable: var<br>
>>>>>>>>>> Type: short<br>
>>>>>>>>>> Total Size: 2880000 bytes<br>
>>>>>>>>>>             1440000 values<br>
>>>>>>>>>> Number of Dimensions: 2<br>
>>>>>>>>>> Dimensions and sizes:    [YDim_grid1km | 1200] x [XDim_grid1km |<br>
>>>>>>>>>> 1200]<br>
>>>>>>>>>> Coordinates:<br>
>>>>>>>>>> Number Of Attributes: 7<br>
>>>>>>>>>>   long_name :    AOT at 0.55 micron<br>
>>>>>>>>>>   scale_factor :    0.001<br>
>>>>>>>>>>   add_offset :       0<br>
>>>>>>>>>>   unit :    None<br>
>>>>>>>>>>   _FillValue :    -28672<br>
>>>>>>>>>>   valid_range :    ( -100, 5000 )<br>
>>>>>>>>>>   hdf_name :    Optical_Depth_055<br>
>>>>>>>>>><br>
>>>>>>>>>> I am trying to read the dimensions of this file as<br>
>>>>>>>>>><br>
>>>>>>>>>> begin<br>
>>>>>>>>>> ;---Read data<br>
>>>>>>>>>>          a = addfile("/media/Local<br>
>>>>>>>>>> Disk/NPL/MODIS_FPC/MAIACTAOT.h<wbr>00v02.20000570505.hdf","r")<br>
>>>>>>>>>><br>
>>>>>>>>>>             var  = a->Optical_Depth_055(:,:)<br>
>>>>>>>>>><br>
>>>>>>>>>>          var&XDim_grid1km = lat<br>
>>>>>>>>>>          var&YDim_grid1km = lon<br>
>>>>>>>>>><br>
>>>>>>>>>>         nlat = 1200<br>
>>>>>>>>>>          nlon = 1200<br>
>>>>>>>>>>             lat = latGlobeFo(nlat, "lat", "latitude",<br>
>>>>>>>>>> "degrees_north")<br>
>>>>>>>>>>         lon = lonGlobeFo(nlon, "lon", "longitude", "degrees_east")<br>
>>>>>>>>>>             lat = lat(::-1)<br>
>>>>>>>>>>             lon = (/ lon - 180. /)  ; subtract 180 from all values<br>
>>>>>>>>>>             lon&lon = lon           ; update coordinates<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>>         var!0 = "lat"<br>
>>>>>>>>>>         var!1 = "lon"<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> But the error appeared as<br>
>>>>>>>>>> fatal:Variable (lat) is undefined<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> So, could anyone please let me know that how to read dimension of<br>
>>>>>>>>>> this file.<br>
>>>>>>>>>><br>
>>>>>>>>>> Thank You<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> Regards<br>
>>>>>>>>>><br>
>>>>>>>>>> Kunal Bali<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> ______________________________<wbr>_________________<br>
>>>>>>>>>> ncl-talk mailing list<br>
>>>>>>>>>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>>>>>>>>>> List instructions, subscriber options, unsubscribe:<br>
>>>>>>>>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> --<br>
>>>>>>>>>><br>
>>>>>>>>>> Adam Phillips<br>
>>>>>>>>>><br>
>>>>>>>>>> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR<br>
>>>>>>>>>><br>
>>>>>>>>>> <a href="http://www.cgd.ucar.edu/staff/asphilli/" rel="noreferrer" target="_blank">www.cgd.ucar.edu/staff/asphill<wbr>i/</a>   303-497-1726<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> ______________________________<wbr>_________________<br>
>>>>>>>>>> ncl-talk mailing list<br>
>>>>>>>>>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>>>>>>>>>> List instructions, subscriber options, unsubscribe:<br>
>>>>>>>>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> --<br>
>>>>>>>>>><br>
>>>>>>>>>> Adam Phillips<br>
>>>>>>>>>><br>
>>>>>>>>>> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR<br>
>>>>>>>>>><br>
>>>>>>>>>> <a href="http://www.cgd.ucar.edu/staff/asphilli/" rel="noreferrer" target="_blank">www.cgd.ucar.edu/staff/asphill<wbr>i/</a>   303-497-1726<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> ______________________________<wbr>_________________<br>
>>>>>>>>>> ncl-talk mailing list<br>
>>>>>>>>>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>>>>>>>>>> List instructions, subscriber options, unsubscribe:<br>
>>>>>>>>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> --<br>
>>>>>>>>>><br>
>>>>>>>>>> Adam Phillips<br>
>>>>>>>>>><br>
>>>>>>>>>> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR<br>
>>>>>>>>>><br>
>>>>>>>>>> <a href="http://www.cgd.ucar.edu/staff/asphilli/" rel="noreferrer" target="_blank">www.cgd.ucar.edu/staff/asphill<wbr>i/</a>   303-497-1726<br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> ______________________________<wbr>_________________<br>
>>>>>>>>>> ncl-talk mailing list<br>
>>>>>>>>>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>>>>>>>>>> List instructions, subscriber options, unsubscribe:<br>
>>>>>>>>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> ______________________________<wbr>_________________<br>
>>>>>>>>>> ncl-talk mailing list<br>
>>>>>>>>>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>>>>>>>>>> List instructions, subscriber options, unsubscribe:<br>
>>>>>>>>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> ______________________________<wbr>_________________ ncl-talk mailing<br>
>>>>>>>>>> list <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a> List instructions, subscriber options, unsubscribe:<br>
>>>>>>>>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
>>>>>>>>>><br>
>>>>>>>>>><br>
>>>>>>>>>> ______________________________<wbr>_________________<br>
>>>>>>>>>> ncl-talk mailing list<br>
>>>>>>>>>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>>>>>>>>>> List instructions, subscriber options, unsubscribe:<br>
>>>>>>>>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
>>>>>>>>>><br>
>>>>>>>>><br>
>>>>>>>><br>
>>>>>>>><br>
>>>>>>>> ______________________________<wbr>_________________<br>
>>>>>>>> ncl-talk mailing list<br>
>>>>>>>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>>>>>>>> List instructions, subscriber options, unsubscribe:<br>
>>>>>>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
>>>>>>>><br>
>>>>>>><br>
>>>>>>><br>
>>>>>>><br>
>>>>>>> --<br>
>>>>>>> Adam Phillips<br>
>>>>>>> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR<br>
>>>>>>> <a href="http://www.cgd.ucar.edu/staff/asphilli/" rel="noreferrer" target="_blank">www.cgd.ucar.edu/staff/asphill<wbr>i/</a>   303-497-1726<br>
>>>>>>><br>
>>>>>><br>
>>>>>><br>
>>>>>> ______________________________<wbr>_________________<br>
>>>>>> ncl-talk mailing list<br>
>>>>>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>>>>>> List instructions, subscriber options, unsubscribe:<br>
>>>>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
>>>>>><br>
>>>>><br>
>>>><br>
>>><br>
>>><br>
>>> ______________________________<wbr>_________________<br>
>>> ncl-talk mailing list<br>
>>> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
>>> List instructions, subscriber options, unsubscribe:<br>
>>> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
>>><br>
>><br>
><br>
><br>
> ______________________________<wbr>_________________<br>
> ncl-talk mailing list<br>
> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
> List instructions, subscriber options, unsubscribe:<br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
><br>
______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/ncl-talk</a><br>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div></div>