@@ -227,6 +227,8 @@ def mt_list(self, value):
227
227
"""
228
228
At this point not implemented, mainly here for inheritance of MTStations
229
229
"""
230
+ if len (self .values ()) != 0 :
231
+ self .logger .warning ("mt_list cannot be set." )
230
232
pass
231
233
232
234
@property
@@ -262,6 +264,7 @@ def add_station(
262
264
survey = None ,
263
265
compute_relative_location = True ,
264
266
interpolate_periods = None ,
267
+ compute_model_error = False ,
265
268
):
266
269
"""
267
270
Add a MT object
@@ -295,6 +298,18 @@ def add_station(
295
298
296
299
m = m .interpolate (interpolate_periods , bounds_error = False )
297
300
301
+ if compute_model_error :
302
+ m .compute_model_z_errors (
303
+ error_value = self .z_model_error .error_value ,
304
+ error_type = self .z_model_error .error_type ,
305
+ floor = self .z_model_error .floor ,
306
+ )
307
+ m .compute_model_t_errors (
308
+ error_value = self .t_model_error .error_value ,
309
+ error_type = self .t_model_error .error_type ,
310
+ floor = self .t_model_error .floor ,
311
+ )
312
+
298
313
self .__setitem__ (f"{ validate_name (m .survey )} .{ m .station } " , m )
299
314
300
315
if compute_relative_location :
@@ -531,9 +546,7 @@ def interpolate(self, new_periods, f_type="period", inplace=True):
531
546
)
532
547
533
548
else :
534
- mt_data .add_station (
535
- new_mt_obj , compute_relative_location = False
536
- )
549
+ mt_data .add_station (new_mt_obj , compute_relative_location = False )
537
550
538
551
if not inplace :
539
552
return mt_data
@@ -557,9 +570,7 @@ def rotate(self, rotation_angle, inplace=True):
557
570
if not inplace :
558
571
rot_mt_obj = mt_obj .copy ()
559
572
rot_mt_obj .rotation_angle = rotation_angle
560
- mt_data .add_station (
561
- rot_mt_obj , compute_relative_location = False
562
- )
573
+ mt_data .add_station (rot_mt_obj , compute_relative_location = False )
563
574
else :
564
575
mt_obj .rotation_angle = rotation_angle
565
576
@@ -657,12 +668,8 @@ def compute_model_errors(
657
668
self .t_model_error .floor = t_floor
658
669
659
670
for mt_obj in self .values ():
660
- mt_obj .compute_model_z_errors (
661
- ** self .z_model_error .error_parameters
662
- )
663
- mt_obj .compute_model_t_errors (
664
- ** self .t_model_error .error_parameters
665
- )
671
+ mt_obj .compute_model_z_errors (** self .z_model_error .error_parameters )
672
+ mt_obj .compute_model_t_errors (** self .t_model_error .error_parameters )
666
673
667
674
def get_nearby_stations (self , station_key , radius , radius_units = "m" ):
668
675
"""
0 commit comments