! Normalised taylor diagrams of global average diatom biomass
set mem/size=100
! Load observational dataset
use "$HOME/Data/PANGAEA/MAREDAT/diatoms/MarEDat20120716Diatoms.nc"
! ! ! ! !  CONVERSION FACTORS  ! ! ! ! ! ! ! 
! convert from per second to per day
let mult = 60*60*24

! convert from per day to per year
let mult1 = mult*365

! convert nitrogen to carbon
let ntc = 6.625

! convert mol C to grams C
let mtg = 12.011
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

! Load new UVic model
use "tavg.01775.01.01.nc"
let uvic = O_diat[d=2,l=25]*ntc !year 2014
let obs = biomass[d=1,gxyz=uvic@nrst,l=@ave]/mtg/1000

canc mode cal

! Load CanESM5-CanOE
use "$HOME/Data/cmip6_data/diat/phydiat_Oyr_CanESM5-CanOE_historical_r1i1p2f1_gn_1850-2014.nc"
let map1 = curv_to_rect_map(longitude[d=3],latitude[d=3],uvic,1)
let can1 = curv_to_rect(phydiat[d=3,l=165],map1)
let canesm = can1[gxyz=uvic@nrst]

! Load EC-Earth3-CC
use "$HOME/Data/cmip6_data/diat/phydiat_Oyr_EC-Earth3-CC_historical_r1i1p1f1_gn_2014-2014.nc"
let map2 = curv_to_rect_map(longitude[d=4],latitude[d=4],uvic,1)
let ec1 = curv_to_rect(phydiat[d=4,l=1],map2)
let ec = ec1[gxyz=uvic@nrst]

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

! Load IPSL
use "$HOME/Data/cmip6_data/diat/phydiat_Oyr_IPSL-CM6A-LR_historical_r32i1p1f1_gn_1850-2014.nc"
let map3 = curv_to_rect_map(nav_lon[d=6],nav_lat[d=6],uvic,1)
let ipsl1 = curv_to_rect(phydiat[d=6,l=165],map3)
let ipsl = ipsl1[gxyz=uvic@nrst]

! Load CMCC-ESM2
use "$HOME/Data/cmip6_data/diat/phydiat_Oyr_CMCC-ESM2_historical_r1i1p1f1_gn_1850-2014.nc"
let map4 = curv_to_rect_map(longitude[d=7],latitude[d=7],uvic,1)
let cmcc1 = curv_to_rect(phydiat[d=7,l=165],map4)
let cmcc = cmcc1[gxyz=uvic@nrst]

! Load CESM
use "$HOME/Data/cmip6_data/diat/phydiat_Oyr_CESM2_historical_r11i1p1f1_gr_2000-2014.nc"
let cesm = phydiat[l=15,d=8,gxyz=uvic@nrst]

!let p = uvic
!let q = uvic
!go regressxyz2

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

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

let p = uvic
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 = uvic
let q = ec
go regressxyz2
go taylor_polymark stdr r /line/fill/nokey/ov/color=black/pal=blue r diamond 2.0

let p = uvic
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 = uvic
let q = canesm
go regressxyz2
go taylor_polymark stdr r /line/fill/nokey/ov/color=black/pal=yellow_dark 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 = uvic
!let q = cmcc
!go regressxyz2
!go taylor_polymark stdr r /line/fill/nokey/ov/color=black/pal=pink r pentagon 2.0

go finalizeFigure_psthicken