if not op.exists(savepath):
os.mkdir(savepath)
# list of stations
slst=[edi[0:-4] for edi in os.listdir(edipath) if edi.find('.edi')>0]
# create an occam data object
ocd = occam2d.Data(edi_path=edipath,
station_list=slst,
interpolate_freq=True,
freq=np.logspace(-3,3,37)
)
ocd.save_path = savepath
ocd.freq_num = 50 # number of frequencies to invert for
##### make data file
# geoelectric strike for rotation, if not specified will calculate from the data
#ocd.geoelectric_strike = 45
# error floors
ocd.res_te_err = 10
ocd.res_tm_err = 10
ocd.phase_te_err = 5
ocd.phase_tm_err = 5
#ocd.model_mode= 4
# -*- coding: utf-8 -*-
"""
Created on Wed Sep 23 13:51:54 2015
@author: jpeacock
"""
import mtpy.modeling.occam2d_rewrite as occam
import numpy as np
s_edi_path = r"d:\Peacock\MTData\SanPabloBay\EDI_Files_dp"
s_list = ["sp{0:02}".format(ii) for ii in range(1, 16)]
ocd = occam.Data(edi_path=s_edi_path, station_list=s_list)
ocd.geoelectric_strike = 0.0
ocd.model_mode = "log_tm_tip"
ocd.phase_tm_err = 2.5
ocd.res_tm_err = 20
ocd.tipper_err = 10
ocd.res_te_err = 20
ocd.phase_te_err = 2.5
ocd.save_path = r"c:\MinGW32-xy\Peacock\occam\SanPabloBay\inv02_tm_tip"
ocd._rotate_to_strike = False
ocd.write_data_file()
ocm = occam.Regularization(station_locations=ocd.station_locations)
ocm.cell_width = 50
ocm.n_layers = 110
ocm.z1_layer = 2
ocm.z_target_depth = 20000
ocm.x_pad_multiplier = 1.5