From d7b3216073c758e9d1ccdfd65e785281cfab31d4 Mon Sep 17 00:00:00 2001 From: Richard Methot Date: Wed, 18 Sep 2024 16:49:53 -0700 Subject: [PATCH] revise calc of spawn_time_season (#620) --- SS_readdata_330.tpl | 5 +++-- SS_write_report.tpl | 2 +- 2 files changed, 4 insertions(+), 3 deletions(-) diff --git a/SS_readdata_330.tpl b/SS_readdata_330.tpl index bfd5905c..a29aca93 100644 --- a/SS_readdata_330.tpl +++ b/SS_readdata_330.tpl @@ -199,9 +199,10 @@ } temp += subseasdur_delta(spawn_seas); } - spawn_time_seas = (temp1 - azero_seas(spawn_seas)) / seasdur(spawn_seas); // remaining fraction of year converted to fraction of season +// spawn_time_seas = (temp1 - azero_seas(spawn_seas)) / seasdur(spawn_seas); // incorrect: remaining fraction of year converted to fraction of season + spawn_time_seas = (temp1 - azero_seas(spawn_seas)); // timing in units of fraction of year such that exp(-Z*spawn_time_seas) will be correct } - echoinput << "SPAWN month: " << spawn_month << "; seas: " << spawn_seas << "; subseas_for_ALK: " << spawn_subseas << "; timing_in_season: " << spawn_time_seas << endl; + echoinput << "SPAWN month: " << spawn_month << "; seas: " << spawn_seas << "; subseas_for_ALK: " << spawn_subseas << "; spawntiming as frac. of year: " << spawn_time_seas << endl; if (spawn_seas > nseas) { warnstream << " spawn_seas index must be <= nseas "; diff --git a/SS_write_report.tpl b/SS_write_report.tpl index 80050111..ed3af3bc 100644 --- a/SS_write_report.tpl +++ b/SS_write_report.tpl @@ -327,7 +327,7 @@ FUNCTION void write_bigoutput() SS2out << "Season_Durations: " << seasdur << endl; SS2out << "Spawn_month: " << spawn_month << endl << "Spawn_seas: " << spawn_seas << endl - << "Spawn_timing_in_season: " << spawn_time_seas << endl; + << "spawntiming as frac. of year: " << spawn_time_seas << endl; SS2out << "N_areas: " << pop << endl; SS2out << "Start_year: " << styr << endl; SS2out << "End_year: " << endyr << endl;