From af10cd990c9316bc3fdbd8aad87371d576221e44 Mon Sep 17 00:00:00 2001 From: daminton Date: Thu, 9 Mar 2017 13:14:48 +0000 Subject: [PATCH] Temporarily removing escape velocity condition from ejecta table in preparation for fixing the mass loss problem --- src/ejecta/ejecta_blanket.f90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/ejecta/ejecta_blanket.f90 b/src/ejecta/ejecta_blanket.f90 index 8362c0bb..8776fc1d 100644 --- a/src/ejecta/ejecta_blanket.f90 +++ b/src/ejecta/ejecta_blanket.f90 @@ -76,10 +76,10 @@ subroutine ejecta_blanket(user,crater,domain,erad,lrad,vejsq,ejang,firstrun) vej1sq = max(ce1sq * abs(user%gaccel * crater%grad) * (crater%grad / erad)**p,0._DP) gterm = ce1sq * (user%gaccel * erad) vejsq = max(vej1sq - gterm - yterm,0._DP) - if (vejsq >= domain%vescsq) then - lrad = VBIG - return - end if + !if (vejsq >= domain%vescsq) then + ! lrad = VBIG + ! return + !end if ! Ejection angle ejang=DEG2RAD * (55._DP - 20 * (erad / crater%grad))