From 225324cc5b6f4c94ce6072ec14942eceaf15a645 Mon Sep 17 00:00:00 2001 From: Nicolas Mounet Date: Mon, 9 Nov 2020 12:39:16 +0100 Subject: [PATCH] Changing default mesh of temperatues to one point every 6 minutes (Geneva temperature) --- cara/models.py | 13 +++++++++++-- 1 file changed, 11 insertions(+), 2 deletions(-) diff --git a/cara/models.py b/cara/models.py index 33776079..471e9739 100644 --- a/cara/models.py +++ b/cara/models.py @@ -137,9 +137,18 @@ class PiecewiseConstant: # Geneva hourly temperatures as piecewise constant function (in Kelvin) -GenevaTemperatures = { +GenevaTemperatures_hourly = { month: PiecewiseConstant(tuple(np.arange(25.)), - tuple(273.15+np.array(temperatures))) + tuple(273.15+np.array(temperatures))) + for month,temperatures in Geneva_hourly_temperatures_celsius_per_hour.items() +} +# same temperatures on a finer temperature mesh +refine_factor = 10 +refined_times = np.linspace(0.,24.,24*refine_factor+1) +GenevaTemperatures = { + month: PiecewiseConstant(tuple(refined_times), + tuple(273.15+np.interp(refined_times[:-1],np.arange(25.), + temperatures+[temperatures[0]]))) for month,temperatures in Geneva_hourly_temperatures_celsius_per_hour.items() }