! Comparing Era Interim Data to Coarse Bling and Aeolus output: lwr up toa, lwr surface down, lwr surface up, lwr surface net
!
! 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

cancel data/all

!spawn cdo -seldate\,0052-01-01T00:00:00\,0104-12-31T23:59:59 $3 LWR_Rad-data.nc
!spawn ncatted -a axis\,level\,c\,c\,Z LWR_Rad-data.nc
!spawn cdo yseasmean LWR_Rad-data.nc LWR_Rad-data_sm.nc
!spawn cdo settaxis,2013-02-28,00:00,3months LWR_Rad-data_sm.nc LWR_Rad-data_sm_taxis.nc 
!spawn ncks -A -v levelb $3 LWR_Rad-data_sm_taxis.nc
!spawn ncatted -a edges\,level\,c\,c\,levelb LWR_Rad-data_sm_taxis.nc
!spawn rm -f LWR_Rad-data.nc LWR_Rad-data_sm.nc
!use "LWR_Rad-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/top_net_ther_rad_sm.nc" !d=4
use "/p/projects/climber3/caesar/Era_interim_data/noleap_data/multi_year_season_mean/ther_rad_sf_down_sm.nc" !d=5
use "/p/projects/climber3/caesar/Era_interim_data/noleap_data/multi_year_season_mean/net_ther_rad_sm.nc" !d=6
use "/p/projects/climber3/caesar/Era_interim_data/noleap_data/multi_year_season_mean/top_net_sol_rad_sm.nc" !d=7

set memory/size=10000

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


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

set win 1
! long wave radiation up(=net) at toa
set viewport UR
shade/nolabels/levels=20 LWR_UP_TOA[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 "Long wave radiation at TOA up (W/m^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) TTR[d=4,l=$1]/43200*-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 "TOA net thermal radiation (W/m^2)"
PPL SHASET RESET

set viewport LL
shade/nolabels/levels=($lev_text) (NETRAD_TOA[d=2,l=60:89@ave]-SWTOA[d=2,l=60:90@ave])*-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 "Netto short wave radiation at TOA - <NL> Total netto radiation at TOA (W/m^2)"
PPL SHASET RESET

set viewport LR
plot/nolabels TTR[d=4,l=$1,i=@ave]/43200*-1,LWR_UP_TOA[d=1,l=$1,i=@ave], (NETRAD_TOA[d=2,l=$1,i=@ave]-SWTOA[d=2,l=$1,i=@ave])*-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.5*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Long wave radiation"
!---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 "Long wave radiation at TOA up - <NL> positive upwards (W/m^2)"
!---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="lwr_radiation_1.ps" Sigh. Postscript export is broken with Qt5
frame/file="lwr_radiation_1.pdf"
spawn pdf2ps lwr_radiation_1.pdf lwr_radiation_1.ps
spawn rm lwr_radiation_1.pdf
endif
!!

set win 2
! Long wave radiation down at the surface

set viewport UR
shade/nolabels/levels=20 LWR_DOWN_SF[d=1,l=$1]*-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 "Long wave radiation at surface down (W/m^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) STRD[d=5,l=$1]/43200; 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 "Surface thermal radiation downwards (W/m^2)"
PPL SHASET RESET

set viewport LL
shade/nolabels/levels=($lev_text) LWDN_SFC[d=2,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 "Long wave radiation at surface down (W/m^2)"
PPL SHASET RESET

set viewport LR
plot/nolabels STRD[d=5,l=$1,i=@ave]/43200, LWR_DOWN_SF[d=1,l=$1,i=@ave]*-1, LWDN_SFC[d=2,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 "Long wave radiation"
!---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 "Long wave radiation at surface down - <NL> positive downwards (W/m^2)"
!---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="lwr_radiation_2.ps" Sigh. Postscript export is broken with Qt5
frame/file="lwr_radiation_2.pdf"
spawn pdf2ps lwr_radiation_2.pdf lwr_radiation_2.ps
spawn rm lwr_radiation_2.pdf
endif
!!

set win 3
! Long wave radiation upwards at the surface
set viewport UR
shade/nolabels/levels=20 LWR_DOWN_SF[d=1,l=$1]*-1-LWR_FLUX_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 "Long wave radiation at surface down - <NL> Netto Long wave radiation at surface (W/m^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) STRD[d=5,l=$1]/43200-STR[d=6,l=$1]/43200*-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 "Surface thermal radiation downwards - <NL> Surface net thermal radiation (W/m^2)"
PPL SHASET RESET

set viewport LL
shade/nolabels/levels=($lev_text) LWDN_SFC[d=2,l=$1]-LWSFC[d=2,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 "Long wave radiation at surface down - <NL> Netto long wave radiation at surface (W/m^2)"
PPL SHASET RESET

set viewport LR
plot/nolabels STRD[d=5,l=$1,i=@ave]/43200-STR[d=6,l=$1,i=@ave]/43200*-1, LWR_DOWN_SF[d=1,l=$1,i=@ave]*-1-LWR_FLUX_SF[d=1,l=$1,i=@ave], LWDN_SFC[d=2,l=$1,i=@ave]-LWSFC[d=2,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 "Long wave radiation"
!---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 "Long wave radiation at surface up - positive upwards <NL>  (W/m^2)"
!---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="lwr_radiation_3.ps" Sigh. Postscript export is broken with Qt5
frame/file="lwr_radiation_3.pdf"
spawn pdf2ps lwr_radiation_3.pdf lwr_radiation_3.ps
spawn rm lwr_radiation_3.pdf
endif
!!

set win 4
! Long wave radiation net at surface
set viewport UR
shade/nolabels/levels=20 LWR_FLUX_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 "Netto long wave radiation at the surface (W/m^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) STR[d=6,l=$1]/43200*-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 "Surface net thermal radiation (W/m^2)"
PPL SHASET RESET

set viewport LL
shade/nolabels/levels=($lev_text) LWSFC[d=2,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 "Netto long wave radiation at the surface (W/m^2)"
PPL SHASET RESET

set viewport LR
plot/nolabels STR[d=6,l=$1,i=@ave]/43200*-1,LWR_FLUX_SF[d=1,l=$1,i=@ave], LWSFC[d=2,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 "Long wave radiation"
!---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 "Netto long wave radiation at the surface - <NL> positive upwards (W/m^2)"
!---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="lwr_radiation_4.ps" Sigh. Postscript export is broken with Qt5
frame/file="lwr_radiation_4.pdf"
spawn pdf2ps lwr_radiation_4.pdf lwr_radiation_4.ps
spawn rm lwr_radiation_4.pdf
endif

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