<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.nc</div><div><div>Scanning self-describing file:  tas_aave_Asia_bandpass30-100.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">http://cola.gmu.edu/grads/gadoc/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/nclscripts/csm/gsn_code.ncl&quot;</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/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>; ***********************************************</div><div>; create the filter weights and apply</div><div>; ***********************************************</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>; ***********************************************</div><div>; create new date array for use on the plot</div><div>; ***********************************************</div><div>  filo = &quot;tas_aave_Asia_bandpass30-100.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>