fixed ejecta_emplace exception bug for quasimc branch

MintonGroup · Sep 15, 2023 · 00c7734 · 00c7734
1 parent 3a6a9fa
commit 00c7734
Showing 1 changed file with 6 additions and 5 deletions.
diff --git a/src/ejecta/ejecta_emplace.f90 b/src/ejecta/ejecta_emplace.f90
@@ -347,7 +347,6 @@ subroutine ejecta_emplace(user,surf,crater,domain,ejb,ejtble,deltaMtot,cumulativ
                ebh = cumulative_elchange(i,j) - crater_profile(user, crater, lrad)
             end if
 
-
             if (user%doregotrack .and. ebh>1.0e-8_DP) then
                call regolith_streamtube(user,surf,crater,domain,ejb,ejtble,xp,yp,xpi,ypi,lrad,ebh,rm,vsq,volm)
                vmelt = vmelt + volm
@@ -356,10 +355,12 @@ subroutine ejecta_emplace(user,surf,crater,domain,ejb,ejtble,deltaMtot,cumulativ
       end do
    end if
 
-   if (totmelt > vmelt) then
-      vmeltsheet = totmelt - vmelt
-   else !give the craters a melt sheet of 1m
-      vmeltsheet = 1.0_DP * user%pix * user%pix * nmeltsheet
+   if (user%doregotrack) then
+      if (totmelt > vmelt) then
+         vmeltsheet = totmelt - vmelt
+      else !give the craters a melt sheet of 1m
+         vmeltsheet = 1.0_DP * user%pix * user%pix * nmeltsheet
+      end if
    end if
 
    ! extra soften calculation