PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800
ElevationChange.hh
Go to the documentation of this file.
1 // Copyright (C) 2011, 2012, 2013, 2014, 2015, 2016, 2017, 2018, 2020, 2021, 2022, 2023 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 namespace pism {
25 namespace atmosphere {
26 
28 {
29 public:
30  ElevationChange(std::shared_ptr<const Grid> g, std::shared_ptr<AtmosphereModel> in);
31  virtual ~ElevationChange() = default;
32 
33 protected:
34  void init_impl(const Geometry &geometry);
35  void update_impl(const Geometry &geometry, double t, double dt);
36 
37  const array::Scalar& precipitation_impl() const;
38  const array::Scalar& air_temperature_impl() const;
39 
40  void begin_pointwise_access_impl() const;
41  void end_pointwise_access_impl() const;
42 
43  void init_timeseries_impl(const std::vector<double> &ts) const;
44  void precip_time_series_impl(int i, int j, std::vector<double> &result) const;
45  void temp_time_series_impl(int i, int j, std::vector<double> &result) const;
46 
47 protected:
48  enum Method {SCALE, SHIFT};
49 
51 
52  // parameter used when m_precip_method == SHIFT
54 
55  // parameters used when m_precip_method == SCALE
58 
59  // surface temperature lapse rate
61 
62  std::shared_ptr<array::Forcing> m_reference_surface;
63 
64  std::shared_ptr<array::Scalar> m_precipitation;
65  std::shared_ptr<array::Scalar> m_temperature;
67 };
68 
69 } // end of namespace atmosphere
70 } // end of namespace pism
71 
72 #endif /* _PAELEVATIONCHANGE_H_ */
A purely virtual class defining the interface of a PISM Atmosphere Model.
void precip_time_series_impl(int i, int j, std::vector< double > &result) const
const array::Scalar & precipitation_impl() const
std::shared_ptr< array::Forcing > m_reference_surface
const array::Scalar & air_temperature_impl() const
std::shared_ptr< array::Scalar > m_precipitation
virtual ~ElevationChange()=default
void init_impl(const Geometry &geometry)
ElevationChange(std::shared_ptr< const Grid > g, std::shared_ptr< AtmosphereModel > in)
std::shared_ptr< array::Scalar > m_temperature
void temp_time_series_impl(int i, int j, std::vector< double > &result) const
void init_timeseries_impl(const std::vector< double > &ts) const
void update_impl(const Geometry &geometry, double t, double dt)
static const double g
Definition: exactTestP.cc:36