[gradsusr] Printing a lot of maps from NetCDF file

Krzysiek krzyba at o2.pl
Thu May 12 14:50:25 EDT 2011


Dear Grads Users,
I have just created a simple script:

'reinit'
'sdfopen /cygdrive/c/ab.nc'
'set lon -20 40'
'set lat 30 70'
'set gxout contour'
'set cint 5'
'set clab forced'
'set ccolor 1'
'set cthick 10'
'set map 1'
'set timelab off'
'set grads off'
'set cterp on'
'set csmooth on'
'd slp/100'
'q time'
res=subwrd(result,3)
ymd=substr(res,4,12)
say ymd
'draw title 'ymd
'printim slp.png x1024 y768 white'
return

NetCDF file contain daily data from 1jan1990 to 31dec1990 and I would 
like to print (at one go) maps which concern all days (365 maps, png 
format). Moreover, this files should be named by the title of particular 
map (I arranged in script date format like 01JAN1990, 02JAN1990 
......etc.).
If anyone knows how to change my scripts in order to do what I want?

Regards

Chris




More information about the gradsusr mailing list