! Comparing Era Interim Data to Coarse Bling and Aeolus output: specific humidity surface, lapse rate
!
! 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
! $4 = dataset flux

cancel data/all

!spawn cdo -seldate\,0052-01-01T00:00:00\,0104-12-31T23:59:59 $3 Humidity-data.nc
!spawn ncatted -a axis\,level\,c\,c\,Z Humidity-data.nc ! set the axis as z-axis
!spawn cdo yseasmean Humidity-data.nc Humidity-data_sm.nc
!spawn cdo settaxis,2013-02-01,12:00,3months Humidity-data_sm.nc Humidity-data_sm_taxis.nc 
!spawn ncks -A -v levelb $3 Humidity-data_sm_taxis.nc
!spawn ncatted -a edges\,level\,c\,c\,levelb Humidity-data_sm_taxis.nc
!spawn rm -f Humidity-data.nc Humidity-data_sm.nc
!use "Humidity-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/specific_humidity_1000mb_sm.nc" !d=4
use "/p/projects/climber3/caesar/lr_sm.cdf" !d=5
use "/p/projects/climber3/caesar/Era_interim_data/multi_year_seasonal_means/specific_humidity_pressure_levels_daily_mean_monthly_means_sm.nc" !d=6
use "/p/projects/climber3/caesar/Era_interim_data/multi_year_seasonal_means/relative_humidity_pressure_levels_sm.nc" !d=7

!spawn cdo -seldate\,0052-01-01T00:00:00\,0104-12-31T23:59:59 $4 Humidity-data_flux.nc
!spawn ncatted -a axis\,level\,c\,c\,Z Humidity-data_flux.nc
!spawn cdo yseasmean Humidity-data_flux.nc Humidity-data_flux_sm.nc
!spawn cdo settaxis,2013-02-01,12:00,3months Humidity-data_flux_sm.nc Humidity-data_flux_sm_taxis.nc 
!spawn ncks -A -v levelb $4 Humidity-data_flux_sm_taxis.nc
!spawn ncatted -a edges\,level\,c\,c\,levelb Humidity-data_flux_sm_taxis.nc
!spawn rm -f Humidity-data_flux.nc Humidity-data_flux_sm.nc
!use "Humidity-data_flux_sm_taxis.nc" !d=8
use $4 !d=8

set memory/size=10000

if `"($PROGRAM_NAME)" eq "Ferret"` then
  set mode metafile: humidity.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
