'reinit'
'open /nfsuser/g03/wx52yx/MK_ELLFB/data/SST/SSTanm_clim.ctl'
'open /nfsuser/g03/wx52yx/MK_ELLFB/data/SL/SLanm.ctl'
'open /nfsuser/g03/wx52yx/MK_ELLFB/data/FSU/tauanm.ctl'
*
'enable print fig.meta'
***
'set rgb 28 191 191 191'
'set rgb 29 127   0   0'
'set rgb 30  0   0  255'
'set rgb 31  0  180 255'
'set rgb 32  63 255 255'
'set rgb 33 127 255 127'
'set rgb 40 148 0 211'
'set rgb 41 0 0 128'
*
*pull white
*pull cbartest
*'run 'white''
'run /nfsuser/g03/wx52yx/grads/lib/white.gs'
**
args = read ("fsst.date")
args = read ("fsst.date")
args = read ("fsst.date")
args = read ("fsst.date")
cdate = sublin(args,2)
dat1=subwrd(cdate,1)
cyr1=substr(dat1,1,4)
cmo1=substr(dat1,5,2)
say cdate
say dat1
say cyr1
say cmo1
'set lat 0'
*'set time Jan80 'xue2(cmo1,cyr1)
'set time jan90 jan2000'
*'q dims'
*temp = sublin(result,5)
*end = subwrd(temp,9)
*beg = end - 120
*'set t 'beg' 'end
**
'set vpage off'
'set parea 0.8 3.5  1. 8.'
'set yflip on'
'set xlint 30'
'set grads off '
'set mproj scaled '
'set grid off'
'set gxout shaded'
*'set clevs  -2 -1.5 -1. -0.5 0 0.5 1. 1.5 2.'
'set clevs  -4 -3 -2 -1 0 1 2 3 4'
'set ccols 41 30   31   32   33 7   8   2   29 40'
'd SST'
'set gxout contour'
'set clab off'
'set cthick 2'
'set ccolor 1'
'set cint 1'
'd SST'
'set strsiz 0.15'
'set string 1 tc 6'
'draw string 2.1 8.3 SST (degree)'
*'run 'cbartest ' 1.3 2.4 0.5 0.13 0'
'run /nfsuser/g03/wx52yx/grads/lib/cbartest.gs 0.5 0.45 0.3 0.13 0'
****
'set vpage off'
'set parea 4.1 6.8  1. 8.'
'set yflip on'
'set xlint 30'
'set strsiz 0.15'
'set string 1 tc 6'
'draw string 5.4 8.3 SL (cm)'
'set grads off'
'set gxout shaded '
'set clevs   -24 -18 -12 -6 0 6 12 18 24'
'set ccols 41 30   31   32   33    7    8   2   29 40'
'd SL.2'
'set gxout contour'
'set clab off'
'set cthick 2'
'set ccolor 1'
'set cint 6'
'd SL.2'
'run /nfsuser/g03/wx52yx/grads/lib/cbartest.gs 3.8 0.45 0.4 0.13 0'
***
'set vpage off'
'set parea 7.4 10.1  1. 8.'
'set yflip on'
'set xlint 30'
'set grads off '
'set grads off'
'set gxout shaded '
*'set clevs   -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8'
'set clevs   -0.6 -0.45 -0.3 -0.15 0 0.15 0.3 0.45 0.6'
'set ccols 41 30   31   32   33    7    8   2   29 40'
*'define var=tloop(ave(taux.3,t-1,t+1))'
'd taux.3'
'set gxout contour'
'set clab off'
'set cthick 2'
'set ccolor 1'
'set cint 0.15'
'd taux.3'
'set strsiz 0.15'
'set string 1 tc 6'
'draw string 8.7 8.3 Zonal Stress (dyn/cm`a2`n)'
'run /nfsuser/g03/wx52yx/grads/lib/cbartest.gs 10.2 2. 0.11 0.5 90'
**
'print '
'disable print'
*'quit'
'! gxeps -c -i fig.meta -o ../eps/sst_sl_wind_90-00.eps'
'! gxgif -i fig.meta -rotate 90 -r -x 825 -y 638 -o ../../pub/web/sst_sl_wind_90-00.gif'

function xue(cmo1,cyr1) 
  if(cmo1='01'); cmo1='January'; endif
  if(cmo1='02'); cmo1='February'; endif
  if(cmo1='03'); cmo1='March'; endif
  if(cmo1='04'); cmo1='April'; endif
  if(cmo1='05'); cmo1='May'; endif
  if(cmo1='06'); cmo1='June'; endif
  if(cmo1='07'); cmo1='July'; endif
  if(cmo1='08'); cmo1='August'; endif
  if(cmo1='09'); cmo1='September'; endif
  if(cmo1='10'); cmo1='October'; endif
  if(cmo1='11'); cmo1='November'; endif
  if(cmo1='12'); cmo1='December'; endif
return (cmo1' 'cyr1)
function xue2(cmo1,cyr1) 
  if(cmo1='01'); cmo1='Jan'; endif
  if(cmo1='02'); cmo1='Feb'; endif
  if(cmo1='03'); cmo1='Mar'; endif
  if(cmo1='04'); cmo1='Apr'; endif
  if(cmo1='05'); cmo1='May'; endif
  if(cmo1='06'); cmo1='Jun'; endif
  if(cmo1='07'); cmo1='Jul'; endif
  if(cmo1='08'); cmo1='Aug'; endif
  if(cmo1='09'); cmo1='Sep'; endif
  if(cmo1='10'); cmo1='Oct'; endif
  if(cmo1='11'); cmo1='Nov'; endif
  if(cmo1='12'); cmo1='Dec'; endif
return (cmo1''cyr1)