Note: datfile run through ss_RewriteDatfile.m on 27-Sep-2006 12:12:50 Header Test how sensitive the results are to the choice of conductivity for lowermost basement. This simulation uses a lower basement value of 2 W/m-K Basal heat flux = 1 W/m^2 100 km thick lower basement Subsidence at 100 m/Ma for 20 Ma to accumulate 2000 m of sediment Sediment-basement-interface at depth of 6 m Variable timing in use : maximum timestep of 4kyr is consistent with node addition criteria of 2 m No Node removal allowed Use porosity function from Davis et al., 1999 JGR ************** sub : type of model, subsidence [sub] or sedimentation [sed] 5 : number of stress periods 6 2.0 3.86e6 0 : layer1 parameters: depth of node at SBI (m), basement conductivity (W/m-K), thermal capacity (J/m^3-K), porosity (decimal) 606 2.0 3.86e6 0 : layer2 parameters: depth of node at layer1/layer2 interface, basement conductivity (W/m-K), thermal capacity (J/m^3-K), porosity (decimal) 4.30e6 : parameter, thermal capacity of water (J/m^3-K) 2.65e6 : parameter, thermal capacity of sediment (J/m^3-K) 0.6 : parameter, thermal conductivity of water (W/m-K) 2.74 : parameter, thermal conductivity of sediment grains (W/m-K) 0.7 : parameter, surface sediment porosity (decimal) 0 0 0 0 0 0.7 -8.333e-4 0 0 : parameter, constants for porosity = f(z). A+Bz+Cz^2+Dz^3+Elnz+Fexp(G*Z)+H^(Iz) m : parameter, porosity = f(z) where z is in [m] or [km] 0.0 : parameter, minimum allowable sediment porosity (decimal) 0 0 : parameter, constants for permeability = f(phi) when pressure term is used to drive seepage: perm = Aexp(B*(porosity/(1-porosity))) 0.5 : parameter, scaling factor theta (for crank-nicholson solution: 0-1) 0=explicit, 1=implicit, 0.5=mixed 2 : parameter, calculate heat flow between surface and this node no : Flag, allow the removal of nodes from upper basement [yes/no] followed by maximum number of nodes to remove. Will remove a basement node when a sediment node is added yes : Flag, write ouput to a text Log file [yes/no] (Will write input-file data regardless) ************** 1.e6 4000 0.1 2000 1.1 : time, length of this stress period (in yrs) followed by maximum time step (in yrs), Variable params: initial timestep (yrs), # timesteps, scaling factor (initial guess) q : boundary condition (lower), ([T] for Temperature (degrees C), [q] for heat flow (W/m^2)) 1 0 0 0 : boundary condition paramaters for lower boundary: [T] or [q] = f(time in yrs): A+B(time)+C/sqrt(D*time) 0 : boundary condition (upper) , temperature in degrees C that boundary is held constant at. 0.000100 : parameter, basement subsidence rate or sedimentation rate during this stress period in m/yr [] [] : parameter, [s] seepage followed by value in m/yr or [p] lower boundary pressure followed by value in (kPa) : (+) down/underpressure 0.0 : parameter, production/sink term for additional nodes during this stress period in (W/m^2) [] : filename containing constants for calculating production/sink (Q=f(t)) values for selected nodes 2 0.5 0 : parameter, subsidence distance (m) for addition of a new node followed by tolerance (0-1). If param 1 ==-1, variable node addition depth scaled by param. 3 3 100 : parameter, number of time steps to increment before storing data in Mat File followed by writing to Log File (Flag to write must be set to [yes] for writing to Log file) 2.e6 4000 4000 100 1.1 : time, length of this stress period (in yrs) followed by maximum time step (in yrs), Variable params: initial timestep (yrs), # timesteps, scaling factor (initial guess) q : boundary condition (lower), ([T] for Temperature (degrees C), [q] for heat flow (W/m^2)) 1 0 0 0 : boundary condition paramaters for lower boundary: [T] or [q] = f(time in yrs): A+B(time)+C/sqrt(D*time) 0 : boundary condition (upper) , temperature in degrees C that boundary is held constant at. 0.000100 : parameter, basement subsidence rate or sedimentation rate during this stress period in m/yr [] [] : parameter, [s] seepage followed by value in m/yr or [p] lower boundary pressure followed by value in (kPa) : (+) down/underpressure 0.0 : parameter, production/sink term for additional nodes during this stress period in (W/m^2) [] : filename containing constants for calculating production/sink (Q=f(t)) values for selected nodes 2 0.5 0 : parameter, subsidence distance (m) for addition of a new node followed by tolerance (0-1). If param 1 ==-1, variable node addition depth scaled by param. 3 3 100 : parameter, number of time steps to increment before storing data in Mat File followed by writing to Log File (Flag to write must be set to [yes] for writing to Log file) 2.e6 4000 4000 100 1.1 : time, length of this stress period (in yrs) followed by maximum time step (in yrs), Variable params: initial timestep (yrs), # timesteps, scaling factor (initial guess) q : boundary condition (lower), ([T] for Temperature (degrees C), [q] for heat flow (W/m^2)) 1 0 0 0 : boundary condition paramaters for lower boundary: [T] or [q] = f(time in yrs): A+B(time)+C/sqrt(D*time) 0 : boundary condition (upper) , temperature in degrees C that boundary is held constant at. 0.000100 : parameter, basement subsidence rate or sedimentation rate during this stress period in m/yr [] [] : parameter, [s] seepage followed by value in m/yr or [p] lower boundary pressure followed by value in (kPa) : (+) down/underpressure 0.0 : parameter, production/sink term for additional nodes during this stress period in (W/m^2) [] : filename containing constants for calculating production/sink (Q=f(t)) values for selected nodes 2 0.5 0 : parameter, subsidence distance (m) for addition of a new node followed by tolerance (0-1). If param 1 ==-1, variable node addition depth scaled by param. 3 3 100 : parameter, number of time steps to increment before storing data in Mat File followed by writing to Log File (Flag to write must be set to [yes] for writing to Log file) 5.e6 4000 4000 100 1.1 : time, length of this stress period (in yrs) followed by maximum time step (in yrs), Variable params: initial timestep (yrs), # timesteps, scaling factor (initial guess) q : boundary condition (lower), ([T] for Temperature (degrees C), [q] for heat flow (W/m^2)) 1 0 0 0 : boundary condition paramaters for lower boundary: [T] or [q] = f(time in yrs): A+B(time)+C/sqrt(D*time) 0 : boundary condition (upper) , temperature in degrees C that boundary is held constant at. 0.000100 : parameter, basement subsidence rate or sedimentation rate during this stress period in m/yr [] [] : parameter, [s] seepage followed by value in m/yr or [p] lower boundary pressure followed by value in (kPa) : (+) down/underpressure 0.0 : parameter, production/sink term for additional nodes during this stress period in (W/m^2) [] : filename containing constants for calculating production/sink (Q=f(t)) values for selected nodes 3 0.5 0 : parameter, subsidence distance (m) for addition of a new node followed by tolerance (0-1). If param 1 ==-1, variable node addition depth scaled by param. 3 10 500 : parameter, number of time steps to increment before storing data in Mat File followed by writing to Log File (Flag to write must be set to [yes] for writing to Log file) 10.e6 8000 4000 100 1.1 : time, length of this stress period (in yrs) followed by maximum time step (in yrs), Variable params: initial timestep (yrs), # timesteps, scaling factor (initial guess) q : boundary condition (lower), ([T] for Temperature (degrees C), [q] for heat flow (W/m^2)) 1 0 0 0 : boundary condition paramaters for lower boundary: [T] or [q] = f(time in yrs): A+B(time)+C/sqrt(D*time) 0 : boundary condition (upper) , temperature in degrees C that boundary is held constant at. 0.000100 : parameter, basement subsidence rate or sedimentation rate during this stress period in m/yr [] [] : parameter, [s] seepage followed by value in m/yr or [p] lower boundary pressure followed by value in (kPa) : (+) down/underpressure 0.0 : parameter, production/sink term for additional nodes during this stress period in (W/m^2) [] : filename containing constants for calculating production/sink (Q=f(t)) values for selected nodes 4 0.5 0 : parameter, subsidence distance (m) for addition of a new node followed by tolerance (0-1). If param 1 ==-1, variable node addition depth scaled by param. 3 10 500 : parameter, number of time steps to increment before storing data in Mat File followed by writing to Log File (Flag to write must be set to [yes] for writing to Log file) ************** 42 : node, number of nodes. Following are node depths (m), initial temps (deg C), radiogenic production/sink (W/m^2), [optional] conductivity (W/m-K) 0 0 0 0 2 2.1116147708694 0 0 4 4.21949761883175 0 0 6 6.32366133741977 0 0 8 7.87481501013001 0 0 10 8.87481501012977 0 0 12 9.8748150101296 0 0 15 11.3748150101295 0 0 20 13.8748150101292 0 0 25 16.3748150101289 0 0 30 18.8748150101287 0 0 40 23.8748150101286 0 0 50 28.8748150101288 0 0 75 41.3748150101286 0 0 106 56.8748150101284 0 0 131 69.3748150101278 0 0 156 81.8748150101298 0 0 206 106.874815010131 0 0 306 156.874815010133 0 0 406 206.874815010132 0 0 606 306.874815010131 0 0 806 406.874815010131 0 0 1006 506.874815010132 0 0 1256 631.874815010132 0 0 1506 756.87481501013 0 0 1756 881.874815010132 0 0 2006 1006.87481501013 0 0 3006 1506.87481501013 0 0 4006 2006.87481501013 0 0 6006 3006.87481501013 0 0 8006 4006.87481501014 0 0 10006 5006.87481501014 0 0 12506 6256.87481501014 0 0 15006 7506.87481501014 0 0 20006 10006.8748150101 0 0 25006 12506.8748150102 0 0 30006 15006.8748150102 0 0 40006 20006.8748150102 0 0 50006 25006.8748150102 0 0 60006 30006.8748150102 0 0 80006 40006.8748150101 0 0 100006 50006.8748150101 0 0