'calado_11\Saltram_L225B366T112D15',
'calado_11\Steena_L182B400T115D15',
'calado_11\Tanker50k_L211B323T117D15',
'calado_11_8\Amazon_L228B400T119D15',
'calado_11_8\BWAustria_L226B366T118D15',
'calado_11_8\DanCisne_L207B322T119D15',
'calado_11_8\Saltram_L225B366T119D15',
'calado_11_8\Tanker50k_L211B323T119D15',
'calado_13\Amazon_L228B400T120D15',
'calado_13\DanCisne_L207B322T130D15',
'calado_13\Saltram_L225B366T120D15',
'calado_13\Tanker50k_L211B323T126D15',
'calado_9_4\Amazon_L228B400T090D15',
'calado_9_4\BWAustria_L226B366T090D15',
'calado_9_4\DanCisne_L207B322T090D15',
'calado_9_4\Saltram_L225B366T090D15',
'calado_9_4\Steena_L182B400T090D15',
'calado_9_4\Tanker50k_L211B323T090D15']
for pst in pastas:
os.chdir(pst)
# m = analysis_wamit.drift_forces_momentum(1)
r = analysis_wamit.raos(1)
os.chdir('..')
nome = pst.split('\\')
plt.savefig(nome[1])
os.chdir('..')
plt.close('all')
multi_fig=True,
T_lim=T_lim,
param_out=p_out)
if wave_forces:
w = analysis_wamit.wave_forces(plota=1,
dof_plot=dof_plot,
inc_plot=inc_plot,
multi_fig=True,
T_lim=T_lim,
param_out=p_out)
if rao:
r = analysis_wamit.raos(plota=1,
dof_plot=dof_plot,
inc_plot=inc_plot,
multi_fig=True,
T_lim=T_lim,
param_out=p_out)
# if damped_rao:
# dr = analysis_wamit.damped_rao(param_out=p_out, mad = a, wf = w)
if drift:
dof_aux = []
for dof in dof_plot:
if dof != 3 and dof != 4 and dof != 5:
dof_aux.append(dof)
for dt in drift_analysis:
d = analysis_wamit.drift_forces(plota=1,
drift_analysis_type=dt,
dof_plot=dof_aux,
p_out = analysis_wamit.output_params()
[params, axis, vol, cb, cg, rest_coef, nome_out, GMt, GMl, M, Bvisc, C,
Cext] = p_out
[g, ulen, rho, water_depth, water_depth_aux, NBODY] = params
if NBODY > 1 and dt == 'm':
raise NameError(
'For multiple bodies the drift analysis must be Pressure or Control Surface!'
)
# Reading RAO
[rao, rao_phase, per, inc, dof, arq4d,
rao_c] = analysis_wamit.raos(param_out=p_out)
# Reading Wave Forces
[wforce, wforce_phase, arq2d] = analysis_wamit.wave_forces(param_out=p_out)
# Reading Drift forces ('m' - Momentum, 'p' - Pressure, 'c' - control surface)
[wdforce, wdforce_phase,
arq8d] = analysis_wamit.drift_forces(drift_analysis_type=dt, param_out=p_out)
# Reading Added Mass and Potential Damping
[added_mass, pot_damp, dof1, arq1d, added_mass_matrix,
pot_damp_matrix] = analysis_wamit.added_mass_pot_damping(param_out=p_out)
dof = np.reshape(dof, (int(np.max(dof) / 6), 6))