WinglinkTools: work with WinGLink MT inversion software

Created on Mon Aug 22 15:19:30 2011

@author: a1185872

MTpy.core.WinglinkTools.plotResponses(outputfile, maxcol=8, plottype='all', **kwargs)

plotResponse will plot the responses modeled from winglink against the observed data.

Inputs:

outputfile = full path and filename to output file maxcol = maximum number of columns for the plot plottype = ‘all’ to plot all on the same plot

‘1’ to plot each respones in a different figure station to plot a single station or enter as a list of stations to plot a few stations [station1,station2]. Does not have to be verbatim but should have similar unique characters input pb01 for pb01cs in outputfile

Outputs:

None

MTpy.core.WinglinkTools.readModelFile(modelfile, profiledirection='ew')

readModelFile reads in the XYZ txt file output by Winglink.

Inputs:

modelfile = fullpath and filename to modelfile profiledirection = ‘ew’ for east-west predominantly, ‘ns’ for

predominantly north-south. This gives column to fix

MTpy.core.WinglinkTools.readOutputFile(outputfile)

readOutputFile will read an output file from winglink and output data in the form of a dictionary.

Input:

outputfile = the full path and filename of outputfile

Output:

idict = dictionary with keys of station name

each idict[station name] is a dictionary with keys corresponding to modeled and observed responses:

‘obsresxy’,’obsphasexy’,’modresxy’,’modphasexy’,’obsresyx’, ‘obsphaseyx’,’modresyx’,’modphaseyx’,’obshzres’, ‘obshzphase’,’modhzres’,’modhzphase’,’period’

rplst = list of dictionaries for each station with keywords:

‘station’ = station name ‘offset’ = relative offset, ‘resxy’ = TE resistivity and error as row 0 and 1 respectively, ‘resyx’= TM resistivity and error as row 0 and 1 respectively, ‘phasexy’= TE phase and error as row 0 and 1 respectively, ‘phaseyx’= Tm phase and error as row 0 and 1 respectively, ‘realtip’= Real Tipper and error as row 0 and 1 respectively, ‘imagtip’= Imaginary Tipper and error as row 0 and 1 respectively

plst = periodlst as the median of all stations. stationlst = list of stations in order from profile title = list of parameters for plotting as [title,profile,inversiontype]