ocr.n_layers = 60 # cell width to aim for, note this is the mesh size (2 mesh blocks per model block) ocr.cell_width = 1250 # controls size of padding ocr.x_pad_multiplier = 1.9 # controls aspect ratio of blocks ocr.trigger= 0.25 #ocr.z_bottom = 200000 # z1 layer and target depth in metres ocr.z1_layer = 50 ocr.z_target_depth = 80000 ocr.save_path=ocd.save_path ocr.build_mesh() ocr.build_regularization() ocr.write_mesh_file() ocr.write_regularization_file() ocr.plot_mesh() #make startup file ocs=occam2d.Startup() ocs.iterations_to_run=40 ocs.data_fn=op.join(ocd.save_path,'OccamDataFile.dat') ocs.resistivity_start=2.0 ocr.get_num_free_params() ocs.param_count=ocr.num_free_param ocs.save_path=ocd.save_path ocs.model_fn=ocr.reg_fn ocs.write_startup_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
ocm.num_x_pad_cells = 9
ocm.num_x_pad_small_cells = 5
ocm.build_mesh()
ocm.save_path = ocd.save_path
ocm.write_mesh_file()
ocm.build_regularization()
ocm.write_regularization_file()
ocm.plot_mesh()
ocs = occam.Startup()
ocs.data_fn = ocd.data_fn
ocs.resistivity_start = 1
ocs.model_fn = ocm.reg_fn
ocs.param_count = ocm.num_free_param
ocs.save_path = ocd.save_path
ocs.write_startup_file()
# mask large tipper values in data file
data_fn = ocd.data_fn
ocd = occam.Data()
ocd.read_data_file(data_fn)
for s_dict in ocd.data:
re_index = np.where(abs(s_dict["re_tip"][0]) > 1.0)
im_index = np.where(abs(s_dict["im_tip"][0]) > 1.0)