pad=1) m_elev[np.where(m_elev > 10)] = 10.0 ### 3.) make a resistivity model that incoorporates topography mod_elev, elev_nodes_z = modem.make_elevation_model(m_elev, m_obj.nodes_z, elevation_cell=elev_cell) ### 4.) write new model file m_obj.nodes_z = elev_nodes_z m_obj.res_model = mod_elev m_obj.write_model_file(save_path=sv_path, model_fn_basename='sm_hydro_topography.rho') #write new data file n_dfn = modem.change_data_elevation(d_obj.data_fn, m_obj.model_fn) # write covariance file cov = modem.Covariance() #cov.grid_dimensions = m_obj.res_model.shape #cov.mask_arr = np.ones_like(m_obj.res_model) #cov.mask_arr[np.where(m_obj.res_model > 1e10)] = 0 cov.smoothing_east = 0.4 cov.smoothing_north = 0.4 cov.smoothing_z = 0.4 cov.save_path = sv_path cov.write_covariance_file(model_fn=m_obj.model_fn) #m_obj.write_vtk_file() #d_obj.write_vtk_station_file()
mod_obj.n_layers = 40 mod_obj.z_target_depth = 40000.0 mod_obj.make_mesh() mod_obj.plot_mesh() mod_obj.write_model_file( model_fn=os.path.join(save_path, r"mshs_modem_sm.rho")) inv_period_list = np.logspace(-np.log10(625.0), 3, num=23) data_obj = modem.Data( edi_list=s_edi_list, station_locations=mod_obj.station_locations, period_list=inv_period_list, ) data_obj.error_type = "egbert" data_obj.error_egbert = 10.0 data_obj.error_tipper = 0.05 data_obj.get_mt_dict() data_obj._fill_data_array() data_obj.data_array["elev"][:] = 0.0 data_obj.write_data_file(save_path=save_path, fn_basename="mshs_modem_data_err10_old.dat") cov = modem.Covariance(grid_dimensions=mod_obj.res_model.shape) cov.smoothing_east = 0.4 cov.smoothing_north = 0.4 cov.smoothing_z = 0.4 cov.write_covariance_file(os.path.join(save_path, "covariance_old.cov"))