! Normalised taylor diagrams of global average opal export
set mem/size=200

! Load observational dataset
go "/gxfs_home/geomar/smomw258/Scripts/David_scripts/readhonjotraps.jnl"
let obs = hfsi[d=1]/1000/365/24/60/60 !mol Si/m2/s

! Load new UVic model
use "tavg.01775.01.01.nc"
let uvic=o_oplexp[l=25,d=2]
let uvic1=samplexy(o_oplexp[l=25,d=2,k=11]*350,hlon[d=1],hlat[d=1]) ! mol Si/m2/y

canc mode cal
 
! Load CESM
use "$HOME/Data/cmip6_data/exp/expsi_Oyr_CESM2_historical_r11i1p1f1_gr_2000-2014.nc"
let cesm1 = expsi[l=15,d=3,gxyz=uvic@nrst]
let cesm = samplexy(cesm1[k=11],hlon[d=1],hlat[d=1])

! Load EC-Earth3-CC
use "$HOME/Data/cmip6_data/exp/expsi_Oyr_EC-Earth3-CC_esm-hist_r1i1p1f1_gn_2014-2014.nc"
let map3 = curv_to_rect_map(longitude[d=4],latitude[d=4],uvic,1)
let ec1 = curv_to_rect(expsi[d=4,l=1],map3)
let ec2 = ec1[gxyz=uvic@nrst]
let ec = samplexy(ec2[k=11],hlon[d=1],hlat[d=1])

! Load GFDL
use "$HOME/Data/cmip6_data/exp/expsi_Oyr_GFDL-ESM4_historical_r1i1p1f1_gr_1950-2014.nc"
let gfdl1 = expsi[l=65,gxyz=uvic@nrst,d=5]
let gfdl = samplexy(gfdl1[k=11],hlon[d=1],hlat[d=1])

! Load IPSL
use "$HOME/Data/cmip6_data/exp/expsi_Oyr_IPSL-CM6A-LR_historical_r32i1p1f1_gn_1850-2014.nc"
let map2 = curv_to_rect_map(nav_lon[d=6],nav_lat[d=6],uvic,1)
let ipsl1 = curv_to_rect(expsi[d=6,l=165],map2)
let ipsl2 = ipsl1[gxyz=uvic@nrst]
let ipsl = samplexy(ipsl2[k=11],hlon[d=1],hlat[d=1])

! Load MPI-ESM-LR
use "$HOME/Data/cmip6_data/exp/MPI_remapped.nc"
let mpi1 = expsi[l=5,gxyz=uvic@nrst,d=7]
let mpi = samplexy(mpi1[k=11],hlon[d=1],hlat[d=1])

! Load CMCC-ESM2
use "$HOME/Data/cmip6_data/exp/expsi_Oyr_CMCC-ESM2_historical_r1i1p1f1_gn_1850-2014.nc"
let map1 = curv_to_rect_map(longitude[d=8],latitude[d=8],uvic,1)
let cmcc1 = curv_to_rect(expsi[d=8,l=165],map1)
let cmcc2 = cmcc1[gxyz=uvic@nrst]
let cmcc = samplexy(cmcc2[k=11],hlon[d=1],hlat[d=1])

let p = obs
let q = uvic1
go regressxyz2

go prepareFigure.jnl taylor_opal_exp
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Plot normalized taylor diagram
set win 1
go taylor_frame 1.6 0.2 " " "Normalised Standard Deviation"
go taylor_rgraticule 1.0 /DASH=(0.2,0.1,0.2,0.1)/COLOR=BLACK/THICK=1
go taylor_polymark 1.0 1.0 /line/fill/nokey/color=black/pal=black r circle 2.0
go taylor_agraticule 1.6

go taylor_polymark stdr r /line/fill/nokey/ov/color=black/pal=red r circle 2.0

let p = obs
let q = gfdl
go regressxyz2
go taylor_polymark stdr r /line/fill/nokey/ov/color=black/pal=green r star 2.0

let p = obs
let q = cesm
go regressxyz2
go taylor_polymark stdr r /line/fill/nokey/ov/color=black/pal=blue_light r square 2.0

let p = obs
let q = ec
go regressxyz2
go taylor_polymark stdr r /line/fill/nokey/ov/color=black/pal=blue r diamond 2.0

let p = obs
let q = ipsl
go regressxyz2
go taylor_polymark stdr r /line/fill/nokey/ov/color=black/pal=purple_dark r hexagon 2.0

let p = obs
let q = cmcc
go regressxyz2
go taylor_polymark stdr r /line/fill/nokey/ov/color=black/pal=orange_dark r triangle 2.0


let p = obs
let q = mpi
go regressxyz2
go taylor_polymark stdr r /line/fill/nokey/ov/color=black/pal=pink r pentagon 2.0

go finalizeFigure_psthicken