<div dir="ltr"><div><div><div>Ditto what Marston posted.<br><br>===<br></div>Please look at NCL's FAQ.<br><br></div>[1] Click "Support" at top of page<br></div>[2] Click FAQ<br><div><div><br></div><div>I think it would be instructive to look at all of the topics.<br><br></div><div>For example: <b>Error messages and other issues<br><br></b></div><div>Look at #7<br><br><br></div><div><b></b></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Sep 19, 2017 at 3:03 AM, Marston Johnston <span dir="ltr"><<a href="mailto:shejo284@gmail.com" target="_blank">shejo284@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div bgcolor="white" link="blue" vlink="purple" lang="EN-GB"><div class="m_-1899063326404824290WordSection1"><p class="MsoNormal"><span>It tells you in the warning message, exactly what to do to remove the warning.<u></u><u></u></span></p><p class="MsoNormal"><span>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.<u></u><u></u></span></p><p class="MsoNormal"><span><u></u> <u></u></span></p><p class="MsoNormal"><span>/M<u></u><u></u></span></p><p class="MsoNormal"><span><u></u> <u></u></span></p><div><p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070c0">~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070c0">Marston S. Ward, PhD<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070c0">Department of Earth Sciences<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070c0">University of Gothenburg, Sweden<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070c0">Email: <a href="mailto:marston.johnston@gu.se" target="_blank"><span style="color:#0070c0">marston.johnston@gu.se</span></a><u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070c0">SkypeID: marston.johnston <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070c0">Phone: <a href="tel:+46%2031%20786%2049%2001" value="+46317864901" target="_blank">+46-31-7864901</a> <u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070c0">Only the fruitful thing is true!<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman",serif;color:#0070c0">~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<u></u><u></u></span></p></div><p class="MsoNormal"><span><u></u> <u></u></span></p><p class="MsoNormal"><span><u></u> <u></u></span></p><div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0cm 0cm 0cm"><p class="MsoNormal"><b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black">ncl-talk <<a href="mailto:ncl-talk-bounces@ucar.edu" target="_blank">ncl-talk-bounces@ucar.edu</a>> on behalf of Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>><br><b>Date: </b>Tuesday, 19 September 2017 at 10:50<br><b>Cc: </b>"<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><b>Subject: </b>Re: [ncl-talk] set_dimension_query<u></u><u></u></span></p></div><div><div class="h5"><div><p class="MsoNormal"><u></u> <u></u></p></div><div><div><p class="MsoNormal" style="margin-bottom:12.0pt">Thank you for providing this information. <u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt">It's working. The script given below is now able to convert multiple files into .netcdf with the desired file name. <br>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? <br><br>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<br>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<br>warning:["Execute.c":8640]:<wbr>Execute: Error occurred at or near line 22 in file netcdf2.ncl<br><br> <br><br>;-----------------------------<wbr>------------------------------<wbr>-----------<br><br>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_csm.ncl"<br>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/contributed.<wbr>ncl"<br><br>begin<br>;---Read data<br><br> diri = "/Users/Pushp/Desktop/test/" ; input directory<br> fili = systemfunc("cd "+diri+" ; ls MAIACTAOT.h00v02*hdf")<br> nfili = dimsizes(fili)<br> print("nfili="+nfili)<br><br> <br> dirnc = "/Users/Pushp/Desktop/test/" ; output (netCDF) 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_<wbr>055(:,:))<br> <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 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> <br> filroot = str_get_cols(fili(nf), 0, 27) ; eg: "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" target="_blank">simple2.nc</a>") <br> ; ncdf = addfile("<a href="http://simple2.nc" target="_blank">simple2.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" 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 applications<br> ;=============================<wbr>==============================<wbr>========<br> filedimdef(ncdf,"time",-1,<wbr>True)<br><br>; ncdf->var&lat = lat<br>; ncdf->var&lon = lon<br> ncdf->var = var<br><br>end<u></u><u></u></p></div><div><p class="MsoNormal"><br clear="all"><u></u><u></u></p><div><div><div><div><div><div><p class="MsoNormal">Kunal Bali<u></u><u></u></p></div><div><p class="MsoNormal">Research Scholar <u></u><u></u></p></div><div><p class="MsoNormal">Radio & Atmospheric Science Division <u></u><u></u></p></div><div><p class="MsoNormal">CSIR - National Physical Laboratory<u></u><u></u></p></div><div><p class="MsoNormal">New Delhi - 110012<u></u><u></u></p></div><div><p class="MsoNormal">India<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p style="margin:0cm;margin-bottom:.0001pt"><span style="font-size:9.0pt;font-family:"Tahoma",sans-serif"><u></u> <u></u></span></p></div></div></div></div></div></div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">On Mon, Sep 18, 2017 at 9:04 PM, Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<u></u><u></u></p><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"><div><div><div><div><div><div><div><div><p class="MsoNormal">Loop over the files ...<u></u><u></u></p></div><div><p class="MsoNormal"><br> diri = "/media/Local Disk/NPL/MODIS_FPC/" ; input directory<u></u><u></u></p></div><p class="MsoNormal"> fili = systemfunc("cd "+diri +; ls <span style="font-size:9.5pt">MAIACTAOT.h00v02*hdf")</span><u></u><u></u></p></div><p class="MsoNormal"><span style="font-size:9.5pt"> nfili = dimsizes(fili)</span><u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:9.5pt"> print("nfili="+nfili)</span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt"> dirnc = "/media/Local Disk/NPL/MODIS_FPC/" ; output (netCDF) directory</span><u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><p class="MsoNormal"><span style="font-size:9.5pt"> do nf=0,nfili-1</span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt"> pthi = diri+fili(nf)</span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt"> f = addfile(pthi,"r")<br><br> ....<br><br> </span>filroot = <strong><span style="font-family:"Calibri",sans-serif">str_get_cols</span></strong>(fili(nf), 0, 27) ; eg: "MAIACTAOT.h00v02.20000570505"<u></u><u></u></p></div><p class="MsoNormal"><span style="font-size:9.5pt"> filnc = filroot+".nc"</span><u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:9.5pt"> pthnc = dirnc + filnc<br></span><br><span style="font-size:9.5pt"> system("/bin/rm -f "+pthnc) <br> ncdf = addfile(pthnc,"c") ; open new netCDF file<br><br> ....<br><br></span><u></u><u></u></p></div><p class="MsoNormal"><span style="font-size:9.5pt">end do ; end 'nf' loop</span><u></u><u></u></p><div><div><p class="MsoNormal"><span style="font-size:9.5pt"><br>=====</span><u></u><u></u></p></div><div><p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:9.5pt">If there is something you do not understand use some print statements.</span><u></u><u></u></p></div><div><p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:9.5pt">Become familiar with the 'string' category. Look at all the functions.<br><br> <a href="https://www.ncl.ucar.edu/Document/Functions/string.shtml" target="_blank">https://www.ncl.ucar.edu/<wbr>Document/Functions/string.<wbr>shtml</a></span><u></u><u></u></p></div><div><p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:9.5pt">Please read the documentation for the function used above<br><br> <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/str_get_cols.shtml" target="_blank">https://www.ncl.ucar.edu/<wbr>Document/Functions/Built-in/<wbr>str_get_cols.shtml</a></span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">Good Luck</span><u></u><u></u></p></div><div><p class="MsoNormal" style="margin-bottom:12.0pt"><u></u> <u></u></p></div></div></div><div><div><div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">On Mon, Sep 18, 2017 at 6:11 AM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<u></u><u></u></p><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"><div><div><p class="MsoNormal">Thnks for the reply. <u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal" style="margin-bottom:12.0pt">that will change the name of only one file. isn't it. but what about the other files <u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt">I have 365 files in one directory and I need to change all the file name one by one.<u></u><u></u></p><div><p class="MsoNormal"><span style="color:#888888"><br clear="all"></span><span class="m_-1899063326404824290m1966040847522573004hoenzb"><span style="color:#888888"><u></u><u></u></span></span></p><div><div><div><div><div><div><p class="MsoNormal"><span style="color:#888888">Kunal Bali</span><u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt"><span style="color:#888888"><u></u> <u></u></span></p><div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p></div><div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p></div><div><p style="margin:0cm;margin-bottom:.0001pt"><span style="font-size:9.0pt;font-family:"Tahoma",sans-serif;color:#888888"><u></u> <u></u></span></p></div></div></div></div></div></div><div><div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">On Mon, Sep 18, 2017 at 5:14 PM, Dennis Shea <<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>> wrote:<u></u><u></u></p><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"><div><p class="MsoNormal">Really? <br><br>Change:<br> system("/bin/rm -f <a href="http://simple.nc" target="_blank">simple.nc</a>") <br> ncdf = addfile("<a href="http://simple.nc" target="_blank">simple.nc</a>" ,"c") <br><br>To:<br> system("/bin/rm -f <a href="http://simple.nc" target="_blank">MAIACTAOT.h00v02.20000570505.<wbr>nc</a>") <br> ncdf = addfile("MAIACTAOT.h<a href="http://00v02.20000570505.nc" target="_blank">00v02.<wbr>20000570505.nc</a>" ,"c") <u></u><u></u></p><div><div><p class="MsoNormal"><u></u> <u></u></p></div></div></div><div><div><div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">On Mon, Sep 18, 2017 at 12:37 AM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<u></u><u></u></p><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"><div><div><div><div><div><div><div><div><div><div><p class="MsoNormal" style="margin-bottom:12.0pt">Thanks, it worked.<u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt">one more question is that.<u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt">I used the code (given below). So it gives <a href="http://simple.nc" target="_blank">simple.nc</a> name as an output file.<br><br>ncdf = addfile("<a href="http://simple.nc" target="_blank">simple.nc</a>" ,"c")<u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt">If I want to keep the original file name with the output file name then what should I do?<u></u><u></u></p></div><p class="MsoNormal">I mean the original file name is MAIACTAOT.h00v02.20003660700.<wbr>hdf<u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt">Now I want to create my netcdf file name as <a href="http://MAIACTAOT.h00v02.20003660700.nc" target="_blank">MAIACTAOT.h00v02.20003660700.<wbr>nc</a> NOT <a href="http://simple.nc" target="_blank">simple.nc</a>. <u></u><u></u></p></div><p class="MsoNormal">I can not simply write <a href="http://simple.nc" target="_blank">simple.nc</a> file name because I have many files. <u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt">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.<u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt">I hope you have got my query.<u></u><u></u></p></div><p class="MsoNormal">please let me know that too.<u></u><u></u></p><div><div><div><div><div><div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">Thank You<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div></div></div></div></div></div></div><div><p class="MsoNormal"><br clear="all"><u></u><u></u></p><div><div><div><div><div><div><p class="MsoNormal">Kunal Bali<u></u><u></u></p></div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p style="margin:0cm;margin-bottom:.0001pt"><span style="font-size:9.0pt;font-family:"Tahoma",sans-serif"><u></u> <u></u></span></p></div></div></div></div></div></div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">On Thu, Sep 14, 2017 at 9:09 PM, Adam Phillips <<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>> wrote:<u></u><u></u></p><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"><div><p class="MsoNormal">Hi Kunal,<u></u><u></u></p><div><p class="MsoNormal">It looks to me like your latitudes are flipped. You are setting up and assigning your coordinate lat/lon variables in this coding:<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">nlat = 1200<br>nlon = 1200<br>lat = latGlobeFo(nlat, "lat", "latitude", "degrees_north")</span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">lon = lonGlobeFo(nlon, "lon", "longitude", "degrees_east") <br>lat = lat(::-1)<br>lon = (/ lon - 180. /) ; subtract 180 from all values </span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">lon&lon = lon ; update coordinates<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 </span><u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">Assuming I'm right on your latitudes being the issue, one of the following two modifications should work:<u></u><u></u></p></div><div><p class="MsoNormal">Delete this line:<u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">lat = lat(::-1)</span><u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">If the above doesn't fix it, try this:</span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">Change this:</span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">lat = lat(::-1)</span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">to this:</span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">lat = lat(::-1)</span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">lat&lat = lat</span><u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">If you continue to have issues please respond to ncl-talk.</span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt;color:#888888">Adam</span><span style="color:#888888"><u></u><u></u></span></p></div></div><div><div><div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">On Wed, Sep 13, 2017 at 12:35 PM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<u></u><u></u></p><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"><div><div><p class="MsoNormal"><u></u> <u></u></p></div><p class="MsoNormal">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. <br><br>;-----------------------------<wbr>------------------------------<wbr>-----------<br><br>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_csm.ncl"<br>load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/contributed.<wbr>ncl"<br><br>begin<br>;---Read data <br> a = addfile("/media/Local Disk/NPL/MODIS_FPC/MAIACTAOT.<wbr>h00v02.20000570505.hdf","r")<br> var = short2flt(a->Optical_Depth_<wbr>055(:,:))<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 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> system("/bin/rm -f <a href="http://simple.nc" target="_blank">simple.nc</a>") <br> ncdf = addfile("<a href="http://simple.nc" target="_blank">simple.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" 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 applications<br> ;=============================<wbr>==============================<wbr>========<br> filedimdef(ncdf,"time",-1,<wbr>True) <br> ncdf->var = var<br><br>end<br> <br> <u></u><u></u></p><div><div><p class="MsoNormal"><u></u> <u></u></p></div></div><div><p class="MsoNormal"><br clear="all"><u></u><u></u></p><div><div><div><div><div><div><p class="MsoNormal">Kunal Bali<u></u><u></u></p></div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p style="margin:0cm;margin-bottom:.0001pt"><span style="font-size:9.0pt;font-family:"Tahoma",sans-serif"><u></u> <u></u></span></p></div></div></div></div></div></div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">On Wed, Sep 13, 2017 at 11:27 PM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<u></u><u></u></p><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"><div><div><p class="MsoNormal" style="margin-bottom:12.0pt">No problem, I sorted out. <u></u><u></u></p></div><p class="MsoNormal">Just changed short to float. <u></u><u></u></p></div><div><p class="MsoNormal"><span style="color:#888888"><br clear="all"></span><span class="m_-1899063326404824290m1966040847522573004m2882238953881155132m4735875726629588789m-4928089736255412493m9186617808376476535m-3688407689937875090hoenzb"><span style="color:#888888"><u></u><u></u></span></span></p><div><div><div><div><div><div><p class="MsoNormal"><span style="color:#888888">Kunal Bali</span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="color:#888888">Research Scholar <u></u><u></u></span></p></div><div><p class="MsoNormal"><span style="color:#888888">Radio & Atmospheric Science Division <u></u><u></u></span></p></div><div><p class="MsoNormal"><span style="color:#888888">CSIR - National Physical Laboratory<u></u><u></u></span></p></div><div><p class="MsoNormal"><span style="color:#888888">New Delhi - 110012<u></u><u></u></span></p></div><div><p class="MsoNormal"><span style="color:#888888">India<u></u><u></u></span></p></div><div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p></div><div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p></div><div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p></div><div><p style="margin:0cm;margin-bottom:.0001pt"><span style="font-size:9.0pt;font-family:"Tahoma",sans-serif;color:#888888"><u></u> <u></u></span></p></div></div></div></div></div></div><div><div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">On Wed, Sep 13, 2017 at 11:24 PM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<u></u><u></u></p><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"><div><div><div><div><p class="MsoNormal" style="margin-bottom:12.0pt">Thanks for providing this information. <u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt">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. <u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt">So, how to correct the values?<u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt">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<span style="color:#888888"><br><br><br></span><u></u><u></u></p><div><p class="MsoNormal"><span style="color:#888888"><br clear="all"></span><span class="m_-1899063326404824290m1966040847522573004m2882238953881155132m4735875726629588789m-4928089736255412493m9186617808376476535m-3688407689937875090m7418340254611092360hoenzb"><span style="color:#888888"><u></u><u></u></span></span></p><div><div><div><div><div><div><p class="MsoNormal"><span style="color:#888888">Kunal Bali</span><u></u><u></u></p></div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p><div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p></div><div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p></div><div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p></div><div><p style="margin:0cm;margin-bottom:.0001pt"><span style="font-size:9.0pt;font-family:"Tahoma",sans-serif;color:#888888"><u></u> <u></u></span></p></div></div></div></div></div></div><div><div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">On Wed, Sep 13, 2017 at 10:28 PM, Adam Phillips <<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>> wrote:<u></u><u></u></p><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"><div><p class="MsoNormal">Hi Kunal,<u></u><u></u></p><div><p class="MsoNormal">Yes, clicking on the output netCDF link from the Applications page: <u></u><u></u></p></div><div><p class="MsoNormal"><a href="http://www.ncl.ucar.edu/Applications/o-netcdf.shtml" target="_blank">http://www.ncl.ucar.edu/<wbr>Applications/o-netcdf.shtml</a><u></u><u></u></p></div><div><p class="MsoNormal">I would recommend following the inefficient method #1. Unless you are writing a file with many large variables, the inefficient method works just fine. <u></u><u></u></p></div><div><p class="MsoNormal"><a href="http://www.ncl.ucar.edu/Applications/method_1.shtml" target="_blank">http://www.ncl.ucar.edu/<wbr>Applications/method_1.shtml</a><u></u><u></u></p></div><div><p class="MsoNormal">Good luck,<u></u><u></u></p></div><div><p class="MsoNormal">Adam <u></u><u></u></p></div></div><div><div><div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">On Wed, Sep 13, 2017 at 10:25 AM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<u></u><u></u></p><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"><div><div><p class="MsoNormal" style="margin-bottom:12.0pt">Thank you so much, it worked<u></u><u></u></p></div><p class="MsoNormal">I want to know one more thing.<br>After rearranging the dimensions can we now convert (or write) this arranged file into netcdf format?<u></u><u></u></p><div><div><div><div><p class="MsoNormal" style="margin-bottom:12.0pt"><br><br><u></u><u></u></p></div></div></div></div><div><div><div><div><div><div><div><p class="MsoNormal">Kunal Bali<u></u><u></u></p></div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p style="margin:0cm;margin-bottom:.0001pt"><span style="font-size:9.0pt;font-family:"Tahoma",sans-serif"><u></u> <u></u></span></p></div></div></div></div></div></div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal">On Wed, Sep 13, 2017 at 8:55 PM, Adam Phillips <<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>> wrote:<u></u><u></u></p><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"><div><p class="MsoNormal">Hi Kunal,<u></u><u></u></p><div><p class="MsoNormal">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:<u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">a = addfile("/media/Local Disk/NPL/MODIS_FPC/MAIACTAOT.<wbr>h00v02.20000570505.hdf","r")<br>var = a->Optical_Depth_055(:,:)<br>nlat = 1200<br>nlon = 1200<br>lat = latGlobeFo(nlat, "lat", "latitude", "degrees_north")</span><u></u><u></u></p></div><div><p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:9.5pt">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</span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">var!0 = "lat" ; you can name these dimensions whatever you want..<br>var!1 = "lon"</span><u></u><u></u></p></div><p class="MsoNormal"><span style="font-size:9.5pt">var&lat = lat ; but make sure you refer to the correct named dimensions</span><u></u><u></u></p><div><div><p class="MsoNormal"><span style="font-size:9.5pt">var&lat = lon </span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">;var&XDim_grid1km = lat<br>;var&YDim_grid1km = lon<br> </span><u></u><u></u></p></div><div><p class="MsoNormal"><span style="font-size:9.5pt">If you have any further questions please respond to the ncl-talk email list.</span><u></u><u></u></p></div><div><p class="MsoNormal" style="margin-bottom:12.0pt"><span style="font-size:9.5pt">Adam </span><u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p></div></div></div><div><p class="MsoNormal"><u></u> <u></u></p><div><div><div><p class="MsoNormal">On Wed, Sep 13, 2017 at 8:38 AM, Kunal Bali <<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>> wrote:<u></u><u></u></p></div></div><blockquote style="border:none;border-left:solid #cccccc 1.0pt;padding:0cm 0cm 0cm 6.0pt;margin-left:4.8pt;margin-right:0cm"><div><div><div><div><div><div><p class="MsoNormal" style="margin-bottom:12.0pt">Dear NCL<u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt">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: [<b>YDim_grid1km | 1200] x [XDim_grid1km | 1200]</b><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><u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt">I am trying to read the dimensions of this file as<br><br>begin<br>;---Read data <br> a = addfile("/media/Local Disk/NPL/MODIS_FPC/MAIACTAOT.<wbr>h00v02.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", "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> <u></u><u></u></p></div><p class="MsoNormal" style="margin-bottom:12.0pt">But the error appeared as <br><b>fatal:Variable (lat) is undefined</b><u></u><u></u></p><div><div><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal" style="margin-bottom:12.0pt">So, could anyone please let me know that how to read dimension of this file. <u></u><u></u></p></div><div><p class="MsoNormal">Thank You<u></u><u></u></p></div><div><p class="MsoNormal"><u></u> <u></u></p><div><p class="MsoNormal"><u></u> <u></u></p></div><div><p class="MsoNormal">Regards<span style="color:#888888"><br clear="all"></span><u></u><u></u></p></div><div><div><div><div><div><div><div><div><p class="MsoNormal"><span style="color:#888888">Kunal Bali<u></u><u></u></span></p></div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p><div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p></div><div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p></div><div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p></div><div><p style="margin:0cm;margin-bottom:.0001pt"><span style="font-size:9.0pt;font-family:"Tahoma",sans-serif;color:#888888"><u></u> <u></u></span></p></div></div></div></div></div></div></div></div></div></div></div></div><p class="MsoNormal"><u></u> <u></u></p></div></div><p class="MsoNormal" style="margin-bottom:12.0pt">______________________________<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" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><u></u><u></u></p></blockquote></div><p class="MsoNormal"><span style="color:#888888"><br><br clear="all"></span><span class="m_-1899063326404824290m1966040847522573004m2882238953881155132m4735875726629588789m-4928089736255412493m9186617808376476535m-3688407689937875090m7418340254611092360m-6183322678709470722m4501774087300120891hoenzb"><span style="color:#888888"><u></u><u></u></span></span></p><div><p class="MsoNormal"><u></u> <u></u></p></div><p class="MsoNormal"><span class="m_-1899063326404824290m1966040847522573004m2882238953881155132m4735875726629588789m-4928089736255412493m9186617808376476535m-3688407689937875090m7418340254611092360m-6183322678709470722m4501774087300120891hoenzb"><span style="color:#888888">-- </span><u></u><u></u></span></p><div><div><div><div><div><div><div><div><div><div><div><p class="MsoNormal"><span style="color:#888888">Adam Phillips </span><u></u><u></u></p></div><p class="MsoNormal"><span style="color:#888888">Associate Scientist, Climate and Global Dynamics Laboratory, NCAR<u></u><u></u></span></p></div></div><div><p class="MsoNormal"><span style="color:#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/<wbr>asphilli/</a> <a href="tel:(303)%20497-1726" target="_blank">303-497-1726</a> <u></u><u></u></span></p></div><div><div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p></div></div></div></div></div></div></div></div></div></div></div></blockquote></div><p class="MsoNormal"><u></u> <u></u></p></div></div><p class="MsoNormal" style="margin-bottom:12.0pt"><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" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><u></u><u></u></p></blockquote></div><p class="MsoNormal"><br><br clear="all"><span class="m_-1899063326404824290m1966040847522573004m2882238953881155132m4735875726629588789m-4928089736255412493m9186617808376476535hoenzb"><span style="color:#888888"><u></u><u></u></span></span></p><div><p class="MsoNormal"><u></u> <u></u></p></div><p class="MsoNormal"><span class="m_-1899063326404824290m1966040847522573004m2882238953881155132m4735875726629588789m-4928089736255412493m9186617808376476535hoenzb"><span style="color:#888888">-- </span><u></u><u></u></span></p><div><div><div><div><div><div><div><div><div><div><div><p class="MsoNormal"><span style="color:#888888">Adam Phillips </span><u></u><u></u></p></div><p class="MsoNormal"><span style="color:#888888">Associate Scientist, Climate and Global Dynamics Laboratory, NCAR<u></u><u></u></span></p></div></div><div><p class="MsoNormal"><span style="color:#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/<wbr>asphilli/</a> <a href="tel:(303)%20497-1726" target="_blank">303-497-1726</a> <u></u><u></u></span></p></div><div><div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p></div></div></div></div></div></div></div></div></div></div></div></div></div></blockquote></div><p class="MsoNormal"><u></u> <u></u></p></div></div></div></div></blockquote></div><p class="MsoNormal"><u></u> <u></u></p></div></div></div></blockquote></div><p class="MsoNormal"><u></u> <u></u></p></div></div><p class="MsoNormal" style="margin-bottom:12.0pt"><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" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><u></u><u></u></p></blockquote></div><p class="MsoNormal"><br><br clear="all"><span class="m_-1899063326404824290m1966040847522573004m2882238953881155132m4735875726629588789hoenzb"><span style="color:#888888"><u></u><u></u></span></span></p><div><p class="MsoNormal"><u></u> <u></u></p></div><p class="MsoNormal"><span class="m_-1899063326404824290m1966040847522573004m2882238953881155132m4735875726629588789hoenzb"><span style="color:#888888">-- </span><u></u><u></u></span></p><div><div><div><div><div><div><div><div><div><div><div><p class="MsoNormal"><span style="color:#888888">Adam Phillips </span><u></u><u></u></p></div><p class="MsoNormal"><span style="color:#888888">Associate Scientist, Climate and Global Dynamics Laboratory, NCAR<u></u><u></u></span></p></div></div><div><p class="MsoNormal"><span style="color:#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/<wbr>asphilli/</a> <a href="tel:(303)%20497-1726" target="_blank">303-497-1726</a> <u></u><u></u></span></p></div><div><div><p class="MsoNormal"><span style="color:#888888"><u></u> <u></u></span></p></div></div></div></div></div></div></div></div></div></div></div></div></div></blockquote></div><p class="MsoNormal"><u></u> <u></u></p></div></div><p class="MsoNormal" style="margin-bottom:12.0pt"><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" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><u></u><u></u></p></blockquote></div><p class="MsoNormal"><u></u> <u></u></p></div></div></div></blockquote></div><p class="MsoNormal"><u></u> <u></u></p></div></div></div></div><p class="MsoNormal" style="margin-bottom:12.0pt"><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" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><u></u><u></u></p></blockquote></div><p class="MsoNormal"><u></u> <u></u></p></div></div></div></blockquote></div><p class="MsoNormal"><u></u> <u></u></p></div><p class="MsoNormal">______________________________<wbr>_________________ ncl-talk mailing list <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a> List instructions, subscriber options, unsubscribe: <a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a> <u></u><u></u></p></div></div></div></div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">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/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>