[ncl-talk] NetCDF output file from NCL not opening in GrADS

[Dennis Shea]

In my view, this question should be sent to the GrADS User forum not
ncl-talk.

Did you look at the created nc file?

%> ncdump  tas_aave_Asia_bandpass30-100.nc  |  less

===========================================
Maybe try making the 'time' dimension unlimited.


  filo = "tas_aave_Asia_bandpass30-100.nc"
  system ("/bin/rm "+filo) ; remove any pre-existing file
  fo  = addfile(filo , "c") ; open output file
  *filedimdef*
<http://www.ncl.ucar.edu/Document/Functions/Built-in/filedimdef.shtml>(fo,"time",-1,True)
; make 'time' unlimited
  tasmax = xBPF


On Wed, Jul 5, 2017 at 10:13 AM, S Br <sbr.climate at gmail.com> wrote:

> Hi,
> 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.
>
> ga-> sdfopen tas_aave_Asia_bandpass30-100.nc
> Scanning self-describing file:  tas_aave_Asia_bandpass30-100.nc
> gadsdf: SDF file has no discernable X coordinate.
>   To open this file with GrADS, use a descriptor file with an XDEF entry.
>   Documentation is at http://cola.gmu.edu/grads/
> gadoc/SDFdescriptorfile.html
>
> Could you please suggest if I am doing somewhere wrong in my NCL script .
> The input file is perfectly fine with GrADS.
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> begin
>
>    fn  = "tas_aave_Asia.nc" ; define filename
>    in  = addfile(fn,"r")                                 ; open netcdf file
>     x  = in->tasmask(:,0,0)                                       ; get
> data
>
> ; ***********************************************
> ; create the filter weights and apply
> ; ***********************************************
>
>   ihp     = 2                             ; band pass
>   sigma   = 1.0                           ; Lanczos sigma
>
>   nWgt    = 201                           ; loose 100 each end
>   fca     = 1./100.                       ; start freq
>   fcb     = 1./30.                        ; last  freq
>   wgt     = filwgts_lanczos (nWgt, ihp, fca, fcb, sigma )
>   xBPF    = wgt_runave ( x, wgt, 0 )      ; 30-100 days
>
>    copy_VarMeta(x,xBPF)
>    dNames= getvardims(xBPF)
>    dSIZES=dimsizes(xBPF)
>    print (dNames)
>    print (dSIZES)
>
> ; ***********************************************
> ; create new date array for use on the plot
> ; ***********************************************
>   filo = "tas_aave_Asia_bandpass30-100.nc"
>   system ("/bin/rm "+filo) ; remove any pre-existing file
>   fo = addfile(filo , "c") ; open output file
>   tasmax=xBPF(:)
>   fo->tasmax = tasmax ; write ud to a file
> end
>
>
> Thanks.
> SB
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170705/8d612799/attachment.html