! Comparing Era Interim Data to Coarse Bling and Aeolus output: convective cloud amount, low- medium, high-level cloud cover
!
! 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 Clouds-data.nc
!spawn ncatted -a axis\,level\,c\,c\,Z Clouds-data.nc
!spawn cdo yseasmean Clouds-data.nc Clouds-data_sm.nc
!spawn cdo settaxis,2013-02-01,12:00,3months Clouds-data_sm.nc Clouds-data_sm_taxis.nc 
!spawn ncks -A -v levelb $3 Clouds-data_sm_taxis.nc
!spawn ncatted -a edges\,level\,c\,c\,levelb Clouds-data_sm_taxis.nc
!spawn rm -f Clouds-data.nc Clouds-data_sm.nc
!use "Clouds-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/ISSCP_data/deep_convective_cloud_amount_sm.nc" !d=4
use "/p/projects/climber3/caesar/Era_interim_data/noleap_data/multi_year_season_mean/cloud_cover_low_level_sm.nc" !d=5
use "/p/projects/climber3/caesar/Era_interim_data/noleap_data/multi_year_season_mean/cloud_cover_medium_level_sm.nc" !d=6
use "/p/projects/climber3/caesar/Era_interim_data/noleap_data/multi_year_season_mean/cloud_cover_high_level_sm.nc" !d=7


set memory/size=10000

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

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

set win 1
! Deep convective cloud amount/cover
set viewport UR
shade/nolabels/levels=20 CLAM_4[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 "Deep convective cloud amout (Percent)"
!---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) CA_DC[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 "ISSCP $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 "Deep convective cloud amount (Percent)"
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 CA_DC[d=4,l=$1,i=@ave], CLAM_4[d=1,l=$1,i=@ave]*100
!---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 "Cloud amount"
!---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 "Deep convective cloud amount (Percent)"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,1,.1 "@P1ISSCP"
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="clouds_1.ps" Sigh. Postscript export is broken with Qt5
frame/file="clouds_1.pdf"
spawn pdf2ps clouds_1.pdf clouds_1.ps
spawn rm clouds_1.pdf
endif
!!

set win 2
! Low level cloud cover
set viewport UR
shade/nolabels/levels=20 CLAM_1[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 "Low level cloud amount (Percent)"
!---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) LCC[d=5,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 "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 "Low level cloud amount (Percent)"
PPL SHASET RESET

set viewport LL
shade/nolabels/levels=($lev_text) LOW_CLD_AMT[d=2,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 "Low level cloud amount (Percent)"
PPL SHASET RESET

set viewport LR
plot/nolabels LCC[d=5,l=$1,i=@ave]*100, CLAM_1[d=1,l=$1,i=@ave]*100,LOW_CLD_AMT[d=2,k=10,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 "Cloud amount"
!---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 "Low level cloud amount (Percent)"
!---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="clouds_2.ps" Sigh. Postscript export is broken with Qt5
frame/file="clouds_2.pdf"
spawn pdf2ps clouds_2.pdf clouds_2.ps
spawn rm clouds_2.pdf
endif
!!

set win 3
! Medium level cloud cover
set viewport UR
shade/nolabels/levels=20 CLAM_2[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 "Medium level cloud amount (Percent)"
!---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) MCC[d=6,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 "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 "Medium level cloud amount (Percent)"
PPL SHASET RESET

set viewport LL
shade/nolabels/levels=($lev_text) MID_CLD_AMT[d=2,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 "Medium level cloud amount (Percent)"
PPL SHASET RESET

set viewport LR
plot/nolabels MCC[d=6,l=$1,i=@ave]*100, CLAM_2[d=1,l=$1,i=@ave]*100,MID_CLD_AMT[d=2,k=10,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 "Cloud amount"
!---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 "Medium level cloud amount (Percent)"
!---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="clouds_3.ps" Sigh. Postscript export is broken with Qt5
frame/file="clouds_3.pdf"
spawn pdf2ps clouds_3.pdf clouds_3.ps
spawn rm clouds_3.pdf
endif
!!

set win 4
! High level cloud cover
set viewport UR
shade/nolabels/levels=20 CLAM_3[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 "High level cloud amount (Percent)"
!---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) HCC[d=7,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 "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 "High level cloud amount (Percent)"
PPL SHASET RESET

set viewport LL
shade/nolabels/levels=($lev_text) HIGH_CLD_AMT[d=2,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 "High level cloud amount (Percent)"
PPL SHASET RESET

set viewport LR
plot/nolabels HCC[d=7,l=$1,i=@ave]*100, CLAM_3[d=1,l=$1,i=@ave]*100,HIGH_CLD_AMT[d=2,k=10,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 "Cloud amount"
!---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 "High level cloud amount (Percent)"
!---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="clouds_4.ps" Sigh. Postscript export is broken with Qt5
frame/file="clouds_4.pdf"
spawn pdf2ps clouds_4.pdf clouds_4.ps
spawn rm clouds_4.pdf
endif

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