!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_default/out/physic"
! use "/Users/hdietze/mom4p1_pubrel_18dec2009/exp/mom_sis_bling_eg/out/physic"
 use "/Users/hdietze/mom4p1_pubrel_18dec2009/exp/mom_sis_bling_kisokgm600/out/physic"


!some generic plot ini 
 define viewport/xlimits=0    ,0.4825   /ylim=0.5,1        t1
 define viewport/xlimits=0.4125,0.895  /ylim=0.5,1        t2
 
 define viewport/xlimits=0    ,0.4825   /ylim=0.17,0.67    t3
 define viewport/xlimits=0.4125,0.895  /ylim=0.17,0.67    t4 

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

 !the actual plot
  se wi 1

!add colorbar
 se viewport t4
 shade/nolabels/set_up/nokey/palette=blue_red_centered/hlimit=70s:35s/vlimit=1034:1037.5/level=(-inf)(-40,40,5)(inf) \
      TY_TRANS_RHO[i=@sum,k=@rsum,d=1,l=1:20@ave]+TY_TRANS_RHO_GM[i=@sum,d=1,l=1:20@ave] 
 ppl shakey 1, 1, 0.15, , 3, 6, 6.5, 6.8, 1.5, 8.9
 ppl shade

   se viewport t1 
    fill/nolables/nokey/hlimit=70s:35s/vlimit=1034:1037.5/level=(-inf)(-40,40,5)(inf)/palette=blue_red_centered \
        TY_TRANS_RHO[i=@sum,k=@rsum,d=1,l=1:20@ave]+TY_TRANS_RHO_GM[i=@sum,d=1,l=1:20@ave]
    contour/nolables/nokey/ov/level=(-40,40,2) TY_TRANS_RHO[i=@sum,k=@rsum,d=1,l=1:20@ave]+TY_TRANS_RHO_GM[i=@sum,d=1,l=1:20@ave]

   se viewport t3
    fill/nolables/nokey/hlimit=70s:35s/vlimit=1034:1037.5/level=(-inf)(-40,40,5)(inf)/palette=blue_red_centered \
        TY_TRANS_RHO[i=@sum,k=@rsum,d=1,l=60:70@ave]+TY_TRANS_RHO_GM[i=@sum,d=1,l=60:70@ave]
    contour/nolables/nokey/ov/level=(-40,40,2) TY_TRANS_RHO[i=@sum,k=@rsum,d=1,l=60:70@ave]+TY_TRANS_RHO_GM[i=@sum,d=1,l=60:70@ave]


   se viewport t2 
    fill/nolables/nokey/hlimit=70s:35s/vlimit=1034:1037.5/level=(-inf)(-40,40,5)(inf)/palette=blue_red_centered \
     TY_TRANS_RHO[i=@sum,k=@rsum,d=1,l=60:70@ave]+TY_TRANS_RHO_GM[i=@sum,d=1,l=60:70@ave] - TY_TRANS_RHO[i=@sum,k=@rsum,d=1,l=1:20@ave] - TY_TRANS_RHO_GM[i=@sum,d=1,l=1:20@ave]
    contour/nolables/nokey/ov/level=(-40,40,2) \
     TY_TRANS_RHO[i=@sum,k=@rsum,d=1,l=60:70@ave]+TY_TRANS_RHO_GM[i=@sum,d=1,l=60:70@ave] - TY_TRANS_RHO[i=@sum,k=@rsum,d=1,l=1:20@ave] - TY_TRANS_RHO_GM[i=@sum,d=1,l=1:20@ave] 

   se viewport t4 
    fill/nolables/nokey/hlimit=70s:35s/vlimit=1034:1037.5/level=(-inf)(-40,40,5)(inf)/palette=blue_red_centered \
        TY_TRANS_RHO_GM[i=@sum,d=1,l=60:70@ave] - TY_TRANS_RHO_GM[i=@sum,d=1,l=1:20@ave]
    contour/nolables/nokey/ov/level=(-40,40,2) \
        TY_TRANS_RHO_GM[i=@sum,d=1,l=60:70@ave] - TY_TRANS_RHO_GM[i=@sum,d=1,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 3 > tmp2.ps 
  spawn ps_no_bleedthrough tmp2.ps > ../ps/moc_southern_sigma_kisokgm600.ps 
  spawn rm tmp.ps tmp2.ps
!  spawn open ../ps/moc_southern_sigma_default.ps