tnode_types.hh - pism - [fork] customized build of PISM, the parallel ice sheet model (tillflux branch)
HTML git clone git://src.adamsgaard.dk/pism
DIR Log
DIR Files
DIR Refs
DIR LICENSE
---
tnode_types.hh (1560B)
---
1 /* Copyright (C) 2016 PISM Authors
2 *
3 * This file is part of PISM.
4 *
5 * PISM is free software; you can redistribute it and/or modify it under the
6 * terms of the GNU General Public License as published by the Free Software
7 * Foundation; either version 3 of the License, or (at your option) any later
8 * version.
9 *
10 * PISM is distributed in the hope that it will be useful, but WITHOUT ANY
11 * WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
12 * FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
13 * details.
14 *
15 * You should have received a copy of the GNU General Public License
16 * along with PISM; if not, write to the Free Software
17 * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
18 */
19
20 #ifndef _NODE_TYPES_H_
21 #define _NODE_TYPES_H_
22
23 namespace pism {
24
25 //! Node types (in the map plane). These are used to implement boundary conditions in the Q1 FEM
26 //! code.
27 /**
28 Note: The condition node_type > 0.5 is also used to detect "exterior" (ice-free) nodes where we
29 prescribe homogeneous Dirichlet B.C. See DirichletData_*::fix_residual_homogeneous.
30
31 This means that interior and boundary types should not use positive values.
32 */
33 enum NodeType {
34 NODE_INTERIOR = -1,
35 NODE_BOUNDARY = 0,
36 NODE_EXTERIOR = 1
37 };
38
39 class IceModelVec2S;
40 class IceModelVec2Int;
41
42 void compute_node_types(const IceModelVec2S &ice_thickness,
43 double thickness_threshold,
44 IceModelVec2Int &result);
45
46 } // end of namespace pism
47
48
49 #endif /* _NODE_TYPES_H_ */