!!! Comparing Era Interim Data to Coarse Bling and Aeolus output: total cloud amount, total precipitation, large-scale precipitation, convective precipitation, cloud
!!!water path, cloud optical depth !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!! 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 ncatted -a axis\,level\,c\,c\,Z $3
spawn cdo -seldate\,0052-01-01T00:00:00\,0104-12-31T23:59:59 $3 data.nc
spawn cdo yseasmean data.nc data_sm.nc
spawn cdo settaxis,2013-02-01,12: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
use "data_sm_taxis.nc"													                          !d=1 AEOLUS DATA
use "/p/projects/climber3/coumou/Atm_Data/CloudData/CloudValidationData.yseasmean.nc"				  				  !d=2
use "/p/projects/climber3/caesar/ISSCP_data/COD_ISCCP_D1_AMPM_19842007.nc"									  !d=3

set memory/size=10000

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

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

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
set win 1
! Total cloud amount/cover
set viewport UR
shade/nolabels/levels=20 CLAM_5[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 "Total 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) TCA[d=2,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 "ISSCP $2 climatology (1979-?)"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Total 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 TCA[d=2,l=$1,i=@ave]*100, CLAM_5[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 (Percent)"
!---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 "Total cloud amount (Percent)"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,0,.1 "@P1ISSCP"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.85*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,0,.1 "@P2Aeolus"
!let xpos=`0.05*($PPL$XLEN)`
!let ypos=`0.8*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,-1,0,.1 "@P3CM2M"
PPL SHASET RESET

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

set win 2
! Precipitation, total
let TOT_PRECIP1 = if TOT_PRECIP gt 0 then TOT_PRECIP else 0
set viewport UR
shade/nolabels/levels=20 PRECTOT[d=1,l=$1]*2.6e6; go land ! going from mm/s to mm/month
!---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 "Total precipitation (mm/month)"
!---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) TOT_PRECIP1[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 "GPCP $2 climatology (1979-2001)" ! Global Precipitation Climatology project
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Total precipitation (mm/month)"
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 TOT_PRECIP1[d=2,l=$1,i=@ave], PRECTOT[d=1,l=$1,i=@ave]*2.6e6
!---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 "Precipitation"
!---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 "Total precipitation (mm/month)"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,0,.1 "@P1Era Interim"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.85*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,0,.1 "@P2Aeolus"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.8*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,-1,0,.1 "@P3CM2M"
PPL SHASET RESET

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

set win 3
! Precipitation, large scale
set viewport UR
shade/nolabels/levels=20 PRECLS[d=1,l=$1]*2.6e6; go land ! going from mm/s to mm/month
!---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 "Large scale precipitation (mm/month)"
!---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) LSC_PRECIP[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 "ERA 40 $2 climatology (1979-2001)"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Large scale precipitation (mm/month)"
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 LSC_PRECIP[d=2,l=$1,i=@ave], PRECLS[d=1,l=$1,i=@ave]*2.6e6
!---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 "Precipitation"
!---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 "Large scale precipitation (mm/month)"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,0,.1 "@P1Era Interim"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.85*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,0,.1 "@P2Aeolus"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.8*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,-1,0,.1 "@P3CM2M"
PPL SHASET RESET

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

set win 4
! Precipitation, convective
set viewport UR
shade/nolabels/levels=20 PRECCONV[d=1,l=$1]*2.6e6; go land ! going from mm/s to mm/month
!---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 "Convective precipitation (mm/month)"
!---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) CON_PRECIP[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 "ERA 40 $2 climatology (1979-2001)"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Convective precipitation (mm/month)"
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 CON_PRECIP[d=2,l=$1,i=@ave], PRECCONV[d=1,l=$1,i=@ave]*2.6e6
!---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 "Precipitation"
!---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 "Convective precipitation (mm/month)"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,0,.1 "@P1Era Interim"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.85*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,0,.1 "@P2Aeolus"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.8*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,-1,0,.1 "@P3CM2M"
PPL SHASET RESET

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

set win 5
! (Total) Cloud water pather ice and liquid
set viewport UR
shade/nolabels/levels=20 CLWATTOT_5[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 "Total clouds water/ice path (kg/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) (CLWP[d=2,l=$1]+CIWP[d=2,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 "CERES $2 climatology (1979-2001)"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Total clouds water/ice path (kg/m^2)"
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 (CLWP[d=2,l=$1,i=@ave]+ CIWP[d=2,l=$1,i=@ave])/1000, CLWATTOT_5[d=1,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 "Clouds water/ice path (kg/m^2)"
!---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 mean total clouds water /ice path (kg/m^2)"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,0,.1 "@P1CERES"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.85*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,0,.1 "@P2Aeolus"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.8*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,-1,0,.1 "@P3CM2M"
PPL SHASET RESET

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

set win 6
! Cloud optical depth
set viewport UR
shade/nolabels/levels=20 T_C[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 "Cloud optical depth ( )"
!---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) A_COD[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 "ISSCP $2 climatology (1984-2007)"
!---titel 2 - below--------------------------------------------
let xpos=`0.5*($PPL$XLEN)`
let ypos=`-0.8*($PPL$YORG)`
label/nouser `xpos`,`ypos`,0,0,.15 "Cloud optical depth ( )"
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 A_COD[d=3,i=@ave,l=$1], T_C[d=1,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.7*($PPL$XORG)`
let ypos=`0.5*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,0,90,.15 "Cloud optical depth ()"
!---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 mean cloud optical depth ( )"
!---key--------------------------------------------------
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.9*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,0,.1 "@P1ISSCP"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.85*($PPL$YLEN)`
label/nouser `xpos`,`ypos`,-1,0,.1 "@P2Aeolus"
let xpos=`0.05*($PPL$XLEN)`
let ypos=`0.8*($PPL$YLEN)`
!label/nouser `xpos`,`ypos`,-1,0,.1 "@P3CM2M"
PPL SHASET RESET

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



!!!!! The end !!!!!
if `"($PROGRAM_NAME)" eq "Ferret"` then
  CANCEL MODE METAFILE
  spawn Fprint -l cps -o clouds2_1.ps -R clouds2.plt
  spawn Fprint -l cps -o clouds2.ps -R clouds2.plt.~*~
  spawn rm clouds2.plt*
  spawn cat clouds2_1.ps >> clouds2.ps
else
  spawn cat clouds2_?.ps > clouds2.ps
endif
spawn rm clouds2_?.ps
spawn ps2pdf clouds2.ps clouds2_seas_$2.pdf
spawn rm clouds2.ps
spawn rm ferret.jnl