diff --git a/models.py b/models.py index 9ffdbc3f5742dc307ba03d7bedaf4d98995a6496..9010ae2b0c9a049c15ad15c979e0bc4f36196b2f 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: