q) + 'lambda' + str(l) + '/'
if not os.path.exists(new_dir): os.makedirs(new_dir)
else:
shutil.rmtree(new_dir)
os.makedirs(new_dir)
shutil.copy(UBC_dir + 'ipinv3d_0' + str(iter_start) + '.chg',
Work_dir)
shutil.copy(UBC_dir + 'ipinv3d.log', Work_dir)
# Read all the input files and extract parameters
meshfile, obsfile, topofile, actv_cells = get_input_files(Home_dir)
# Load mesh and create dimension vectors
numcell, origin, dX, dY, dZ = get_UBC_mesh(Home_dir + meshfile)
nX = numcell[0]
nY = numcell[1]
nZ = numcell[2]
mcell = nX * nY * nZ
actv_cells = array(actv_cells)
#actv_cells.resize((nY,nX,nZ))
# Read *.log and extract information
ndata, phid, Lxyz, beta_in, ref_model, str_ref_model = read_UBC_log(
Work_dir, Home_dir, iter_start, mcell)
# Load reference model
if not os.path.exists(new_dir): os.makedirs(new_dir)
else:
shutil.rmtree(new_dir)
os.makedirs(new_dir)
shutil.copy(UBC_dir + 'maginv3d_00' + str(iter_start) + '.sus',
Work_dir)
#os.rename(work_dir + 'maginv3d_0' + str(iter_start) + '.con','maginv3d_01.con')
shutil.copy(UBC_dir + 'maginv3d.log', Work_dir)
# Read all the input files and extract parameters
meshfile, obsfile, topo_check = get_input_files(Home_dir)
# Load mesh and create dimension vectors
nX, nY, nZ, origin, dX, dY, dZ = get_UBC_mesh(Home_dir + meshfile)
print dX.size
print dY.size
print dZ.size
mcell = nX * nY * nZ
#actv_cells.resize((nY,nX,nZ))
# Read *.log and extract information
ndata, phid, Lxyz, beta_in, ref_model = read_UBC_log_primer(
Work_dir, Home_dir, iter_start, mcell)
## Set-up parameters
# Length scales