PIHMv2.2
 All Classes Files Functions Variables Typedefs Macros
pihm.h
Go to the documentation of this file.
1 /**********************************************************************************
2  * File : pihm.h *
3  * Function : Declaration and Definition of global variables and data structure*
4  * Developer of PIHM 2.0: Mukesh Kumar (muk139@psu.edu) *
5  * Developer of PIHM 1.0: Yizhong Qu (quyizhong@gmail.com) *
6  * Version : Nov, 2007 (2.0) *
7  *--------------------------------------------------------------------------------*
8  * *
9  *..............MODIFICATIONS/ADDITIONS in PIHM 2.0...............................*
10  * a) Definition of new variables for ELEMENT data model related to *
11  * i. Subsurface: KsatH, KsatV, infKsatV,infD, RzD, macD, macKsatH, macKsatV *
12  * vAreaF, hAreaF *
13  * ii. ET: LAImax, vegFrac, Albedo, Rs_ref, Rmin *
14  * iii. Snow: meltF *
15  * iv. Surface: dhBydx, dhBYdy, surfH *
16  * b) Definition of new variables for RIVER data model *
17  * i. KsatH, KsatV, bedThick *
18  * c) Definition of new structures: *
19  * i. geology *
20  * ii. Land Cover *
21  * iii. Calibration *
22  * d) Definition of New Control Parameters for each file *
23  *--------------------------------------------------------------------------------*
24  * For questions or comments, please contact *
25  * --> Mukesh Kumar (muk139@psu.edu) *
26  * --> Prof. Chris Duffy (cxd11@psu.edu) *
27  * This code is free for research purpose only. *
28  * Please provide relevant references if you use this code in your research work *
29  *--------------------------------------------------------------------------------*
30  **********************************************************************************/
31 
32 #include "sundials_types.h"
33 #include "nvector_serial.h"
34 
35 /* Definition of Global Variable Types */
36 
38 typedef struct element_type { /* Data model for a triangular element */
39  int index; /* Element No. */
40  int node[3];/* anti-clock-wise */
41  int nabr[3];/* neighbor i shares edge i (0: on boundary) */
42 
43  realtype edge[3];/* edge i is from node i to node i+1 */
44  realtype area; /* area of element */
45 
46  realtype x; /* x of centroid */
47  realtype y; /* y of centroid */
48  realtype zmin; /* z_min of centroid */
49  realtype zmax; /* z_max of centroid */
50 
51  realtype KsatH; /* horizontal geologic saturated hydraulic
52  * conductivity */
53  realtype KsatV; /* vertical geologic saturated hydraulic
54  * conductivity */
55  realtype infKsatV; /* vertical surface saturated
56  * hydraulic conductivity */
57  realtype Porosity;
58  realtype infD; /* depth from ground surface accross which
59  * head is calculated during infiltration */
60  realtype Alpha; /* Alpha from van-genuchten eqn which is
61  * given by satn =
62  * 1/pow(1+pow(abs(Alpha*psi),Beta),1-1/Beta) */
63  realtype Beta;
64  realtype RzD; /* Root zone depth */
65  realtype macD; /* macropore Depth */
66  realtype macKsatH; /* macropore horizontal saturated
67  * hydraulic conductivity */
68  realtype macKsatV; /* macropore vertical saturated
69  * hydraulic conductivity */
70  realtype vAreaF; /* macropore area fraction on a vertical
71  * cross-section */
72  realtype hAreaF; /* macropore area fraction on a horizontal
73  * cross-section */
74  int Macropore; /* 1: macropore; 0: regular soil */
75 
76  realtype LAImax; /* maxm. LAI accross all seasons for a
77  * vegetation type */
78  realtype VegFrac;/* areal vegetation fraction in a triangular
79  * element */
80  realtype Albedo; /* albedo of a triangular element */
81  realtype Rs_ref; /* reference incoming solar flux for
82  * photosynthetically active canopy */
83  realtype Rmin; /* minimum canopy resistance */
84  realtype Rough; /* surface roughness of an element */
85 
86  realtype windH; /* wind measurement height */
87 
88 
89  int soil; /* soil type */
90  int geol; /* geology type */
91  int LC; /* Land Cover type */
92  int IC; /* initial condition type */
93  int BC[3]; /* boundary type. 0:natural bc (no flow);
94  * 1:Dirichlet BC; 2:Neumann BC */
95  int prep; /* precipitation (forcing) type */
96  int temp; /* temperature (forcing) type */
97  int humidity; /* humidity type */
98  int WindVel;/* wind velocity type */
99  int Rn; /* net radiation input */
100  int G; /* radiation into ground */
101  int pressure; /* pressure type */
102  int source; /* source (well) type */
103  int meltF; /* meltFactor */
104  /* for calculation of dh/ds */
105  realtype surfH[3]; /* Total head in neighboring cells */
106  realtype surfX[3]; /* Center X location of neighboring
107  * cells */
108  realtype surfY[3]; /* Center Y location of neighboring
109  * cells */
110  realtype dhBYdx; /* Head gradient in x dirn. */
111  realtype dhBYdy; /* Head gradient in y dirn. */
112 } element;
113 
114 typedef struct nodes_type { /* Data model for a node */
115  int index; /* Node no. */
116 
117  realtype x; /* x coordinate */
118  realtype y; /* y coordinate */
119  realtype zmin; /* z bed rock elevation */
120  realtype zmax; /* z surface elevation */
121 
122 } nodes;
123 
124 typedef struct element_IC_type {/* Initial state variable conditions on each
125  * element */
126  int index;
127 
128  realtype interception; /* Interception storage (Note all
129  * these variables have dimension of
130  * L */
131  realtype snow; /* Snow depth */
132  realtype surf; /* Overland flow depth */
133  realtype unsat; /* unsaturated zone depth */
134  realtype sat; /* saturated zone depth */
135 
136 } element_IC;
137 
138 typedef struct soils_type {
139  int index; /* index */
140 
141  realtype KsatV; /* vertical saturated soil conductivity */
142  realtype ThetaS; /* soil porosity */
143  realtype ThetaR; /* soil moisture residual */
144  realtype Alpha; /* soil curve parameter 1 */
145  realtype Beta; /* soil curve parameter 2 */
146 
147  realtype hAreaF; /* macroporous area fraction on horizontal
148  * section */
149  realtype macKsatV; /* macroporous saturated vertical
150  * conductivity */
151 
152  realtype infD; /* depth from ground surface accross which
153  * head is calculated during infiltration */
154 
155 } soils;
156 
157 typedef struct geol_type {
158  int index; /* index */
159 
160  realtype KsatH; /* horizontal saturated geology conductivity */
161  realtype KsatV; /* vertical saturated geology conductivity */
162  realtype ThetaS; /* geology porosity */
163  realtype ThetaR; /* residual porosity */
164  realtype Alpha; /* van genuchten parameter */
165  realtype Beta; /* van genuchten parameter */
166 
167  realtype vAreaF; /* macroporous area fraction on vertical
168  * section */
169  realtype macKsatH; /* macroporous saturated horizontal
170  * conductivity */
171  realtype macD;
172 
173 } geol;
174 
175 typedef struct lc_type {
176  int index; /* index */
177 
178  realtype LAImax; /* max LAI */
179  realtype VegFrac;/* Canopy Fracn */
180  realtype Albedo; /* Albedo */
181  realtype Rs_ref;
182  realtype Rmin; /* Minimum stomatal resistance */
183  realtype Rough; /* surface roughness factor */
184  realtype RzD; /* rootZone Depth */
185 } LC;
186 
187 typedef struct river_segment_type {
188  int index;
189 
190  realtype x; /* x of river segment */
191  realtype y;
192  realtype zmin; /* bed elevation */
193  realtype zmax; /* bank elevation */
194  realtype depth; /* max depth */
195  realtype Length; /* Riv segment Length */
196  realtype Rough; /* Manning's roughness coeff */
197  realtype KsatH; /* Side conductivity */
198  realtype KsatV; /* Bed conductivity */
199  realtype bedThick;
200  realtype coeff; /* Coefficient c in D = c*pow(B/2,interpOrd)
201  * where D is depth */
202  int FromNode; /* Upstream Node no. */
203  int ToNode; /* Dnstream Node no. */
204  int down; /* down stream segment */
205  int LeftEle;/* Left neighboring element */
206  int RightEle; /* Right neighboring element */
207  int shape; /* shape type */
208  int material; /* material type */
209  int IC; /* IC type */
210  int BC; /* BC type */
212 
213 } river_segment;
214 
215 typedef struct river_shape_type {
216  int index;
217  realtype depth; /* depth */
218  int interpOrd; /* Interpolation order for river
219  * shape: Order =1 (rectangle),
220  * 2(triangle), 3(quadratic) and
221  * 4(cubic) */
222  realtype coeff; /* Coefficient c in D = c*pow(B/2,interpOrd) */
223 
224 } river_shape;
225 
226 typedef struct river_material_type {
227  int index;
228  realtype Rough;
229  realtype Cwr; /* Weir Discharge Coefficient */
230  realtype KsatH; /* Conductivity of river banks */
231  realtype KsatV; /* Conductivity of river bed */
232  realtype bedThick; /* thickeness of conductive river bed */
234 
235 typedef struct river_IC_type {
236  int index;
237  realtype value; /* initial flow depth */
238 
239 } river_IC;
240 
241 typedef struct TSD_type {
242  char name[5];
243  int index;
244  int length; /* length of time series */
245  int iCounter; /* interpolation counter */
246  realtype **TS; /* 2D time series data */
247 
248 } TSD;
249 
250 typedef struct global_calib {
251  realtype KsatH; /* For explanation of each calibration
252  * variable, look for corresponding variables
253  * above */
254  realtype KsatV;
255  realtype infKsatV;
256  realtype macKsatH;
257  realtype macKsatV;
258  realtype infD;
259  realtype RzD;
260  realtype macD;
261  realtype Porosity;
262  realtype Alpha;
263  realtype Beta;
264  realtype vAreaF;
265  realtype hAreaF;
266  realtype Temp;
267  realtype Prep;
268  realtype VegFrac;
269  realtype Albedo;
270  realtype Rough;
271 
272  realtype rivRough;
273  realtype rivKsatH;
274  realtype rivKsatV;
275  realtype rivbedThick;
276  realtype rivDepth;
277  realtype rivShapeCoeff;
278 } globalCal;
279 
280 typedef struct process_control {
281  realtype Et0;
282  realtype Et1;
283  realtype Et2;
284 } processCal;
285 
286 typedef struct model_data_structure { /* Model_data definition */
287  int UnsatMode; /* Unsat Mode */
288  int SurfMode; /* Surface Overland Flow Mode */
289  int RivMode;/* River Routing Mode */
290 
291  int NumEle; /* Number of Elements */
292  int NumNode;/* Number of Nodes */
293  int NumRiv; /* Number of Rivere Segments */
294 
295  int NumPrep;/* Number of Precipatation time series types */
296  int NumTemp;/* Number of Temperature time series types */
297  int NumHumidity; /* Number of Humidity time series
298  * types */
299  int NumWindVel; /* Number of Wind Velocity time
300  * series types */
301  int NumRn; /* Number of Net Radiation time series types */
302  int NumG; /* Number of Ground Heat time series types */
303  int NumP; /* Number of Pressure time series types */
304  int NumSource; /* Number of Source time series types */
305 
306  int NumSoil;/* Number of Soils */
307  int NumGeol;/* Number of Geologies */
308  int NumRes; /* Number of Reservoir */
309  int NumLC; /* Number of Land Cover Index Data */
310 
311  int NumMeltF; /* Number of Melt Factor Time series */
312 
313  int Num1BC; /* Number of Dirichlet BC */
314  int Num2BC; /* Number of Numann BC */
315  int NumEleIC; /* Number of Element Initial Condtion */
316 
317  int NumRivShape; /* Number of River Shape */
318  int NumRivMaterial; /* Number of River Bank/Bed Material */
319  int NumRivIC; /* Number of River Initial Condition */
320  int NumRivBC; /* Number of River Boundary Condition */
321 
322  element *Ele; /* Store Element Information */
323  nodes *Node; /* Store Node Information */
324  element_IC *Ele_IC; /* Store Element Initial Condtion */
325  soils *Soil; /* Store Soil Information */
326  geol *Geol; /* Store Geology Information */
327  LC *LandC; /* Store Land Cover Information */
328 
329  river_segment *Riv; /* Store River Segment Information */
330  river_shape *Riv_Shape; /* Store River Shape Information */
331  river_material *Riv_Mat;/* Store River Bank Material Information */
332  river_IC *Riv_IC; /* Store River Initial Condition */
333 
334  TSD *TSD_Inc;/* Infiltration Capacity Time Series Data */
335  TSD *TSD_LAI;/* Leaves Area Index Time Series Data */
336  //TSD * TSD_DH; /* Zero plane Displacement Height */
337  TSD *TSD_RL; /* Roughness Length */
338  realtype *ISFactor; /* ISFactor is used to calculate
339  * ISMax from LAI */
340  realtype *windH; /* Height at which wind velocity is measured */
341  TSD *TSD_MeltF; /* Monthly Varying Melt Factor for
342  * Temperature Index model */
343 
344  TSD *TSD_EleBC; /* Element Boundary Condition Time
345  * Series Data */
346  TSD *TSD_Riv;/* River Related Time Series Data */
347  TSD *TSD_Prep; /* RainFall Time Series Data */
348  TSD *TSD_Temp; /* Temperature Time Series Data */
349  TSD *TSD_Humidity; /* Humidity Time Series Data */
350  TSD *TSD_WindVel; /* Wind Velocity Time Series Data */
351  TSD *TSD_Rn; /* Net Radiation Time Series Data */
352  TSD *TSD_G; /* Radiation into Ground Time Series Data */
353  TSD *TSD_Pressure; /* Vapor Pressure Time Series data */
354  TSD *TSD_Source; /* Source (well) Time Series data */
355 
356  realtype **FluxSurf; /* Overland Flux */
357  realtype **FluxSub;/* Subsurface Flux */
358  realtype **FluxRiv;/* River Segement Flux */
359 
360  realtype *ElePrep;/* Precep. on each element */
361  realtype *EleETloss;
362  realtype *EleNetPrep; /* Net precep. on each elment */
363  realtype *EleViR; /* Variable infiltration rate */
364  realtype *Recharge; /* Recharge rate to GW */
365  realtype *EleSnow;/* Snow depth on each element */
366  realtype *EleSnowGrnd; /* Snow depth on ground element */
367  realtype *EleSnowCanopy; /* Snow depth on canopy element */
368  realtype *EleIS; /* Interception storage */
369  realtype *EleISmax; /* Maximum interception storage
370  * (liquid precep) */
371  realtype *EleISsnowmax; /* Maximum interception storage
372  * (snow) */
373  realtype *EleTF; /* Through Fall */
374  realtype **EleET; /* Evapo-transpiration (from canopy, ground,
375  * subsurface, transpiration) */
376  realtype Q;
377  realtype *DummyY;
378  realtype *PrintVar[24];
380 } *Model_Data;
381 
382 typedef struct control_data_structure {
383  int Verbose;
384  int Debug;
385 
386  int Solver; /* Solver type */
387  int NumSteps; /* Number of external time steps
388  * (when results can be printed) for
389  * the whole simulation */
390 
391  int gwD; /* File boolean, Choose 1 if your want to
392  * print ground water */
393  int surfD; /* File boolean, Choose 1 if your want to
394  * print overland flow */
395  int snowD; /* File boolean, Choose 1 if your want to
396  * print snow Depth */
397  int rivStg; /* File boolean, Choose 1 if your want to
398  * print river Stage */
399  int Rech; /* File boolean, Choose 1 if your want to
400  * print recharge to ground water */
401  int IsD; /* File boolean, Choose 1 if your want to
402  * print interception depth */
403  int usD; /* File boolean, Choose 1 if your want to
404  * print unsaturated depth */
405  int et[3]; /* File boolean, Choose 1 if your want to
406  * print individual evapo-transpiration
407  * components */
408  int rivFlx[10]; /* File boolean, Choose 1 if your
409  * want to print river/river bed
410  * fluxes */
411 
412  int gwDInt; /* Time interval to output average val of
413  * variables */
414  int surfDInt;
415  int snowDInt;
417  int RechInt;
418  int IsDInt;
419  int usDInt;
420  int etInt;
422 
423  int init_type; /* initialization mode */
424 
425  realtype abstol; /* absolute tolerance */
426  realtype reltol; /* relative tolerance */
427  realtype InitStep; /* initial step size */
428  realtype MaxStep;/* Maximum step size */
429  realtype ETStep; /* Step for et from interception */
430 
431  int GSType, MaxK; /* Maximum Krylov order */
432  realtype delt;
433 
434  realtype StartTime; /* Start time of simulation */
435  realtype EndTime;/* End time of simulation */
436 
437 
438  int outtype;
439  realtype a; /* External time stepping controls */
440  realtype b;
441 
442  realtype *Tout;
443 
444  globalCal Cal; /* Convert this to pointer for localized
445  * calibration */
446 } Control_Data;
447 
448 
449 void FPrintFinalStats(FILE *, long int iopt[], realtype ropt[]);
450 void PrintData(FILE **, Control_Data *, Model_Data, N_Vector, realtype);
realtype x
Definition: pihm.h:190
Definition: pihm.h:241
int rivStgInt
Definition: pihm.h:416
realtype rivShapeCoeff
Definition: pihm.h:277
realtype zmin
Definition: pihm.h:48
realtype Beta
Definition: pihm.h:263
struct element_IC_type element_IC
realtype * EleTF
Definition: pihm.h:373
realtype edge[3]
Definition: pihm.h:43
realtype KsatV
Definition: pihm.h:53
realtype macKsatV
Definition: pihm.h:68
realtype infKsatV
Definition: pihm.h:55
Definition: pihm.h:124
realtype Albedo
Definition: pihm.h:180
realtype * EleNetPrep
Definition: pihm.h:362
realtype y
Definition: pihm.h:118
int gwD
Definition: pihm.h:391
realtype macD
Definition: pihm.h:171
realtype macKsatH
Definition: pihm.h:169
realtype snow
Definition: pihm.h:131
realtype VegFrac
Definition: pihm.h:179
realtype zmin
Definition: pihm.h:119
int humidity
Definition: pihm.h:97
realtype ** EleET
Definition: pihm.h:374
int RechInt
Definition: pihm.h:417
realtype rivRough
Definition: pihm.h:272
int index
Definition: pihm.h:236
realtype Et0
Definition: pihm.h:281
int rivFlx[10]
Definition: pihm.h:408
globalCal Cal
Definition: pihm.h:444
int index
Definition: pihm.h:216
int rivFlxInt
Definition: pihm.h:421
Definition: pihm.h:215
realtype infD
Definition: pihm.h:258
int Verbose
Definition: pihm.h:383
int surfDInt
Definition: pihm.h:414
int ToNode
Definition: pihm.h:203
int GSType
Definition: pihm.h:431
realtype * Tout
Definition: pihm.h:442
realtype Alpha
Definition: pihm.h:144
TSD * TSD_Temp
Definition: pihm.h:348
int NumRn
Definition: pihm.h:301
river_IC * Riv_IC
Definition: pihm.h:332
struct river_material_type river_material
TSD * TSD_EleBC
Definition: pihm.h:344
int NumWindVel
Definition: pihm.h:299
realtype Alpha
Definition: pihm.h:164
realtype ** FluxSurf
Definition: pihm.h:356
int nabr[3]
Definition: pihm.h:41
int Solver
Definition: pihm.h:386
int NumSteps
Definition: pihm.h:387
geol * Geol
Definition: pihm.h:326
int NumRes
Definition: pihm.h:308
realtype Rough
Definition: pihm.h:84
realtype KsatV
Definition: pihm.h:141
int IC
Definition: pihm.h:209
realtype * ElePrep
Definition: pihm.h:360
int shape
Definition: pihm.h:207
int NumRivIC
Definition: pihm.h:319
realtype InitStep
Definition: pihm.h:427
char name[5]
Definition: pihm.h:242
int Rech
Definition: pihm.h:399
int index
Definition: pihm.h:176
int iCounter
Definition: pihm.h:245
realtype zmax
Definition: pihm.h:49
int outtype
Definition: pihm.h:438
int LC
Definition: pihm.h:91
TSD * TSD_Pressure
Definition: pihm.h:353
realtype VegFrac
Definition: pihm.h:78
nodes * Node
Definition: pihm.h:323
realtype interception
Definition: pihm.h:128
realtype zmin
Definition: pihm.h:192
int index
Definition: pihm.h:126
realtype ** FluxSub
Definition: pihm.h:357
realtype a
Definition: pihm.h:439
realtype bedThick
Definition: pihm.h:232
realtype Cwr
Definition: pihm.h:229
realtype infD
Definition: pihm.h:58
realtype macD
Definition: pihm.h:260
int NumSource
Definition: pihm.h:304
int length
Definition: pihm.h:244
struct global_calib globalCal
int NumG
Definition: pihm.h:302
int pressure
Definition: pihm.h:101
realtype * EleViR
Definition: pihm.h:363
realtype hAreaF
Definition: pihm.h:147
int node[3]
Definition: pihm.h:40
realtype abstol
Definition: pihm.h:425
realtype Beta
Definition: pihm.h:145
struct river_segment_type river_segment
realtype Rough
Definition: pihm.h:270
realtype b
Definition: pihm.h:440
int WindVel
Definition: pihm.h:98
realtype ThetaS
Definition: pihm.h:162
int index
Definition: pihm.h:158
processCal pcCal
Definition: pihm.h:379
int NumLC
Definition: pihm.h:309
int NumNode
Definition: pihm.h:292
int NumHumidity
Definition: pihm.h:297
int usD
Definition: pihm.h:403
realtype x
Definition: pihm.h:117
realtype Temp
Definition: pihm.h:266
int index
Definition: pihm.h:188
void PrintData(FILE **, Control_Data *, Model_Data, N_Vector, realtype)
Definition: print.c:122
int LeftEle
Definition: pihm.h:205
LC * LandC
Definition: pihm.h:327
realtype Alpha
Definition: pihm.h:262
realtype depth
Definition: pihm.h:194
int MaxK
Definition: pihm.h:431
int index
Definition: pihm.h:243
Definition: pihm.h:226
int RivMode
Definition: pihm.h:289
int NumRivShape
Definition: pihm.h:317
soils * Soil
Definition: pihm.h:325
realtype Prep
Definition: pihm.h:267
int NumSoil
Definition: pihm.h:306
int Macropore
Definition: pihm.h:74
struct element_type element
void FPrintFinalStats(FILE *, long int iopt[], realtype ropt[])
int interpOrd
Definition: pihm.h:218
realtype ETStep
Definition: pihm.h:429
int Num1BC
Definition: pihm.h:313
int gwDInt
Definition: pihm.h:412
int temp
Definition: pihm.h:96
realtype y
Definition: pihm.h:191
realtype Alpha
Definition: pihm.h:60
TSD * TSD_Prep
Definition: pihm.h:347
int meltF
Definition: pihm.h:103
int soil
Definition: pihm.h:89
realtype * EleISmax
Definition: pihm.h:369
realtype * EleSnowGrnd
Definition: pihm.h:366
int etInt
Definition: pihm.h:420
int NumPrep
Definition: pihm.h:295
realtype ThetaR
Definition: pihm.h:163
int rivStg
Definition: pihm.h:397
int usDInt
Definition: pihm.h:419
realtype LAImax
Definition: pihm.h:76
realtype * windH
Definition: pihm.h:340
int IsDInt
Definition: pihm.h:418
realtype RzD
Definition: pihm.h:64
realtype value
Definition: pihm.h:237
realtype ThetaS
Definition: pihm.h:142
realtype rivbedThick
Definition: pihm.h:275
TSD * TSD_Rn
Definition: pihm.h:351
realtype Beta
Definition: pihm.h:63
int NumP
Definition: pihm.h:303
int SurfMode
Definition: pihm.h:288
int et[3]
Definition: pihm.h:405
realtype windH
Definition: pihm.h:86
realtype y
Definition: pihm.h:47
realtype Porosity
Definition: pihm.h:57
realtype KsatH
Definition: pihm.h:230
realtype Porosity
Definition: pihm.h:261
realtype * ISFactor
Definition: pihm.h:338
int surfD
Definition: pihm.h:393
int NumGeol
Definition: pihm.h:307
realtype Albedo
Definition: pihm.h:269
int NumRivMaterial
Definition: pihm.h:318
realtype rivKsatH
Definition: pihm.h:273
realtype macKsatH
Definition: pihm.h:66
TSD * TSD_WindVel
Definition: pihm.h:350
realtype Rough
Definition: pihm.h:183
int index
Definition: pihm.h:115
struct process_control processCal
struct geol_type geol
realtype EndTime
Definition: pihm.h:435
Definition: pihm.h:175
FILE * riv_state_file
Definition: pihm.h:37
realtype vAreaF
Definition: pihm.h:264
realtype RzD
Definition: pihm.h:259
realtype KsatH
Definition: pihm.h:197
int Debug
Definition: pihm.h:384
realtype ThetaR
Definition: pihm.h:143
realtype * EleETloss
Definition: pihm.h:361
int NumRivBC
Definition: pihm.h:320
struct lc_type LC
int NumRiv
Definition: pihm.h:293
realtype * EleSnow
Definition: pihm.h:365
realtype surf
Definition: pihm.h:132
int material
Definition: pihm.h:208
struct model_data_structure * Model_Data
int RightEle
Definition: pihm.h:206
int Rn
Definition: pihm.h:99
realtype KsatV
Definition: pihm.h:161
int index
Definition: pihm.h:39
Definition: pihm.h:187
int source
Definition: pihm.h:102
realtype StartTime
Definition: pihm.h:434
realtype KsatV
Definition: pihm.h:231
struct nodes_type nodes
realtype coeff
Definition: pihm.h:222
realtype macKsatH
Definition: pihm.h:256
realtype infKsatV
Definition: pihm.h:255
int prep
Definition: pihm.h:95
element_IC * Ele_IC
Definition: pihm.h:324
struct control_data_structure Control_Data
realtype * PrintVar[24]
Definition: pihm.h:378
int reservoir
Definition: pihm.h:211
struct river_shape_type river_shape
realtype rivKsatV
Definition: pihm.h:274
realtype macKsatV
Definition: pihm.h:149
Definition: pihm.h:38
int BC[3]
Definition: pihm.h:93
TSD * TSD_Humidity
Definition: pihm.h:349
int FromNode
Definition: pihm.h:202
realtype LAImax
Definition: pihm.h:178
TSD * TSD_Inc
Definition: pihm.h:334
int IC
Definition: pihm.h:92
river_shape * Riv_Shape
Definition: pihm.h:330
realtype macKsatV
Definition: pihm.h:257
Definition: pihm.h:382
realtype Length
Definition: pihm.h:195
struct soils_type soils
realtype KsatH
Definition: pihm.h:160
realtype RzD
Definition: pihm.h:184
realtype * EleSnowCanopy
Definition: pihm.h:367
realtype surfX[3]
Definition: pihm.h:106
Definition: pihm.h:250
realtype Rs_ref
Definition: pihm.h:181
realtype dhBYdy
Definition: pihm.h:111
TSD * TSD_LAI
Definition: pihm.h:335
realtype surfY[3]
Definition: pihm.h:108
int geol
Definition: pihm.h:90
int G
Definition: pihm.h:100
Definition: pihm.h:138
realtype * DummyY
Definition: pihm.h:377
realtype Et1
Definition: pihm.h:282
int NumMeltF
Definition: pihm.h:311
Definition: pihm.h:280
realtype Rmin
Definition: pihm.h:182
TSD * TSD_G
Definition: pihm.h:352
realtype Et2
Definition: pihm.h:283
int NumTemp
Definition: pihm.h:296
Definition: pihm.h:286
TSD * TSD_MeltF
Definition: pihm.h:341
realtype Rs_ref
Definition: pihm.h:81
realtype sat
Definition: pihm.h:134
realtype x
Definition: pihm.h:46
realtype VegFrac
Definition: pihm.h:268
int snowD
Definition: pihm.h:395
realtype hAreaF
Definition: pihm.h:72
realtype KsatV
Definition: pihm.h:198
realtype Rough
Definition: pihm.h:228
realtype delt
Definition: pihm.h:432
realtype hAreaF
Definition: pihm.h:265
Definition: pihm.h:114
int IsD
Definition: pihm.h:401
int down
Definition: pihm.h:204
realtype * EleIS
Definition: pihm.h:368
realtype * Recharge
Definition: pihm.h:364
int init_type
Definition: pihm.h:423
struct river_IC_type river_IC
TSD * TSD_Source
Definition: pihm.h:354
realtype MaxStep
Definition: pihm.h:428
realtype Q
Definition: pihm.h:376
realtype unsat
Definition: pihm.h:133
int BC
Definition: pihm.h:210
river_material * Riv_Mat
Definition: pihm.h:331
realtype ** FluxRiv
Definition: pihm.h:358
int NumEle
Definition: pihm.h:291
realtype dhBYdx
Definition: pihm.h:110
realtype area
Definition: pihm.h:44
int snowDInt
Definition: pihm.h:415
realtype * EleISsnowmax
Definition: pihm.h:371
realtype zmax
Definition: pihm.h:120
realtype KsatH
Definition: pihm.h:51
realtype zmax
Definition: pihm.h:193
int index
Definition: pihm.h:139
int Num2BC
Definition: pihm.h:314
int UnsatMode
Definition: pihm.h:287
realtype Rough
Definition: pihm.h:196
river_segment * Riv
Definition: pihm.h:329
realtype vAreaF
Definition: pihm.h:167
int NumEleIC
Definition: pihm.h:315
Definition: pihm.h:157
realtype Beta
Definition: pihm.h:165
int index
Definition: pihm.h:227
realtype coeff
Definition: pihm.h:200
realtype reltol
Definition: pihm.h:426
struct TSD_type TSD
realtype macD
Definition: pihm.h:65
realtype ** TS
Definition: pihm.h:246
realtype Rmin
Definition: pihm.h:83
realtype infD
Definition: pihm.h:152
realtype KsatV
Definition: pihm.h:254
realtype KsatH
Definition: pihm.h:251
TSD * TSD_RL
Definition: pihm.h:337
Definition: pihm.h:235
realtype surfH[3]
Definition: pihm.h:105
realtype depth
Definition: pihm.h:217
realtype vAreaF
Definition: pihm.h:70
TSD * TSD_Riv
Definition: pihm.h:346
realtype Albedo
Definition: pihm.h:80
element * Ele
Definition: pihm.h:322
realtype bedThick
Definition: pihm.h:199
realtype rivDepth
Definition: pihm.h:276