Test file created June 23, 2006 Non-uniform spacing The region -L0 ************** sub : type of model, subsidence [sub] or sedimentation [sed] 1 : number of stress periods [] [] [] [] : layer1 parameters: depth of node at SBI (m), basement conductivity (W/m-K), thermal capacity (J/m^3-K), porosity (decimal) [] [] [] [] : layer2 parameters: depth of node at layer1/layer2 interface, basement conductivity (W/m-K), thermal capacity (J/m^3-K), porosity (decimal) 2.65e6 : parameter, thermal capacity of water (J/m^3-K) 2.65e6 : parameter, thermal capacity of sediment (J/m^3-K) 1 : parameter, thermal conductivity of water (W/m-K) 1 : parameter, thermal conductivity of sediment grains (W/m-K) 0. : parameter, surface sediment porosity (decimal) 0 0 0 0 0 0 0 0 0 0 0 : parameter, constants for porosity = f(z). A+Bz+Cz^2+Dz^3+Elnz+Fexp(G*Z)+H^(Iz)+Jz^K km : parameter, porosity = f(z) where z is in [m] or [km] 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 0 : 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) ************** 500 1 0.1 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) T : boundary condition (lower), ([T] for Temperature (degrees C), [q] for heat flow (W/m^2)) 0 0 0 0 0 0 : boundary condition paramaters for lower boundary: [T] or [q] = f(time in yrs): A+B(time)+C(time^2)+D(time^3)+E/sqrt(F*time) 0 : boundary condition (upper) , temperature in degrees C that boundary is held constant at. 0.00000 : 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 0 0 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 1 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 ************** 101 0 0 0 0 1 0 0.1 0 2 0 0.1 0 3 0 0.1 0 4 0 0.1 0 5 0 0.1 0 6 0 0.1 0 7 0 0.1 0 8 0 0.1 0 9 0 0.1 0 10 0 0.1 0 11 0 0.1 0 12 0 0.1 0 13 0 0.1 0 14 0 0.1 0 15 0 0.1 0 16 0 0.1 0 17 0 0.1 0 18 0 0.1 0 19 0 0.1 0 20 0 0.1 0 21 0 0.1 0 22 0 0.1 0 23 0 0.1 0 24 0 0.1 0 25 0 0.1 0 26 0 0.1 0 27 0 0.1 0 28 0 0.1 0 29 0 0.1 0 30 0 0.1 0 31 0 0.1 0 32 0 0.1 0 33 0 0.1 0 34 0 0.1 0 35 0 0.1 0 36 0 0.1 0 37 0 0.1 0 38 0 0.1 0 39 0 0.1 0 40 0 0.1 0 41 0 0.1 0 42 0 0.1 0 43 0 0.1 0 44 0 0.1 0 45 0 0.1 0 46 0 0.1 0 47 0 0.1 0 48 0 0.1 0 49 0 0.1 0 50 0 0.1 0 51 0 0.1 0 52 0 0.1 0 53 0 0.1 0 54 0 0.1 0 55 0 0.1 0 56 0 0.1 0 57 0 0.1 0 58 0 0.1 0 59 0 0.1 0 60 0 0.1 0 61 0 0.1 0 62 0 0.1 0 63 0 0.1 0 64 0 0.1 0 65 0 0.1 0 66 0 0.1 0 67 0 0.1 0 68 0 0.1 0 69 0 0.1 0 70 0 0.1 0 71 0 0.1 0 72 0 0.1 0 73 0 0.1 0 74 0 0.1 0 75 0 0.1 0 76 0 0.1 0 77 0 0.1 0 78 0 0.1 0 79 0 0.1 0 80 0 0.1 0 81 0 0.1 0 82 0 0.1 0 83 0 0.1 0 84 0 0.1 0 85 0 0.1 0 86 0 0.1 0 87 0 0.1 0 88 0 0.1 0 89 0 0.1 0 90 0 0.1 0 91 0 0.1 0 92 0 0.1 0 93 0 0.1 0 94 0 0.1 0 95 0 0.1 0 96 0 0.1 0 97 0 0.1 0 98 0 0.1 0 99 0 0.1 0 100 0 0 0