From 9c5b013a667b7a04cd9df061be9670c6cebe4452 Mon Sep 17 00:00:00 2001 From: Stefano Serafin <serafin@jet01.img.univie.ac.at> Date: Thu, 27 Feb 2025 11:46:13 +0100 Subject: [PATCH] fixed a bug for no-PE runs --- models.py | 19 +++++++++++-------- 1 file changed, 11 insertions(+), 8 deletions(-) diff --git a/models.py b/models.py index 9ffdbc3..9010ae2 100644 --- a/models.py +++ b/models.py @@ -42,15 +42,18 @@ class CBL: self.parameter_true = np.array([self.parameter_true]) self.parameter_ensemble_min = np.array([self.parameter_ensemble_min]) self.parameter_ensemble_max = np.array([self.parameter_ensemble_max]) - self.parameter_inflation_rtps_alpha = np.array([self.parameter_inflation_rtps_alpha]) - if self.parameter_transform == 'pit': - self.parameter_transform = [ transform_pit ] - elif self.parameter_transform == 'log': - self.parameter_transform = [ transform_log ] - elif self.parameter_transform == 'none': - self.parameter_transform = [ transform_none ] + #self.parameter_inflation_rtps_alpha = np.array([self.parameter_inflation_rtps_alpha]) if not hasattr(self,'do_parameter_estimation'): self.do_parameter_estimation == False + if self.do_parameter_estimation: + if self.parameter_transform == 'pit': + self.parameter_transform = [ transform_pit ] + elif self.parameter_transform == 'log': + self.parameter_transform = [ transform_log ] + elif self.parameter_transform == 'none': + self.parameter_transform = [ transform_none ] + else: + self.parameter_transform = [ transform_none ] def initialize(self,argin,coordinates=None): @@ -154,7 +157,7 @@ class CBL: assert (self.parameter_ensemble_min[k]<=self.parameter_ensemble_max[k]),\ f'parameter min must be less or equal to max,\ got {self.parameter_ensemble_min[k],self.parameter_ensemble_max[k]}' - x0[-self.parameter_number:,:] = self.initial_perturbed_parameters + x0[-self.parameter_number:,:] = self.initial_perturbed_parameters # Squeeze ensemble dimension if needed if self.nens==1: -- GitLab