diff --git a/Compare-multiple-ecospace-outputs.R b/Compare-multiple-ecospace-outputs.R index e61928b..17f5d1e 100644 --- a/Compare-multiple-ecospace-outputs.R +++ b/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) diff --git a/Scenario_comps/Compare_PP_Drivers/PDF_plots/BxY_scaled_10y-multispa.PDF b/Scenario_comps/Compare_PP_Drivers/PDF_plots/BxY_scaled_10y-multispa.PDF new file mode 100644 index 0000000..3dc4e3c Binary files /dev/null and b/Scenario_comps/Compare_PP_Drivers/PDF_plots/BxY_scaled_10y-multispa.PDF differ diff --git a/Scenario_comps/Compare_PP_Drivers/PDF_plots/BxY_scaled_1y-multispa.PDF b/Scenario_comps/Compare_PP_Drivers/PDF_plots/BxY_scaled_1y-multispa.PDF new file mode 100644 index 0000000..a074c68 Binary files /dev/null and b/Scenario_comps/Compare_PP_Drivers/PDF_plots/BxY_scaled_1y-multispa.PDF differ diff --git a/Scenario_comps/Compare_PP_Drivers/PDF_plots/BxY_scaled_36y-multispa.PDF b/Scenario_comps/Compare_PP_Drivers/PDF_plots/BxY_scaled_36y-multispa.PDF new file mode 100644 index 0000000..332bc6c Binary files /dev/null and b/Scenario_comps/Compare_PP_Drivers/PDF_plots/BxY_scaled_36y-multispa.PDF differ diff --git a/Scenario_comps/Compare_PP_Drivers/PDF_plots/BxY_scaled_5y-multispa.PDF b/Scenario_comps/Compare_PP_Drivers/PDF_plots/BxY_scaled_5y-multispa.PDF new file mode 100644 index 0000000..126b209 Binary files /dev/null and b/Scenario_comps/Compare_PP_Drivers/PDF_plots/BxY_scaled_5y-multispa.PDF differ diff --git a/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_10y/Ecospace-fits-summed-10y.csv b/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_10y/Ecospace-fits-summed-10y.csv new file mode 100644 index 0000000..31416c5 --- /dev/null +++ b/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" diff --git a/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_10y/Fit_Metrics_10y.xlsx b/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_10y/Fit_Metrics_10y.xlsx new file mode 100644 index 0000000..af39893 Binary files /dev/null and b/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_10y/Fit_Metrics_10y.xlsx differ diff --git a/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_1y/Ecospace-fits-summed-1y.csv b/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_1y/Ecospace-fits-summed-1y.csv new file mode 100644 index 0000000..c7a752a --- /dev/null +++ b/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 diff --git a/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_1y/Ecospace-fits-summed-1y.xlsx b/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_1y/Ecospace-fits-summed-1y.xlsx new file mode 100644 index 0000000..95c58d9 Binary files /dev/null and b/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_1y/Ecospace-fits-summed-1y.xlsx differ diff --git a/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_1y/Fit_Metrics_1y.xlsx b/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_1y/Fit_Metrics_1y.xlsx new file mode 100644 index 0000000..bd366ab Binary files /dev/null and b/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_1y/Fit_Metrics_1y.xlsx differ diff --git a/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_36y/Ecospace-fits-summed-36y.csv b/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_36y/Ecospace-fits-summed-36y.csv new file mode 100644 index 0000000..2befdfc --- /dev/null +++ b/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" diff --git a/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_36y/Fit_Metrics_36y.xlsx b/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_36y/Fit_Metrics_36y.xlsx new file mode 100644 index 0000000..e4b2cb0 Binary files /dev/null and b/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_36y/Fit_Metrics_36y.xlsx differ diff --git a/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_5y/Ecospace-fits-summed-5y.csv b/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_5y/Ecospace-fits-summed-5y.csv new file mode 100644 index 0000000..9b643b7 --- /dev/null +++ b/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" diff --git a/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_5y/Fit_Metrics_5y.xlsx b/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_5y/Fit_Metrics_5y.xlsx new file mode 100644 index 0000000..e75c09b Binary files /dev/null and b/Scenario_comps/Compare_PP_Drivers/Tables/Fit_metrics_scaled_5y/Fit_Metrics_5y.xlsx differ diff --git a/Scenario_comps/Habitat foraging usage.xlsx b/Scenario_comps/Habitat foraging usage.xlsx new file mode 100644 index 0000000..82c9901 Binary files /dev/null and b/Scenario_comps/Habitat foraging usage.xlsx differ diff --git a/Scenario_comps/Test_EnvDrivers/BxM_scaled_6y-test-STEdrivers.PDF b/Scenario_comps/Test_EnvDrivers/BxM_scaled_6y-test-STEdrivers.PDF new file mode 100644 index 0000000..dea9256 Binary files /dev/null and b/Scenario_comps/Test_EnvDrivers/BxM_scaled_6y-test-STEdrivers.PDF differ diff --git a/Scenario_comps/Test_EnvDrivers/BxY_scaled_6y-test-STEdrivers.PDF b/Scenario_comps/Test_EnvDrivers/BxY_scaled_6y-test-STEdrivers.PDF new file mode 100644 index 0000000..ac18c89 Binary files /dev/null and b/Scenario_comps/Test_EnvDrivers/BxY_scaled_6y-test-STEdrivers.PDF differ diff --git a/Scenario_comps/Test_EnvDrivers/Ecospace-fits-summed-test-STEdrivers.csv b/Scenario_comps/Test_EnvDrivers/Ecospace-fits-summed-test-STEdrivers.csv new file mode 100644 index 0000000..f54db84 --- /dev/null +++ b/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" diff --git a/Scenario_comps/Test_EnvDrivers/Fit_metrics_test-STEdrivers.xlsx b/Scenario_comps/Test_EnvDrivers/Fit_metrics_test-STEdrivers.xlsx new file mode 100644 index 0000000..5eb451a Binary files /dev/null and b/Scenario_comps/Test_EnvDrivers/Fit_metrics_test-STEdrivers.xlsx differ