def test_modular_MPI_NLCG_004(self):
wd = op.normpath(op.join(SAMPLE_DIR, 'ModEM'))
filestem = 'Modular_MPI_NLCG_004'
datafn = 'ModEM_Data.dat'
station = 'pb23'
plot_z = False
ro = PlotResponse(data_fn=op.join(wd, datafn),
resp_fn=op.join(wd, filestem + '.dat'),
plot_type=[station],
plot_z=plot_z)
ro.plot()
if __name__== "__main__":
#### Inputs ####
wd = r'E:\Githubz\Alison_Bugs\data'
savepath = r'E:\Githubz\Alison_Bugs\output'
filestem = 'Modular_MPI_NLCG_108'
datafn = 'ModEM_Data.dat'
Resist_Only = False
# True to plot impedance,
# False for plotting resistivity and phase
respfn = filestem+'.dat'
station = ['GB08','GB09']
ro = PlotResponse(data_fn=os.path.join(wd,datafn),
resp_fn=os.path.join(wd,respfn),
plot_type=station,
plot_style=1,
plot_z= Resist_Only,
save_plots=True,
ctem='b',ctmm='r',
# mtem=
# plot_yn = False,
# fig_size=[3,2],
# font_size=4
)
ro.plot()
"""
import os
import os.path as op
os.chdir(r'C:\Git\mtpy')
from mtpy.modeling.modem import Data, PlotResponse
#### Inputs ####
wd = r'C:\Git\mtpy\examples\model_files\ModEM'
filestem = 'Modular_MPI_NLCG_004'
datafn = 'ModEM_Data.dat'
station = 'pb23'
plot_z = False
################
respfn = filestem + '.dat'
# plot responses at a station
ro = PlotResponse(
data_fn=op.join(wd, datafn),
resp_fn=op.join(wd, respfn),
plot_type=[station],
plot_z=plot_z,
mtmm='-',
mtem='-',
mtmd='.',
mted='.',
# res_limits=(.01,1000)
)
ro.plot()
"""
import os.path as op
#from mtpy.imaging.plot_response import PlotResponse
from mtpy.modeling.modem import PlotResponse
#### Inputs ####
wd = r'C:\mtpywin\mtpy\examples\model_files\ModEM_2'
savepath = r'U:\Software\mtpy\example_plots'
filestem = 'Modular_MPI_NLCG_004'
datafn = 'ModEM_Data.dat'
plot_z = False
respfn = filestem + '.dat'
station = 'Synth02'
ro = PlotResponse(
data_fn=op.join(wd, datafn),
resp_fn=op.join(wd, respfn),
plot_type=[station],
plot_style=3,
plot_z=False,
# fig_size=[3,2],
# font_size=4
)
ro.plot()
#ro.save_figure(r'U:\Software\mtpy\example_plots')
File "mtpy\modeling\modem.py", line 5015, in plot
z_obj._compute_res_phase()
AttributeError: 'Z' object has no attribute '_compute_res_phase'
"""
import os
import os.path as op
os.chdir(r'C:\Git\mtpy')
from mtpy.modeling.modem import PlotResponse
from mtpy.modeling.modem import Data
#### Inputs ####
wd = r'C:\Git\mtpy\examples\model_files\ModEM'
filestem = 'Modular_MPI_NLCG_004'
datafn = 'ModEM_Data.dat'
station = 'pb23'
plot_z = False
################
respfn = filestem + '.dat'
# plot responses at a station
ro = PlotResponse(
data_fn=op.join(wd, datafn),
resp_fn=op.join(wd, respfn),
plot_type=[station],
plot_z=plot_z,
# ctmm='r',ctem='b',
# res_limits=(.01,1000)
)
ro.plot()
#from mtpy.imaging.plot_response import PlotResponse
from mtpy.modeling.modem import PlotResponse
#### Inputs ####
wd = r'C:\mtpywin\mtpy\examples\model_files\ModEM_2'
#wd = r'E:\Githubz\mtpy\examples\model_files\ModEM_2'
savepath = r'C:/tmp'
filestem = 'Modular_MPI_NLCG_004'
datafn = 'ModEM_Data.dat'
plot_z = False
respfn = filestem+'.dat'
station = 'Synth02'
ro = PlotResponse(data_fn=os.path.join(wd,datafn),
resp_fn=os.path.join(wd,respfn),
plot_type=[station],
plot_style=2, # 1 for 4-colums; 2 for 2-columns
plot_z=False,
# fig_size=[3,2],
# font_size=4
)
ro.plot()
ro.save_figure(os.path.join(savepath,'response.png'),
fig_dpi=400) # change fig_dpi to your desired resolution
#### Inputs ####
wd = r'C:\mtpywin\mtpy\examples\model_files\ModEM_2'
#wd = r'E:\Githubz\mtpy\examples\model_files\ModEM_2'
savepath = r'C:/tmp'
filestem = 'Modular_MPI_NLCG_004'
datafn = 'ModEM_Data.dat'
plot_z = False
respfn = filestem+'.dat'
station = 'Synth02'
ro = PlotResponse(data_fn=os.path.join(wd,datafn),
resp_fn=os.path.join(wd,respfn),
plot_type=[station],
plot_style=3, # 1 for 4-columns; 2 for 2-columns, 3 for
# 1-column with diagonals semi-transparent
plot_z=False,
res_limits=[1e-2,1e4],
phase_limits=[-180,180],
shift_yx_phase=False,
fig_size=[3,6],
font_size=10
)
ro.plot()
#ro.save_figure(os.path.join(savepath,'response.png'),
# fig_dpi=400) # change fig_dpi to your desired resolution
from mtpy.modeling.modem import PlotResponse
#### Inputs ####
wd = r'U:\RegionalSurvey\MT046_GeorginaArunta\Modelling\ModEM\GBinv52tip2'
#wd = r'E:\Githubz\mtpy\zprivate\modem_plot_response_issues\GBinv52tip2'
savepath = wd # r'U:\Software\mtpy\example_plots'
filestem = 'Modular_MPI_NLCG_108'
datafn = 'ModEM_Data.dat'
plot_z = False
respfn = filestem + '.dat'
#station = 'GB08' # no tipper data, was not working for plot() with style=1 which is intended to produce 4-column figure
station = 'GB09' # has tipper data
ro = PlotResponse(
data_fn=op.join(wd, datafn),
resp_fn=op.join(wd, respfn),
plot_type=[station],
plot_style=2, # 1 for 4-colums; 2 for 2-columns
plot_z=False,
# fig_size=[3,2],
# font_size=4
)
# print("calling plot() ....")
ro.plot()
#ro.save_figure(op.join(wd,))