fixed the depth model bug once and for all

MintonGroup · Feb 10, 2023 · e795580 · e795580
1 parent 9253b08
commit e795580
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/src/regolith/regolith_depth_model.f90 b/src/regolith/regolith_depth_model.f90
@@ -97,8 +97,9 @@ subroutine regolith_depth_model(user,domain,finterval,nflux,p)
          f    = 0.5_DP * dr * ( fmin + fmax )
          psum = psum + f   
       end do
-      psumfunc = exp(-1.0_DP * PI * psum * t)
-      if (psumfunc > epsilon(psum)) then
+      !psumfunc = exp(-1.0_DP * PI * psum * t)
+      psumfunc = -1.0_DP * PI * psum * t
+      if (psumfunc > log(epsilon(psum))) then
          p(2,i) = 1.0_DP - exp(-1.0_DP * PI * psum * t)
       else
          p(2,i) = 1.0_DP