Plotting

::::::::::::::

mld2.ctl

::::::::::::::

dset dat/mld.%y4%m2.bin

undef -999.9

options template

xdef 101 linear 24.5 1.0

ydef 126 levels

-30.5 -29.5 -28.5 -27.5 -26.5 -25.5 -24.502 -23.513 -22.537

-21.581 -20.651 -19.75 -18.883 -18.054 -17.266 -16.522 -15.823

-15.169 -14.562 -14 -13.482 -13.005 -12.566 -12.163 -11.79

-11.442 -11.115 -10.803 -10.5 -10.2 -9.9 -9.6 -9.3 -9 -8.7

-8.4 -8.1 -7.8 -7.5 -7.2 -6.9 -6.6 -6.3 -6 -5.7 -5.4 -5.1

-4.8 -4.5 -4.2 -3.9 -3.6 -3.3 -3 -2.7 -2.4 -2.1 -1.8 -1.5

-1.2 -0.9 -0.6 -0.3 0 0.3 0.6 0.9 1.2 1.5 1.8 2.1 2.4 2.7

3 3.3 3.6 3.9 4.2 4.5 4.8 5.1 5.4 5.7 6 6.3 6.6 6.9 7.2

7.5 7.8 8.1 8.4 8.7 9 9.3 9.6 9.9 10.2 10.5 10.803 11.115

11.442 11.789 12.163 12.566 13.005 13.482 14 14.562 15.169

15.823 16.522 17.266 18.054 18.883 19.75 20.651 21.581 22.537

23.513 24.502 25.5 26.5 27.5 28.5 29.5

zdef 21 levels

5 15 25 35 45 55 65 75 85 95 105 115 125 135 145

155.565 168.295 185.205 208.675 239.915 279.91

tdef 264 linear 01jan1993 1mo

vars 51

tm 0 99 t,y,x wind speed

sm 0 99 t,y,x wind speed

mld 0 99 t,y,x wind speed

ild 0 99 t,y,x wind speed

blt 0 99 t,y,x wind speed

d20 0 99 t,y,x wind speed

um 0 99 t,y,x wind speed

vm 0 99 t,y,x wind speed

uadv1 0 99 t,y,x adv

vadv1 0 99 t,y,x adv

wadv1 0 99 t,y,x adv

xdif1 0 99 t,y,x adv

ydif1 0 99 t,y,x adv

zdif1 0 99 t,y,x adv

gmdif1 0 99 t,y,x adv

kppmix1 0 99 t,y,x adv

kf_cor1 0 99 t,y,x adv

sflux1 0 99 t,y,x adv

sdamp1 0 99 t,y,x adv

total1 0 99 t,y,x adv

dtdt 0 99 t,y,x adv

ent1 0 99 t,y,x adv

uadv2 0 99 t,y,x adv

vadv2 0 99 t,y,x adv

wadv2 0 99 t,y,x adv

xdif2 0 99 t,y,x adv

ydif2 0 99 t,y,x adv

zdif2 0 99 t,y,x adv

gmdif2 0 99 t,y,x adv

kppmix2 0 99 t,y,x adv

kf_cor2 0 99 t,y,x adv

sflux2 0 99 t,y,x adv

sdamp2 0 99 t,y,x adv

total2 0 99 t,y,x adv

dsdt 0 99 t,y,x adv

ent2 0 99 t,y,x adv

advx1 0 99 t,y,x adv

advy1 0 99 t,y,x adv

advz1 0 99 t,y,x adv

subx1 0 99 t,y,x adv

suby1 0 99 t,y,x adv

we1 0 99 t,y,x adv

advx2 0 99 t,y,x adv

advy2 0 99 t,y,x adv

advz2 0 99 t,y,x adv

subx2 0 99 t,y,x adv

suby2 0 99 t,y,x adv

we2 0 99 t,y,x adv

t 21 99 t,y,x wind speed

s 21 99 t,y,x wind speed

r 21 99 t,y,x wind speed

endvars

::::::::::::::

map.blt.02.gs

::::::::::::::

*

'open mld2.ctl'

*

'set t 1 12'

'define bltclm = ave(blt, t+0, t=264, 12)'

'modify bltclm seasonal'

*

var=BLT

title.0=var' (ECCO)'

title.1='(a) Jan'

title.2='(b) Mar'

title.3='(c) May'

title.4='(d) Jul'

title.5='(e) Sep'

title.6='(f) Nov'

'q dims'

say result

*

xmax = 2

ymax = 3

xwid = 5.0/xmax ;*9.0/xmax

ywid = 7.5/ymax ;*6.0/ymax

say 'xwid ' xwid

say 'ywid ' ywid

say

'cc'

nmap = 1

ymap = 1

while (ymap <= ymax)

xmap = 1

while (xmap <= xmax)

*

xs = 1.0 + (xwid+0.75)*(xmap-1)

xe = xs + xwid

#ye = 7.5 - (ywid+0.75)*(ymap-1)

ye = 10.0 - (ywid+0.50)*(ymap-1)

ys = ye - ywid

say xs ' ' xe

say ys ' ' ye

say

month = (nmap-1)*2 + 1

*

if (ymap = ymax)

'set xlopts 1 4 0.12'

else

'set xlopts 1 4 0.00'

'set xlopts 1 4 0.12'

endif

if (xmap = 1)

'set ylopts 1 4 0.12'

else

'set ylopts 1 4 0.00'

'set ylopts 1 4 0.12'

endif

'set vpage 0.0 11.0 0.0 8.5'

'set parea 'xs' 'xe' 'ys' 'ye

'set grads off'

'set font 4'

*

'set xlint 20'

'set ylint 10'

*

'set t 'month

'set gxout shaded'

