Set to 3 scenarios.

Update PDFs stored on GitHub

Compare-multiple-ecospace-outputs.R

@@ -16,11 +16,9 @@ sim_scenario = "sim-spa_01"
 obs_TS_name  = "TS_updated_IB13"
 srt_year     = 1980
 spa_scenarios  = c("spa_ST00_base-no-drivers", "spa_ST01a_surf-sal", 
-                   "spa_ST01b_temp", "spa_ST01c_PP-MODIS", 
-                   "spa_ST01c_PP-MOM6")
+                   "spa_ST01b_temp", "spa_ST01c_PP-MODIS")
 spa_scen_names = c("01 No drivers",  "02 Salinity", 
-                   "03 Temperature", "04 PP (MODIS)",
-                   "05 PP (MOM6)")
+                   "03 Temperature", "04 PP (MODIS)")
 
 
 ## User-defined output parameters ----------------------------------------------
@@ -462,10 +460,12 @@ pdf_file_name_xM = paste0(dir_pdf_out, plot_name_xM)
   dev.off()    
 
 
-
   ## -----------------------------------------------------------------------------
   ##
   ## Plot by month (xM)
+
+  ## Make Ecosim line thinner and slightly darker
+  col_sim = rgb(0.2, 0.7, .1, alpha = 0.8); sim_lwd = 1
 
   print(paste("Writing", pdf_file_name_xM))
   pdf(pdf_file_name_xM, onefile = TRUE)

Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_10y/Ecospace-fits-summed-10y.csv

@@ -0,0 +1,13 @@
+"","scenario","nll_spa_obs","nll_spa_sim","nll_sim_obs","mae_spa_obs","mae_spa_sim","mae_sim_obs","weighting"
+"1","01 Base",1528.42,278.21,1780.01,27.65,18.64,28.62,"none"
+"2","01 Base",10574.291307125,5148.832830511,12001.342015046,193.523766425,120.171832494,206.322433078,"weighted"
+"3","01 Base",3566.341120738,1322.23572869733,4128.73772971725,65.0801581957053,41.2482146142378,68.9495539955147,"root-weighted"
+"4","02 MODIS-ChlA",1467.83,592.08,1780.01,26.22,20.19,28.62,"none"
+"5","02 MODIS-ChlA",10324.474429731,7523.78382826,12001.342015046,186.669077377,143.017539965,206.322433078,"weighted"
+"6","02 MODIS-ChlA",3453.20532351707,2074.3050814097,4128.73772971725,62.4015392489925,47.0652913469262,68.9495539955147,"root-weighted"
+"7","03 MOM6-ChlA",1820.66,857.77,1780.01,29.98,21.84,28.62,"none"
+"8","03 MOM6-ChlA",12721.993169623,7343.571399654,12001.342015046,211.766264941,137.540479711,206.322433078,"weighted"
+"9","03 MOM6-ChlA",4278.61525088078,2341.90143231264,4128.73772971725,71.1466079874312,48.013831544879,68.9495539955147,"root-weighted"
+"10","04 MOM6-Vint",1697.76,632.56,1780.01,28.41,19.8,28.62,"none"
+"11","04 MOM6-Vint",11887.474305725,6206.748921238,12001.342015046,204.794261454,128.596143316,206.322433078,"weighted"
+"12","04 MOM6-Vint",3990.63135504594,1891.98435207886,4128.73772971725,68.2266459919954,44.1697323961675,68.9495539955147,"root-weighted"

Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_1y/Ecospace-fits-summed-1y.csv

@@ -0,0 +1,13 @@
+,scenario,nll_spa_obs,nll_spa_sim,nll_sim_obs,mae_spa_obs,mae_spa_sim,mae_sim_obs,weighting
+1,01 Base,3005,1187,3307,325,23,327,none
+4,02 MODIS-ChlA,2977,1799,3307,321,24,327,none
+7,03 MOM6-ChlA,3234,1531,3307,326,26,327,none
+10,04 MOM6-Vint,3052,1418,3307,324,24,327,none
+3,01 Base,6521,2875,7765,488,51,499,root-weighted
+6,02 MODIS-ChlA,6174,4150,7765,482,56,499,root-weighted
+9,03 MOM6-ChlA,7100,3567,7765,496,58,499,root-weighted
+12,04 MOM6-Vint,6770,3233,7765,491,53,499,root-weighted
+2,01 Base,18401,8920,23284,842,148,882,weighted
+5,02 MODIS-ChlA,16781,12735,23284,827,172,882,weighted
+8,03 MOM6-ChlA,19961,10336,23284,868,167,882,weighted
+11,04 MOM6-Vint,19286,9453,23284,855,156,882,weighted

Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_36y/Ecospace-fits-summed-36y.csv

