'reinit'
  'open PRCP_CU_GAUGE_V1.0GLB_0.50deg.lnx.RT.ctl'
  'set dfile 1'
  'set t 1'
  'set z 1'
  'set lat -60 90'
  'set lon -30 330'
*
  'set time 19jan2008' 
  'set rgb 91 220 220 220'
  'run /export/lnx56/mingyue/grads/scripts/schemm'
  'set gxout shaded'
*  'set mpdset mres'
  'set grads off'
  'set grid off'
*
xmin=1.0
xmax=8.0
ymin=1.0
ymax=10.0
b=1
dy=(ymax-ymin-b)/2
*
nt=75
plot=1
while (plot=1)  
ntime=nt
  'set t 'ntime
xstr=8.5/2
ystr=10.8
  'set strsiz 0.20'
  'set string 1 c 6 0'
  'draw string 'xstr' 'ystr' Daily precip [mm/day] Day 'nt' from 20070101'
*
x1=xmin
x2=xmax
y2=ymax
y1=y2-dy
  'set gxout shaded'
  'set parea 'x1' 'x2' 'y1' 'y2
  'set ccols  91    33   35   37 39 42 44  47  22  23  25  27  29  81  82'
  'set clevs  0.01  0.2  0.5  1  2  5  10  15  20  30  40  50  100'
  'd rain.1/10'   
  'run /export/lnx56/mingyue/grads/scripts/cbar.gs' 
  'draw title CPC Unified Analysis'
y2=y1-b
y1=y2-dy
  'set gxout shaded'
  'set parea 'x1' 'x2' 'y1' 'y2
  'set ccols  91    2   3    4   5'
  'set clevs  0     1    2    3 '
  'd gnum.1'   
  'run /export/lnx56/mingyue/grads/scripts/cbar.gs' 
  'draw title stn'
say 'more plot?'
pull con
con=subwrd(con,1)
if (con != n)
nt=nt+1
'clear'
else
plot=0
endif
endwhile
**************************************