sphere_cg = ((0, 0, -2.0), ) # center of gravity
sphere_name = ('sphere_cpt', ) # body name
sphere_w = np.linspace(0.02, 8.4, 3) # wave frequencies. 420 for full run
sphere_headings = np.linspace(0, 0, 1) # wave heading
sphere_depth = 50.0 # water depth
sphere_ncFile = os.getcwd(
) + os.path.sep + 'test.nc' # path for output .nc file
# ----------------------------------------------------------------------------#
# check that old output is not being overwritten (runs take awhile)
if os.path.isfile(sphere_ncFile):
print(
f'Output ({sphere_ncFile}) file already exists and will be overwritten. '
'Do you wish to proceed? (y/n)')
ans = input()
if ans.lower() != 'y':
print('\nEnding simulation. file not overwritten')
sys.exit(0)
# Run Capytaine
cd, p = cc.call_capy(meshFName=sphere_file,
wCapy=sphere_w,
CoG=sphere_cg,
headings=sphere_headings,
ncFName=sphere_ncFile,
body_name=sphere_name,
depth=sphere_depth,
density=1000.0)
) # mesh files, .dat nemoh, .gdf wamit
rm3_cg = ((0, 0, -0.72), (0, 0, -21.29)) # centers of gravity
rm3_name = ('rm3_float', 'rm3_spar') # body names
rm3_w = np.linspace(0.02, 5.2, 3) # wave frequencies. 260 for full run
rm3_headings = np.linspace(0, 0, 1) # wave headings
rm3_depth = np.infty # water depth
rm3_ncFile = os.getcwd() + os.path.sep + 'test.nc' # path for output .nc file
# ----------------------------------------------------------------------------#
# check that old output is not being overwritten (runs take awhile)
if os.path.isfile(rm3_ncFile):
print(
f'Output ({rm3_ncFile}) file already exists and will be overwritten. '
'Do you wish to proceed? (y/n)')
ans = input()
if ans.lower() != 'y':
print('\nEnding simulation. file not overwritten')
sys.exit(0)
# Run Capytaine
cd, p = cc.call_capy(meshFName=rm3_file,
wCapy=rm3_w,
CoG=rm3_cg,
headings=rm3_headings,
ncFName=rm3_ncFile,
body_name=rm3_name,
depth=rm3_depth,
density=1000.0)
) # mesh file, .dat nemoh, .gdf wamit
bem_cg = ((0, 0, -2.0), ) # center of gravity
bem_name = ('sphere_cpt', ) # body name
bem_w = np.linspace(0.02, 8.4, 3) # wave frequencies. 420 for full run
bem_headings = np.linspace(0, 0, 1) # wave heading
bem_depth = 50.0 # water depth
bem_ncFile = os.getcwd() + os.path.sep + 'test.nc' # path for output .nc file
# ----------------------------------------------------------------------------#
# check that old output is not being overwritten (runs take awhile)
if os.path.isfile(bem_ncFile):
print(
f'Output ({bem_ncFile}) file already exists and will be overwritten. '
'Do you wish to proceed? (y/n)')
ans = input()
if ans.lower() != 'y':
print('\nEnding simulation. file not overwritten')
sys.exit(0)
# Run Capytaine
cc.call_capy(meshFName=bem_file,
wCapy=bem_w,
CoG=bem_cg,
headings=bem_headings,
ncFName=bem_ncFile,
body_name=bem_name,
depth=bem_depth,
density=1000.0)
cubes_cg = ((0, 0, -2.5), (100, 0, -2.5)) # centers of gravity
cubes_name = ('r_cube_capytaine', 't_cube_capytaine') # body names
cubes_w = np.linspace(0.03, 15.0, 3) # wave frequencies. 500 for full run
cubes_headings = np.linspace(0, 90, 10) # wave headings
cubes_depth = 20.0 # water depth
cubes_ncFile = os.getcwd(
) + os.path.sep + 'test.nc' # path for output .nc file
# ----------------------------------------------------------------------------#
# check that old output is not being overwritten (runs take awhile)
if os.path.isfile(cubes_ncFile):
print(
f'Output ({cubes_ncFile}) file already exists and will be overwritten. '
'Do you wish to proceed? (y/n)')
ans = input()
if ans.lower() != 'y':
print('\nEnding simulation. file not overwritten')
sys.exit(0)
# Run Capytaine
cd = cc.call_capy(meshFName=cubes_file,
wCapy=cubes_w,
CoG=cubes_cg,
headings=cubes_headings,
ncFName=cubes_ncFile,
body_name=cubes_name,
depth=cubes_depth,
density=1000.0)
cylinder_cg = ((0, 0, 0), ) # center of gravity
cylinder_name = ('cylinder_cpt', ) # body name
cylinder_w = np.linspace(0.04, 21.0, 3) # wave frequencies. 525 for full run
cylinder_headings = np.linspace(0, 0, 1) # wave heading
cylinder_depth = 3.0 # water depth
cylinder_ncFile = os.getcwd(
) + os.path.sep + 'test.nc' # path for output .nc file
# ----------------------------------------------------------------------------#
# check that old output is not being overwritten (runs take awhile)
if os.path.isfile(cylinder_ncFile):
print(
f'Output ({cylinder_ncFile}) file already exists and will be overwritten. '
'Do you wish to proceed? (y/n)')
ans = input()
if ans.lower() != 'y':
print('\nEnding simulation. file not overwritten')
sys.exit(0)
# Run Capytaine
cc.call_capy(meshFName=cylinder_file,
wCapy=cylinder_w,
CoG=cylinder_cg,
headings=cylinder_headings,
ncFName=cylinder_ncFile,
body_name=cylinder_name,
depth=cylinder_depth,
density=1000.0)
ellipsoid_cg = ((0, 0, 0), ) # center of gravity
ellipsoid_name = ('ellipsoid_cpt', ) # body name
ellipsoid_w = np.linspace(0.03, 9.24, 3) # wave frequencies. 308 for full run
ellipsoid_headings = np.linspace(0, 0, 1) # wave heading
ellipsoid_depth = np.infty # water depth
ellipsoid_ncFile = os.getcwd(
) + os.path.sep + 'test.nc' # path for output .nc file
# ----------------------------------------------------------------------------#
# check that old output is not being overwritten (runs take awhile)
if os.path.isfile(ellipsoid_ncFile):
print(
f'Output ({ellipsoid_ncFile}) file already exists and will be overwritten. '
'Do you wish to proceed? (y/n)')
ans = input()
if ans.lower() != 'y':
print('\nEnding simulation. file not overwritten')
sys.exit(0)
# Run Capytaine
cc.call_capy(meshFName=ellipsoid_file,
wCapy=ellipsoid_w,
CoG=ellipsoid_cg,
headings=ellipsoid_headings,
ncFName=ellipsoid_ncFile,
body_name=ellipsoid_name,
depth=ellipsoid_depth,
density=1000.0)
coercomp_cg = ((0, 0, -0.2), ) # center of gravity
coercomp_name = ('coercomp_cpt', ) # body name
coercomp_w = np.linspace(0.1, 80, 3) # wave frequencies. 800 for full run
coercomp_headings = np.linspace(0, 90, 10) # wave headings
coercomp_depth = 2.20 # water depth
coercomp_ncFile = os.getcwd(
) + os.path.sep + 'test.nc' # path for output .nc file
# ----------------------------------------------------------------------------#
# check that old output is not being overwritten (runs take awhile)
if os.path.isfile(coercomp_ncFile):
print(
f'Output ({coercomp_ncFile}) file already exists and will be overwritten. '
'Do you wish to proceed? (y/n)')
ans = input()
if ans.lower() != 'y':
print('\nEnding simulation. file not overwritten')
sys.exit(0)
# Run Capytaine
cc.call_capy(meshFName=coercomp_file,
wCapy=coercomp_w,
CoG=coercomp_cg,
headings=coercomp_headings,
ncFName=coercomp_ncFile,
body_name=coercomp_name,
depth=coercomp_depth,
density=1000.0)
oswec_cg = ((0, 0, -3.90), (0, 0, -10.90)) # centers of gravity
oswec_name = ('oswec_flap', 'oswec_base') # body names
oswec_w = np.linspace(0.04, 20.0, 3) # wave frequencies. 500 for full run
oswec_headings = np.linspace(0, 90, 10) # wave headings
oswec_depth = 10.90 # water depth
oswec_ncFile = os.getcwd(
) + os.path.sep + 'test.nc' # path for output .nc file
# ----------------------------------------------------------------------------#
# check that old output is not being overwritten (runs take awhile)
if os.path.isfile(oswec_ncFile):
print(
f'Output ({oswec_ncFile}) file already exists and will be overwritten. '
'Do you wish to proceed? (y/n)')
ans = input()
if ans.lower() != 'y':
print('\nEnding simulation. file not overwritten')
sys.exit(0)
# Run Capytaine
cc.call_capy(meshFName=oswec_file,
wCapy=oswec_w,
CoG=oswec_cg,
headings=oswec_headings,
ncFName=oswec_ncFile,
body_name=oswec_name,
depth=oswec_depth,
density=1000.0)
