tIBSurfaceModel.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
---
tIBSurfaceModel.hh (2266B)
---
1 // Copyright (C) 2011, 2012, 2013, 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 #ifndef _IBSURFACEMODEL_H_
20 #define _IBSURFACEMODEL_H_
21
22 #include "pism/util/iceModelVec.hh"
23 #include "pism/coupler/AtmosphereModel.hh"
24 #include "pism/coupler/SurfaceModel.hh"
25
26 namespace pism {
27 namespace icebin {
28
29 //! \brief A class implementing a constant-in-time surface model for the surface mass balance.
30 //!
31 //! Reads data from a PISM input file.
32 //!
33 //! Ice surface temperature is parameterized as in PISM-IBSurfaceModel, using a latitude
34 //! and surface elevation-dependent formula.
35
36 class IBSurfaceModel : public pism::surface::SurfaceModel {
37 public:
38 IBSurfaceModel(IceGrid::ConstPtr grid);
39
40 protected:
41 virtual void init_impl(const Geometry &geometry);
42 virtual void update_impl(const Geometry &geometry, double t, double dt);
43 virtual MaxTimestep max_timestep_impl(double t) const;
44
45 virtual void define_model_state_impl(const File &output) const;
46 virtual void write_model_state_impl(const File &output) const;
47
48 virtual const IceModelVec2S& mass_flux_impl() const;
49 virtual const IceModelVec2S& temperature_impl() const;
50
51 protected:
52 bool _initialized;
53 std::string m_input_file;
54
55 public:
56 // Inputs from IceBin
57 pism::IceModelVec2S icebin_wflux;
58 pism::IceModelVec2S icebin_deltah;
59 pism::IceModelVec2S icebin_massxfer; // [kg m-2 s-1]
60 pism::IceModelVec2S icebin_enthxfer; // [J m-2 s-1]
61 // Calculated
62 pism::IceModelVec2S surface_temp;
63 };
64
65 } // end of namespace surface
66 } // end of namespace pism
67
68 #endif /* _IBSURFACEMODEL_H_ */