Test file created July 3, 2006 Test seepage against the Bredehoeft & Papadopulos 1965 analytical solution Water Resources Research, Vol. 1 No. 2, equation 4, page 326 Seepage = -0.1 m/yr (up) Use 100 meters of sediment with constant properties ************** 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 10 0.1 500 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)) 100 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 s -0.1 : 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 : node, number of nodes. Following are node depths, initial temps, and production/sink term 0 0 0 0 1 1 0 0 2 2 0 0 3 3 0 0 4 4 0 0 5 5 0 0 6 6 0 0 7 7 0 0 8 8 0 0 9 9 0 0 10 10 0 0 11 11 0 0 12 12 0 0 13 13 0 0 14 14 0 0 15 15 0 0 16 16 0 0 17 17 0 0 18 18 0 0 19 19 0 0 20 20 0 0 21 21 0 0 22 22 0 0 23 23 0 0 24 24 0 0 25 25 0 0 26 26 0 0 27 27 0 0 28 28 0 0 29 29 0 0 30 30 0 0 31 31 0 0 32 32 0 0 33 33 0 0 34 34 0 0 35 35 0 0 36 36 0 0 37 37 0 0 38 38 0 0 39 39 0 0 40 40 0 0 41 41 0 0 42 42 0 0 43 43 0 0 44 44 0 0 45 45 0 0 46 46 0 0 47 47 0 0 48 48 0 0 49 49 0 0 50 50 0 0 51 51 0 0 52 52 0 0 53 53 0 0 54 54 0 0 55 55 0 0 56 56 0 0 57 57 0 0 58 58 0 0 59 59 0 0 60 60 0 0 61 61 0 0 62 62 0 0 63 63 0 0 64 64 0 0 65 65 0 0 66 66 0 0 67 67 0 0 68 68 0 0 69 69 0 0 70 70 0 0 71 71 0 0 72 72 0 0 73 73 0 0 74 74 0 0 75 75 0 0 76 76 0 0 77 77 0 0 78 78 0 0 79 79 0 0 80 80 0 0 81 81 0 0 82 82 0 0 83 83 0 0 84 84 0 0 85 85 0 0 86 86 0 0 87 87 0 0 88 88 0 0 89 89 0 0 90 90 0 0 91 91 0 0 92 92 0 0 93 93 0 0 94 94 0 0 95 95 0 0 96 96 0 0 97 97 0 0 98 98 0 0 99 99 0 0 100 100 0 0