! Comparing Aeolus to Era Interim observational data and CM2M: u surface, v surface, u_za, v_za, vertical, zonal and meridional synoptic kinetic energy, w_za
!
! Arguments
! $1 = Season index: 1=DJF, 2=MAM, 3=JJA, 4=SON, ANNUAL=1:4@ave
! $2 = Season: DJF, MAM, JJA, SON, ANNUAL
! $3 = dataset aeolus
! $4 = dataset atmos model down

cancel data/all

!spawn cdo -seldate\,0052-01-01T00:00:00\,0104-12-31T23:59:59 $3 data.nc
!spawn ncatted -a axis\,level\,c\,c\,Z data.nc
!spawn cdo yseasmean data.nc data_sm.nc
!spawn cdo settaxis,2013-02-01,00:00,3months data_sm.nc data_sm_taxis.nc 
!spawn ncks -A -v levelb $3 data_sm_taxis.nc
!spawn ncatted -a edges\,level\,c\,c\,levelb data_sm_taxis.nc
!spawn rm -f data.nc data_sm.nc
!use "data_sm_taxis.nc" !d=1 AEOLUS DATA
use $3 !d=1 AEOLUS DATA
use "/p/projects/climber3/caesar/CM2M/multi_year_seasonal_means/00530101.atmos-dly_sm.nc" !d=2 CM2M DATA
use "/p/projects/climber3/caesar/CM2M/multi_year_seasonal_means/00530101.flux-dly_sm.nc" !d=3 CM2M DATA
use "/p/projects/climber3/caesar/Era_interim_data/noleap_data/multi_year_season_mean/10_m_u_wind_sm.nc" !d=4
use "/p/projects/climber3/caesar/Era_interim_data/noleap_data/multi_year_season_mean/10_m_v_wind_sm.nc" !d=5
use "/p/projects/climber3/caesar/Era_interim_data/noleap_data/multi_year_season_mean/U_wind_pressure_levels_daily_means_monthly_means_sm.nc" !d=6
use "/p/projects/climber3/caesar/Era_interim_data/noleap_data/multi_year_season_mean/V_wind_pressure_levels_daily_means_monthly_means_sm.nc" !d=7

!spawn cdo -seldate\,0052-01-01T00:00:00\,0104-12-31T23:59:59 $4 data_down.nc
!spawn ncatted -a axis\,level\,c\,c\,Z data_down.nc
!spawn cdo yseasmean data_down.nc data_down_sm.nc
!spawn cdo settaxis,2013-02-01,00:00,3months data_down_sm.nc data_down_sm_taxis.nc 
!spawn ncks -A -v levelb $4 data_down_sm_taxis.nc
!spawn ncatted -a edges\,level\,c\,c\,levelb data_down_sm_taxis.nc
!spawn rm -f data_down.nc data_down_sm.nc
!use "data_down_sm_taxis.nc" !d=8 AEOLUS DATA
use $4 !d=8 atmos model down

use "/p/projects/climber3/coumou/Atm_Data/AlexijeData/input/ERA40_synwind.yseasmean.withstratosphere.cdf" !d=9 w'^2
use "/p/projects/climber3/caesar/Era_interim_data/multi_year_seasonal_means/vertical_velocity_monthly_means_daily_means_sm.nc" !d=10 vertical velocity w


set memory/size=10000

if `"($PROGRAM_NAME)" eq "Ferret"` then
  set mode metafile: wind.plt
endif

!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!! Make the plots! !!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!

!DEF AXIS/Z/UNITS=mbar/BOUNDS/DEPTH zax = {24,325,417,535,687,882}, {0,287, 287,368, 368,472, 472,607, 607,779, 779,1000}

