texactTestsFG.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
---
texactTestsFG.hh (2093B)
---
1 /*
2 Copyright (C) 2004-2006, 2014, 2016 Jed Brown and Ed Bueler and Constantine Khroulev
3
4 This file is part of PISM.
5
6 PISM is free software; you can redistribute it and/or modify it under the
7 terms of the GNU General Public License as published by the Free Software
8 Foundation; either version 3 of the License, or (at your option) any later
9 version.
10
11 PISM is distributed in the hope that it will be useful, but WITHOUT ANY
12 WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
13 FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
14 details.
15
16 You should have received a copy of the GNU General Public License
17 along with PISM; if not, write to the Free Software
18 Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
19 */
20
21 #ifndef EXACTTESTSFG_H
22 #define EXACTTESTSFG_H
23
24 #include <vector>
25
26 namespace pism {
27
28 /*
29 ELB 9/12/05; 05/12/06; 10/14/06; 5/30/08
30 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
31 ! exactTestsFG is a C++ implementation of the exact solutions Test F & G for a
32 ! thermocoupled ice sheet. References:
33 !
34 ! Ed Bueler, Jed Brown, and Craig Lingle, "Exact solutions to the
35 ! thermomechanically coupled shallow ice approximation: effective
36 ! tools for verification," J. Glaciol. 53 (182), 499--516.
37 !
38 ! Ed Bueler and Jed Brown, "On exact solutions for cold, shallow, and
39 ! thermocoupled ice sheets," preprint arXiv:physics/0610106, 2006
40 !
41 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
42 */
43
44 struct TestFGParameters {
45 TestFGParameters(int Mz)
46 : T(Mz), U(Mz), w(Mz), Sig(Mz), Sigc(Mz) {
47 // empty
48 }
49
50 double H, M;
51 std::vector<double> T, U, w, Sig, Sigc;
52 };
53
54 TestFGParameters exactFG(double t, double r, const std::vector<double> &z, double Cp);
55
56 /*
57 * NOTE: Units returned for Sig and Sigc are K/s (i.e. temperature) not J/s.
58 * This matches the published sources above but requires conversion in
59 * PISM as of revision 311.
60 */
61
62 } // end of namespace pism
63
64 #endif /* EXACTTESTSFG_H */