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

aaa=1
while(aaa<=4)

'reset'

clim=1
while(clim<=1)
zzz=1
while(zzz<=3)
'reset'
'set strsiz 0.27 0.3'
if(aaa= 1) ; varnum='Prec GPCCv6 and GPCCv6       ' ; endif
if(aaa= 2) ; varnum='Prec GPCPv1.2 and GPCCv6     ' ; endif
if(aaa= 3) ; varnum='Prec APHRODITEv110 and GPCCv6' ; endif
if(aaa= 4) ; varnum='Prec GSMAP MWRv484 and GPCCv6' ; endif

if(aaa= 1) ; anum='gpcc' ; endif
if(aaa= 2) ; anum='gpcp' ; endif
if(aaa= 3) ; anum='aphro' ; endif
if(aaa= 4) ; anum='gsmap' ; endif

'draw string 0.6 8.1 'varnum''

if(zzz= 1) ; 'draw string 9.0 8.1 Prec' ; endif
if(zzz= 2) ; 'draw string 9.0 8.1 Evap' ; endif
if(zzz= 3) ; 'draw string 9.0 8.1 Roff' ; endif

if(zzz= 1) ; znum='prec' ; endif
if(zzz= 2) ; znum='evap' ; endif
if(zzz= 3) ; znum='roff' ; endif

'set mpdset hires'
'set gxout shaded'
'define_colors.gs'
'set grads off'
'set lon 90 110'
'set lat  0  35'
'set xlint 5'
'set ylint 5'
'presen.gs -alabel 0.2'
*-----------------------------------
area=1
while(area<=1)


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

if(zzz=1) ; 'colorcb_k 9 c r' ; endif
if(zzz=2) ; 'colorcb_k 9 c r' ; endif
if(zzz=3) ; 'colorcb_k 9 c r' ; endif

'set clevs -500 -250 -100 -50 50 100 250 500'
*if(zzz=1) ; 'set clevs -500 -250 -200 -100 -50 50 100 200 250 500' ; endif
*if(zzz=2) ; 'set clevs -500 -250 -200 -100 -50 50 100 200 250 500' ; endif
*if(zzz=3) ; 'set clevs -500 -250 -200 -100 -50 50 100 200 250 500' ; endif

*if(zzz=2) ; 'set clevs 100 200 300 400 500 600 700 800 900 1000' ; endif
*if(zzz=3) ; 'set clevs 100 200 300 400 500 600 700 800 900 1000' ; endif

t1 = 1
t2 = aaa
if(zzz=1) ; 'd (prec(t='t2')-prec(t='t1'))*86400' ; endif
if(zzz=2) ; 'd (evap(t='t2')-evap(t='t1'))*86400' ; endif
if(zzz=3) ; 'd (roff(t='t2')-roff(t='t1'))*86400' ; endif

endif


'cbar.gs'
area=area+1
endwhile

if(month<10) ; 'printim GIFfile/comp1yr_'znum''anum'.gif' ; endif
if(month> 9) ; 'printim GIFfile/comp1yr_'znum''anum'.gif' ; endif


*'q pos'
zzz=zzz+1
endwhile

*'! convert -loop 1 -delay 30 TEMPfile/im* ANIMfile/anim'clim''varnum'd.gif' 

clim=clim+1
endwhile

aaa=aaa+1
endwhile

return