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

aaa=7
while(aaa<=7)

'reset'

clim=4
while(clim<=5)
month=13
while(month<=13)
'reset'
'set strsiz 0.27 0.3'
if(aaa= 1) ; varnum='Yield' ; endif
if(aaa= 2) ; varnum='Plant' ; endif
if(aaa= 3) ; varnum='Head' ; endif
if(aaa= 4) ; varnum='Harvest' ; endif
if(aaa= 5) ; varnum='SSSS' ; endif
if(aaa= 6) ; varnum='ghigh' ; endif
if(aaa= 7) ; varnum='glow' ; endif
if(aaa= 8) ; varnum='stress' ; endif

*'draw string 0.6 8.1 'varnum''
if(aaa= 1) ; 'draw string 0.6 8.1 Yield (g/m2)' ; endif
if(aaa= 2) ; 'draw string 0.6 8.1 Plant (doy)' ; endif
if(aaa= 3) ; 'draw string 0.6 8.1 Head (doy)' ; endif
if(aaa= 4) ; 'draw string 0.6 8.1 Harvest (doy)' ; endif
if(aaa= 5) ; 'draw string 0.6 8.1 SSSS (MJ)' ; endif
if(aaa= 6) ; 'draw string 0.6 8.1 ghigh (stress)' ; endif
if(aaa= 7) ; 'draw string 0.6 8.1 glow (stress)' ; endif
if(aaa= 8) ; 'draw string 0.6 8.1 glow (stress)' ; endif

if(clim=1) ; 'draw string 5.0 8.1 Present Climate' ; endif
if(clim=2) ; 'draw string 5.0 8.1 Future1 Climate' ; endif
if(clim=3) ; 'draw string 5.0 8.1 Future2 Climate' ; endif
if(clim=4) ; 'draw string 5.0 8.1 Future1-Present(diff)' ; endif
if(clim=5) ; 'draw string 5.0 8.1 Future2-Present(diff)' ; endif

'set mpdset hires'
'set gxout shaded'
'set gxout grfill'
'define_colors.gs'
'set grads off'

'set parea 0.5 9.0 0.5 8.0'
'set lat 30 46'
'set lon 128 146'

if(clim<=3)
'set ccols 15 14 9 4 5 3 10 7 8 2 6' 
if(aaa= 1) ; 'set ccols 15 4 3 7 8 2' ; endif

*if(aaa= 1) ; 'set clevs 300 350 400 450 500 550 600 650 700 750' ; endif
if(aaa= 1) ; 'set clevs 300 400 500 600 700' ; endif
if(aaa= 2) ; 'set clevs 115 140 165 190 215 240 265 290 315 340' ; endif
if(aaa= 3) ; 'set clevs 115 140 165 190 215 240 265 290 315 340' ; endif
if(aaa= 4) ; 'set clevs 115 140 165 190 215 240 265 290 315 340' ; endif
if(aaa= 6) ; 'set clevs 0   0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9' ; endif
if(aaa= 7) ; 'set clevs 0   0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9' ; endif
'set t 'clim'' 
if(aaa= 1) ; 'd yield' ; endif
if(aaa= 2) ; 'd plant' ; endif
if(aaa= 3) ; 'd head' ; endif
if(aaa= 4) ; 'd harvest' ; endif
if(aaa= 5) ; 'd ssss' ; endif
if(aaa= 6) ; 'd ghigh' ; endif
if(aaa= 7) ; 'd glow' ; endif
'cbar.gs'
endif

if(clim>=4)
*'set ccols 15 14 9 4 5 3 10 7 8 2 6' 
*if(aaa= 1) ; 'set ccols 15 4 3 7 8 2' ; endif

'set gxout grfill'
if(aaa= 1) ; 'color -100 100 1 -kind mediumvioletred->khaki->lime' ; endif
if(aaa= 5) ; 'color -20 20 1 -kind deepskyblue->white->orange' ; endif
if(aaa= 6) ; 'color -0.1 0.1 0.05 -kind deepskyblue->white->orange' ; endif
if(aaa= 7) ; 'color -10 10 1 -kind fuchsia->white->lawngreen' ; endif

ttt1=1
ttt2=clim-2
if(aaa= 1) ; 'd (yield(t='ttt2')-yield(t='ttt1'))' ; endif
if(aaa= 2) ; 'd (plant(t='ttt2')-plant(t='ttt1'))' ; endif
if(aaa= 3) ; 'd (head(t='ttt2')-head(t='ttt1'))' ; endif
if(aaa= 4) ; 'd (harvest(t='ttt2')-harvest(t='ttt1'))' ; endif
if(aaa= 5) ; 'd (ssss(t='ttt2')-ssss(t='ttt1'))' ; endif
if(aaa= 6) ; 'd (ghigh(t='ttt2')*100-ghigh(t='ttt1')*100)' ; endif
if(aaa= 7) ; 'd (glow(t='ttt2')*100-glow(t='ttt1')*100)' ; endif

if(aaa= 1) ; 'xcbar.gs -fstep 25' ; endif
if(aaa= 5) ; 'xcbar.gs -fstep 10' ; endif
if(aaa= 6) ; 'xcbar.gs -fstep 0.2' ; endif
if(aaa= 7) ; 'xcbar.gs -fstep 5' ; endif


'set gxout contour'
'set clevs 0 400'
'set t 1'
'set clab off'
'd yield'
'hatch yield 0 300 -int 0.1 -angle 45,135'
endif



if(month<10) ; 'printim TEMPfile/im'varnum'clim'clim'0'month'.gif' ; endif
if(month> 9) ; 'printim TEMPfile/im'varnum'clim'clim''month'.gif' ; endif

*'q pos'
month=month+1
endwhile

*'! convert -loop 1 -delay 30 TEMPfile/im* ANIMfile/anim'varnum'climall.gif' 
'!mv TEMPfile/im* GIFfile/'

clim=clim+1
endwhile

aaa=aaa+1
endwhile

return