script = "run_script_"+str(run_no)+".sh"
timeseries_files = []#["fort.63"]
nontimeseries_files = ["tinun.63", "maxvel.63", "maxele.63"]#, "timemax63"]
# NoNx12/TpN where NoN is number of nodes and TpN is tasks per node, 12 is the
# number of cores per node See -pe line in submission_script <TpN>way<NoN x
# 12>
nprocs = 8 # number of processors per PADCIRC run
ppnode = 12
NoN = 32
num_of_parallel_runs = (ppnode*NoN)/nprocs # procs_pnode * NoN / nproc
domain = dom.domain(grid_dir)
domain.update()
main_run = rmw.runSet(grid_dir, save_dir, basis_dir, num_of_parallel_runs,
base_dir=adcirc_dir, script_name=script)
main_run.initialize_random_field_directories(num_procs=nprocs)
# Set samples
lam_domain = np.array([[.07, .15], [.1, .2]])
lam1 = np.linspace(lam_domain[0, 0], lam_domain[0, 1], 52)
lam2 = np.linspace(lam_domain[1, 0], lam_domain[1, 1], 50)
lam4 = 0.012
lam1, lam2, lam4 = np.meshgrid(lam1, lam2, lam4)
lam_samples = np.column_stack((lam1.ravel(), lam2.ravel(), lam4.ravel()))
mann_pts = lam_samples.transpose()
num_walls = 21
a = run_no*walls_per_run
b = (run_no+1)*walls_per_run
script = "run_script_"+str(run_no)+".sh"
timeseries_files = []#["fort.63"]
nontimeseries_files = ["tinun.63", "maxvel.63", "maxele.63"]#, "timemax63"]
# NoNx12/TpN where NoN is number of nodes and TpN is tasks per node, 12 is the
# number of cores per node See -pe line in submission_script <TpN>way<NoN x
# 12>
nprocs = 8 # number of processors per PADCIRC run
ppnode = 12
NoN = 32
num_of_parallel_runs = (ppnode*NoN)/nprocs # procs_pnode * NoN / nproc
domain = dom.domain(grid_dir)
domain.update()
main_run = rmw.runSet(grid_dir, save_dir, basis_dir, num_of_parallel_runs,
base_dir=adcirc_dir, script_name=script)
main_run.initialize_random_field_directories(num_procs=nprocs)
# Set samples
lam_domain = np.array([[.07, .15], [.1, .2]])
lam1 = np.linspace(lam_domain[0, 0], lam_domain[0, 1], 52)
lam2 = np.linspace(lam_domain[1, 0], lam_domain[1, 1], 50)
lam4 = 0.012
lam1, lam2, lam4 = np.meshgrid(lam1, lam2, lam4)
lam_samples = np.column_stack((lam1.ravel(), lam2.ravel(), lam4.ravel()))
mann_pts = lam_samples.transpose()
num_walls = 21
a = run_no*walls_per_run
b = (run_no+1)*walls_per_run