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