reinit
*THIS SCRIPT IS TO DISPLAY 3 MAPS OF SOIL MOISTURE 1m, 0.2m and ratio 0.2m/1m
* &0 = s or w 
* &1 = SPRING or WINTER
open  grd.&0100.ctl
open  grd.&020.ctl
open  grd.&0r.ctl
open mask.ctl
define ocean=mask.4(t=1)
set display color white
c
* &2 = time step number (2-73)
* &3 = January,..., December
* &4 = 8,18 or 28
* &5 = year 87 or 88
*enable print pr3.&2.gx
set vpage 0 8.5 7 11
set parea 1 8 0.6 3
exec setp100
set mpdset hires
set grads off
set lon 20 120
set lat 40 62
set t &2 
set gxout shaded
set cterp on
d &0100
set gxout contour
set cmark 0
set cthick 2
set dignum 0
set ccolor 1
set cstyle 1
set cthick 4
set cint 5
set cterp on
set csmooth on
set clab forced
d &0100
exec barp100
*************************************
*maskout ocean points
set gxout fgrid
set fgvals 1 0
d ocean
draw title PLANT AVAILABLE SOIL MOISTURE (CM)\&1 cereal crops. &3 &4, 19&5\Upper 100 cm soil layer
** ** **
set vpage 0 8.5 3.75 11
set parea 1 8 0.6 3
set mpdset hires
set grads off
set lon 20 120
set lat 40 62
set t &2
exec setp20
set gxout shaded
set cterp on
d &020.2
set gxout contour
set cmark 0
set cthick 2
set dignum 0
set ccolor 1
set cstyle 1
set cthick 4
set cint 1
set cterp on
set csmooth on
set clab forced
d &020.2
exec barp20
*************************************
*maskout ocean points
set gxout fgrid
set fgvals 1 0
d ocean
draw title \\Upper 20 cm soil layer
set vpage 0 8.5 0.25 11
set parea 1 8 0.6 3
set mpdset hires
set grads off
set lon 20 120
set lat 40 62
set t &2
exec setpr
set gxout shaded
set cterp on
d (100*&0r.3)
set gxout contour
set cmark 0
set cthick 2
set dignum 0
set ccolor 1
set cstyle 1
set cthick 4
set cint 5
set cterp on
set csmooth on
set clab forced
d (100*&0r.3)
exec barpr
*************************************
*maskout ocean points
set gxout fgrid
set fgvals 1 0
d ocean
draw title \Ratio of plant available water content (%)\Upper 20 cm and 100 cm layers
*print
*disable print
**********************************************************


*exec mpc3 w Winter 10 March 28 87
*exec mpc3 w Winter 11 April 8 87
*exec mpc3 w Winter 12 April 18 87
*exec mpc3 w Winter 13 April 28 87
*exec mpc3 w Winter 14 May 8 87
*exec mpc3 w Winter 15 May 18 87
*exec mpc3 w Winter 16 May 28 87
*exec mpc3 w Winter 17 June 8 87
*exec mpc3 w Winter 18 June 18 87
*exec mpc3 w Winter 19 June 28 87
*exec mpc3 w Winter 20 July 8 87
*exec mpc3 w Winter 21 July 18 87
*exec mpc3 w Winter 22 July 28 87
*exec mpc3 w Winter 23 August 8 87
*exec mpc3 w Winter 24 August 18 87
*exec mpc3 w Winter 25 August 28 87
*exec mpc3 w Winter 26 September 8 87
*exec mpc3 w Winter 27 September 18 87
*exec mpc3 w Winter 28 September 28 87
*exec mpc3 w Winter 29 October 8 87
*exec mpc3 w Winter 30 October 18 87
*exec mpc3 w Winter 31 October 28 87

*exec mpc3 w Winter 47 April 8 88
*exec mpc3 w Winter 48 April 18 88
*exec mpc3 w Winter 49 April 28 88
*exec mpc3 w Winter 50 May 8 88
*exec mpc3 w Winter 51 May 18 88
*exec mpc3 w Winter 52 May 28 88
*exec mpc3 w Winter 52 June 8 88
*exec mpc3 w Winter 54 June 18 88
*exec mpc3 w Winter 55 June 28 88
*exec mpc3 w Winter 56 July 8 88
*exec mpc3 w Winter 57 July 18 88
*exec mpc3 w Winter 58 July 28 87
*exec mpc3 w Winter 59 August 8 88
*exec mpc3 w Winter 60 August 18 88
*exec mpc3 w Winter 61 August 28 88
*exec mpc3 w Winter 62 September 8 88
*exec mpc3 w Winter 63 September 18 88
*exec mpc3 w Winter 64 September 28 88
*exec mpc3 w Winter 65 October 8 88
*exec mpc3 w Winter 66 October 18 88
*exec mpc3 w Winter 67 October 28 88
**********************************************************


*exec mpc3 s Summer 11 April 8 87
*exec mpc3 s Summer 12 April 18 87
*exec mpc3 s Summer 13 April 28 87
*exec mpc3 s Summer 14 May 8 87
*exec mpc3 s Summer 15 May 18 87
*exec mpc3 s Summer 16 May 28 87
*exec mpc3 s Summer 17 June 8 87
*exec mpc3 s Summer 18 June 18 87
*exec mpc3 s Summer 19 June 28 87
*exec mpc3 s Summer 20 July 8 87
*exec mpc3 s Summer 21 July 18 87
*exec mpc3 s Summer 22 July 28 87
*exec mpc3 s Summer 23 August 8 87
*exec mpc3 s Summer 24 August 18 87
*exec mpc3 s Summer 25 August 28 87
*exec mpc3 s Summer 26 September 8 87
*exec mpc3 s Summer 27 September 18 87
*exec mpc3 s Summer 28 September 28 87


*exec mpc3 s Summer 47 April 8 88
*exec mpc3 s Summer 48 April 18 88
*exec mpc3 s Summer 49 April 28 88
*exec mpc3 s Summer 50 May 8 88
*exec mpc3 s Summer 51 May 18 88
*exec mpc3 s Summer 52 May 28 88
*exec mpc3 s Summer 53 June 8 88
*exec mpc3 s Summer 54 June 18 88
*exec mpc3 s Summer 55 June 28 88
*exec mpc3 s Summer 56 July 8 88
*exec mpc3 s Summer 57 July 18 88
*exec mpc3 s Summer 58 July 28 88
*exec mpc3 s Summer 59 August 8 88
*exec mpc3 s Summer 60 August 18 88
*exec mpc3 s Summer 61 August 28 88
*exec mpc3 s Summer 62 September 8 88
*exec mpc3 s Summer 63 September 18 88
*exec mpc3 s Summer 64 September 28 88