diff --git a/examples/global-lunar-bombardment/ctem.in b/examples/global-lunar-bombardment/ctem.in
new file mode 100755
index 00000000..9467c18b
--- /dev/null
+++ b/examples/global-lunar-bombardment/ctem.in
@@ -0,0 +1,60 @@
+! CTEM Input file: Compatible with v. 1.1
+
+! IDL driver inputs
+interval     5.0  
+numintervals 100                                ! Total number of intervals
+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 
+saveshaded   F                                  ! Output shaded relief images  
+saverego     F                                  ! Output regolith map images 
+savepres     F                                  ! Output simplified console display images (presentation-compatible images)
+savetruelist T                                  ! Save the true cumulative crater distribution for each interval (large file size)
+sfdcompare   LOLASethCraterCatalogv8gt20-binned.dat         ! File name for the SFD comparison file used in the console display
+shadedminh   -85.0                              ! Minimum height for shaded relief map (m) (Default - automatic)
+shadedmaxh   85.0                               ! Maximum height for shaded relief map (m) (Default - automatic)
+runtype      single                             ! Run type: options are normal / statistical
+                                                !    single:      craters accumulate in successive intervals
+                                                !    statistical: surface is reset between intervals
+
+! CTEM required inputs
+seed         11082                    ! Random number generator seed
+gridsize     2000                          ! Size of grid in pixels
+numlayers    10                                 ! Number of perched layers
+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 
+kv_b         0.20e0                             ! Experimentally derived parameter for bedrock crater scaling law 
+trho_b       2250.0e0                           ! Target density (bedrock) (kg/m**3)
+ybar_b       0.0e6                               ! Bedrock strength (Pa)
+! Regolith scaling parameters
+mu_r         0.55e0                             ! Experimentally derived parameter for regolith crater scaling law 
+kv_r         0.20e0                             ! Experimentally derived parameter for regolith crater scaling law
+trho_r       2250.0e0                           ! Target density (regolith) (kg/m**3)
+ybar_r       0.00e6                              ! Regolith strength (Pa)
+! Body parameters
+gaccel       1.62e0                             ! Gravitational acceleration at target (m/s**2)
+trad         1737.35e3                          ! Target radius (m)
+prho         2500.0e0                           ! Projectile density (kg/m**3)
+sfdfile      production.dat                     ! Impactor SFD file (col 1: Dimp (m), col 2: ! impactors > D
+velfile      lunar-MBA-impactor-velocities.dat  ! Impactor velocity dist file
+
+! Seismic shaking input (required if seismic shaking is set to T, otherwise ignored)
+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.15e0                          ! 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.0001
+psi             2.000
+fe              10.0
+ejecta_truncation 20.0
+dorays  T 
+superdomain F
+