'reinit' 'open prec_monthly_clim_v1.0.sgi.ctl' 'enable print fig.meta' 'set gxout grfill' 'set grid off' 'set lat -65 65' 'set lon 0 390' 'set z 1' * xmin=1.0 xmax=10.0 ymin=1.0 ymax=7.0 * mon=1 while(mon<13) x1=xmin x2=xmax y2=ymax y1=ymin 'set t 'mon 'set grads off' 'set parea 'x1' 'x2' 'y1' 'y2 'set rgb 91 220 220 220' 'run scripts_color.gs' 'set ccols 91 31 33 35 37 39 43 45 47 49 23 25' 'set clevs 0 0.2 0.5 1.0 2.0 3.0 5.0 6.0 7.0 9.0 11.0' 'd clim' 'draw title PREC Climatology 1951-90 (mm/day) month= 'mon 'run scripts_cbar.gs' say 'more plot (y/n)?' pull con if (con != n) 'clear' mon = mon + 1 else mon = 13 endif endwhile ***************************************************