sdfopen /grove3/iesuw/foam/ho.F1.0k.seae486605.nc
  sdfopen /grove3/iesuw/foam/ho.F1.PBC_TP_CLEAN.seae486605.nc

  set lat -90 90
  set lev -1000 -20

  enable print yz_crosssect.prt

* turn grads grid off
  set grid off

* set time - order for seae is DJF,MAM,JJA,SON,ANN
  set t 5

  define miss = 9e+20
  maskt1 = maskout(t.1, -t.1+70)
  maskt2 = maskout(t.2, -t.2+70)

  avet1 = ave(maskt1,lon=300,lon=360)
  avet2 = ave(maskt2,lon=300,lon=360)

* draw strings on page to label plots
* draw string xpos ypos string
  draw string 7.5 8.25 ANN

* set dimensions on page of first map
* top center of page
* set vpage xmin xmax ymin ymax
  set vpage 0 7.5 6.0 10.5

* turn grads label off (must do this before every plot is drawn)
  set grads off

  set xflip on

* put this 'set x 1' command AFTER averaging (avet1,etc) is done
  set x 1

* default cint ANN range -0.2 to 0.4, cint 0.05
  set cint 0.05

  d avet2-avet1

* label y axis
  draw ylab levels (m)

* draw title
* putting in back slash makes text after the slash go onto a new line
  draw title FOAM T(C) PBC_TP_CLEAN-FC Atlantic 60W-0

* print contour label at bottom
* could also use 'draw string' for this, but
* draw xlab gives nice big type and centers at bottom of plot
  draw xlab cint 0.05

  print
  disable print