!meridional oberturning circualtion in the souther ocean on sigma  

!clean
 cancel da/all
 cancel re/all

!define input
 use "/Users/hdietze/mom4p1_pubrel_18dec2009/exp/mom_sis_bling_kisokgm600/out/physic"
 use "/Users/hdietze/mom4p1_pubrel_18dec2009/exp/mom_sis_bling_default/out/physic"
 use "/Users/hdietze/mom4p1_pubrel_18dec2009/exp/mom_sis_bling_eg/out/physic"

!define region
 set region/y=-60:-45/z=0:4000

!some generic plot ini 
 define viewport/xlimits=0  ,0.257    /ylim=0.25,0.5        t1
 define viewport/xlimits=0.2,0.4570   /ylim=0.25,0.5        t2
 define viewport/xlimits=0.4,0.6570   /ylim=0.25,0.5        t3

!some specific plot ini
 go mytransparency.jnl
 cancel mode logo
 set mode meta


!some calculations
 let diff_kisokgm600 = TY_TRANS[i=@sum,k=@rsum,d=1,l=69]+TY_TRANS_GM[i=@sum,d=1,l=69] - \
        (TY_TRANS[i=@sum,k=@rsum,d=1,l=1:20@ave]+TY_TRANS_GM[i=@sum,d=1,l=1:20@ave])

 let diff_default    = TY_TRANS[i=@sum,k=@rsum,d=2,l=69]+TY_TRANS_GM[i=@sum,d=2,l=69] - \
        (TY_TRANS[i=@sum,k=@rsum,d=2,l=1:20@ave]+TY_TRANS_GM[i=@sum,d=2,l=1:20@ave])

 let diff_eg         =  TY_TRANS[i=@sum,k=@rsum,d=3,l=69]+TY_TRANS_GM[i=@sum,d=3,l=69] -       \
       (TY_TRANS[i=@sum,k=@rsum,d=3,l=1:20@ave]+TY_TRANS_GM[i=@sum,d=3,l=1:20@ave])

 !the actual plot
  se wi 1


   se viewpor t1 
    shade/nolables/nokey/level=(-1,12,1)         \
         diff_kisokgm600
    contour/nolables/nokey/ov/level=(-1,12,1)                                          \
         diff_kisokgm600 

   se viewport t2
    shade/nolables/nokey/level=(-1,12,1)         \
         diff_default
    contour/nolables/nokey/ov/level=(-1,12,1)                                          \
         diff_default

   se viewport t3 
    !add colorbar
     shade/nolables/set_up/nokey/level=(-1,12,1)       \
          diff_eg
     ppl shakey 1, 0, 0.15, , 5, 10, -5.8, 3, 0.5,0.65 
     ppl shade
    contour/nolables/nokey/ov/level=(-1,12,1)                                               \
           diff_eg

 !bless and print the plot
 
  go remove_logo.jnl
  cancel mode meta
  ppl clsplt 
  spawn Fprint -l cps -o tmp.ps -R metafile.plt
  spawn rm metafile.plt
  spawn ps_thicken tmp.ps 2 > tmp2.ps 
  spawn ps_no_bleedthrough tmp2.ps > ../ps/moc_southern_z_winddiff.ps 
  spawn rm tmp.ps tmp2.ps
  spawn open ../ps/moc_southern_z_winddiff.ps

!some calculations
 list diff_kisokgm600[j=@max,k=@max]
 list diff_default[j=@max,k=@max]
 list diff_eg[j=@max,k=@max]

 list diff_kisokgm600[j=@ave,k=@ave]
 list diff_default[j=@ave,k=@ave]
 list diff_eg[j=@ave,k=@ave]