tTemperatureModel.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
---
tTemperatureModel.hh (2293B)
---
1 /* Copyright (C) 2016, 2017 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 TEMPERATUREMODEL_H
21 #define TEMPERATUREMODEL_H
22
23 #include "EnergyModel.hh"
24
25 namespace pism {
26 namespace energy {
27
28 class TemperatureModel : public EnergyModel {
29 public:
30 TemperatureModel(IceGrid::ConstPtr grid, stressbalance::StressBalance *stress_balance);
31
32 const IceModelVec3 & temperature() const;
33
34 protected:
35 void restart_impl(const File &input_file, int record);
36
37 void bootstrap_impl(const File &input_file,
38 const IceModelVec2S &ice_thickness,
39 const IceModelVec2S &surface_temperature,
40 const IceModelVec2S &climatic_mass_balance,
41 const IceModelVec2S &basal_heat_flux);
42
43 void initialize_impl(const IceModelVec2S &basal_melt_rate,
44 const IceModelVec2S &ice_thickness,
45 const IceModelVec2S &surface_temperature,
46 const IceModelVec2S &climatic_mass_balance,
47 const IceModelVec2S &basal_heat_flux);
48
49 using EnergyModel::update_impl;
50 void update_impl(double t, double dt, const Inputs &inputs);
51
52 void define_model_state_impl(const File &output) const;
53 void write_model_state_impl(const File &output) const;
54
55 void column_drainage(const double rho, const double c, const double L,
56 const double z, const double dz,
57 double *Texcess, double *bwat) const;
58
59 IceModelVec3 m_ice_temperature;
60 };
61
62 } // end of namespace energy
63 } // end of namespace pism
64
65 #endif /* TEMPERATUREMODEL_H */