! Siren Ruehs, November 2017
!!! About this PyFerret script:
! This is a script to reproduce Figures 3-5 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 filename1=Plot_Ruehs2017_Figure3
define symbol filename2=Plot_Ruehs2017_Figure4
define symbol filename3=Plot_Ruehs2017_Figure5
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! Loading data
use "INALT01_topo_k2.nc"
use "Obs_Diffu_5degbins.nc"
use "SIMeddy5d_Diffu_5degbins.nc"
use "Obs_EKE_5degbins.nc"
use "SIMeddy5d_EKE_5degbins.nc"
use "Diffu_topo_mask.nc"
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! Defining Variables
define axis/y/units="degrees_N" ylat=lat[d=2,i=1]
define axis/x/units="degrees_E" xlon=lon[d=2,j=1]
let TopoMask = TOPOMASK50[d=6,gx=xlon@asn,gy=ylat@asn]
! Kmax and kasy have been calculated appr. 2° land inwards (because of 5° square) -> have to mask it!
let yp2={-49, -43, -43, -49, -49}
let xp2={-10, -10, -3.5, -3.5, -3.5}
!!!!!!!!!!!!!!!!!!!!!
!!! OBS
!! diffusivities in m^2/s
let Kmax_obs = Kmax[D=2,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(Kmax_obs)
let Kmax_disp_obs= Kmax_dis[D=2,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(Kmax_disp_obs)
let Kmax_dav_obs = Kmax_dav[D=2,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(Kmax_dav_obs)
let Kinf_obs = Kasymptotic[D=2,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(Kinf_obs)
let Kinf_disp_obs= Kasymptotic_dis[D=2,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(Kinf_disp_obs)
let Kinf_dav_obs = Kasymptotic_dav[D=2,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(Kinf_dav_obs)
let tmp=Kmax_obs/Kinf_obs*TopoMask
stats ignore0(tmp)
let tmp_disp=Kmax_disp_obs/Kinf_disp_obs*TopoMask
stats ignore0(tmp_disp)
let tmp_dav=Kmax_dav_obs/Kinf_dav_obs*TopoMask
stats ignore0(tmp_dav)
!! velocity in m/s, EKE in m^2/s^2
let speed_obs = CURRENTVELOCITY[d=4,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(speed_obs)
let angle_obs = CURRENTANGLE[d=4,gx=xlon@asn,gy=ylat@asn]*TopoMask
let mask_lowspeed_obs = if speed_obs le 0.1 then 1 else 0
let mask_highspeed_obs = if speed_obs gt 0.1 then 1 else 0
let vel_v_obs = sin(angle_sim*(3.14159/180))*speed_obs
set var/bad=0 vel_v_obs
let vel_u_obs = cos(angle_obs*(3.14159/180))*speed_obs
set var/bad=0 vel_u_obs
let eke_obs = EKE[D=4,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(eke_obs)
let Ncountstime_obs = DRIFTERTIMESTATS[D=4,gx=xlon@asn,gy=ylat@asn]
let Ncounts_obs = DRIFTERSTATS[D=4,gx=xlon@asn,gy=ylat@asn]
!!!!!!!!!!!!!!!!!!!!!
!!! INALT01
!! diffusivities in m^2/s
let Kmax_sim = Kmax[D=3,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(Kmax_sim)
let Kmax_dav_sim = Kmax_Davis[D=3,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(Kmax_dav_sim)
let Kmax_disp_sim = Kmax_Cloud[D=3,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(Kmax_disp_sim)
let Kinf_sim = Kasymptotic[D=3,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(Kinf_sim)
let Kinf_dav_sim = Kasymptotic_Davis[D=3,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(Kinf_dav_sim)
let Kinf_disp_sim = Kasymptotic_Cloud[D=3,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(Kinf_disp_sim)
let tmp_sim=Kmax_sim/Kinf_sim*TopoMask
stats ignore0(tmp_sim)
let tmp_disp_sim=Kmax_disp_sim/Kinf_disp_sim*TopoMask
stats ignore0(tmp_disp_sim)
let tmp_dav_sim=Kmax_dav_sim/Kinf_dav_sim*TopoMask
stats ignore0(tmp_sim)
!! velocity in m/s, EKE in m^2/s^2
let speed_sim = CURRENTVELOCITY[d=5,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(speed_sim)
let angle_sim = CURRENTANGLE[d=5,gx=xlon@asn,gy=ylat@asn]*TopoMask
let mask_lowspeed_sim = if speed_sim le 0.1 then 1 else 0
let mask_highspeed_sim = if speed_sim gt 0.1 then 1 else 0
let vel_v_sim = sin(angle_sim*(3.14159/180))*speed_sim*TopoMask
set var/bad=0 vel_v_sim
let vel_u_sim = cos(angle_sim*(3.14159/180))*speed_sim*TopoMask
set var/bad=0 vel_u_sim
let eke_sim = EKE[D=5,gx=xlon@asn,gy=ylat@asn]*TopoMask
stats ignore0(eke_sim)
let Ncountstime_sim = DRIFTERTIMESTATS[D=5,gx=xlon@asn,gy=ylat@asn]
let Ncounts_sim = DRIFTERSTATS[D=5,gx=xlon@asn,gy=ylat@asn]
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! Plotting
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! Manuscript Figure3: asymptotic diffusivities !!!!!!!!!!!!!!!!!!!!!!!!!
set wi/aspect=0.75/thicken=0.5 1
DEFINE VIEWPORT/XLIM=0.0,0.34/YLIM=0.5,0.9 V1
DEFINE VIEWPORT/XLIM=0.27,0.61/YLIM=0.5,0.9 V2
DEFINE VIEWPORT/XLIM=0.54,0.88/YLIM=0.5,0.9 V3
DEFINE VIEWPORT/XLIM=0.0,0.34/YLIM=0.27,0.67 V4
DEFINE VIEWPORT/XLIM=0.27,0.61/YLIM=0.27,0.67 V5
DEFINE VIEWPORT/XLIM=0.54,0.88/YLIM=0.27,0.67 V6
DEFINE VIEWPORT/XLIM=0.0,0.34/YLIM=0,0.4 V7
DEFINE VIEWPORT/XLIM=0.27,0.61/YLIM=0,0.4 V8
DEFINE VIEWPORT/XLIM=0.54,0.88/YLIM=0,0.4 V9
set text/font=arial
set text/font=arial cbar
set text/font=arial/bold moveable
set text/font=arial contour
cancel mode ascii
ppl axlint 3,2 !label intervals of large tics
ppl axnmtc,2,2
ppl tics 0.1,0.2,0.1,0.2,1,1
ppl axlsze 0.18, 0.18
define symbol fsize1=0.18
define symbol fsize=0.18
set viewport V1
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=broad/levels=(0)(0.18)(0.25)(0.35)(0.5)(0.7)(1.0)(1.25)(2)(2.83)(4)(5.66)(8)(11.31)(16)(22.63)(32)(inf)/set Kinf_dav_sim
ppl axlabp,-1,1
ppl shakey 0
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(1)(2.5)(5)(10)(15) Kinf_dav_sim
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "a)"
label 20,-17,0,0,($fsize2) "Kinfdavis"
label -15,-34,0,90,($fsize2) "SIMeddy-5d"
set viewport V2
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=broad/levels=(0)(0.18)(0.25)(0.35)(0.5)(0.7)(1.0)(1.25)(2)(2.83)(4)(5.66)(8)(11.31)(16)(22.63)(32)(inf)/set Kinf_disp_sim
ppl axlabp,-1,1
ppl shakey 0
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(1)(2.5)(5)(10)(15) Kinf_disp_sim
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "b)"
label 20,-17.0,0,0,($fsize2) "Kinfdisp"
set viewport V3
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=broad/levels=(0)(0.18)(0.25)(0.35)(0.5)(0.7)(1.0)(1.25)(2)(2.83)(4)(5.66)(8)(11.31)(16)(22.63)(32)(inf)/set Kinf_sim
ppl axlabp,-1,0
ppl shakey 0
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(1)(2.5)(5)(10)(15) Kinf_sim
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "c)"
label 20,-17.0,0,0,($fsize2) "Kinf"
set viewport V4
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=broad/levels=(0)(0.18)(0.25)(0.35)(0.5)(0.7)(1.0)(1.25)(2)(2.83)(4)(5.66)(8)(11.31)(16)(22.63)(32)(inf)/set Kinf_dav_obs
ppl axlabp,0,1
ppl shakey 0
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(1)(2.5)(5)(10)(15) Kinf_dav_obs
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "d)"
label -15,-34.0,0,90,($fsize2) "OBS"
set viewport V5
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=broad/levels=(0)(0.18)(0.25)(0.35)(0.5)(0.7)(1.0)(1.25)(2)(2.83)(4)(5.66)(8)(11.31)(16)(22.63)(32)(inf)/set Kinf_disp_obs
ppl axlabp,0,1
ppl shakey 0
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(1)(2.5)(5)(10)(15) Kinf_disp_obs
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "e)"
set viewport V6
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=broad/levels=(0)(0.18)(0.25)(0.35)(0.5)(0.7)(1.0)(1.25)(2)(2.83)(4)(5.66)(8)(11.31)(16)(22.63)(32)(inf)/set Kinf_obs
ppl axlabp,0,0
ppl shakey 1,1,($fsize1),2,4,4,5.35,5.6,1.4,6.4
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(1)(2.5)(5)(10)(15) Kinf_obs
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "f)"
label 60,14,0,0,($fsize1) "103"
label 60,10,0,0,($fsize1) "m2 s-1"
set viewport V7
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=blue_darkred/levels=(-inf)(-1)(-0.8)(-0.6)(-0.4)(-0.2)(0)(0.2)(0.4)(0.6)(0.8)(1)(inf)/set (Kinf_dav_sim-Kinf_dav_obs)/(Kinf_dav_obs)
ppl axlabp,-1,1
ppl shakey 0
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(-1)(-0.5)(0)(0.5)(1) (Kinf_dav_sim-Kinf_dav_obs)/(Kinf_dav_obs)
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
label -18,-34.0,0,90,($fsize2) "(SIMeddy-5d - OBS) /"
label -13,-34.0,0,90,($fsize2) "OBS"
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "g)"
set viewport V8
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=blue_darkred/levels=(-inf)(-1)(-0.8)(-0.6)(-0.4)(-0.2)(0)(0.2)(0.4)(0.6)(0.8)(1)(inf)/set (Kinf_disp_sim-Kinf_disp_obs)/(Kinf_disp_obs)
ppl axlabp,-1,1
ppl shakey 0
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(-1)(-0.5)(0)(0.5)(1) (Kinf_disp_sim-Kinf_disp_obs)/(Kinf_disp_obs)
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "h)"
set viewport V9
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=blue_darkred/levels=(-inf)(-1)(-0.8)(-0.6)(-0.4)(-0.2)(0)(0.2)(0.4)(0.6)(0.8)(1)(inf)/set (Kinf_sim-Kinf_obs)/(Kinf_obs)
ppl axlabp,-1,0
ppl shakey 1,1,($fsize1),2,4,4,5.35,5.6,1.4,4
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(-1)(-0.5)(0)(0.5)(1) (Kinf_sim-Kinf_obs)/(Kinf_obs)
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "i)"
FRAME/FILE=($filename1).pdf
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! Manuscript Figure 4: Maximum Diffusivities !!!!!!!!!!!!!!!!!!!!!!!!!!!
set wi/aspect=0.75/thicken=0.5 2
DEFINE VIEWPORT/XLIM=0.0,0.34/YLIM=0.5,0.9 V1
DEFINE VIEWPORT/XLIM=0.27,0.61/YLIM=0.5,0.9 V2
DEFINE VIEWPORT/XLIM=0.58,0.92/YLIM=0.5,0.9 V3
DEFINE VIEWPORT/XLIM=0.0,0.34/YLIM=0.27,0.67 V4
DEFINE VIEWPORT/XLIM=0.27,0.61/YLIM=0.27,0.67 V5
DEFINE VIEWPORT/XLIM=0.58,0.92/YLIM=0.27,0.67 V6
DEFINE VIEWPORT/XLIM=0.0,0.34/YLIM=0,0.4 V7
DEFINE VIEWPORT/XLIM=0.27,0.61/YLIM=0,0.4 V8
DEFINE VIEWPORT/XLIM=0.58,0.92/YLIM=0,0.4 V9
set text/font=arial
set text/font=arial cbar
set text/font=arial/bold moveable
set text/font=arial contour
cancel mode ascii
ppl axlint 3,2 !label intervals of large tics
ppl axnmtc,2,2
ppl tics 0.1,0.2,0.1,0.2,1,1
ppl axlsze 0.18, 0.18
define symbol fsize1=0.18
define symbol fsize2=0.18
set viewport V2
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=broad/levels=(0)(0.18)(0.25)(0.35)(0.5)(0.7)(1.0)(1.25)(2)(2.83)(4)(5.66)(8)(11.31)(16)(22.63)(32)(inf)/set Kmax_sim
ppl axlabp,-1,0
ppl shakey 0
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(1)(2.5)(5)(10)(15) Kmax_sim
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
label 20,-17.0,0,0,($fsize2) "Kmax"
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "a)"
label -15,-34,0,90,($fsize2) "SIMeddy-5d"
set viewport V3
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=orange_red/levels=(1,3,0.1)(inf)/set Kmax_sim/Kinf_sim
ppl axlabp,-1,0
ppl shakey 0
ppl shade
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
vector/aspect/ov/nolab/length=0.15/xskip=2/yskip=1/thick=1 vel_u_sim*TopoMask*mask_lowspeed_sim,vel_v_sim*TopoMask*mask_lowspeed_sim
vector/aspect/ov/nolab/length=0.3/xskip=2/yskip=1/thick=3 vel_u_sim*TopoMask*mask_highspeed_sim,vel_v_sim*TopoMask*mask_highspeed_sim
label 20,-17.0,0,0,($fsize2) "Kmax / Kinf"
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "b)"
set viewport V5
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=broad/levels=(0)(0.18)(0.25)(0.35)(0.5)(0.7)(1.0)(1.25)(2)(2.83)(4)(5.66)(8)(11.31)(16)(22.63)(32)(inf)/set Kmax_obs
ppl axlabp,0,0
ppl shakey 1,1,($fsize1),2,4,4,5.35,5.6,1.4,6.4
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(1)(2.5)(5)(10)(15) Kmax_obs
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "c)"
label -15,-34,0,90,($fsize2) "OBS"
label 60,14,0,0,($fsize1) "103"
label 60,10,0,0,($fsize1) "m2 s-1"
set viewport V6
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=orange_red/levels=(1,3,0.1)(inf)/set Kmax_obs/Kinf_obs
ppl axlabp,0,0
ppl shakey 1,1,($fsize1),4,4,4,5.35,5.6,1.4,6.4
ppl shade
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
vector/aspect/ov/nolab/length=0.15/xskip=2/yskip=1/thick=1 vel_u_obs*TopoMask*mask_lowspeed_obs,vel_v_obs*TopoMask*mask_lowspeed_obs
vector/aspect/ov/nolab/length=0.3/xskip=2/yskip=1/thick=3 vel_u_obs*TopoMask*mask_highspeed_obs,vel_v_obs*TopoMask*mask_highspeed_obs
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "d)"
set viewport V8
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=blue_darkred/levels=(-inf)(-1)(-0.8)(-0.6)(-0.4)(-0.2)(0)(0.2)(0.4)(0.6)(0.8)(1)(inf)/set (Kmax_sim-Kmax_obs)/(Kmax_obs)
ppl axlabp,-1,1
ppl shakey 0
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(-1)(-0.5)(0)(0.5)(1) (Kmax_sim-Kmax_obs)/(Kmax_obs)
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
label -18,-34.0,0,90,($fsize2) "(SIMeddy-5d - OBS) /"
label -13,-34.0,0,90,($fsize2) "OBS"
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "e)"
set viewport V9
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=blue_darkred/levels=(-inf)(-1)(-0.8)(-0.6)(-0.4)(-0.2)(0)(0.2)(0.4)(0.6)(0.8)(1)(inf)/set ((Kmax_sim/Kinf_sim)-(Kmax_obs/Kinf_obs))/(Kmax_obs/Kinf_obs)
ppl axlabp,-1,0
ppl shakey 1,1,($fsize1),2,4,4,5.35,5.6,1.4,4
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(-1)(-0.5)(0)(0.5)(1) ((Kmax_sim/Kinf_sim)-(Kmax_obs/Kinf_obs))/(Kmax_obs/Kinf_obs)
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "f)"
FRAME/FILE=($filename2).pdf
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! Manuscript Figure 5: Speed and EKE derived from Lagrangian data !!!!!!
set wi/aspect=0.75/thicken=0.5 3
set text/font=arial
set text/font=arial cbar
set text/font=arial/bold moveable
cancel mode ascii
ppl axlint 3,2 !label intervals of large tics
ppl axnmtc,2,2
ppl tics 0.1,0.2,0.1,0.2,1,1
ppl axlsze 0.18, 0.18
define symbol fsize1=0.18
define symbol fsize2=0.18
set viewport V2
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/levels=(0,20,1)(inf)/pal=mycolbar/set speed_sim*100
ppl axlabp,-1,0
ppl shakey 0
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(2.5)(0,30,5) speed_sim*100
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
vector/aspect/ov/nolab/length=0.15/xskip=2/yskip=1/thick=1 vel_u_sim*TopoMask*mask_lowspeed_sim,vel_v_sim*TopoMask*mask_lowspeed_sim
vector/aspect/ov/nolab/length=0.3/xskip=2/yskip=1/thick=3 vel_u_sim*TopoMask*mask_highspeed_sim,vel_v_sim*TopoMask*mask_highspeed_sim
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "a)"
label 20,-17,0,0,($fsize2) "velocity"
label -15,-34.0,0,90,($fsize2) "SIMeddy-5d"
set viewport V3
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=yellow_orange_brown/levels=(0)(2.5)(5)(7.5)(10)(30)(50)(75)(100)(200)(300)(500)(700)(1000)(1300)(inf)/set eke_sim*10000
ppl axlabp,-1,0
ppl shakey 0
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(0)(5)(10)(50)(100)(300)(700)(1300) eke_sim*10000
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "b)"
label 20,-17,0,0,($fsize2) "EKE"
set viewport V5
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/levels=(0,20,1)(inf)/pal=mycolbar/set speed_obs*100
ppl axlabp,0,0
ppl shakey 1,1,($fsize1),4,4,4,5.3,5.55,1.4,6.4
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(2.5)(0,30,5) speed_obs*100
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
vector/aspect/ov/nolab/length=0.15/xskip=2/yskip=1/thick=1 vel_u_obs*TopoMask*mask_lowspeed_obs,vel_v_obs*TopoMask*mask_lowspeed_obs
vector/aspect/ov/nolab/length=0.3/xskip=2/yskip=1/thick=3 vel_u_obs*TopoMask*mask_highspeed_obs,vel_v_obs*TopoMask*mask_highspeed_obs
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "c)"
label -15,-34.0,0,90,($fsize2) "OBS"
label 60,10,0,0,($fsize1) "cm s-1"
set viewport V6
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=yellow_orange_brown/levels=(0)(2.5)(5)(7.5)(10)(30)(50)(75)(100)(200)(300)(500)(700)(1000)(1300)(inf)/set eke_obs*10000
ppl axlabp,0,0
ppl shakey 1,1,($fsize1),4,4,4,5.4,5.65,1.4,6.4
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(0)(5)(10)(50)(100)(300)(700)(1300) eke_obs*10000
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "d)"
label 60,10,0,0,($fsize1) "cm2 s-2"
set viewport V8
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=blue_darkred/levels=(-inf)(-1)(-0.8)(-0.6)(-0.4)(-0.2)(0)(0.2)(0.4)(0.6)(0.8)(1)(inf)/set (speed_sim-speed_obs)/(speed_obs)
ppl axlabp,-1,1
ppl shakey 0
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(-1)(-0.5)(0)(0.5)(1) (speed_sim-speed_obs)/(speed_obs)
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "e)"
label -18,-34.0,0,90,($fsize2) "(SIMeddy-5d - OBS) /"
label -13,-34.0,0,90,($fsize2) "OBS"
set viewport V9
shade/nolab/hlim=-10:50:6/vlim=-49:-19:6/pal=blue_darkred/levels=(-inf)(-1)(-0.8)(-0.6)(-0.4)(-0.2)(0)(0.2)(0.4)(0.6)(0.8)(1)(inf)/set (eke_sim-eke_obs)/(eke_obs)
ppl axlabp,-1,0
ppl shakey 1,1,($fsize1),2,4,4,5.3,5.55,1.3,4.05
ppl shade
contour/ov/pen=7/size=($fsize1)/nolab/lev=(-1)(-0.5)(0)(0.5)(1) (eke_sim-eke_obs)/(eke_obs)
shade/ov/nolab/pal=gray_light/lev=(0) TopoMask
shade/ov/nolab/pal=gray/lev=(0) topo[d=1]
polygon/ov/nolab/pal=white/color=black/thick=2/opacity=100 xp2,yp2
label -9,-47.5,-1,0,($fsize2) "f)"
FRAME/FILE=($filename3).pdf