updates to statevector for paper draft

da/sf6/statevector.py

@@ -168,7 +168,7 @@ class SF6StateVector(StateVector):
         """    
 
         if newmean == None: 
-            newmean = np.ones(self.nparams)*0.1   # emma 05-02
+            newmean = np.ones(self.nparams)*0.15   # emma 05-02
 
         if covariancematrix == None: 
             covariancematrix = np.identity(self.nparams)
@@ -241,7 +241,8 @@ class SF6StateVector(StateVector):
             dimgrid = ncf.add_latlon_dim()
             logging.warning('Parameters for TM5 are NOT taken with exponential')
 
-            data = mem.param_values
+            # Explicitly set maximum allowable value to 0.0
+            data = np.where(mem.param_values < 0, 0.0, mem.param_values)
 
             savedict = io.std_savedict.copy()
             savedict['name'] = "parametervalues"