set win 1
! u wind surface
set viewport UR
shade/nolabels/levels=20 U_SF[d=1,l=$1]; go land
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Longitude"
!---ylabel-------------------------------------------------
let xpos=`-0.5*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Latitude"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Aeolus $2 mean"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Zonal surface wind (m/s)"
!---timestemp--------------------------------------------------
let xpos=`($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.1 "@ASPIK/L.Caesar<nl>($SESSION_DATE)"
PPL SHASET RESET

set viewport UL
shade/nolabels/levels=(-inf)(`($lev_min)+($lev_del)`,`($lev_max)-($lev_del)`,($lev_del))(inf) U10[d=4,l=$1]; go land
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Longitude"
!---ylabel-------------------------------------------------
let xpos=`-0.5*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Latitude"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Era Interim $2 climatology (1979-2013)"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "10 metre U wind component (m/s)"
PPL SHASET RESET

set viewport LL
shade/nolabels/levels=($lev_text) U_REF[d=3,l=$1]; go land
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Longitude"
!---ylabel-------------------------------------------------
let xpos=`-0.5*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Latitude"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "CM2M $2 mean"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Zonal wind component at 10 m (m/s)"
PPL SHASET RESET

set viewport LR
plot/nolabels U10[d=4,l=$1,i=@ave],U_SF[d=1,l=$1,i=@ave], U_REF[d=3,l=$1,i=@ave]
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.5*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "U wind"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "$2 Zonal Means"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Zonal wind at approx. 10 m (m/s)"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P1Era Interim"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.85*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P2Aeolus"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.8*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P3CM2M"
PPL SHASET RESET

if `"($PROGRAM_NAME)" eq "PyFerret"` then
!frame/file="wind_1.ps" Sigh. Postscript export is broken with Qt5
frame/file="wind_1.pdf"
spawn pdf2ps wind_1.pdf wind_1.ps
spawn rm wind_1.pdf
endif
!!

set win 2
! v wind surface
set viewport UR
shade/nolabels/levels=20 V_SF[d=1,l=$1]; go land
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Longitude"
!---ylabel-------------------------------------------------
let xpos=`-0.5*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Latitude"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Aeolus $2 mean"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Meridional surface wind (m/s)"
!---timestemp--------------------------------------------------
let xpos=`($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.1 "@ASPIK/L.Caesar<nl>($SESSION_DATE)"
PPL SHASET RESET

set viewport UL
shade/nolabels/levels=(-inf)(`($lev_min)+($lev_del)`,`($lev_max)-($lev_del)`,($lev_del))(inf) V10[d=5,l=$1]; go land
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Longitude"
!---ylabel-------------------------------------------------
let xpos=`-0.5*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Latitude"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Era Interim $2 climatology (1979-2013)"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "10 metre V wind component (m/s)"
PPL SHASET RESET


set viewport LL
shade/nolabels/levels=($lev_text) V_REF[d=3,l=$1]; go land
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Longitude"
!---ylabel-------------------------------------------------
let xpos=`-0.5*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Latitude"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "CM2M $2 mean"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Meridional wind component at 10 m (m/s)"
PPL SHASET RESET

set viewport LR
plot/nolabels V10[d=5,l=$1,i=@ave],V_SF[d=1,l=$1,i=@ave], V_REF[d=3,l=$1,i=@ave]
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.5*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "V wind"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "$2 Zonal Means"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Meridional wind at approx. 10 m (m/s)"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P1Era Interim"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.85*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P2Aeolus"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.8*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P3CM2M"
PPL SHASET RESET

if `"($PROGRAM_NAME)" eq "PyFerret"` then
!frame/file="wind_2.ps" Sigh. Postscript export is broken with Qt5
frame/file="wind_2.pdf"
spawn pdf2ps wind_2.pdf wind_2.ps
spawn rm wind_2.pdf
endif
!!

set win 3
! 3D U wind field
set viewport UR
shade/nolabels/levels=20 U_ZA[d=1,l=$1,Z=1000:287,gz=zax] ! 0 - 10 000 m
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.6*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Aeolus $2 mean"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Zonal-mean zonal wind (m/s)"
!---timestemp--------------------------------------------------
let xpos=`($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.1 "@ASPIK/L.Caesar<nl>($SESSION_DATE)"
PPL SHASET RESET

set viewport UL
shade/nolabels/levels=(-inf)(`($lev_min)+($lev_del)`,`($lev_max)-($lev_del)`,($lev_del))(inf) 'U'[d=6,i=@ave,l=$1,k=16:37] ! 1000 - 200 mb 
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.6*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Era Interim $2 climatology (1979-2013)"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Zonal-mean U component wind velocity (m/s)"
PPL SHASET RESET

!set viewport LL
!shade/nolabels/levels=(-10,40,2) U_REF[d=3,l=60:90@ave,i=@ave]; go land
!!---xlabel-------------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`-0.5*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!!---ylabel-------------------------------------------------
!let xpos=`-0.5*($PPL$XORG)`
!let ypos=`0.5*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,0,90,.15 "Height"
!!---titel 1 - above--------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "CM2M coarse bling March average"
!!---titel 2 - below--------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`-0.8*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "Zonal-mean of zonal wind component ar 10 m (m/s)"
!PPL SHASET RESET

set viewport LR
plot/nolabels 'U'[d=6,l=$1,i=@ave,k=16:37@ave], U_ZA[d=1,l=$1,Z=1000:287@ave,gz=zax] 
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.6*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Zonal-mean zonal wind"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Vertical Means"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "$2 Zonal-mean of zonal wind component (m/s)"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P1Era Interim"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.85*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P2Aeolus"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.8*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,-1,1,.1 "@P3CM2M"
PPL SHASET RESET

if `"($PROGRAM_NAME)" eq "PyFerret"` then
!frame/file="wind_3.ps" Sigh. Postscript export is broken with Qt5
frame/file="wind_3.pdf"
spawn pdf2ps wind_3.pdf wind_3.ps
spawn rm wind_3.pdf
endif
!!

set win 4
! 3D V wind field
set viewport UR
shade/nolabels/levels=20 V_ZA[d=1,l=$1,Z=1000:287,gz=zax] ! 0 - 10 000 m
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.6*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Aeolus $2 mean"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Zonal-mean meridional wind (m/s)"
!---timestemp--------------------------------------------------
let xpos=`($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.1 "@ASPIK/L.Caesar<nl>($SESSION_DATE)"
PPL SHASET RESET

set viewport UL
shade/nolabels/levels=(-inf)(`($lev_min)+($lev_del)`,`($lev_max)-($lev_del)`,($lev_del))(inf) 'V'[d=7,i=@ave,l=$1,k=16:37] ! 1000 - 200 mb 
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.6*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Era Interim $2 climatology (1979-2013)"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Zonal-mean V component wind velocity (m/s)"
PPL SHASET RESET

!set viewport LL
!shade/nolabels/levels=(-3,3,0.2) V_REF[d=3,l=60:90@ave,i=@ave]; go land
!!---xlabel-------------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`-0.5*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!!---ylabel-------------------------------------------------
!let xpos=`-0.5*($PPL$XORG)`
!let ypos=`0.5*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,0,90,.15 "Height"
!!---titel 1 - above--------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "CM2M coarse bling March average"
!!---titel 2 - below--------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`-0.8*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "Zonal-mean of meridional wind component ar 10 m (m/s)"
!PPL SHASET RESET

set viewport LR
plot/nolabels 'V'[d=7,l=$1,i=@ave,k=16:37@ave], V_ZA[d=1,l=$1,Z=1000:287@ave,gz=zax] 
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.8*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Zonal-mean meridional wind"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Vertical Means"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "$2 Zonal-mean of meridional wind component (m/s)"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P1Era Interim"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.85*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P2Aeolus"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.8*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,-1,1,.1 "@P3CM2M"
PPL SHASET RESET

if `"($PROGRAM_NAME)" eq "PyFerret"` then
!frame/file="wind_4.ps" Sigh. Postscript export is broken with Qt5
frame/file="wind_4.pdf"
spawn pdf2ps wind_4.pdf wind_4.ps
spawn rm wind_4.pdf
endif
!!

let W2=-1*WW/(9.81*1.3*Z)*1000

set win 5
! Synoptic component vertical velocity squared: w'^2
set viewport UR
shade/nolabels/levels=20 WW[d=1,l=$1,Z=1000:287,gz=zax,i=@ave] ! 0 - 10 000 m
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.6*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Aeolus $2 mean"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Vertical synoptic kinetic energy (m^2/s^2)"
!---timestemp--------------------------------------------------
let xpos=`($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.1 "@ASPIK/L.Caesar<nl>($SESSION_DATE)"
PPL SHASET RESET

set viewport UL
shade/nolabels/levels=(-inf)(`($lev_min)+($lev_del)`,`($lev_max)-($lev_del)`,($lev_del))(inf) W2[d=9,i=@ave,l=$1]^2 ! 1000 - 200 mb 
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.6*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Era 40 $2 climatology (1976-2001)"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Vertical synoptic kinetic energy (m^2/s^2)"
PPL SHASET RESET

!set viewport LL
!shade/nolabels/levels=20 W2[d=9,i=@ave,l=$1]^2 ! 1000 - 200 mb 
!!---xlabel-------------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`-0.5*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!!---ylabel-------------------------------------------------
!let xpos=`-0.6*($PPL$XORG)`
!let ypos=`0.5*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!!---titel 1 - above--------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "Era 40 $2 climatology (1976-2001)/DIFFERENT COLORBAR"
!!---titel 2 - below--------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`-0.8*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "Vertical synoptic kinetic energy (m^2/s^2)"
!PPL SHASET RESET

set viewport LR
plot/nolabels W2[d=9,l=$1,i=@ave,k=@ave]^2, WW[d=1,l=$1,Z=1000:287@ave,gz=zax,i=@ave] 
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.8*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Vertical-mean energy"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Vertical Means"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "$2 Vertical-mean of vert. synoptic kinetic energy (m^2/s^2)"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P1Era 40"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.85*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P2Aeolus"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.8*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,-1,1,.1 "@P3CM2M"
PPL SHASET RESET

if `"($PROGRAM_NAME)" eq "PyFerret"` then
!frame/file="wind_5.ps" Sigh. Postscript export is broken with Qt5
frame/file="wind_5.pdf"
spawn pdf2ps wind_5.pdf wind_5.ps
spawn rm wind_5.pdf
endif
!!

set win 6
! Zonal synoptic kinetic energy: uu (but also u'^2)
set viewport UR
shade/nolabels/levels=20 UU[d=1,l=$1,Z=1000:287,gz=zax,i=@ave] ! 0 - 10 000 m
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.6*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Aeolus $2 mean"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Zonal synoptic kinetic energy (m^2/s^2)"
!---timestemp--------------------------------------------------
let xpos=`($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.1 "@ASPIK/L.Caesar<nl>($SESSION_DATE)"
PPL SHASET RESET

set viewport UL
shade/nolabels/levels=(-inf)(`($lev_min)+($lev_del)`,`($lev_max)-($lev_del)`,($lev_del))(inf) UU[d=9,i=@ave,l=$1]^2 ! 1000 - 200 mb 
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.6*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Era 40 $2 climatology (1976-2001)"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Zonal synoptic kinetic energy (m^2/s^2)"
PPL SHASET RESET

!set viewport LL
!shade/nolabels/levels=20 UU[d=9,i=@ave,l=$1]^2 ! 1000 - 200 mb 
!!---xlabel-------------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`-0.5*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!!---ylabel-------------------------------------------------
!let xpos=`-0.6*($PPL$XORG)`
!let ypos=`0.5*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!!---titel 1 - above--------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "Era 40 $2 climatology (1976-2001)/DIFFERENT COLORBAR"
!!---titel 2 - below--------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`-0.8*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "Zonal synoptic kinetic energy (m^2/s^2)"
!PPL SHASET RESET

set viewport LR
plot/nolabels UU[d=9,l=$1,i=@ave,k=@ave]^2, UU[d=1,l=$1,Z=1000:287@ave,gz=zax,i=@ave] 
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.8*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Vertical-mean energy"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Vertical Means"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "$2 Vertical-mean of zonal sysnoptic kinetic energy (m^2/s^2)"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P1Era 40"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.85*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P2Aeolus"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.8*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,-1,1,.1 "@P3CM2M"
PPL SHASET RESET

if `"($PROGRAM_NAME)" eq "PyFerret"` then
!frame/file="wind_6.ps" Sigh. Postscript export is broken with Qt5
frame/file="wind_6.pdf"
spawn pdf2ps wind_6.pdf wind_6.ps
spawn rm wind_6.pdf
endif
!!

set win 7
! Meridional synoptic kinetic energy: vv (but also v'^2)
set viewport UR
shade/nolabels/levels=20 VV[d=1,l=$1,Z=1000:287,gz=zax,i=@ave] ! 0 - 10 000 m
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.6*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Aeolus $2 mean"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Meridional synoptic kinetic energy (m^2/s^2)"
!---timestemp--------------------------------------------------
let xpos=`($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.1 "@ASPIK/L.Caesar<nl>($SESSION_DATE)"
PPL SHASET RESET

set viewport UL
shade/nolabels/levels=(-inf)(`($lev_min)+($lev_del)`,`($lev_max)-($lev_del)`,($lev_del))(inf) VV[d=9,i=@ave,l=$1]^2 ! 1000 - 200 mb 
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.6*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Era 40 $2 climatology (1976-2001)"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Meridional synoptic kinetic energy (m^2/s^2)"
PPL SHASET RESET

!set viewport LL
!shade/nolabels/levels=20 VV[d=9,i=@ave,l=$1]^2 ! 1000 - 200 mb 
!!---xlabel-------------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`-0.5*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!!---ylabel-------------------------------------------------
!let xpos=`-0.6*($PPL$XORG)`
!let ypos=`0.5*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!!---titel 1 - above--------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "Era 40 $2 climatology (1976-2001)/DIFFERENT COLORBAR"
!!---titel 2 - below--------------------------------------------
!let xpos=`0.5*($PPL$XLEN)`
!let ypos=`-0.8*($PPL$YORG)`
!label/nouser `xpos`,`ypos`,0,0,.15 "Meridional synoptic kinetic energy (m^2/s^2)"
!PPL SHASET RESET

set viewport LR
plot/nolabels VV[d=9,l=$1,i=@ave,k=@ave]^2, VV[d=1,l=$1,Z=1000:287@ave,gz=zax,i=@ave] 
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.8*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Vertical-mean energy"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Vertical Means"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "$2 Vertical-mean of meridional sysnoptic kinetic energy (m^2/s^2)"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P1Era 40"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.85*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P2Aeolus"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.8*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,-1,1,.1 "@P3CM2M"
PPL SHASET RESET

if `"($PROGRAM_NAME)" eq "PyFerret"` then
!frame/file="wind_7.ps" Sigh. Postscript export is broken with Qt5
frame/file="wind_7.pdf"
spawn pdf2ps wind_7.pdf wind_7.ps
spawn rm wind_7.pdf
endif
!!

let W_m=-1*W/(9.81*1.3*Z)*1000

set win 8
! 3D W wind field
set viewport UR
shade/nolabels/levels=20 W[d=1,l=$1,i=@ave,Z=1000:287,gz=zax] ! 0 - 10 000 m
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.6*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Aeolus $2 mean"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Zonal-mean vertical wind (m/s)"
!---timestemp--------------------------------------------------
let xpos=`($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.1 "@ASPIK/L.Caesar<nl>($SESSION_DATE)"
PPL SHASET RESET

set viewport UL
shade/nolabels/levels=(-inf)(`($lev_min)+($lev_del)`,`($lev_max)-($lev_del)`,($lev_del))(inf) W_m[d=10,i=@ave,l=$1,k=16:37] ! 1000 - 200 mb 
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.6*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Era Interim $2 climatology (1979-2014)"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Zonal-mean W component wind velocity (m/s)"
PPL SHASET RESET

set viewport LL
shade/nolabels/levels=20 W_m[d=10,i=@ave,l=$1] ! 1000 - 200 mb 
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.6*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Height (mbar)"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Era Interim $2 climatology (1979-2014) - diff colorbar, h"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Zonal-mean W component wind velocity (m/s)"
PPL SHASET RESET

set viewport LR
plot/nolabels W_m[d=10,l=$1,i=@ave,k=16:37@ave], W[d=1,l=$1,i=@ave,Z=1000:287@ave,gz=zax] 
!---xlabel-------------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.5*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Latitude"
!---ylabel-------------------------------------------------
let xpos=`-0.8*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Zonal-mean vertical wind"
!---titel 1 - above--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`($PPL$YLEN)+0.2*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Vertical Means"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "$2 Zonal-mean of vertical wind component (m/s)"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P1Era Interim"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.85*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P2Aeolus"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.8*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,-1,1,.1 "@P3CM2M"
PPL SHASET RESET

if `"($PROGRAM_NAME)" eq "PyFerret"` then
!frame/file="wind_8.ps" Sigh. Postscript export is broken with Qt5
frame/file="wind_8.pdf"
spawn pdf2ps wind_8.pdf wind_8.ps
spawn rm wind_8.pdf
endif
!!


!!!!! The end !!!!!
if `"($PROGRAM_NAME)" eq "Ferret"` then
  CANCEL MODE METAFILE
  spawn Fprint -l cps -o wind_1.ps -R wind.plt
  spawn Fprint -l cps -o wind.ps -R wind.plt.~*~
  spawn rm wind.plt*
  spawn cat wind_1.ps >> wind.ps
else
  spawn cat wind_?.ps > wind.ps
endif
spawn rm wind_?.ps
spawn ps2pdf wind.ps wind_seas_$2.pdf
spawn rm ferret.jnl
spawn rm wind.ps
cancel data/all
!spawn rm -f data_sm_taxis.nc data_down_sm_taxis.nc