Test file created June21, 2006 Variable timing Validate the model against the analytical solution from Carslaw and Jaeger, 1967, page 96, equation 6. Create a slab (0-100 m). Ends at zero, initial temperature of 100 degrees m-file analytical1.m creates the analytical solution for comparison ************** 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.001 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 : node, number of nodes, depths, initial temps, and production/sink term 0 0 0 0 1 100 0 0 2 100 0 0 3 100 0 0 4 100 0 0 5 100 0 0 6 100 0 0 7 100 0 0 8 100 0 0 9 100 0 0 10 100 0 0 11 100 0 0 12 100 0 0 13 100 0 0 14 100 0 0 15 100 0 0 16 100 0 0 17 100 0 0 18 100 0 0 19 100 0 0 20 100 0 0 21 100 0 0 22 100 0 0 23 100 0 0 24 100 0 0 25 100 0 0 26 100 0 0 27 100 0 0 28 100 0 0 29 100 0 0 30 100 0 0 31 100 0 0 32 100 0 0 33 100 0 0 34 100 0 0 35 100 0 0 36 100 0 0 37 100 0 0 38 100 0 0 39 100 0 0 40 100 0 0 41 100 0 0 42 100 0 0 43 100 0 0 44 100 0 0 45 100 0 0 46 100 0 0 47 100 0 0 48 100 0 0 49 100 0 0 50 100 0 0 51 100 0 0 52 100 0 0 53 100 0 0 54 100 0 0 55 100 0 0 56 100 0 0 57 100 0 0 58 100 0 0 59 100 0 0 60 100 0 0 61 100 0 0 62 100 0 0 63 100 0 0 64 100 0 0 65 100 0 0 66 100 0 0 67 100 0 0 68 100 0 0 69 100 0 0 70 100 0 0 71 100 0 0 72 100 0 0 73 100 0 0 74 100 0 0 75 100 0 0 76 100 0 0 77 100 0 0 78 100 0 0 79 100 0 0 80 100 0 0 81 100 0 0 82 100 0 0 83 100 0 0 84 100 0 0 85 100 0 0 86 100 0 0 87 100 0 0 88 100 0 0 89 100 0 0 90 100 0 0 91 100 0 0 92 100 0 0 93 100 0 0 94 100 0 0 95 100 0 0 96 100 0 0 97 100 0 0 98 100 0 0 99 100 0 0 100 0 0 0