! Normalised taylor diagrams of global average Si concentration
set mem/size=100
! Load observational dataset
use "data/O_silica.nc"

! Load new UVic model
use "tavg.01775.01.01.nc"
let uvic = O_sil[d=2,l=25] !year 2014
let obs = o_silica[d=1,gxyz=uvic@nrst]
canc mode cal

!define grid/like=o_sil[d=2]/t=o_sil[d=2] mg 

! Load CESM
use "$HOME/Data/cmip6_data/dSi/si_Oyr_CESM2_historical_r11i1p1f1_gr_2000-2014.nc"
let cesm = si[l=15,d=3,gxyz=uvic@nrst]

! Load EC-Earth3-CC
use "$HOME/Data/cmip6_data/dSi/si_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(si[d=4,l=1],map3)
let ec = ec1[gxyz=uvic@nrst]

! Load GFDL
use "$HOME/Data/cmip6_data/dSi/si_Oyr_GFDL-ESM4_historical_r1i1p1f1_gr_1950-2014.nc"
let gfdl = si[l=65,gxyz=uvic@nrst,d=5]

! Load IPSL
!use "$HOME/Data/cmip6_data/dSi/si_Oyr_IPSL-CM5A2-INCA_historical_r1i1p1f1_gn_1850-2014.nc"
!let ipsl1 = si[l=165,gxyz=uvic@nrst,d=6]

use "$HOME/Data/cmip6_data/dSi/si_Oyr_IPSL-CM6A-LR_historical_r22i1p1f1_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(si[d=6,l=165],map2)
let ipsl = ipsl1[gxyz=uvic@nrst]

! Load MPI-ESM-LR
use "$HOME/Data/cmip6_data/dSi/MPI_remapped.nc"
let mpi = si[l=5,gxyz=uvic@nrst,d=7]

! Load CMCC-ESM2
use "$HOME/Data/cmip6_data/dSi/si_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(si[d=8,l=165],map1)
let cmcc = cmcc1[gxyz=uvic@nrst]

! Load KIOST-ESM!not using- units issues!
! Load NorESM
use "$HOME/Data/cmip6_data/dSi/si_Oyr_NorESM2-LM_historical_r2i1p1f1_gr_2010-2014.nc"
let map4 = curv_to_rect_map(longitude[d=9],latitude[d=9],uvic,1)
let noresm1 = curv_to_rect(si[d=9,l=5],map4)
let noresm = noresm1[gxyz=uvic@nrst]

! Load NorCPM
!use "$HOME/Data/cmip6_data/dSi/si_Oyr_NorCPM1_historical_r10i1p1f1_gr_2010-2014.nc"
!let map4 = curv_to_rect_map(longitude[d=9],latitude[d=9],uvic,1)
!let noresm1 = curv_to_rect(si[d=9,l=5],map4)
!let noresm = noresm1[gxyz=uvic@nrst]

let p = obs
let q = uvic
go regressxyz2

go prepareFigure.jnl Fig15
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! 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

let p = obs
let q = noresm
go regressxyz2
go taylor_polymark stdr r /line/fill/nokey/ov/color=black/pal=brown r plus 2.0
go finalizeFigure_psthicken