<div dir="ltr"><div>Thank you for providing this information. <br><br></div>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]:Execute: Error occurred at or near line 22 in file netcdf2.ncl<br><br> <br><br>;----------------------------------------------------------------------<br><br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.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_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">simple2.nc</a>") <br>   ;     ncdf = addfile("<a href="http://simple2.nc">simple2.nc</a>" ,"c")  ; open output netCDF file<br><br>    ;===================================================================<br>    ; create global attributes of the file (optional)<br>    ;===================================================================<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">original-file.nc</a>"<br>       fAtt@Conventions   = "None"<br>       fAtt@creation_date = systemfunc ("date")<br>       fileattdef( ncdf, fAtt )            ; copy file attributes<br><br>    ;===================================================================<br>    ; make time an UNLIMITED dimension; recommended  for most applications<br>    ;===================================================================<br>       filedimdef(ncdf,"time",-1,True)<br><br>;       ncdf->var&lat = lat<br>;       ncdf->var&lon = lon<br>       ncdf->var = var<br><br>end<br><br></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Kunal Bali<br></div><div>Research Scholar <br></div><div>Radio & Atmospheric Science Division <br></div><div>CSIR - National Physical Laboratory<br></div><div>New Delhi - 110012<br></div><div>India<br></div><div><br></div><div><br></div><div><br></div><div><p style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left: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 Mon, Sep 18, 2017 at 9:04 PM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div>Loop over the files ...<br></div><div><br> diri = "/media/Local Disk/NPL/MODIS_FPC/"   ; input directory<br></div> fili  = systemfunc("cd "+diri +; ls <span><span style="font-size:12.8px">MAIACTAOT.</span><span style="font-size:12.8px">h00v02*hdf")<br></span></span></div><span><span style="font-size:12.8px"> nfili = dimsizes(fili)<br></span></span></div><span><span style="font-size:12.8px"> print("nfili="+nfili)<br><br></span></span></div><div><span><span style="font-size:12.8px"> dirnc = </span></span><span><span style="font-size:12.8px">"/media/Local Disk/NPL/MODIS_FPC/"   ; output (netCDF) directory</span></span></div><div><span><span style="font-size:12.8px"><br></span></span></div><span><span style="font-size:12.8px"> do nf=0,nfili-1<br></span></span></div><div><span><span style="font-size:12.8px">      pthi = diri+fili(nf)<br></span></span></div><div><span><span style="font-size:12.8px">      f      = addfile(pthi,"r")<br><br>      ....<br><br>     </span></span>filroot = <strong>str_get_cols</strong>(fili(nf), 0, 27)  ; eg: "MAIACTAOT.h00v02.20000570505"<span><span style="font-size:12.8px"></span></span></div><span><span style="font-size:12.8px">     filnc    = filroot+".nc"<br></span></span></div><span><span style="font-size:12.8px">     pthnc  = dirnc + filnc<br></span></span><br><span><span style="font-size:12.8px">     system("/bin/rm -f "+pthnc)  <br>     ncdf = addfile(pthnc,"c")  ; open new netCDF file<br><br>     ....<br><br><br></span></span></div><span><span style="font-size:12.8px">end do    ; end 'nf' loop<br></span></span><div><div><span><span style="font-size:12.8px"><br>=====<br></span></span></div><div><span><span style="font-size:12.8px">If there is something you do not understand use some print statements.<br><br></span></span></div><div><span><span style="font-size:12.8px">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><br><br></span></span></div><div><span><span style="font-size:12.8px">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><br><br></span></span></div><div><span><span style="font-size:12.8px">Good Luck<br></span></span></div><div><span><span style="font-size:12.8px"><br><br></span></span></div></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Sep 18, 2017 at 6:11 AM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@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 dir="ltr"><div>Thnks for the reply. <br></div><div><br></div><div>that will change the name of only one file. isn't it. but what about the other files <br><br></div>I have 365 files in one directory and I need to change all the file name one by one.<span class="m_1966040847522573004HOEnZb"><font color="#888888"><br><br></font></span><div class="gmail_extra"><span class="m_1966040847522573004HOEnZb"><font color="#888888"><br clear="all"><div><div class="m_1966040847522573004m_2882238953881155132gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Kunal Bali<br></div><br><br><div><br></div><div><br></div><div><p style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></p></div></div></div></div></div></div></font></span><div><div class="m_1966040847522573004h5">
<br><div class="gmail_quote">On Mon, Sep 18, 2017 at 5:14 PM, Dennis Shea <span dir="ltr"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Really? <br><br>Change:<span><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></span>To:<br>      system("/bin/rm -f <a href="http://simple.nc" target="_blank">MAIACTAOT.h00v02.20000570505.n<wbr>c</a>")  <br>      ncdf = addfile("MAIACTAOT.h<a href="http://00v02.20000570505.nc" target="_blank">00v02.2000<wbr>0570505.nc</a>" ,"c") <br><div><div><br></div></div></div><div class="m_1966040847522573004m_2882238953881155132HOEnZb"><div class="m_1966040847522573004m_2882238953881155132h5"><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Sep 18, 2017 at 12:37 AM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@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 dir="ltr"><div><div><div><div><div><div><div><div><div>Thanks, it worked.<br><br></div>one more question is that.<br><br></div>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")<br><br></div>If I want to keep the original file name with the output file name then what should I do?<br><br></div>I mean the original file name is MAIACTAOT.h00v02.20003660700.h<wbr>df<br></div>Now I want to create my netcdf file name  as <a href="http://MAIACTAOT.h00v02.20003660700.nc" target="_blank">MAIACTAOT.h00v02.20003660700.n<wbr>c</a> NOT <a href="http://simple.nc" target="_blank">simple.nc</a>. <br><br></div>I can not simply write <a href="http://simple.nc" target="_blank">simple.nc</a> file name because I have many files. <br></div>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.<br><br></div>I hope you have got my query.<br><br></div>please let me know that too.<br><div><div><div><div><div><div><div><br></div><div>Thank You<br></div><div><br></div></div></div></div></div></div></div><div class="gmail_extra"><br clear="all"><div><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493gmail_signature" data-smartmail="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-top:0px;margin-right:0px;margin-bottom:0px;margin-left: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, Sep 14, 2017 at 9:09 PM, Adam Phillips <span dir="ltr"><<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Kunal,<div>It looks to me like your latitudes are flipped. You are setting up and assigning your coordinate lat/lon variables in this coding:</div><span><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">nlat = 1200</span><br style="font-size:12.8px"><span style="font-size:12.8px">nlon = 1200</span><br style="font-size:12.8px"><span style="font-size:12.8px">lat = latGlobeFo(nlat, "lat", "latitude", "degrees_north")</span></div><div><span style="font-size:12.8px">lon = lonGlobeFo(nlon, "lon", "longitude", "degrees_east")   </span><br style="font-size:12.8px"><span style="font-size:12.8px">lat = lat(::-1)</span><br style="font-size:12.8px"><span style="font-size:12.8px">lon = (/ lon - 180. /)  ; subtract 180 from all values </span></div><div><span style="font-size:12.8px">lon&lon = lon           ; update coordinates</span><br style="font-size:12.8px"><br style="font-size:12.8px"><span style="font-size:12.8px">var!0 = "lat"    ; you can name these dimensions whatever you want..</span><br style="font-size:12.8px"><span style="font-size:12.8px">var!1 = "lon"</span><br style="font-size:12.8px"><span style="font-size:12.8px">var&lat = lat     ; but make sure you refer to the correct named dimensions</span><br style="font-size:12.8px"><span style="font-size:12.8px">var&lon = lon </span><br></div><div><br></div></span><div>Assuming I'm right on your latitudes being the issue, one of the following two modifications should work:</div><div>Delete this line:</div><div><span style="font-size:12.8px">lat = lat(::-1)</span><br></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">If the above doesn't fix it, try this:</span></div><div><span style="font-size:12.8px">Change this:</span></div><div><span style="font-size:12.8px">lat = lat(::-1)</span><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">to this:</span></div><div><span style="font-size:12.8px">lat = lat(::-1)</span><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">lat&lat = lat</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">If you continue to have issues please respond to ncl-talk.</span></div><span class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493HOEnZb"><font color="#888888"><div><span style="font-size:12.8px">Adam</span></div></font></span></div><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493HOEnZb"><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Sep 13, 2017 at 12:35 PM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@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 dir="ltr"><div><br></div>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/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 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 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,True<wbr>) <br>       ncdf->var = var<br><br>end<br>       <br>     <br><div><div><br></div></div><div class="gmail_extra"><br clear="all"><div><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090gmail_signature" data-smartmail="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-top:0px;margin-right:0px;margin-bottom:0px;margin-left: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 Wed, Sep 13, 2017 at 11:27 PM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@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 dir="ltr"><div>No problem, I sorted out. <br><br></div>Just changed short to float. <br></div><div class="gmail_extra"><span class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090HOEnZb"><font color="#888888"><br clear="all"><div><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Kunal Bali<br></div><div>Research Scholar <br></div><div>Radio & Atmospheric Science Division <br></div><div>CSIR - National Physical Laboratory<br></div><div>New Delhi - 110012<br></div><div>India<br></div><div><br></div><div><br></div><div><br></div><div><p style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></p></div></div></div></div></div></div></font></span><div><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090h5">
<br><div class="gmail_quote">On Wed, Sep 13, 2017 at 11:24 PM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@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 dir="ltr"><div><div><div>Thanks for providing this information. <br><br></div>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. <br><br></div>So, how to correct the values?<br><br></div>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 class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360HOEnZb"><font color="#888888"><br><br><br><br></font></span><div class="gmail_extra"><span class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360HOEnZb"><font color="#888888"><br clear="all"><div><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360m_-6183322678709470722gmail_signature" data-smartmail="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-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;border-collapse:collapse;font-family:Tahoma,Verdana;font-size:12px"><font color="#1F497D"><br></font></p></div></div></div></div></div></div></font></span><div><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360h5">
<br><div class="gmail_quote">On Wed, Sep 13, 2017 at 10:28 PM, Adam Phillips <span dir="ltr"><<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Kunal,<div>Yes, clicking on the output netCDF link from the Applications page: </div><div><a href="http://www.ncl.ucar.edu/Applications/o-netcdf.shtml" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/o-netcdf.shtml</a><br></div><div>I would recommend following the inefficient method #1. Unless you are writing a file with many large variables, the inefficient method works just fine. <br></div><div><a href="http://www.ncl.ucar.edu/Applications/method_1.shtml" target="_blank">http://www.ncl.ucar.edu/Applic<wbr>ations/method_1.shtml</a></div><div>Good luck,</div><div>Adam  </div></div><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360m_-6183322678709470722HOEnZb"><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360m_-6183322678709470722h5"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Sep 13, 2017 at 10:25 AM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@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 dir="ltr"><div>Thank you so much, it worked<br><br></div>I want to know one more thing.<br>After rearranging the dimensions can we now convert (or write) this arranged file into netcdf format?<br><div><div><div><div><br><br><br></div></div></div></div><div class="gmail_extra"><div><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360m_-6183322678709470722m_4501774087300120891m_7649397337960907787gmail_signature" data-smartmail="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-top:0px;margin-right:0px;margin-bottom:0px;margin-left: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 Wed, Sep 13, 2017 at 8:55 PM, Adam Phillips <span dir="ltr"><<a href="mailto:asphilli@ucar.edu" target="_blank">asphilli@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Kunal,<div>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:</div><div><span><span style="font-size:12.8px">a = addfile("/media/Local Disk/NPL/MODIS_FPC/MAIACTAOT.</span><span style="font-size:12.8px">h<wbr>00v02.20000570505.hdf","r")</span><br style="font-size:12.8px"><span style="font-size:12.8px">var  = a->Optical_Depth_055(:,:)</span><br style="font-size:12.8px"></span><span><span style="font-size:12.8px">nlat = 1200</span><br style="font-size:12.8px"><span style="font-size:12.8px">nlon = 1200</span><br style="font-size:12.8px"><span style="font-size:12.8px">lat = latGlobeFo(nlat, "lat", "latitude", "degrees_north")</span></span></div><span><div><span style="font-size:12.8px">lon = lonGlobeFo(nlon, "lon", "longitude", "degrees_east")   </span><br style="font-size:12.8px"><span style="font-size:12.8px">lat = lat(::-1)</span><br style="font-size:12.8px"><span style="font-size:12.8px">lon = (/ lon - 180. /)  ; subtract 180 from all values </span><br style="font-size:12.8px"><span style="font-size:12.8px">lon&lon = lon           ; update coordinates</span><br style="font-size:12.8px"><span style="font-size:12.8px"><br></span></div></span><div><span style="font-size:12.8px">var!0 = "lat"    ; you can name these dimensions whatever you want..</span><br style="font-size:12.8px"><span style="font-size:12.8px">var!1 = "lon"</span><span style="font-size:12.8px"><br></span></div><span style="font-size:12.8px">var&lat = lat     ; but make sure you refer to the correct named dimensions</span><div><div><span style="font-size:12.8px">var&lat = lon</span><span style="font-size:12.8px"> </span></div><div><span style="font-size:12.8px">;var&XDim_grid1km = lat</span><br style="font-size:12.8px"><span style="font-size:12.8px">;var&YDim_grid1km = lon</span><br style="font-size:12.8px"><span style="font-size:12.8px"> </span></div><div><span style="font-size:12.8px">If you have any further questions please respond to the ncl-talk email list.</span></div><div><span style="font-size:12.8px">Adam </span><br style="font-size:12.8px"><br></div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360m_-6183322678709470722m_4501774087300120891m_7649397337960907787h5">On Wed, Sep 13, 2017 at 8:38 AM, Kunal Bali <span dir="ltr"><<a href="mailto:kunal.bali9@gmail.com" target="_blank">kunal.bali9@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360m_-6183322678709470722m_4501774087300120891m_7649397337960907787h5"><div dir="ltr"><div><div><div>Dear NCL<br><br></div>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><br></div>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.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", "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></div>But the error appeared as  <br><b>fatal:Variable (lat) is undefined</b><br><br><div><div><div><br></div><div>So, could anyone please let me know that how to read dimension of this file. <br><br></div><div>Thank You<br></div><div><br><div><br></div><div>Regards<span class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360m_-6183322678709470722m_4501774087300120891m_7649397337960907787m_-8297990099082400272HOEnZb"><font color="#888888"><br clear="all"></font></span></div><span class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360m_-6183322678709470722m_4501774087300120891m_7649397337960907787m_-8297990099082400272HOEnZb"><font color="#888888"><div><div><div><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360m_-6183322678709470722m_4501774087300120891m_7649397337960907787m_-8297990099082400272m_2392430217485589217gmail_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>
</div></div></font></span></div></div></div></div>
<br></div></div>______________________________<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></blockquote></div><span class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360m_-6183322678709470722m_4501774087300120891m_7649397337960907787HOEnZb"><font color="#888888"><br><br clear="all"><span class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360m_-6183322678709470722m_4501774087300120891HOEnZb"><font color="#888888"><div><br></div>-- <br><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360m_-6183322678709470722m_4501774087300120891m_7649397337960907787m_-8297990099082400272gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist,  </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphill<wbr>i/</a>   </font></span><span><font color="#888888"><a href="tel:(303)%20497-1726" value="+13034971726" target="_blank">303-497-1726</a> </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</font></span></font></span></div>
</blockquote></div><br></div></div>
<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></blockquote></div><br><br clear="all"><span class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535HOEnZb"><font color="#888888"><div><br></div>-- <br><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535m_-3688407689937875090m_7418340254611092360m_-6183322678709470722m_4501774087300120891gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist,  </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphill<wbr>i/</a>   </font></span><span><font color="#888888"><a href="tel:(303)%20497-1726" value="+13034971726" target="_blank">303-497-1726</a> </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</font></span></div><span class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535HOEnZb"><font color="#888888">
</font></span></div></div></blockquote></div><span class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535HOEnZb"><font color="#888888"><br></font></span></div></div></div></div><span class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535HOEnZb"><font color="#888888">
</font></span></blockquote></div><span class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535HOEnZb"><font color="#888888"><br></font></span></div></div></div>
</blockquote></div><br></div></div>
<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></blockquote></div><br><br clear="all"><span class="m_1966040847522573004m_2882238953881155132m_4735875726629588789HOEnZb"><font color="#888888"><div><br></div>-- <br><div class="m_1966040847522573004m_2882238953881155132m_4735875726629588789m_-4928089736255412493m_9186617808376476535gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div><div><span><font color="#888888">Adam Phillips <br></font></span></div><span><font color="#888888">Associate Scientist,  </font></span><span><font color="#888888">Climate and Global Dynamics Laboratory, NCAR<br></font></span></div></div><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli/" target="_blank">www.cgd.ucar.edu/staff/asphill<wbr>i/</a>   </font></span><span><font color="#888888"><a href="tel:(303)%20497-1726" value="+13034971726" target="_blank">303-497-1726</a> </font></span></div><span><font color="#888888"></font></span><div><div><span><font color="#888888"><br></font></span><div><span><font color="#888888"><a href="http://www.cgd.ucar.edu/staff/asphilli" target="_blank"></a></font></span></div></div></div></div></div></div></div></div></div></div></div>
</font></span></div>
</div></div></blockquote></div><br></div></div>
<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></blockquote></div><br></div>
</div></div></blockquote></div><br></div></div></div></div>
<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></blockquote></div><br></div>
</div></div></blockquote></div><br></div>