PISM, A Parallel Ice Sheet Model  stable v2.1-1-g6902d5502 committed by Ed Bueler on 2023-12-20 08:38:27 -0800
Uniform.cc
Go to the documentation of this file.
1 /* Copyright (C) 2018, 2019, 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 
20 #include "pism/coupler/atmosphere/Uniform.hh"
21 
22 #include "pism/geometry/Geometry.hh"
23 
24 namespace pism {
25 namespace atmosphere {
26 
27 Uniform::Uniform(std::shared_ptr<const Grid> grid)
28  : AtmosphereModel(grid, std::shared_ptr<AtmosphereModel>()) {
31 }
32 
33 void Uniform::init_impl(const Geometry &geometry) {
34  (void) geometry;
35 
36  m_log->message(2, "* Initializing the test atmosphere model...\n");
37 
38  m_temperature->set(m_config->get_number("atmosphere.uniform.temperature", "Kelvin"));
39  m_precipitation->set(m_config->get_number("atmosphere.uniform.precipitation", "kg m-2 s-1"));
40 }
41 
42 void Uniform::update_impl(const Geometry &geometry, double t, double dt) {
43  (void) geometry;
44  (void) t;
45  (void) dt;
46 }
47 
49  return *m_precipitation;
50 }
51 
53  return *m_temperature;
54 }
55 
57  m_precipitation->begin_access();
58  m_temperature->begin_access();
59 }
60 
62  m_precipitation->end_access();
63  m_temperature->end_access();
64 }
65 
66 void Uniform::init_timeseries_impl(const std::vector<double> &ts) const {
67  m_ts_times = ts;
68 }
69 
70 void Uniform::temp_time_series_impl(int i, int j, std::vector<double> &values) const {
71  for (size_t k = 0; k < m_ts_times.size(); ++k) {
72  values[k] = (*m_temperature)(i, j);
73  }
74 }
75 
76 void Uniform::precip_time_series_impl(int i, int j, std::vector<double> &values) const {
77  for (size_t k = 0; k < m_ts_times.size(); ++k) {
78  values[k] = (*m_precipitation)(i, j);
79  }
80 }
81 
82 } // end of namespace atmosphere
83 } // end of namespace pism
std::shared_ptr< const Grid > grid() const
Definition: Component.cc:105
const Config::ConstPtr m_config
configuration database used by this component
Definition: Component.hh:158
const Logger::ConstPtr m_log
logger (for easy access)
Definition: Component.hh:162
static std::shared_ptr< array::Scalar > allocate_temperature(std::shared_ptr< const Grid > grid)
static std::shared_ptr< array::Scalar > allocate_precipitation(std::shared_ptr< const Grid > grid)
std::vector< double > m_ts_times
A purely virtual class defining the interface of a PISM Atmosphere Model.
Uniform(std::shared_ptr< const Grid > g)
Definition: Uniform.cc:27
const array::Scalar & air_temperature_impl() const
Definition: Uniform.cc:52
void init_timeseries_impl(const std::vector< double > &ts) const
Definition: Uniform.cc:66
void init_impl(const Geometry &geometry)
Definition: Uniform.cc:33
void begin_pointwise_access_impl() const
Definition: Uniform.cc:56
void temp_time_series_impl(int i, int j, std::vector< double > &values) const
Definition: Uniform.cc:70
std::shared_ptr< array::Scalar > m_precipitation
Definition: Uniform.hh:47
std::shared_ptr< array::Scalar > m_temperature
Definition: Uniform.hh:47
void update_impl(const Geometry &geometry, double t, double dt)
Definition: Uniform.cc:42
const array::Scalar & precipitation_impl() const
Definition: Uniform.cc:48
void precip_time_series_impl(int i, int j, std::vector< double > &values) const
Definition: Uniform.cc:76
void end_pointwise_access_impl() const
Definition: Uniform.cc:61
static const double k
Definition: exactTestP.cc:42