From b6c7c5f9fb4f1ec714a083dd5e755e710af72560 Mon Sep 17 00:00:00 2001 From: Austin Blevins Date: Thu, 31 Aug 2023 18:52:40 -0400 Subject: [PATCH] hacky fix so ages work properly in multi-interval runs, but multi-interval runs are slower because of this --- src/crater/crater_populate.f90 | 12 ++++++++++-- 1 file changed, 10 insertions(+), 2 deletions(-) diff --git a/src/crater/crater_populate.f90 b/src/crater/crater_populate.f90 index 7cf773e1..d0f6734a 100644 --- 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