tElevationChange.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
---
tElevationChange.hh (2109B)
---
1 // Copyright (C) 2011, 2012, 2013, 2014, 2015, 2016, 2017, 2018, 2020 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 _PAELEVATIONCHANGE_H_
20 #define _PAELEVATIONCHANGE_H_
21
22 #include "pism/coupler/AtmosphereModel.hh"
23
24 #include "pism/util/iceModelVec2T.hh"
25
26 namespace pism {
27 namespace atmosphere {
28
29 class ElevationChange : public AtmosphereModel
30 {
31 public:
32 ElevationChange(IceGrid::ConstPtr g, std::shared_ptr<AtmosphereModel> in);
33 virtual ~ElevationChange();
34
35 protected:
36 void init_impl(const Geometry &geometry);
37 void update_impl(const Geometry &geometry, double t, double dt);
38
39 const IceModelVec2S& mean_precipitation_impl() const;
40 const IceModelVec2S& mean_annual_temp_impl() const;
41
42 void begin_pointwise_access_impl() const;
43 void end_pointwise_access_impl() const;
44
45 void init_timeseries_impl(const std::vector<double> &ts) const;
46 void precip_time_series_impl(int i, int j, std::vector<double> &result) const;
47 void temp_time_series_impl(int i, int j, std::vector<double> &result) const;
48
49 protected:
50 enum Method {SCALE, SHIFT};
51
52 Method m_precip_method;
53 double m_precip_lapse_rate;
54 double m_precip_exp_factor;
55 double m_temp_lapse_rate;
56
57 IceModelVec2T::Ptr m_reference_surface;
58
59 IceModelVec2S::Ptr m_precipitation;
60 IceModelVec2S::Ptr m_temperature;
61 IceModelVec2S m_surface;
62 };
63
64 } // end of namespace atmosphere
65 } // end of namespace pism
66
67 #endif /* _PAELEVATIONCHANGE_H_ */