'set clevs 0 5 10 15 20 25 30 35 40 45'

'd bltclm'

#'d maskout(bltclm,bltclm-5)'

*

xt = xs+(xe-xs)/2.0 ;* xs+0.2

yt = ye

'set string 1 c 4 0'

'set strsiz 0.15 0.15'

'draw string 'xt' 'yt' 'title.nmap

*

if (xmap = xmax & ymap = 1)

'run xcbar.gs 'xe+0.2' 'xe+0.3' 'ys' 'ye' -direction vertical -line on'

endif

*

'set parea off'

'set vpage off'

*

nmap = nmap + 1

xmap = xmap + 1

endwhile

ymap = ymap + 1

endwhile

xt = 3.8

yt = 10.5

'set string 1 c 4 0'

'set strsiz 0.25 0.25'

'draw string 'xt' 'yt' 'title.0

note='Stepwise_MLD'

xt = 3.8; yt = 0.5

'set string 1 c 4 0'

'set strsiz 0.2 0.2'

'draw string 'xt' 'yt' 'note

ofile=var'.ECCO.dr080.bimonthly.'note'.eps'

'gxprint 'ofile

'allclose'

::::::::::::::

map.hb2

::::::::::::::

*

'open mld2.ctl'

*

'set t 1 12'

'define mldclm = ave(mld, t+0, t=264, 12)'

'modify mldclm seasonal'

'define bltclm = ave(blt, t+0, t=264, 12)'

'modify bltclm seasonal'

'define ildclm = ave(ild, t+0, t=264, 12)'

'modify ildclm seasonal'

*'define d20clm = ave(d20, t+0, t=264, 12)'

*'modify d20clm seasonal'

'define dtdtc = ave(dtdt, t+0, t=264, 12)'

'modify dtdtc seasonal'

'define totalc = ave(total1, t+0, t=264, 12)'

'modify totalc seasonal'

'define advc = ave(uadv1+vadv1+wadv1, t+0, t=264, 12)'

'modify advc seasonal'

'define difc = ave(xdif1+ydif1+zdif1+gmdif1+kppmix1, t+0, t=264, 12)'

'modify difc seasonal'

'define fluxc = ave(sflux1+sdamp1, t+0, t=264, 12)'

'modify fluxc seasonal'

'define tmc = ave(tm, t+0, t=264, 12)'

'modify tmc seasonal'

'define adv2c = ave(advx1+advy1+advz1, t+0, t=264, 12)'

'modify adv2c seasonal'

'define entc = ave(we1, t+0, t=264, 12)'

'modify entc seasonal'

*

'set lon 95 '

'set lat -2'

'set t 1 12'

*

title.1='(a)BLT,MLD,ILD'

title.2='(b)dTm/dt,total'

title.3='(c)adv,dif,flux+sdamp'

title.4='(d)Tm'

*

xmax = 1

ymax = 2

xwid = 8.5/xmax

ywid = 5.5/ymax

*

nmap = 1

ymap = 1

while (ymap <= ymax)

xmap = 1

while (xmap <= xmax)

*

xs = 1.0 + (xwid+0.80)*(xmap-1)

xe = xs + xwid

ye = 7.5 - (ywid+0.75)*(ymap-1)

*ye = 10.5 - (ywid+0.25)*(ymap-1)

ys = ye - ywid

*

month = (nmap-1)*3 + 1

*

if (ymap = ymax)

'set xlopts 1 4 0.15'

else

'set xlopts 1 4 0.00'

'set xlopts 1 4 0.15'

endif

if (xmap = 1)

'set ylopts 1 4 0.15'

else

'set ylopts 1 4 0.00'

'set ylopts 1 4 0.15'

endif

'set vpage 0.0 11.0 0.0 8.5'

'set parea 'xs' 'xe' 'ys' 'ye

'set grads off'

'set font 4'

*

'set xlint 10'

'set ylint 5'

*

if (nmap = 1)

'set axlim -0.6 0.6'

'set ylint 0.1'

'd ave(ave(dtdtc,lon=90,lon=100),lat=-1,lat=1)'

'd ave(ave(totalc,lon=90,lon=100),lat=-1,lat=1)'

*

'run cbar_l -x 'xs+0.3' -y 'ys+0.8' -n 2 -t "dTdt" "total"'

endif

*

if (nmap = 2)

'set axlim -0.6 0.6'

'set ylint 0.1'

'd ave(ave(advc,lon=90,lon=100),lat=-1,lat=1)'

'd ave(ave(difc,lon=90,lon=100),lat=-1,lat=1)'

'd ave(ave(fluxc,lon=90,lon=100),lat=-1,lat=1)'

'd ave(ave(adv2c,lon=90,lon=100),lat=-1,lat=1)'

'd ave(ave(entc,lon=90,lon=100),lat=-1,lat=1)'

*

'run cbar_l2 -x 'xs+1.3' -y 'ye-0.2' -n 5 -t "Flux" "Dif" "surf" "adv" "ent"'

endif

*

if (nmap = 4)

'set axlim 29 31'

'set ylint 0.5'

'd ave(ave(tmc,lon=90,lon=100),lat=-1,lat=1)'

*

'run cbar_l -x 'xs+0.3' -y 'ys+0.8' -n 1 -t "Tm"'

endif

*

xt = xs+0.2

yt = ye+0.2

'set string 1 l 4 0'

'set strsiz 0.18 0.18'

'draw string 'xt' 'yt' 'title.nmap

*

if (xmap = xmax)

*'run xcbar.gs 'xe+0.1' 'xe+0.2' 'ys' 'ye' -direction vertical -line on'

endif

*

'set parea off'

'set vpage off'

*

nmap = nmap + 1

xmap = xmap + 1

endwhile

ymap = ymap + 1

endwhile