'reinit'
'open wstress_project.ctl'
'set display color white'
'reset'

aaa=3
while(aaa<=3)

clim=8
while(clim<=17)

'reset'
'set strsiz 0.27 0.3'
if(aaa= 1) ; varnum='Resources' ; endif
if(aaa= 2) ; varnum='Demand' ; endif
if(aaa= 3) ; varnum='Stress' ; endif

'draw string 0.6 8.1 'varnum''

if(aaa=1) ; 'draw string 3.0 8.1 (mm)' ; endif
if(aaa=2) ; 'draw string 3.0 8.1 (mm)' ; endif
if(aaa=3) ; 'draw string 3.0 8.1 (CWD)' ; endif

if(clim=1) ; 'draw string 7.4 8.1 PRESENT' ; endif
if(clim=2) ; 'draw string 7.4 8.1 FUTURE1' ; endif
if(clim=3) ; 'draw string 7.4 8.1 FUTURE2' ; endif
if(clim=4) ; 'draw string 7.4 8.1 ADVANCE' ; endif
if(clim=5) ; 'draw string 7.4 8.1 DELAYED' ; endif
if(clim=6) ; 'draw string 7.4 8.1 FORHEAD' ; endif
if(clim=7) ; 'draw string 7.4 8.1 WAITBIT' ; endif

if(clim=8) ; 'draw string 5.0 8.1 diff FUTURE1-PRESENT' ; endif
if(clim=9) ; 'draw string 5.0 8.1 diff FUTURE2-PRESENT' ; endif
if(clim=10) ; 'draw string 5.0 8.1 diff ADVANCE-PRESENT' ; endif
if(clim=11) ; 'draw string 5.0 8.1 diff DELAYED-PRESENT' ; endif
if(clim=12) ; 'draw string 5.0 8.1 diff FORHEAD-PRESENT' ; endif
if(clim=13) ; 'draw string 5.0 8.1 diff WAITBIT-PRESENT' ; endif

if(clim=14) ; 'draw string 5.0 8.1 diff ADVANCE-FUTURE2' ; endif
if(clim=15) ; 'draw string 5.0 8.1 diff DELAYED-FUTURE2' ; endif
if(clim=16) ; 'draw string 5.0 8.1 diff FORHEAD-FUTURE2' ; endif
if(clim=17) ; 'draw string 5.0 8.1 diff WAITBIT-FUTURE2' ; endif


*'draw string 8.6 8.1 GDBD' 


'set grads off'

area=1
while(area<=2)
'set mpdraw off'
'set gxout shaded'
if(area=1)
'set parea 0.5 9.0 0.5 8.0'
'set lat 30 46'
'set lon 128 146'
'set xlopts 1 0 0'
'set ylopts 1 0 0'
endif
if(area=2)
'set parea 6.07 8.9 0.5 2.84'
'set xlab off'
'set ylab off'
'set lat 24 29'
'set lon 126 130'
endif

'set ccols 0 15 14 9 4 5 3 10 7 8 2 6' 
* ------setting color bar--------
*light yellow to dark red
'set rgb 21 255 250 170'
'set rgb 22 255 232 120'
'set rgb 23 255 192  60'
'set rgb 24 255 160   0'
'set rgb 25 255  96   0'
'set rgb 26 255  50   0'
'set rgb 27 225  20   0'
'set rgb 28 192   0   0'
'set rgb 29 165   0   0'
*light blue to dark blue
'set rgb 41 225 255 255'
'set rgb 42 180 240 250'
'set rgb 43 150 210 250'
'set rgb 44 120 185 250'
'set rgb 45  80 165 245'
'set rgb 46  60 150 245'
'set rgb 47  40 130 240'
'set rgb 48  30 110 235'
'set rgb 49  20 100 210'
* ------setting color bar--------
if(clim<=7)
if(aaa=1) ; 'set ccols 20 41 42 43 44 45 46 47 48 49' ; endif
if(aaa=2) ; 'set ccols 20 21 22 23 24 25 26 27 28 29' ; endif
if(aaa=3) ; 'set ccols 20 29 28 27 26 25 24 23 22 21' ; endif

if(aaa=1) ; 'set clevs 0 10 25 50 100 250 500 1000 2000' ; endif
if(aaa=2) ; 'set clevs 0 10 25 50 100 250 500 1000 2000' ; endif
if(aaa=3) ; 'set clevs 0. 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9' ; endif

ttt=clim
'set t 'ttt''
if(aaa=1) ; 'd disc' ; endif
if(aaa=2) ; 'd demand' ; endif
if(aaa=3) ; 'd stress' ; endif
endif

if(clim>=8)

if(clim<=13)
'set ccols 47 45 43 41 21 23 25 27 29'
if(aaa=1) ; 'set ccols 29 27 25 23 21 41 43 45 47' ; endif
if(aaa=2) ; 'set ccols 47 45 43 41 21 23 25 27 29' ; endif
***if(aaa=3) ; 'set ccols 29 27 25 23 21 41 43 45 47' ; endif


if(aaa=1) ; 'set clevs -500 -250 -100 -50 50 100 250 500' ; endif
if(aaa=2) ; 'set clevs -500 -250 -100 -50 50 100 250 500' ; endif
***if(aaa=3) ; 'set clevs -0.3 -0.1 -0.05 -0.02 0.02 0.05 0.1 0.3' ; endif

if(aaa=3) ; 'color -10 10 1 -kind fuchsia->gold->deepskyblue' ;endif

ttt1=1
ttt2=clim-6
if(aaa=1) ; 'd disc(t='ttt2')-disc(t='ttt1')' ; endif
if(aaa=2) ; 'd demand(t='ttt2')-demand(t='ttt1')' ; endif
if(aaa=3) ; 'd stress(t='ttt2')*100-stress(t='ttt1')*100' ; endif
endif


if(clim>=14)

if(aaa=1) ; 'set ccols 29 27 25 23 21 41 43 45 47' ; endif
if(aaa=2) ; 'set ccols 47 45 43 41 21 23 25 27 29' ; endif
***if(aaa=3) ; 'set ccols 29 27 25 23 21 41 43 45 47' ; endif

if(aaa=1) ; 'set clevs -500 -250 -100 -50 50 100 250 500' ; endif
if(aaa=2) ; 'set clevs -500 -250 -100 -50 50 100 250 500' ; endif
***if(aaa=3) ; 'set clevs -0.3 -0.1 -0.05 -0.02 0.02 0.05 0.1 0.3' ; endif
if(aaa=3) ; 'color -10 10 1 -kind fuchsia->gold->deepskyblue' ;endif

ttt1=3
ttt2=clim-10
if(aaa=1) ; 'd disc(t='ttt2')-disc(t='ttt1')' ; endif
if(aaa=2) ; 'd demand(t='ttt2')-demand(t='ttt1')' ; endif
if(aaa=3) ; 'd stress(t='ttt2')*100-stress(t='ttt1')*100' ; endif
endif

endif

if(aaa=1) ; 'cbar.gs' ; endif
if(aaa=2) ; 'cbar.gs' ; endif
if(aaa=3) ; 'xcbar.gs -fstep 5' ; endif

area=area+1
endwhile

if(clim<10)
'printim GIFfile/sts'varnum'0'clim'.gif' 
endif
if(clim>9)
'printim GIFfile/sts'varnum''clim'.gif' 
endif

*'q pos'

clim=clim+1
endwhile

aaa=aaa+1
endwhile

return