tForceThickness.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
---
tForceThickness.hh (2258B)
---
1 // Copyright (C) 2011, 2012, 2013, 2014, 2015, 2016, 2017, 2018, 2019 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 _PSFORCETHICKNESS_H_
20 #define _PSFORCETHICKNESS_H_
21
22 #include "pism/coupler/SurfaceModel.hh"
23
24 namespace pism {
25
26 class IceModelVec2CellType;
27
28 namespace surface {
29
30 //! A class implementing a modified surface mass balance which forces
31 //! ice thickness to a given target by the end of the run.
32 class ForceThickness : public SurfaceModel {
33 public:
34 ForceThickness(IceGrid::ConstPtr g, std::shared_ptr<SurfaceModel> input);
35 virtual ~ForceThickness();
36 protected:
37 void init_impl(const Geometry &geometry);
38 void update_impl(const Geometry &geometry, double t, double dt);
39
40 void define_model_state_impl(const File &output) const;
41 void write_model_state_impl(const File &output) const;
42
43 const IceModelVec2S& mass_flux_impl() const;
44
45 const IceModelVec2S& accumulation_impl() const;
46 const IceModelVec2S& melt_impl() const;
47 const IceModelVec2S& runoff_impl() const;
48
49 MaxTimestep max_timestep_impl(double t) const;
50 private:
51 void adjust_mass_flux(double time,
52 const IceModelVec2S &ice_thickness,
53 const IceModelVec2CellType &cell_type,
54 IceModelVec2S &result) const;
55
56 double m_alpha, m_alpha_ice_free_factor, m_ice_free_thickness_threshold;
57 double m_start_time;
58 IceModelVec2S m_target_thickness;
59 IceModelVec2Int m_ftt_mask;
60
61 IceModelVec2S::Ptr m_mass_flux;
62 };
63
64 } // end of namespace surface
65 } // end of namespace pism
66
67 #endif /* _PSFORCETHICKNESS_H_ */