@@ -169,13 +169,13 @@ def physical_params(self, **kwargs):
169169 known = [getattr (self , i ) for i in self .physical_parameters ]
170170 return {i .name : kwargs .get (i .name , i ) or i for i in known }
171171
172- def _gen_phys_param (self , field , name , space_order , is_param = True ,
173- default_value = 0 , avg_mode = 'arithmetic' , ** kwargs ):
172+ def _gen_phys_param (self , field , name , space_order , default_value = 0 ,
173+ avg_mode = 'arithmetic' , ** kwargs ):
174174 if field is None :
175175 return default_value
176176 if isinstance (field , np .ndarray ):
177177 function = Function (name = name , grid = self .grid , space_order = space_order ,
178- parameter = is_param , avg_mode = avg_mode )
178+ avg_mode = avg_mode )
179179 initialize_function (function , field , self .padsizes )
180180 else :
181181 function = Constant (name = name , value = field , dtype = self .grid .dtype )
@@ -306,9 +306,8 @@ def _initialize_physics(self, vp, space_order, **kwargs):
306306 # Initialize elastic with Lame parametrization
307307 if 'vs' in kwargs :
308308 vs = kwargs .pop ('vs' )
309- self .lam = self ._gen_phys_param ((vp ** 2 - 2. * vs ** 2 )/ b , 'lam' , space_order ,
310- is_param = True )
311- self .mu = self ._gen_phys_param (vs ** 2 / b , 'mu' , space_order , is_param = True ,
309+ self .lam = self ._gen_phys_param ((vp ** 2 - 2. * vs ** 2 )/ b , 'lam' , space_order )
310+ self .mu = self ._gen_phys_param (vs ** 2 / b , 'mu' , space_order ,
312311 avg_mode = 'harmonic' )
313312 else :
314313 # All other seismic models have at least a velocity
0 commit comments