problem using sdfopen

Arindam Chakraborty arch at IO.MET.FSU.EDU
Mon Jan 23 16:19:40 EST 2006


hi Mukul,
        try making a descriptor file for those netcdf files with
template option and open using xdfopen:

% cat prec.des
dset ^precip.%y4.nc
options template
tdef time 18615 linear 1jan1948 1dy

* assuming NetCDF file time dimnsion name is 'time'

ga-> xdfopen prec.des


----------------------------
 ARINDAM CHAKRABORTY
 Department of Meteorology
 Florida State University
 Tallahassee, FL-32306, USA
 Tel: +001-850-6443524 (Off)
      +001-850-5758550 (Res)
 Fax: +001-850-6449642
____________________________

On Mon, 23 Jan 2006, Mukul Tewari wrote:

> Hi,
>
> I have 50 yrs of precip data in netcdf format. I am using the following
> script to plot the difference in precip but it gives segmentation fault. If
> I try to plot the difference for the first year it does it correctly. I see
> the problem when I try to plot the differnce in precip for 2 different
> years??? I am using the following script
>
>
> ------------------------------------------------
> 'sdfopen precip.1948.nc precip.%y4.nc 18615'
>
> 'run rgbset.gs'
> 'set display color white'
> 'c'
>
> 'set grads off'
> 'set grid off'
> 'set gxout shaded'
> 'set mpdset hires'
> 'set clevs  -1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.1 0.1 0.5 0.6 0.7 0.8 0.9 1.0 1.2
> 1.4 1.6 1.8 2'
> 'set ccols 49 47 46 45 44 43 42 0 33 35 37 39 21 22 23 24 25 27 29'
>
> 'd precip (t=368) - precip (t=360)'
>
> 'run cbarn.gs'
> -------------------------------------------------
>
> thanks for the help,
> Mukul
>



More information about the gradsusr mailing list