delta = 1e-10 file_list = [] new_dir = Home_dir + 'lp' + str(p) + 'lq' + str( 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(
file_list = [] new_dir = Home_dir + 'lp' + str(p) + 'lq' + str( 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 + '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)