GeoscienceAustralia / mtpy2

This Repository has been moved to https://github.com/MTgeophysics/mtpy
GNU General Public License v3.0
5 stars 1 forks source link

mtpy.modeling.modem_new.Data.write_data_file() wrong error calculation? #25

Closed YingzhiGou closed 7 years ago

YingzhiGou commented 7 years ago

in mtpy.modeling.modem_new.Data.write_data_file() the code below that calculate the error of vertical component

# compute relative error
if comp.find('t') == 0:
    if 'floor' in self.error_type:
        abs_err = max(self.error_tipper, self.data_array[ss]['tip_err'][ff, 0, z_ii])
    else:
        abs_err = self.error_tipper

where z_ii is always 0.

should I change change to [ff, z_ii, z_jj] or [ff, 0, z_jj] ? @zhang01GA

NOTE: I checked the code from the geophysics/mtpy, the same issue exists