'reinit'
'open /export-8/cacsrv1/cpc/anal/cdas/daily/flux/y1979_cur/PRESmsl.daily.79_cur.ctl'
'open /export-8/cacsrv1/cpc/anal/cdas/daily/flux/clim/PRESmsl.daily.clim.ctl'
'set vpage 0 8.5 0 11.0'
'set parea 1.0 4.0 8.5 10.0'
'set grads off'
'set map 1 1 8'
'set looping off' 
'set csmooth on'
'run /export/lnx286/wd52vk/grads/util/rgbset.gs'
'set clevs -16 -14 -12 -10 -8 -6 -4 -2 0 2 4 6 8 10 12 14 16'
'set ccols 29 28 27 26 25 24 23 22 21 41 42 43 44 45 46 47 48 49'
'set string 1 tc 7'
'set strsiz 0.2 0.2'
'draw string 4.25 10.5 SLP Anomaly t=1 July 1980'
'set strsiz 0.1 0.1'
'draw string 2.5 10.15 t-5'
'set lat -60 0'
'set lon 210 340'
'set time 1jul1980'
cdate=mydate()
cdy1=subwrd(cdate,1)
cmo1=subwrd(cdate,2)
cyr1=subwrd(cdate,3)
'set time 'cdy1''cmo1''cyr1
'define u1=0.01*presmsl.1(t-5)'
'define ua1=0.01*presmsl.1(t-5)-0.01*pres.2(t-5)'
'set gxout shaded'
'set ylint 10'
'set xlab off'
'd ua1'
'set gxout contour'
'set cthick 2'
'set ccolor 1'
'set cint 2'
'set clab off'
'd ua1'
'set vpage 0 8.5 0 11.0'
'set parea 1.0 4.0 6.7 8.2'
'set grads off'
'set map 1 1 8'
'set looping off' 
'set csmooth on'
'set string 1 tc 7'
'set strsiz 0.1 0.1'
'draw string 2.5 8.35 t-4'
'set time 'cdy1''cmo1''cyr1
'define u1=0.01*presmsl.1(t-4)'
'define ua1=0.01*presmsl.1(t-4)-0.01*pres.2(t-4)'
'set clevs -16 -14 -12 -10 -8 -6 -4 -2 0 2 4 6 8 10 12 14 16'
'set ccols 29 28 27 26 25 24 23 22 21 41 42 43 44 45 46 47 48 49'
'set gxout shaded'
'set ylint 10'
'set xlab off'
'd ua1'
'set gxout contour'
'set cthick 2'
'set ccolor 1'
'set cint 2'
'set clab off'
'd ua1'
'set vpage 0 8.5 0 11.0'
'set parea 1.0 4.0 4.9 6.4'
'set grads off'
'set map 1 1 8'
'set looping off' 
'set csmooth on'
'set string 1 tc 7'
'set strsiz 0.1 0.1'
'draw string 2.5 6.55 t-3'
'set time 'cdy1''cmo1''cyr1
'define u1=0.01*presmsl.1(t-3)'
'define ua1=0.01*presmsl.1(t-3)-0.01*pres.2(t-3)'
'set clevs -16 -14 -12 -10 -8 -6 -4 -2 0 2 4 6 8 10 12 14 16'
'set ccols 29 28 27 26 25 24 23 22 21 41 42 43 44 45 46 47 48 49'
'set gxout shaded'
'set ylint 10'
'set xlab off'
'd ua1'
'set gxout contour'
'set cthick 2'
'set ccolor 1'
'set cint 2'
'set clab off'
'd ua1'
'set vpage 0 8.5 0 11.0'
'set parea 1.0 4.0 3.1 4.6'
'set grads off'
'set map 1 1 8'
'set looping off' 
'set csmooth on'
'set string 1 tc 7'
'set strsiz 0.1 0.1'
'draw string 2.5 4.75 t-2'
'set time 'cdy1''cmo1''cyr1
'define u1=0.01*presmsl.1(t-2)'
'define ua1=0.01*presmsl.1(t-2)-0.01*pres.2(t-2)'
'set clevs -16 -14 -12 -10 -8 -6 -4 -2 0 2 4 6 8 10 12 14 16'
'set ccols 29 28 27 26 25 24 23 22 21 41 42 43 44 45 46 47 48 49'
'set gxout shaded'
'set ylint 10'
'set xlab off'
'd ua1'
'set gxout contour'
'set cthick 2'
'set ccolor 1'
'set cint 2'
'set clab off'
'd ua1'
'set vpage 0 8.5 0 11.0'
'set parea 1.0 4.0 1.3 2.8'
'set grads off'
'set map 1 1 8'
'set looping off' 
'set csmooth on'
'set string 1 tc 7'
'set strsiz 0.1 0.1'
'draw string 2.5 2.95 t-1'
'set time 'cdy1''cmo1''cyr1
'define u1=0.01*presmsl.1(t-1)'
'define ua1=0.01*presmsl.1(t-1)-0.01*pres.2(t-1)'
'set clevs -16 -14 -12 -10 -8 -6 -4 -2 0 2 4 6 8 10 12 14 16'
'set ccols 29 28 27 26 25 24 23 22 21 41 42 43 44 45 46 47 48 49'
'set gxout shaded'
'set ylint 10'
'set xlint 20'
'set xlab on'
'd ua1'
'set gxout contour'
'set cthick 2'
'set ccolor 1'
'set cint 2'
'set clab off'
'd ua1'
'set vpage 0 8.5 0 11.0'
'set parea 4.5 7.5 8.5 10.0'
'set grads off'
'set map 1 1 8'
'set looping off' 
'set csmooth on'
'set string 1 tc 7'
'set strsiz 0.1 0.1'
'draw string 6.0 10.15 t'
'set time 'cdy1''cmo1''cyr1
'define u1=0.01*presmsl.1'
'define ua1=0.01*presmsl.1-0.01*pres.2'
'set clevs -16 -14 -12 -10 -8 -6 -4 -2 0 2 4 6 8 10 12 14 16'
'set ccols 29 28 27 26 25 24 23 22 21 41 42 43 44 45 46 47 48 49'
'set gxout shaded'
'set ylab off'
'set xlab off'
'd ua1'
'set gxout contour'
'set cthick 2'
'set ccolor 1'
'set cint 2'
'set clab off'
'd ua1'
'set vpage 0 8.5 0 11.0'
'set parea 4.5 7.5 6.7 8.2'
'set grads off'
'set map 1 1 8'
'set looping off' 
'set csmooth on'
'set string 1 tc 7'
'set strsiz 0.1 0.1'
'draw string 6.0 8.35 t+1'
'set time 'cdy1''cmo1''cyr1
'define u1=0.01*presmsl.1(t+1)'
'define ua1=0.01*presmsl.1(t+1)-0.01*pres.2(t+1)'
'set clevs -16 -14 -12 -10 -8 -6 -4 -2 0 2 4 6 8 10 12 14 16'
'set ccols 29 28 27 26 25 24 23 22 21 41 42 43 44 45 46 47 48 49'
'set gxout shaded'
'set ylab off'
'set xlab off'
'd ua1'
'set gxout contour'
'set cthick 2'
'set ccolor 1'
'set cint 2'
'set clab off'
'd ua1'
'set vpage 0 8.5 0 11.0'
'set parea 4.5 7.5 4.9 6.4'
'set grads off'
'set map 1 1 8'
'set looping off' 
'set csmooth on'
'set string 1 tc 7'
'set strsiz 0.1 0.1'
'draw string 6.0 6.55 t+2'
'set time 'cdy1''cmo1''cyr1
'define u1=0.01*presmsl.1(t+2)'
'define ua1=0.01*presmsl.1(t+2)-0.01*pres.2(t+2)'
'set clevs -16 -14 -12 -10 -8 -6 -4 -2 0 2 4 6 8 10 12 14 16'
'set ccols 29 28 27 26 25 24 23 22 21 41 42 43 44 45 46 47 48 49'
'set gxout shaded'
'set ylab off'
'set xlab off'
'd ua1'
'set gxout contour'
'set cthick 2'
'set ccolor 1'
'set cint 2'
'set clab off'
'd ua1'
'set vpage 0 8.5 0 11.0'
'set parea 4.5 7.5 3.1 4.6'
'set grads off'
'set map 1 1 8'
'set looping off' 
'set csmooth on'
'set string 1 tc 7'
'set strsiz 0.1 0.1'
'draw string 6.0 4.75 t+3'
'set time 'cdy1''cmo1''cyr1
'define u1=0.01*presmsl.1(t+3)'
'define ua1=0.01*presmsl.1(t+3)-0.01*pres.2(t+3)'
'set clevs -16 -14 -12 -10 -8 -6 -4 -2 0 2 4 6 8 10 12 14 16'
'set ccols 29 28 27 26 25 24 23 22 21 41 42 43 44 45 46 47 48 49'
'set gxout shaded'
'set ylab off'
'set xlab off'
'd ua1'
'set gxout contour'
'set cthick 2'
'set ccolor 1'
'set cint 2'
'set clab off'
'd ua1'
'set vpage 0 8.5 0 11.0'
'set parea 4.5 7.5 1.3 2.8'
'set grads off'
'set map 1 1 8'
'set looping off' 
'set csmooth on'
'set string 1 tc 7'
'set strsiz 0.1 0.1'
'draw string 6.0 2.95 t+4'
'set time 'cdy1''cmo1''cyr1
'define u1=0.01*presmsl.1(t+4)'
'define ua1=0.01*presmsl.1(t+4)-0.01*pres.2(t+4)'
'set clevs -16 -14 -12 -10 -8 -6 -4 -2 0 2 4 6 8 10 12 14 16'
'set ccols 29 28 27 26 25 24 23 22 21 41 42 43 44 45 46 47 48 49'
'set gxout shaded'
'set ylab off'
'set xlab on'
'set xlint 20'
'd ua1'
'set gxout contour'
'set cthick 2'
'set ccolor 1'
'set cint 2'
'set clab off'
'd ua1'
'run /export/lnx286/wd52vk/grads/util/cbarmerc.gs'
'printim slp-lead-lag-cf-26jun-5jul1980.gif white'
function mydate
  'query time'
  sres = subwrd(result,3)
  i = 1
  while (substr(sres,i,1)!='Z')
    i = i + 1
  endwhile
  hour = substr(sres,1,i)
  isav = i
  i = i + 1
  while (substr(sres,i,1)>='0' & substr(sres,i,1)<='9')
    i = i + 1
  endwhile
  day = substr(sres,isav+1,i-isav-1)
  month = substr(sres,i,3)
  year = substr(sres,i+3,4)
  return (day' 'month' 'year)