Note: datfile run through ss_RewriteDatfile.m on 29-Sep-2006 09:45:10 Header This simulation uses two basement layers, with the upper one corresponding to a 600 m-thick basement aquifer (K=1700 W/m-K) and the lower one having a MIXED conductivity: 2.0 W/m-K extends from the base of the aquifer for to a depth of 2000 m, and the remainder has a conductivity of K=2.9 W/m-K (the default) 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 timing in use -- maximum timestep = 1000 yr No Node removal allowed Use porosity function from Davis et al., 1999 JGR ************** sub : type of model, subsidence [sub] or sedimentation [sed] 6 : number of stress periods 6 1700 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.9 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 5 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 20000 1.e6 : 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 50 5 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 2000 1.e6 : 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 500 50 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 200 : 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 1000 500 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 20 1.e6 : 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 5000 1000 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 20 1.e6 : 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 5000 5000 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 20 1.e6 : 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.11348206522113 0 0 4 4.22696038463974 0 0 6 6.34043798825707 0 0 8 7.39252532738655 0 0 10 7.40577585003221 0 0 12 7.3996546362042 0 0 15 7.4012255269089 0 0 20 7.40493161589553 0 0 25 7.40775678431119 0 0 30 7.41069795712701 0 0 40 7.41658334178056 0 0 50 7.42246542161839 0 0 75 7.43717131316502 0 0 106 7.45540660723008 0 0 131 7.47011248958318 0 0 156 7.48481837193639 0 0 206 7.51423013664278 0 0 306 7.57305366605543 0 0 406 7.63187719546781 0 0 606 7.74952425429189 0 0 806 57.8083477836404 0 2 1006 157.808347783609 0 2 1256 282.8083477836 0 2 1506 407.808347783597 0 2 1756 532.808347783597 0 2 2006 657.808347783598 0 2 3006 1080.22214088705 0 0 4006 1425.04972709394 0 0 6006 2114.70489950773 0 0 8006 2804.36007192152 0 0 10006 3494.01524433531 0 0 12506 4356.08420985255 0 0 15006 5218.15317536979 0 0 20006 6942.29110640427 0 0 25006 8666.42903743876 0 0 30006 10390.5669684732 0 0 40006 13838.8428305422 0 0 50006 17287.1186926112 0 0 60006 20735.3945546802 0 0 80006 27631.9462788181 0 0 100006 34528.4980029562 0 0