diff --git a/examples/global-lunar-bombardment/ctem.in b/examples/global-lunar-bombardment/ctem.in index 986b5672..2c8f6c75 100755 --- 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 index b1863ebe..a51c843c 100644 --- 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 index dc233529..7214c3cc 100644 --- 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)