revise calc of spawn_time_season

nmfs-ost · Sep 18, 2024 · da4b8ca · da4b8ca
1 parent 605a4c2
commit da4b8ca
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 3 deletions.
diff --git 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
@@ -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;