!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=-80:0

!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

 !the actual plot
  se wi 1


   se viewpor t1 
    shade/nolables/nokey/level=(-inf)(-40,40,5)(inf)/palette=blue_red_centered      \
        TY_TRANS[i=@sum,k=@rsum,d=1,l=1:20@ave]+TY_TRANS_GM[i=@sum,d=1,l=1:20@ave]
    contour/nolables/nokey/ov/level=(-40,40,5)                                      \
        TY_TRANS[i=@sum,k=@rsum,d=1,l=1:20@ave]+TY_TRANS_GM[i=@sum,d=1,l=1:20@ave] 

   se viewport t2
    shade/nolables/nokey/level=(-inf)(-40,40,5)(inf)/palette=blue_red_centered       \
        TY_TRANS[i=@sum,k=@rsum,d=2,l=1:20@ave]+TY_TRANS_GM[i=@sum,d=2,l=1:20@ave]
    contour/nolables/nokey/ov/level=(-40,40,5)                                       \
        TY_TRANS[i=@sum,k=@rsum,d=2,l=1:20@ave]+TY_TRANS_GM[i=@sum,d=2,l=1:20@ave]

   se viewport t3 
    !add colorbar
     shade/nolables/set_up/nokey/level=(-inf)(-40,40,5)(inf)/palette=blue_red_centered       \
        TY_TRANS[i=@sum,k=@rsum,d=3,l=1:20@ave]+TY_TRANS_GM[i=@sum,d=3,l=1:20@ave]
     ppl shakey 1, 0, 0.15, , 5, 10, -5.8, 3, 0.5,0.65 
     ppl shade
    contour/nolables/nokey/ov/level=(-40,40,5)                                               \
        TY_TRANS[i=@sum,k=@rsum,d=3,l=1:20@ave]+TY_TRANS_GM[i=@sum,d=3,l=1:20@ave]


 !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_z.ps 
  spawn rm tmp.ps tmp2.ps
  spawn open ../ps/moc_z.ps