tInitialization.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
---
tInitialization.hh (2147B)
---
1 /* Copyright (C) 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 POINITIALIZATION_H
21 #define POINITIALIZATION_H
22
23 #include "pism/coupler/OceanModel.hh"
24
25 namespace pism {
26 namespace ocean {
27
28 /*! Ocean model "modifier" that helps with initialization.
29 *
30 * This modifier saves *all* fields a ocean model provides as a part of the model state and re-loads
31 * them during initialization so that they are available *before* the first time step in a
32 * re-started run.
33 *
34 * It is
35 *
36 * - not visible to the user,
37 * - is added automatically, and
38 * - does not have a corresponding "keyword" in ocean::Factory.
39 */
40 class InitializationHelper : public OceanModel {
41 public:
42 InitializationHelper(IceGrid::ConstPtr g, std::shared_ptr<OceanModel> in);
43
44 private:
45 void define_model_state_impl(const File &output) const;
46 void write_model_state_impl(const File &output) const;
47
48 void update_impl(const Geometry &geometry, double t, double dt);
49 void init_impl(const Geometry &geometry);
50
51 const IceModelVec2S& shelf_base_temperature_impl() const;
52 const IceModelVec2S& shelf_base_mass_flux_impl() const;
53 const IceModelVec2S& melange_back_pressure_fraction_impl() const;
54
55 // storage for melange_back_pressure_fraction is inherited from OceanModel
56 IceModelVec2S::Ptr m_shelf_base_temperature;
57 IceModelVec2S::Ptr m_shelf_base_mass_flux;
58 };
59
60 } // end of namespace ocean
61 } // end of namespace pism
62
63 #endif /* POINITIALIZATION_H */