<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">The first step I told is to reshape &quot;tasmask&quot; using reshape function. </div><div class="gmail_quote"><br></div><div class="gmail_quote">Secondly what I told is to put some arbitrary lat-lon coordinate value.  </div><div class="gmail_quote"><br></div><div class="gmail_quote">On 6 July 2017 at 16:29, S Br <span dir="ltr">&lt;<a href="mailto:sbr.climate@gmail.com" target="_blank">sbr.climate@gmail.com</a>&gt;</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 Gibies,<div>Thanks for your reply. I have added as you suggested. I get some error in the code. Do you think, I have done correctly.</div><div><br></div><div><span class=""><div>; ******************************<wbr>*****************</div><div>; create new date array for use on the plot</div><div>; ******************************<wbr>*****************</div><div>  filo = &quot;tas_aave_Asia_bandpass30-100.<wbr>nc&quot;</div><div>  system (&quot;/bin/rm &quot;+filo) ; remove any pre-existing file</div><div>  fo = addfile(filo , &quot;c&quot;) ; open output file</div></span><span class=""><div>  filedimdef(fo,&quot;time&quot;,-1,True)       ; make &#39;time&#39; unlimited</div></span><div>  tasmask=xBPF</div><div>  tasmask!0=&quot;time&quot;</div><div>  tasmask!1=&quot;lat&quot;</div><div>  tasmask&amp;lat=<u><font color="#ff0000">&quot;lat&quot; </font></u></div><div>  tasmask!2=&quot;lon&quot;</div><div>  tasmask&amp;lon=<u><font color="#ff0000">&quot;lon&quot;</font></u></div><div>  fo-&gt;tasmask = tasmask ; write ud to a file</div><div>end</div></div><div><br></div><div>I get the error as,</div><div><br></div><div><div>fatal:Variable (tasmask) has (1) dimensions can not write to dimension (1)</div><div>fatal:[&quot;Execute.c&quot;:8575]:<wbr>Execute: Error occurred at or near line 38 in file bandpass.ncl</div></div><div><br></div><div>Thanks.</div><span class="HOEnZb"><font color="#888888"><div>SB</div><div><br></div><div><br></div><div><br></div></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jul 6, 2017 at 10:13 AM, gibies george <span dir="ltr">&lt;<a href="mailto:gibiesgeorge@gmail.com" target="_blank">gibiesgeorge@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="auto"><div>I have seen those mails from Dr. Dennis and Dr. Adam.<div dir="auto"><br></div><div dir="auto">You can use reshape function to add additional coordinate to a variable. Define some single value coordinate for latitude and longitude.</div><div dir="auto"><br></div><div dir="auto">varname!0 =&quot;time&quot;</div><div dir="auto">varname!1 = &quot;lat&quot;</div><div dir="auto">varname&amp;lat = .......</div><div dir="auto">varname!2 = &quot;lon&quot;</div><div dir="auto">varname&amp;lon = ......</div><div><div class="m_3752893363194926388h5"><div dir="auto"><br></div><div dir="auto"><br></div><br><div class="gmail_extra"><br><div class="gmail_quote">On Jul 6, 2017 2:08 PM, &quot;S Br&quot; &lt;<a href="mailto:sbr.climate@gmail.com" target="_blank">sbr.climate@gmail.com</a>&gt; wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Gibies,<div>NCL also provided the same solution as you have mentioned.<br><div><div>Could you please tell me how can I add single value coordinates (time, lat, lon) to the variable.</div></div></div><div>Hope, you can write me the command for me.</div><div><br></div><div>Thanks.</div><div>SB</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jul 5, 2017 at 5:58 PM, gibies george <span dir="ltr">&lt;<a href="mailto:gibiesgeorge@gmail.com" target="_blank">gibiesgeorge@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="auto">Yes, you can add single value coordinates (time lat lon) to the variable before saving it to nc file.</div><div class="m_3752893363194926388m_-3993977731678294765m_-2811224259625025701m_7970586787079507880HOEnZb"><div class="m_3752893363194926388m_-3993977731678294765m_-2811224259625025701m_7970586787079507880h5"><div class="gmail_extra"><br><div class="gmail_quote">On Jul 5, 2017 9:58 PM, &quot;S Br&quot; &lt;<a href="mailto:sbr.climate@gmail.com" target="_blank">sbr.climate@gmail.com</a>&gt; wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Is it possible to do something inside the NCL script so that it can be opened by GrADS.</div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jul 5, 2017 at 5:25 PM, gibies george <span dir="ltr">&lt;<a href="mailto:gibiesgeorge@gmail.com" target="_blank">gibiesgeorge@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="auto">This is because grads required standard file coordinates structure. Time lat lon etc.</div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="m_3752893363194926388m_-3993977731678294765m_-2811224259625025701m_7970586787079507880m_-138813461870546186m_2678604644904053558h5">On Jul 5, 2017 9:43 PM, &quot;S Br&quot; &lt;<a href="mailto:sbr.climate@gmail.com" target="_blank">sbr.climate@gmail.com</a>&gt; wrote:<br type="attribution"></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="m_3752893363194926388m_-3993977731678294765m_-2811224259625025701m_7970586787079507880m_-138813461870546186m_2678604644904053558h5"><div dir="ltr">Hi,<div>I have applied bandpass filter to my one dimentional NetCDF file. I have written the output to a NetCDF file but this output file is not opening in GrADS. I get the following error.</div><div><br></div><div>ga-&gt; sdfopen tas_aave_Asia_bandpass30-100.n<wbr>c</div><div><div>Scanning self-describing file:  tas_aave_Asia_bandpass30-100.<wbr>nc</div><div>gadsdf: SDF file has no discernable X coordinate.</div><div>  To open this file with GrADS, use a descriptor file with an XDEF entry.</div><div>  Documentation is at <a href="http://cola.gmu.edu/grads/gadoc/SDFdescriptorfile.html" target="_blank">http://cola.gmu.edu/grads/gado<wbr>c/SDFdescriptorfile.html</a></div></div><div><br></div><div>Could you please suggest if I am doing somewhere wrong in my NCL script . The input file is perfectly fine with GrADS.<br></div><div><br></div><div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/gsn_code.ncl&quot;</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/gsn_csm.ncl&quot;</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/contributed.ncl&quot;</div><div>begin</div><div><br></div><div>   fn  = &quot;tas_aave_Asia.nc&quot; ; define filename</div><div>   in  = addfile(fn,&quot;r&quot;)                                 ; open netcdf file</div><div>    x  = in-&gt;tasmask(:,0,0)                                       ; get data</div><div><br></div><div>; ******************************<wbr>*****************</div><div>; create the filter weights and apply</div><div>; ******************************<wbr>*****************</div><div><br></div><div>  ihp     = 2                             ; band pass</div><div>  sigma   = 1.0                           ; Lanczos sigma</div><div><br></div><div>  nWgt    = 201                           ; loose 100 each end</div><div>  fca     = 1./100.                       ; start freq</div><div>  fcb     = 1./30.                        ; last  freq</div><div>  wgt     = filwgts_lanczos (nWgt, ihp, fca, fcb, sigma )</div><div>  xBPF    = wgt_runave ( x, wgt, 0 )      ; 30-100 days</div><div><br></div><div>   copy_VarMeta(x,xBPF)</div><div>   dNames= getvardims(xBPF)</div><div>   dSIZES=dimsizes(xBPF)</div><div>   print (dNames)</div><div>   print (dSIZES)</div><div><br></div><div>; ******************************<wbr>*****************</div><div>; create new date array for use on the plot</div><div>; ******************************<wbr>*****************</div><div>  filo = &quot;tas_aave_Asia_bandpass30-100.<wbr>nc&quot;</div><div>  system (&quot;/bin/rm &quot;+filo) ; remove any pre-existing file</div><div>  fo = addfile(filo , &quot;c&quot;) ; open output file</div><div>  tasmax=xBPF(:)</div><div>  fo-&gt;tasmax = tasmax ; write ud to a file</div><div>end</div></div><div><br></div><div><br></div><div>Thanks.</div><div>SB</div><div><br></div><div><br></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></div>
</blockquote></div><br></div>
</blockquote></div></div>
</div></div></blockquote></div><br></div>
</blockquote></div></div>
</div></div></div></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div></div>