reinit
*THIS SCRIPT IS TO DISPLAY 3 MAPS OF SOIL MOISTURE 1m, 0.5m and 0.2m layers
* &0 = s or w 
* &1 = Spring or WINTER
open  grd.&0100.ctl
open  grd.&050.ctl
open  grd.&020.ctl
open mask.ctl
define ocean=mask.4(t=1)
set display color white
c
* &2 = time step number (2-37)
* &3 = January,..., December
* &4 = 8,18 or 28

*enable print pr3.&2.gx
set vpage 0 8.5 6.8 11
set parea 1 8 0.85 3.2
exec setp100
set mpdset hires
set grads off
set lon 20 140
set lat 34 70
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 MULTI-YEAR AVERAGE OF PLANT AVAILABLE SOIL MOISTURE (CM)\&1 cereal crops. &3 &4\Upper 100 cm soil layer
** ** **
set vpage 0 8.5 3.4 11
set parea 1 8 0.85 3.2
set mpdset hires
set grads off
set lon 20 140
set lat 34 70
set t &2
exec setp50
set gxout shaded
set cterp on
d &050.2
set gxout contour
set cmark 0
set cthick 2
set dignum 0
set ccolor 1
set cstyle 1
set cthick 4
set cint 2.5
set cterp on
set csmooth on
set clab forced
d &050.2
exec barp50
*************************************
*maskout ocean points
set gxout fgrid
set fgvals 1 0
d ocean
draw title \\Upper 50 cm soil layer
************
set vpage 0 8.5 0 11
set parea 1 8 0.85 3.2
set mpdset hires
set grads off
set lon 20 140
set lat 34 70
set t &2
exec setp20
set gxout shaded
set cterp on
d &020.3
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.3
exec barp20
*************************************
*maskout ocean points
set gxout fgrid
set fgvals 1 0
d ocean
draw title \\Upper 20 cm soil layer

*print
*disable print
**********************************************************
*exec mpc3 w Winter 2 January 8
*exec mpc3 w Winter 3 January 18
*exec mpc3 w Winter 4 January 28
*exec mpc3 w Winter 5 February 8
*exec mpc3 w Winter 6 February 18
*exec mpc3 w Winter 7 February 28
*exec mpc3 w Winter 8 March 8
*exec mpc3 w Winter 9 March 18
*exec mpc3 w Winter 10 March 28
*exec mpc3 w Winter 11 April 8
*exec mpc3 w Winter 12 April 18
*exec mpc3 w Winter 13 April 28
*exec mpc3 w Winter 14 May 8
*exec mpc3 w Winter 15 May 18
*exec mpc3 w Winter 16 May 28
*exec mpc3 w Winter 17 June 8
*exec mpc3 w Winter 18 June 18
*exec mpc3 w Winter 19 June 28
*exec mpc3 w Winter 20 July 8
*exec mpc3 w Winter 21 July 18
*exec mpc3 w Winter 22 July 28
*exec mpc3 w Winter 23 August 8
*exec mpc3 w Winter 24 August 18
*exec mpc3 w Winter 25 August 28
*exec mpc3 w Winter 26 September 8
*exec mpc3 w Winter 27 September 18
*exec mpc3 w Winter 28 September 28
*exec mpc3 w Winter 29 October 8
*exec mpc3 w Winter 30 October 18
*exec mpc3 w Winter 31 October 28
*exec mpc3 w Winter 32 November 8
*exec mpc3 w Winter 33 November 18
*exec mpc3 w Winter 34 November 28
*exec mpc3 w Winter 35 December 8
*exec mpc3 w Winter 36 December 18
*exec mpc3 w Winter 37 December 28
**********************************************************

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