diff --git a/examples/global-lunar-bombardment/ctem.in b/examples/global-lunar-bombardment/ctem.in index e3ecff98..e3eb94fa 100755 --- a/examples/global-lunar-bombardment/ctem.in +++ b/examples/global-lunar-bombardment/ctem.in @@ -1,8 +1,23 @@ -! CTEM Input file: Compatible with v. 1.1 +! CTEM Input file -! IDL driver inputs -interval 1.0 -numintervals 300 ! Total number of intervals + +! Testing input. These are used to perform non-Monte Carlo tests. +testflag F ! Set to T to create a single crater with user-defined impactor properties +testimp 15000.0 ! Diameter of test impactor (m) +testvel 18.3e3 ! Velocity of test crater (m/s) +testang 45.0 ! Impact angle of test crater (deg) - Default 90.0 +testxoffset 0.0 ! x-axis offset of crater center from grid center (m) - Default 0.0 +testyoffset 0.0 ! y-axis offset of crater center from grid center (m) - Default 0.0 +tallyonly F ! Tally the craters without generating any craters +testtally F +quasimc T ! MC run constrained by non-MC 'real' craters given in a list +realcraterlist craterlist.in ! list of 'real' craters for Quasi-MC runs + + + +! IDL driver in uts +interval 6.280 +numintervals 100 ! Total number of intervals (total time = interval * numintervals) <--when runtype is 'single' restart F ! Restart a previous run impfile NPFextrap.dat ! Impactor SFD rate file (col 1: Dimp (m), col 2: ! impactors > D (m**(-2) y**(-1)) popupconsole F ! Pop up console window every output interval @@ -18,10 +33,10 @@ runtype single ! Run type: options are normal / ! statistical: surface is reset between intervals ! CTEM required inputs -seed 11082 ! Random number generator seed -gridsize 200 ! Size of grid in pixels +seed 76535 ! Random number generator seed +gridsize 2000 ! Size of grid in pixels numlayers 10 ! Number of perched layers -pix 3.08e4 ! 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 @@ -44,18 +59,17 @@ velfile lunar-MBA-impactor-velocities.dat ! Impactor velocity dist file doseismic F ! Perform seismic shaking calculations with each impact - Default F ! Optional inputF These have internally set default values that work reasonable well. Comment them out with -deplimit 9d99 ! Depth limit for craters (m) - Default is to ignore. -maxcrat 0.20e0 ! Fraction of gridsize that maximum crater can be - Default 1.0 +deplimit 9e99 ! Depth limit for craters (m) - Default is to ignore. +maxcrat 3.15e-2 ! 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 dosoftening T ! Do ejecta softening - Default T doangle T ! Vary the impact angle. Set to F to have only vertical impacts - Default T -Kd1 0.0315151258095091 +Kd1 0.0001 psi 2.000 -psi2 1.250 -rbreak 100.0 -fe 5.0 -ejecta_truncation 5.0 -dorays F +fe 4.00 +ejecta_truncation 4.0 +dorays T superdomain F + diff --git a/src/Makefile.am b/src/Makefile.am index 1eaeb9e2..196c87f8 100755 --- a/src/Makefile.am +++ b/src/Makefile.am @@ -1,9 +1,9 @@ bin_PROGRAMS = CTEM #ifort optimized flags -AM_FCFLAGS = -O3 -qopenmp -parallel -xHost -assume byterecl -m64 -heap-arrays -FR +#AM_FCFLAGS = -O3 -qopenmp -parallel -xHost -assume byterecl -m64 -heap-arrays -FR #AM_FCFLAGS = -O3 -qopenmp -parallel -assume byterecl -m64 -heap-arrays -FR #ifort debug flags -#AM_FCFLAGS = -O3 -p -g -openmp -debug all -traceback -CB -assume byterecl -m64 -heap-arrays -FR +AM_FCFLAGS = -O0 -p -g -qopenmp -debug all -traceback -CB -assume byterecl -m64 -heap-arrays -FR #gfortran optimized flags #AM_FCFLAGS = -O3 -fopenmp -ffree-form -g -fbounds-check -fbacktrace