PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800
exactTestsFG.hh
Go to the documentation of this file.
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 
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 */
TestFGParameters exactFG(double t, double r, const std::vector< double > &z, double Cp)
Definition: exactTestsFG.cc:38
std::vector< double > U
Definition: exactTestsFG.hh:51
std::vector< double > Sig
Definition: exactTestsFG.hh:51
std::vector< double > Sigc
Definition: exactTestsFG.hh:51
std::vector< double > w
Definition: exactTestsFG.hh:51
std::vector< double > T
Definition: exactTestsFG.hh:51