hacky fix so ages work properly in multi-interval runs, but multi-int…

…erval runs are slower because of this
MintonGroup · Aug 31, 2023 · b6c7c5f · b6c7c5f
1 parent 06dfe51
commit b6c7c5f
Showing 1 changed file with 10 additions and 2 deletions.
diff --git a/src/crater/crater_populate.f90 b/src/crater/crater_populate.f90
@@ -179,9 +179,17 @@ subroutine crater_populate(user,surf,crater,domain,prod,production_list,vdist,nt
       end if
       if (crater%timestamp < 2330._DP) then
          if (oldGa > 0._DP) then 
-            crater%timestampGa = util_t_from_scale(maxage-crater%timestamp,agemin,oldGa)
+            if (user%numintervals .eq. 1) then
+               crater%timestampGa = util_t_from_scale(maxage-crater%timestamp,agemin,oldGa)
+            else
+               crater%timestampGa = util_t_from_scale(maxage-crater%timestamp,1e-10_DP,oldGa)
+            end if
          else
-            crater%timestampGa = util_t_from_scale(maxage-crater%timestamp,agemin,maxageGa)
+            if (user%numintervals .eq. 1) then
+               crater%timestampGa = util_t_from_scale(maxage-crater%timestamp,agemin,maxageGa)
+            else
+               crater%timestampGa = util_t_from_scale(maxage-crater%timestamp,1e-10_DP,maxageGa)
+            end if
          end if
       else
          crater%timestampGa = 4.5_DP