tPSVerification.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
---
tPSVerification.hh (2192B)
---
1 /* Copyright (C) 2014, 2015, 2016, 2017, 2018 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 _PSVERIFICATION_H_
21 #define _PSVERIFICATION_H_
22
23 #include "pism/coupler/surface/Formulas.hh"
24 #include "pism/util/EnthalpyConverter.hh"
25
26 namespace pism {
27 namespace surface {
28
29 //! Climate inputs for verification tests.
30 class Verification : public PSFormulas {
31 public:
32 Verification(IceGrid::ConstPtr g, EnthalpyConverter::Ptr EC, int test);
33 ~Verification();
34 private:
35 void init_impl(const Geometry &geometry);
36 void update_impl(const Geometry &geometry, double t, double dt);
37
38 void define_model_state_impl(const File &output) const;
39 void write_model_state_impl(const File &output) const;
40
41 MaxTimestep max_timestep_impl(double t) const;
42
43 int m_testname;
44 EnthalpyConverter::Ptr m_EC;
45 void update_ABCDH(double t);
46 void update_FG(double t);
47 void update_KO();
48 void update_L();
49 void update_V();
50
51 // FIXME: get rid of code duplication (these constants are defined
52 // both here and in IceCompModel).
53 static const double secpera, ablationRateOutside;
54
55 static const double ST; // K m^-1; surface temperature gradient: T_s = ST * r + Tmin
56 static const double Tmin; // K; minimum temperature (at center)
57 static const double LforFG; // m; exact radius of tests F&G ice sheet
58 static const double ApforG; // m; magnitude A_p of annular perturbation for test G;
59
60 };
61
62 } // end of namespace surface
63 } // end of namespace pism
64
65 #endif /* _PSVERIFICATION_H_ */