Loading [MathJax]/extensions/tex2jax.js
PISM, A Parallel Ice Sheet Model 2.2.2-d6b3a29ca committed by Constantine Khrulev on 2025-03-28
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Uniform.cc
Go to the documentation of this file.
1/* Copyright (C) 2018, 2019, 2023, 2024 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
24namespace pism {
25namespace atmosphere {
26
27Uniform::Uniform(std::shared_ptr<const Grid> grid)
28 : AtmosphereModel(grid, std::shared_ptr<AtmosphereModel>()) {
31}
32
33void 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
42void Uniform::update_impl(const Geometry &geometry, double t, double dt) {
43 (void) geometry;
44 (void) t;
45 (void) dt;
46}
47
51
55
60
62 m_precipitation->end_access();
63 m_temperature->end_access();
64}
65
66void Uniform::init_timeseries_impl(const std::vector<double> &ts) const {
67 m_ts_times = ts;
68}
69
70void 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
76void 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
virtual void begin_access() const
Checks if an Array is allocated and calls DAVecGetArray.
Definition Array.cc:568
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)
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