! Siren Ruehs, November 2017
!!! About this PyFerret script:
! This is a script to reproduce Figure 2 of Ruehs et al. (2017): Eddy diffusivity estimates
! from lagrangian trajectories simulated with ocean models and surface drifter data -
! a case study for the greater Agulhas system
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! Note for NESH users only:
!To use PyFerret at NESH PyFerret has to be activated from bash shell (c-shell does not work):
!bash
!source /gpfs/fs6/home-geomar/smomw122/TM/software/miniconda3_latest/bin/activate pyferret
!pyferret
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! General settings
cancel region; cancel view; cancel win/all; cancel var/all; cancel d/all
set mem/size=2000
ppl conset,2
define symbol filename=Plot_Ruehs2017_Figure2
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! Loading data
use "SIMeddy5d_Diffu_Disp_5degbin_eSAG.nc"
use "SIMeddy5d_Diffu_Disp_5degbin_AC.nc"
use "SIMeddy5d_Diffu_Disp_5degbin_AR.nc"
use "SIMeddy5d_Diffu_Disp_5degbin_ARC.nc"
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! Plotting
DEFINE VIEWPORT/XLIM=0.07,0.59/YLIM=0.6,0.97 V1
DEFINE VIEWPORT/XLIM=0.5,0.99/YLIM=0.6,0.97 V2
DEFINE VIEWPORT/XLIM=0.07,0.59/YLIM=0.4,0.77 V3
DEFINE VIEWPORT/XLIM=0.5,0.99/YLIM=0.4,0.77 V4
DEFINE VIEWPORT/XLIM=0.07,0.59/YLIM=0.2,0.57 V5
DEFINE VIEWPORT/XLIM=0.5,0.99/YLIM=0.2,0.57 V6
DEFINE VIEWPORT/XLIM=0.07,0.59/YLIM=0,0.37 V7
DEFINE VIEWPORT/XLIM=0.5,0.99/YLIM=0,0.37 V8
set wi/aspect=1.1/thick=0.6 1
let xp2={36, 36, 40, 40, 36}
set text/font=arial
set text/font=arial/bold cbar
set text/font=arial/bold moveable
set text/font=arial hlab
set text/font=arial/bold vlab
cancel mode ascii
ppl axlint 2,2
ppl axnsig,2,0
ppl axnmtc,0,0
ppl tics 0.1,0.2,0.1,0.2,1,1
ppl axlsze 0.18, 0.18
ppl labset 0.18,0.18,0.18,0.18
define symbol fsize1=0.18
ppl color 5 35 55 75
ppl color 4 10 20 40
ppl color 2 80 10 20
set viewport V1
let yp2={12, 15, 15, 12, 12}
plot/nolab/hlim=1:40:5/vlim=0:15:3/color=black/thick=2/dashed=(0.2,0.15,0.2,0.15)/set DP2[d=1]/1000000000
! plot dispersion in 10^3 km^2
ppl ylab eSAG
ppl axlabp 0,-1
ppl plot
plot/ov/color=4/thick=3/nolab DXX[d=1]/1000000000
plot/ov/color=5/thick=3/nolab DYY[d=1]/1000000000
plot/ov/color=black/thick=3/dashed=(0.2,0.15,0.2,0.15)/nolab DP2[d=1]/1000000000
plot/vs/ov/color=4/thick=3/nolab/line {3,7},{13.5,13.5}
plot/ov/color=5/thick=3/nolab/vs/line {3,7},{11.7,11.7}
plot/ov/color=black/thick=3/dashed=(0.2,0.15,0.2,0.15)/nolab/vs/line {3,7},{9.9,9.9}
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label 37,13.2,-1,0,($fsize1) "a)"
label 21,16.5,0,0,($fsize1) "eddy dispersion in 103 km2"
label 56,16.5,0,0,($fsize1) "eddy diffusivity in 103 m2 s-1"
label 8,13.25,-1,0,($fsize1) "sxx"
label 8,11.1,-1,0,($fsize1) "syy"
label 8,9.3,-1,0,($fsize1) "sp2"
cancel var yp2
set viewport V2
let Klat_asy = KLAT[d=1,i=15:20@ave]
let yp2={2, 2.5, 2.5, 2, 2}
let xp={15, 15, 20, 20, 15}
let yp={0, 2.5, 2.5, 0, 0}
plot/nolab/hlim=1:40:5/vlim=0:2.5:0.5/color=black/thick=3/set KLAT[d=1]
ppl axlabp 0,-1
ppl plot
poly/ov/nolab/pal=myred/opacity=50 xp,yp
plot/ov/nolab/vs/line/color=(80,15,25)/thick=3 {15,40},{`Klat_asy`,`Klat_asy`}
plot/ov/nolab/color=black/thick=3 KLAT[d=1]
plot/ov/color=black/thick=3/dashed=(0.05,0.15,0.05,0.15)/nolab KDAVIS[d=1]
plot/ov/color=black/thick=3/dashed=(0.2,0.15,0.2,0.15)/nolab KDISPERS[d=1]
plot/ov/color=black/thick=3/dashed=(0.05,0.15,0.05,0.15)/nolab/vs/line {3,7},{2.25,2.25}
plot/ov/color=black/thick=3/dashed=(0.2,0.15,0.2,0.15)/nolab/vs/line {3,7},{1.95,1.95}
plot/ov/nolab/color=black/thick=3/vs/line {3,7},{1.65,1.65}
plot/ov/nolab/color=(80,15,25)/thick=3/vs/line {33,37},{0.25,0.25}
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label 37,2.2,-1,0,($fsize1) "b)"
label 8,2.15,-1,0,0.16 "kp2davis"
label 8,1.85,-1,0,0.16 "kp2disp"
label 8,1.55,-1,0,0.16 "K"
label 28,0.25,-1,0,0.16 "Kinf"
cancel var Klat_asy,xp,yp,yp2
set viewport V3
let yp2={32, 40, 40, 32, 32}
plot/nolab/vlim=0:40:8/hlim=1:40:5/color=black/thick=2/dashed=(0.2,0.15,0.2,0.15)/set DP2[d=4]/1000000000 ! in km^2
ppl axlabp 0,-1
ppl ylab ARC
ppl plot
plot/ov/color=4/thick=3/nolab DXX[d=4]/1000000000
plot/ov/color=5/thick=3/nolab DYY[d=4]/1000000000
plot/ov/color=black/thick=3/dashed=(0.2,0.15,0.2,0.15)/nolab DP2[d=4]/1000000000
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label 37,35.2,-1,0,($fsize1) "c)"
cancel var yp2
set viewport V4
let Klat_asy = KLAT[d=4,i=15:20@ave]
let yp2={12, 15, 15, 12, 12}
let xp={15, 15, 20, 20, 15}
let yp={0, 15, 15, 0, 0}
plot/nolab/hlim=1:40:5/vlim=0:15:3/color=black/thick=3/set KLAT[d=4]
ppl axlabp 0,-1
ppl plot
poly/ov/nolab/pal=myred/opacity=50 xp,yp
plot/ov/nolab/vs/line/color=(80,15,25)/thick=3 {15,40},{`Klat_asy`,`Klat_asy`}
plot/ov/color=black/thick=3/dashed=(0.05,0.15,0.05,0.15)/nolab KDAVIS[d=4]
plot/ov/color=black/thick=3/dashed=(0.2,0.15,0.2,0.15)/nolab KDISPERS[d=4]
plot/ov/nolab/color=black/thick=3 KLAT[d=4]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label 37,13.22,-1,0,($fsize1) "d)"
cancel var Klat_asy,yp,xp,yp2
set viewport V5
let yp2={80, 100, 100, 80, 80}
plot/nolab/VLIM=0:100:20/hlim=1:40:5/color=black/thick=2/dashed=(0.2,0.15,0.2,0.15)/set DP2[d=3]/1000000000 ! in km^2
ppl axlabp 0,-1
ppl ylab AR
ppl plot
plot/ov/color=4/thick=3/nolab DXX[d=3]/1000000000
plot/ov/color=5/thick=3/nolab DYY[d=3]/1000000000
plot/ov/color=black/thick=3/dashed=(0.2,0.15,0.2,0.15)/nolab DP2[d=3]/1000000000
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label 37,88,-1,0,($fsize1) "e)"
cancel var yp2
set viewport V6
let Klat_asy = KLAT[d=3,i=15:20@ave]
let yp2={20, 25, 25, 20, 20}
let xp={15, 15, 20, 20, 15}
let yp={0, 25, 25, 0, 0}
plot/nolab/hlim=1:40:5/vlim=0:25:5/color=black/thick=3/set KLAT[d=3]
ppl axlabp 0,-1
ppl plot
poly/ov/nolab/pal=myred/opacity=50 xp,yp
plot/ov/nolab/vs/line/color=(80,15,25)/thick=3 {15,40},{`Klat_asy`,`Klat_asy`}
plot/ov/color=black/thick=3/dashed=(0.05,0.15,0.05,0.15)/nolab KDAVIS[d=3]
plot/ov/color=black/thick=3/dashed=(0.2,0.15,0.2,0.15)/nolab KDISPERS[d=3]
plot/ov/nolab/color=black/thick=3 KLAT[d=3]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label 37,22,-1,0,($fsize1) "f)"
cancel var Klat_asy,yp,xp
set viewport V7
let yp2={80, 100, 100, 80, 80}
plot/nolab/VLIM=0:100:20/hlim=1:40:5/color=black/thick=2/dashed=(0.2,0.15,0.2,0.15)/set DP2[d=2]/1000000000 ! in km^2
ppl axlabp -1,-1
ppl xlab "time lag t in days"
ppl ylab AC
ppl plot
plot/ov/color=4/thick=3/nolab DXX[d=2]/1000000000
plot/ov/color=5/thick=3/nolab DYY[d=2]/1000000000
plot/ov/color=black/thick=3/dashed=(0.2,0.15,0.2,0.15)/nolab DP2[d=2]/1000000000
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label 37,88,-1,0,($fsize1) "g)"
cancel var yp2
set viewport V8
let Klat_asy = KLAT[d=2,i=15:20@ave]
let yp2={8, 10, 10, 8, 8}
let xp={15, 15, 20, 20, 15}
let yp={0, 10, 10, 0, 0}
plot/nolab/hlim=1:40:5/vlim=0:10:2/color=black/thick=3/set KLAT[d=2]
ppl axlabp -1,-1
ppl xlab "time lag t in days"
ppl plot
poly/ov/nolab/pal=myred/opacity=50 xp,yp
plot/ov/nolab/vs/line/color=(80,15,25)/thick=3 {15,40},{`Klat_asy`,`Klat_asy`}
plot/ov/color=black/thick=3/dashed=(0.05,0.15,0.05,0.15)/nolab KDAVIS[d=2]
plot/ov/color=black/thick=3/dashed=(0.2,0.15,0.2,0.15)/nolab KDISPERS[d=2]
plot/ov/nolab/color=black/thick=3 KLAT[d=2]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label 37,8.8,-1,0,($fsize1) "h)"
cancel var Klat_asy,yp,xp
FRAME/FILE=($filename).pdf