fixed a bug for no-PE runs

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: