Test file created June21, 2006 Use Variable timing to match Benfield analytical solution Start with 50 meters of sediment (doesn't matter because properties are uniform) Subsidence at 550 m/Ma Match with analytical solution from Benfield, 1949: Qs=Qb(4i^2erfc[(Vt^1/2)/(2kappa^1/2)]) conductivity is 1 thermal capacity is 2.65e6 measure heat flux over top meter Qbase = 1 W/m^2 ************** sub : type of model, subsidence [sub] or sedimentation [sed] 1 : number of stress periods 50 1 2.65e6 0 : 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 : parameter, constants for porosity = f(z). A+Bz+Cz^2+Dz^3+Elnz+Fexp(G*Z)+H^(Iz) km : parameter, porosity = f(z) where z is in [m] or [km] 0.0 : parameter, minimum allowable sediment porosity (decimal) [] [] : 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 3 : 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) ************** 1000000 1000 0.001 5000 1.10 : 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 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.000550 : 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 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) ************** 205 : node, number of nodes. Following are node depths, initial temps, and production/sink term 0 0 0 0 0.5 0.5 0 0 1 1 0 0 1.5 1.5 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 110 110 0 0 120 120 0 0 130 130 0 0 140 140 0 0 150 150 0 0 160 160 0 0 170 170 0 0 180 180 0 0 190 190 0 0 200 200 0 0 210 210 0 0 220 220 0 0 230 230 0 0 240 240 0 0 250 250 0 0 260 260 0 0 270 270 0 0 280 280 0 0 290 290 0 0 300 300 0 0 310 310 0 0 320 320 0 0 330 330 0 0 340 340 0 0 350 350 0 0 360 360 0 0 370 370 0 0 380 380 0 0 390 390 0 0 400 400 0 0 410 410 0 0 420 420 0 0 430 430 0 0 440 440 0 0 450 450 0 0 460 460 0 0 470 470 0 0 480 480 0 0 490 490 0 0 500 500 0 0 600 600 0 0 700 700 0 0 800 800 0 0 900 900 0 0 1000 1000 0 0 1100 1100 0 0 1200 1200 0 0 1300 1300 0 0 1400 1400 0 0 1500 1500 0 0 1600 1600 0 0 1700 1700 0 0 1800 1800 0 0 1900 1900 0 0 2000 2000 0 0 2100 2100 0 0 2200 2200 0 0 2300 2300 0 0 2400 2400 0 0 2500 2500 0 0 2600 2600 0 0 2700 2700 0 0 2800 2800 0 0 2900 2900 0 0 3000 3000 0 0 3100 3100 0 0 3200 3200 0 0 3300 3300 0 0 3400 3400 0 0 3500 3500 0 0 3600 3600 0 0 3700 3700 0 0 3800 3800 0 0 3900 3900 0 0 4000 4000 0 0 4100 4100 0 0 4200 4200 0 0 4300 4300 0 0 4400 4400 0 0 4500 4500 0 0 4600 4600 0 0 4700 4700 0 0 4800 4800 0 0 4900 4900 0 0 5000 5000 0 0 5200 5200 0 0 5400 5400 0 0 5600 5600 0 0 5800 5800 0 0 6000 6000 0 0 6500 6500 0 0 7000 7000 0 0 7500 7500 0 0 8000 8000 0 0 9000 9000 0 0 10000 10000 0 0 12000 12000 0 0 14000 14000 0 0 16000 16000 0 0 20000 20000 0 0 25000 25000 0 0 30000 30000 0 0