Temporarily removing escape velocity condition from ejecta table in p…

…reparation for fixing the mass loss problem
MintonGroup · Mar 9, 2017 · af10cd9 · af10cd9
1 parent 0cd1e8c
commit af10cd9
Showing 1 changed file with 4 additions and 4 deletions.
diff --git 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))