! Specific humidity at the surface
set viewport UR
shade/nolabels/levels=20 QS[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 "Specific humidity at surface (no topography) (kg/kg)"
!---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)-3*($lev_del)`,($lev_del))(inf) Q[d=4,l=$1]; go land ! Watch out: times 3 due to ferret correction
!---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 "Specific humidity at 1000mb (kg/kg)"
PPL SHASET RESET

set viewport LL
shade/nolabels/levels=($lev_text) Q_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 "Specific humidity at surface (no topography) (kg/kg)"
PPL SHASET RESET

set viewport LR
plot/nolabels Q[d=4,l=$1,i=@ave], QS[d=1,l=$1,i=@ave], Q_REF[d=3,i=@ave,l=$1]
!---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.7*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Specific humidity"
!---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 "Specific humidity at approximately 1000mb (kg/kg)"
!---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="humidity_1.ps" Sigh. Postscript export is broken with Qt5
frame/file="humidity_1.pdf"
spawn pdf2ps humidity_1.pdf humidity_1.ps
spawn rm humidity_1.pdf
endif
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

set win 2
! Specific humidity 3d
set viewport UR
shade/nolabels/levels=20 Q3D[d=1,i=@ave,Z=1000:287,gz=zax,l=$1]; 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 (mb)"
!---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 specific humidity (kg/kg)"
!---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))(inf) Q[d=6,i=@ave,Z=1000:287,l=$1]; 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 (mb)"
!---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 specific humidity (kg/kg)"
PPL SHASET RESET

set viewport LL
shade/nolabels/levels=($lev_text) SPHUM[d=2,i=@ave,Z=1000:287,l=$1]; 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 (mb)"
!---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 mean specific humidity (kg/kg)"
PPL SHASET RESET

set viewport LR
plot/nolabels  Q[d=6,l=$1,i=@ave,Z=1000:287@ave], Q3D[d=1,l=$1,i=@ave,Z=1000:287@ave,gz=zax], SPHUM[d=2,i=@ave,Z=1000:287@ave,l=$1]
!---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.7*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Specific humidity"
!---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 "1000mb - 287 mb vertical mean specific humidity (kg/kg)"
!---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="humidity_2.ps" Sigh. Postscript export is broken with Qt5
frame/file="humidity_2.pdf"
spawn pdf2ps humidity_2.pdf humidity_2.ps
spawn rm humidity_2.pdf
endif
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

set win 3
! Relative humidity at the surface
set viewport UR
shade/nolabels/levels=20 QS[d=1,l=$1]/Q_SAT_TSA[d=1,l=$1]*100; 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 "Relative humidity at the surface (@AS%)"
!---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) R[d=7,Z=1000,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 "Relative humidity at 1000mb (@AS%)"
PPL SHASET RESET

set viewport LL
shade/nolabels/levels=($lev_text) RH[d=2,l=$1,Z=1000]; 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 "Relative humidity at 1000mb (@AS%)"
PPL SHASET RESET

set viewport LR
plot/nolabels  R[d=7,Z=1000,l=$1,i=@ave], QS[d=1,l=$1,i=@ave]/Q_SAT_TSA[d=1,l=$1,i=@ave]*100, RH[d=2,i=@ave,l=$1,Z=1000]
!---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.7*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Relative humidity"
!---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 "Relative humidity at approximately 1000mb (@AS%)"
!---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="humidity_3.ps" Sigh. Postscript export is broken with Qt5
frame/file="humidity_3.pdf"
spawn pdf2ps humidity_3.pdf humidity_3.ps
spawn rm humidity_3.pdf
endif
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

set win 4
! Relative humidity at lowest atmospheric level
set viewport UR
shade/nolabels/levels=20 RHUM[d=1,l=$1,Z=1000:779@ave]*100; 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 "Relative humidity at approx. 1000m (@AS%)"
!---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) R[d=7,Z=1000:779@ave,l=$1]; go land ! Watch out: times 3 due to ferret correction
!---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 "Relative humidity at 1000mb:779mb ave (@AS%)"
PPL SHASET RESET

set viewport LL
shade/nolabels/levels=($lev_text) RH[d=2,l=$1,Z=1000:779@ave]; 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 "Relative humidity at 1000mb:779mb ave (@AS%)"
PPL SHASET RESET

set viewport LR
plot/nolabels  R[d=7,Z=1000:779@ave,l=$1,i=@ave], RHUM[d=1,l=$1,i=@ave,Z=1000:779@ave]*100, RH[d=2,i=@ave,l=$1,Z=1000:779@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.7*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Relative humidity"
!---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 "Relative humidity at approximately 1000m (@AS%)"
!---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="humidity_4.ps" Sigh. Postscript export is broken with Qt5
frame/file="humidity_4.pdf"
spawn pdf2ps humidity_4.pdf humidity_4.ps
spawn rm humidity_4.pdf
endif
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

set win 5
! Relative humidity 3d
set viewport UR
shade/nolabels/levels=(0,100,5) RHUM[d=1,i=@ave,Z=1000:287,gz=zax,l=$1]*100; 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 (mb)"
!---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 relative humidity (@AS%)"
!---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) R[d=7,i=@ave,Z=1000:287,l=$1]; go land ! Watch out: times 3 due to ferret correction
!---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 (mb)"
!---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 relative humidity (@AS%)"
PPL SHASET RESET

set viewport LL
shade/nolabels/levels=($lev_text) RH[d=2,i=@ave,Z=1000:287,l=$1]; 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 (mb)"
!---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 mean relative humidity (@AS%)"
PPL SHASET RESET

set viewport LR
plot/nolabels  R[d=7,l=$1,i=@ave,Z=1000:287@ave],RHUM[d=1,l=$1,i=@ave,Z=1000:287@ave,gz=zax]*100, RH[d=2,i=@ave,Z=1000:287@ave,l=$1]
!---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.7*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Relative humidity"
!---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 "1000mb - 287 mb vertical mean relative humidity (@AS%)"
!---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="humidity_5.ps" Sigh. Postscript export is broken with Qt5
frame/file="humidity_5.pdf"
spawn pdf2ps humidity_5.pdf humidity_5.ps
spawn rm humidity_5.pdf
endif
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

set win 6
! Lapse rate
set viewport UR
shade/nolabels/levels=20 LR[d=1,Z=1000,l=$1]*1000; 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 "Lapse rate (K/km)"
!---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) LR[d=5,l=$1]*1000; 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 "Lapse rate (K/km)"
PPL SHASET RESET

!set viewport LL
!shade/nolabels/levels=($lev_text) SHFLX[d=3,k=10,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 "Sensible heat (W/m^2)"
!PPL SHASET RESET

set viewport LR
plot/nolabels LR[d=5,l=$1,i=@ave]*1000, LR[d=1,k=3,l=$1,i=@ave]*1000
!---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 "Lapse rate"
!---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 "Lapse rate (K/km)"
!---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="humidity_6.ps" Sigh. Postscript export is broken with Qt5
frame/file="humidity_6.pdf"
spawn pdf2ps humidity_6.pdf humidity_6.ps
spawn rm humidity_6.pdf
endif
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!! The end !!!!!
if `"($PROGRAM_NAME)" eq "Ferret"` then
  CANCEL MODE METAFILE
  spawn Fprint -p landscape -l cps -o humidity_1.ps -R humidity.plt
  spawn Fprint -p landscape -l cps -o humidity.ps -R humidity.plt.~*~
  !spawn Fprint -l cps -o humidity_1.ps -R humidity.plt
  !spawn Fprint -l cps -o humidity.ps -R humidity.plt.~*~
  spawn rm humidity.plt*
  spawn cat humidity_1.ps >> humidity.ps
else
  spawn cat humidity_?.ps > humidity.ps
endif
spawn rm humidity_?.ps
spawn ps2pdf humidity.ps humidity_seas_$2.pdf
spawn rm humidity.ps
spawn rm ferret.jnl
cancel data/all
!spawn rm -f Humidity-data_sm_taxis.nc Humidity-data_flux_sm_taxis.nc