'reinit'
'open ens_GFS_00Z_z700.ctl'

*First fine the average for 700mb, day 4

'set t 9'
'define zave=ave(h,z=2,z=12)'

'query file'
ln=sublin(result,5)
lm=subwrd(ln,9)
say lm
l=1
while(l <= lm)
'define p2=lev'
'set z 'l
say l
'define mem'l'= h'
l=l+1
endwhile
*
'define dev02=abs(zave-mem2)'
'define dev03=abs(zave-mem3)'
'define dev04=abs(zave-mem4)'
'define dev05=abs(zave-mem5)'
'define dev06=abs(zave-mem6)'
'define dev07=abs(zave-mem7)'
'define dev08=abs(zave-mem8)'
'define dev09=abs(zave-mem9)'
'define dev10=abs(zave-mem10)'
'define dev11=abs(zave-mem11)'
*Take area average for northern hemisphere:
'define dev02=aave(dev02,lon=0,lon=360,lat=0,lat=90)'
'define dev03=aave(dev03,lon=0,lon=360,lat=0,lat=90)'
'define dev04=aave(dev04,lon=0,lon=360,lat=0,lat=90)'
'define dev05=aave(dev05,lon=0,lon=360,lat=0,lat=90)'
'define dev06=aave(dev06,lon=0,lon=360,lat=0,lat=90)'
'define dev07=aave(dev07,lon=0,lon=360,lat=0,lat=90)'
'define dev08=aave(dev08,lon=0,lon=360,lat=0,lat=90)'
'define dev09=aave(dev09,lon=0,lon=360,lat=0,lat=90)'
'define dev10=aave(dev10,lon=0,lon=360,lat=0,lat=90)'
'define dev11=aave(dev11,lon=0,lon=360,lat=0,lat=90)'
'd dev02'
ln=sublin(result,1)
lm2=subwrd(ln,4)
say lm2
'd dev03'
ln=sublin(result,1)
lm3=subwrd(ln,4)
say lm3
'd dev04'
ln=sublin(result,1)
lm4=subwrd(ln,4)
say lm4
'd dev05'
ln=sublin(result,1)
lm5=subwrd(ln,4)
say lm5
'd dev06'
ln=sublin(result,1)
lm6=subwrd(ln,4)
say lm6
'd dev07'
ln=sublin(result,1)
lm7=subwrd(ln,4)
say lm7
'd dev08'
ln=sublin(result,1)
lm8=subwrd(ln,4)
say lm8
'd dev09'
ln=sublin(result,1)
lm9=subwrd(ln,4)
say lm9
'd dev10'
ln=sublin(result,1)
lm10=subwrd(ln,4)
say lm10
'd dev11'
ln=sublin(result,1)
lm11=subwrd(ln,4)
say lm11

*Now find lowest area-averaged value, or that value
*which is closest to area-averaged mean value
mtot=lm2
if (lm3 < mtot)
mtot=lm3
endif
if (lm4 < mtot)
mtot=lm4
endif
if (lm5 < mtot)
mtot=lm5
endif
if (lm6 < mtot)
mtot=lm6
endif
if (lm7 < mtot)
mtot=lm7
endif
if (lm8 < mtot)
mtot=lm8
endif
if (lm9 < mtot)
mtot=lm9
endif
if (lm10 < mtot)
mtot=lm10
endif
if (lm11 < mtot)
mtot=lm11
endif
say mtot
var=0
*

if (mtot = lm2)
tt=02
say tt
'define var'tt'=0'
'q define'
res=subwrd(result,49)
final=substr(res,4,2)
write('consmem.txt',final)
endif

if (mtot = lm3)
tt=03
say tt
'define var'tt'=0'
'q define'
res=subwrd(result,49)
final=substr(res,4,2)
write('consmem.txt',final)
endif

if (mtot = lm4)
tt=04
say tt
'define var'tt'=0'
'q define'
res=subwrd(result,49)
final=substr(res,4,2)
write('consmem.txt',final)
endif

if (mtot = lm5)
tt=05
say tt
'define var'tt'=0'
'q define'
res=subwrd(result,49)
final=substr(res,4,2)
write('consmem.txt',final)
endif

if (mtot = lm6)
tt=06
say tt
'define var'tt'=0'
'q define'
res=subwrd(result,49)
final=substr(res,4,2)
write('consmem.txt',final)
endif


if (mtot = lm7)
tt=07
say tt
'define var'tt'=0'
'q define'
res=subwrd(result,49)
final=substr(res,4,2)
write('consmem.txt',final)
endif

if (mtot = lm8)
tt=08
say tt
'define var'tt'=0'
'q define'
res=subwrd(result,49)
final=substr(res,4,2)
write('consmem.txt',final)
endif

if (mtot = lm9)
tt=09
say tt
'define var'tt'=0'
'q define'
res=subwrd(result,49)
final=substr(res,4,2)
write('consmem.txt',final)
endif

if (mtot = lm10)
tt=10
say tt
'define var'tt'=0'
'q define'
res=subwrd(result,49)
final=substr(res,4,2)
write('consmem.txt',final)
endif

if (mtot = lm11)
tt=11
say tt
'define var'tt'=0'
'q define'
res=subwrd(result,49)
final=substr(res,4,2)
write('consmem.txt',final)
endif

*'d dm'var''