'reinit'

'open hugo-14JAN2016.hepi.ctl '

'set fwrite fig15-1.bin'
'set gxout fwrite'

'set t 9'
'set lon 0 359.5'
'set lat -89.75 89.75 '
* from 565.6 to 0.01033 hPa: 

iz = 12
while  ( iz <= 132 )
  'set z 'iz
  'd t-tl'
  iz = iz + 1
endwhile

iz = 12
while  ( iz <= 132 )
  'set z 'iz
  'd psi '
  iz = iz + 1
endwhile

'disable fwrite'

'reinit'

'open hugo.14JAN2016.hegw '

'set fwrite fig15-2.bin'
'set gxout fwrite'

'set t 3'
'set lon 0 356 '
'set lat -88 88 '
* pick the levels z=12 to z=132:
'set lev 565.6  0.01033 '  

"define nbe = ave(fcor*wir - goeL - ug*fcory - 2*vyg*vyg - 2*(wir+uyg)*uyg ,t=2,t=4)"
"define mks = ave( (vv - uu)*vyg + uu*(wir+2*uyg)-uo*dugdp-vo*dvgdp , t=2,t=4)"
"define mpc  = ave( -(ugox + vgoy) + 287.04*To/p ,t=2,t=4) "
"define mksv = ave( -uo*dugdp-vo*dvgdp , t=2,t=4)"
"define mpcv = ave( gdiv + 287.04*To/p ,t=2,t=4) "

iz = 12
while  ( iz <= 132 )
  'set z 'iz
  'd nbe '
  iz = iz + 1
endwhile

iz = 12
while  ( iz <= 132 )
  'set z 'iz
  'd mks '
  iz = iz + 1
endwhile

iz = 12
while  ( iz <= 132 )
  'set z 'iz
  'd mpc '
  iz = iz + 1
endwhile

iz = 12
while  ( iz <= 132 )
  'set z 'iz
  'd mksv '
  iz = iz + 1
endwhile

iz = 12
while  ( iz <= 132 )
  'set z 'iz
  'd mpcv '
  iz = iz + 1
endwhile

iz = 12
while  ( iz <= 132 )
  'set z 'iz
  'd psi  '
  iz = iz + 1
endwhile

iz = 12
while  ( iz <= 132 )
  'set z 'iz
  'd ug   '
  iz = iz + 1
endwhile

iz = 12
while  ( iz <= 132 )
  'set z 'iz
  'd vg   '
  iz = iz + 1
endwhile

'disable fwrite'

exit