def _read_edi_file(self):
"""
read in edi file and set attributes accordingly
"""
self.edi_object = MTedi.Edi(self.fn, datatype=self._data_type)
self.lat = self.edi_object.lat
self.lon = self.edi_object.lon
self.elev = self.edi_object.elev
self.Z = self.edi_object.Z
self.Tipper = self.edi_object.Tipper
self.station = self.edi_object.station
#--> get utm coordinates from lat and lon
self._get_utm()
#--> make sure things are ordered from high frequency to low
self._check_freq_order()
#--> compute phase tensor
self.pt = MTpt.PhaseTensor(z_object=self.Z, freq=self.Z.freq)
#--> compute invariants
self.zinv = MTinv.Zinvariants(z_object=self.Z)
def _set_Z(self, z_object):
"""
set z_object
recalculate phase tensor and invariants, which shouldn't change except
for strike angle
"""
self._Z = z_object
self._Z._compute_res_phase()
#--> compute phase tensor
self.pt = MTpt.PhaseTensor(z_object=self._Z, freq=self._Z.freq)
#--> compute invariants
self.zinv = MTinv.Zinvariants(z_object=self._Z)