@@ -0,0 +1,13 @@
+"","scenario","nll_spa_obs","nll_spa_sim","nll_sim_obs","mae_spa_obs","mae_spa_sim","mae_sim_obs","weighting"
+"1","01 Base",187.96,-296.94,256.66,4.82,15.15,4.88,"none"
+"2","01 Base",2983.4140787,1222.809301761,3146.951432,56.9436994,97.597720245,55.3353167,"weighted"
+"3","01 Base",776.149780813616,26.3192835452572,858.254970125903,15.5687387187602,33.5381787203903,15.2554491405406,"root-weighted"
+"4","02 MODIS-ChlA",203.56,-218.14,256.66,5.24,16.19,4.88,"none"
+"5","02 MODIS-ChlA",3053.8223136,2606.311796629,3146.951432,58.2092391,113.916960098,55.3353167,"weighted"
+"6","02 MODIS-ChlA",803.317400753428,364.564758444052,858.254970125903,16.1688652889165,37.582209868035,15.2554491405406,"root-weighted"
+"7","03 MOM6-ChlA",248.37,89.76,256.66,5.33,17.47,4.88,"none"
+"8","03 MOM6-ChlA",3471.9962728,2866.378214278,3146.951432,62.0706409,111.631852659,55.3353167,"weighted"
+"9","03 MOM6-ChlA",920.883139702998,720.23981913814,858.254970125903,17.0170704076916,38.5677126940699,15.2554491405406,"root-weighted"
+"10","04 MOM6-Vint",256,-160.39,256.66,5.38,15.99,4.88,"none"
+"11","04 MOM6-Vint",3338.9583635,1635.833140696,3146.951432,61.204098,105.125771274,55.3353167,"weighted"
+"12","04 MOM6-Vint",902.183003268733,245.499826728828,858.254970125903,16.8998927955922,35.886166650139,15.2554491405406,"root-weighted"

Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_5y/Ecospace-fits-summed-5y.csv

@@ -0,0 +1,13 @@
+"","scenario","nll_spa_obs","nll_spa_sim","nll_sim_obs","mae_spa_obs","mae_spa_sim","mae_sim_obs","weighting"
+"1","01 Base",1590.83,468.4,1945.1,27.72,19.68,29.86,"none"
+"2","01 Base",11059.316539584,5827.803308788,14298.377916095,179.775891386,123.128974413,204.311417055,"weighted"
+"3","01 Base",3683.90667832266,1602.62588093292,4673.60056433216,62.0096078250895,42.6679296218221,69.1749259298803,"root-weighted"
+"4","02 MODIS-ChlA",1551.02,945.07,1945.1,26.2,20.78,29.86,"none"
+"5","02 MODIS-ChlA",10215.662559155,9145.225225705,14298.377916095,171.273263761,146.579960509,204.311417055,"weighted"
+"6","02 MODIS-ChlA",3479.34797613811,2672.66233981866,4673.60056433216,58.9875922366575,48.2549107108908,69.1749259298803,"root-weighted"
+"7","03 MOM6-ChlA",1949.72,1123.73,1945.1,30.79,23.42,29.86,"none"
+"8","03 MOM6-ChlA",13861.499524024,8478.410398694,14298.377916095,209.518973211,143.827881292,204.311417055,"weighted"
+"9","03 MOM6-ChlA",4612.86461258168,2802.49774079135,4673.60056433216,71.352843235564,50.7644929610327,69.1749259298803,"root-weighted"
+"10","04 MOM6-Vint",1817.04,922.5,1945.1,28.82,20.99,29.86,"none"
+"11","04 MOM6-Vint",12975.747030271,7319.883328647,14298.377916095,197.830166109,133.655368233,204.311417055,"weighted"
+"12","04 MOM6-Vint",4292.47454680922,2357.13078747501,4673.60056433216,66.9425092965426,46.2240178769225,69.1749259298803,"root-weighted"

Scenario_comps/Test_EnvDrivers/Ecospace-fits-summed-test-STEdrivers.csv

@@ -0,0 +1,16 @@
+"scenario","ssr_spa_obs","ssr_spa_sim","mae_spa_obs","mae_spa_sim","weighting"
+"00 Ecosim",955,NA,68,NA,"none"
+"01 No drivers",1068,827,72,251,"none"
+"02 Salinity",1060,824,71,250,"none"
+"03 Temperature",1006,728,66,240,"none"
+"04 PP (MODIS)",1007,739,65,225,"none"
+"00 Ecosim",2015,NA,144,NA,"root-weighted"
+"01 No drivers",2822,1603,190,308,"root-weighted"
+"02 Salinity",2797,1592,188,306,"root-weighted"
+"03 Temperature",2600,1260,170,276,"root-weighted"
+"04 PP (MODIS)",2546,1330,165,267,"root-weighted"
+"00 Ecosim",5972,NA,427,NA,"weighted"
+"01 No drivers",8990,4237,606,508,"weighted"
+"02 Salinity",8906,4196,598,504,"weighted"
+"03 Temperature",8256,3062,539,404,"weighted"
+"04 PP (MODIS)",7935,3325,514,418,"weighted"