Note: datfile run through ss_RewriteDatfile.m on 29-Sep-2006 14:45:05 Header This simulation uses two basement layers, with the upper one corresponding to a 600 m-thick basement aquifer (K=1700 W/m-K : Nu=1000) 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) Heat sinks are evenly distributed in the aquifer 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 Heat sinks are active for the first 5 Ma (extracting 80%), then rampdown linearly over the following 5 Ma 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] 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 2.e4 : 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) 1.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) 3.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 1000 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) Sub100_20MaQcHSRD5Ma-80.sinks-off_600.hflx: 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 0.422322954173274 0 0 4 0.843899523765146 0 0 6 1.26473226749584 0 0 8 1.47508061317902 -0.004 0 10 1.47532072089854 -0.004 0 12 1.47556531895887 -0.004 0 15 1.47593947749933 -0.006 0 20 1.47658064879038 -0.01 0 25 1.4772512373358 -0.01 0 30 1.47795123731761 -0.01 0 40 1.47941006084764 -0.02 0 50 1.48098653143582 -0.02 0 75 1.48522182555337 -0.05 0 106 1.49138535496506 -0.062 0 131 1.49726770790619 -0.05 0 156 1.50388535496499 -0.05 0 206 1.51859123731792 -0.1 0 306 1.55388535496495 -0.2 0 406 1.60094417849432 -0.2 0 606 1.71859123731777 0 0 806 51.7774147667301 0 2 1006 151.777414766731 0 2 1256 276.777414766731 0 2 1506 401.777414766732 0 2 1756 526.777414766733 0 2 2006 651.777414766734 0 2 3006 1074.19120787018 0 0 4006 1419.01879407707 0 0 6006 2108.67396649087 0 0 8006 2798.32913890469 0 0 10006 3487.98431131845 0 0 12506 4350.05327683573 0 0 15006 5212.12224235313 0 0 20006 6936.26017338811 0 0 25006 8660.39810442261 0 0 30006 10384.5360354576 0 0 40006 13832.8118975273 0 0 50006 17281.0877596008 0 0 60006 20729.3636216604 0 0 80006 27625.9153457978 0 0 100006 34522.4670699267 0 0