blob: ae824fe7c3579f9042187885a4ec571f9332d9c4 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
|
// -*- mode:c++ -*-
#ifndef CONSTANTS_HPP
#define CONSTANTS_HPP
#include <iostream>
#include <vector>
#include <array>
#include <string> // string class
#include <algorithm> // min
#include <float.h> // DBL_MAX
#include <math.h> // fabs, sqrt, ceil
#include <fstream> // ifstream
#include <sstream> // stringstream
#include <string.h> // strcmp
#include <random> // random number generator
#include <omp.h>
#include <H5Cpp.h>
#include <H5File.h>
#define DIM 3
#define SURFACE_NCOL 4
#define SUBSURFACE_NCOL 7
#define INFILTRATE_TRUE 1
#define INFILTRATE_FALSE 2
#define SURFACE_MESH_MIXED_ELEMENTS "/0/Mesh/MixedElements"
#define SURFACE_MESH_NODES "/0/Mesh/Nodes"
#define SUBSURFACE_MESH_MIXED_ELEMENTS "/77/Mesh/MixedElements"
#define SUBSURFACE_MESH_NODES "/77/Mesh/Nodes"
#define SURFACE_EXCHANGE_FLUX "surface-surface_subsurface_flux.cell.0"
#define SURFACE_WATER_DEPTH "surface-ponded_depth.cell.0"
#define SURFACE_VELOCITY_X "surface-velocity.cell.0"
#define SURFACE_VELOCITY_Y "surface-velocity.cell.1"
#define SURFACE_EVAPOTRANSPIRATION_FLUX "surface-potential_transpiration.cell.0"
#define SURFACE_AREA "surface-cell_volume.cell.0"
#define DARCY_VELOCITY_X "darcy_velocity.cell.0"
#define DARCY_VELOCITY_Y "darcy_velocity.cell.1"
#define DARCY_VELOCITY_Z "darcy_velocity.cell.2"
#define SUBSURFACE_FLAG 1
#define SURFACE_FLAG 2
#define OUTLET_FLAG 3
#define EVAPORATION_FLAG 4
#define TIME_FUZZ 1.0e-3
#define WATERDEPTH_DRY 1.0e-8
#define ZERO_VELOCITY 1.0e-10
#define DT_MAX 432000.0
#define MOLARVOLUME 1.8e-5
#define CFL 1.0
typedef std::array<double, DIM> point;
typedef std::vector<std::vector<std::vector<int>>> index3d;
typedef std::vector<std::vector<std::vector<double>>> grid3d;
#endif
|
