quasi-MC mode works in test run

MintonGroup · Apr 2, 2021 · 354734b · 354734b
1 parent 32d9262
commit 354734b
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 5 deletions.
diff --git a/examples/global-lunar-bombardment/ctem.in b/examples/global-lunar-bombardment/ctem.in
@@ -36,7 +36,7 @@ runtype      single                             ! Run type: options are normal /
 seed         33790                    ! Random number generator seed
 gridsize     2000                          ! Size of grid in pixels
 numlayers    10                                 ! Number of perched layers
-pix          2.0e2                           ! Pixel size (m)
+pix          3.08e3                           ! Pixel size (m)
 mat          rock                               ! Material (rock or ice)
 ! Bedrock scaling parameters
 mu_b         0.55e0                             ! Experimentally derived parameter for bedrock crater scaling law 
@@ -60,7 +60,7 @@ doseismic    F                                  ! Perform seismic shaking calcul
 
 ! Optional inputF These have internally set default values that work reasonable well. Comment them out with
 deplimit        9e99                            ! Depth limit for craters (m) - Default is to ignore.
-maxcrat         1.00e0                          ! Fraction of gridsize that maximum crater can be - Default 1.0
+maxcrat         0.20e0                          ! Fraction of gridsize that maximum crater can be - Default 1.0
 killatmaxcrater F                               ! Stop the run if a crater larger than the maximum is produced - Default F      
 basinimp        35.0e3                          ! Size of impactor to switch to lunar basin scaling law - Default is to ignore
 docollapse      T                               ! Do slope collapse - Default T

diff --git a/src/crater/crater_populate.f90 b/src/crater/crater_populate.f90
@@ -144,7 +144,7 @@ subroutine crater_populate(user,surf,crater,domain,prod,production_list,vdist,nt
       !if in quasiMC mode: check to see if it's time for a real crater
       if (user%doquasimc) then
          if (crater%timestamp > user%rctime) then
-            write(*,*) "real crater at time ", rclist(6, rccount)
+            write(*,*) "Real crater at time ", rclist(6, rccount)
             user%testflag = .true.
             user%testimp = rclist(1, rccount)
             user%testvel = rclist(2, rccount)
@@ -162,7 +162,7 @@ subroutine crater_populate(user,surf,crater,domain,prod,production_list,vdist,nt
             user%testflag = .false.
             rccount = rccount + 1
             if (rccount > domain%rcnum) then
-               write(*,*) "real crater list complete."
+               write(*,*) "Real crater list complete."
                user%rctime = 1e30
             else
                user%rctime = rclist(6,rccount)

diff --git a/src/main/CTEM.f90 b/src/main/CTEM.f90
@@ -89,7 +89,7 @@ program CTEM
 
 ! Read in real crater list for quasi-MC run
 call io_read_craterlist(rclist,user,domain)
-write(*,*) rclist
+!write(*,*) rclist
 
 write(*,*) "Initializing simulation domain and determining minimum impactor size"
 call init_domain(user,crater,domain,prod,pdist,vdist,crtscl,nflux)