# OrcaWave runs via OrcFxAPI (same as OrcaFlex)
python3 -c "
import OrcFxAPI
model = OrcFxAPI.Model()
model.LoadData('diffraction_model.owd')
model.CalculateStatics()
model.RunSimulation()
model.SaveSimulation('results.owr')
"

# AQWA runs via ANSYS Workbench or command line
# Typical AQWA-LINE execution:
aqwa_line input.dat output.lis

# AQWA-DRIFT for QTFs:
aqwa_drift input.dat output.lis

import OrcFxAPI

def extract_diffraction_results(owd_path):
    """Extract hydrodynamic data from OrcaWave for OrcaFlex.

    OrcaWave results are accessed via VesselType diffraction data,
    not the Vessel instance directly.
    """
    model = OrcFxAPI.Model()
    model.LoadSimulation(owd_path.replace('.owd', '.owr'))

    # Access vessel type (diffraction data lives on VesselType, not Vessel)
    vessel_type = model['Vessel Type1']

    results = {
        "raos": {},
        "added_mass": {},
        "damping": {},
    }

    # Extract RAOs via VesselType hydrodynamic database
    # RAO data is loaded via VesselType.LoadHydrodynamicData() or
    # accessed after OrcaWave populates the vessel type
    dof_names = ['Surge', 'Sway', 'Heave', 'Roll', 'Pitch', 'Yaw']
    for i, dof in enumerate(dof_names):
        results["raos"][dof] = {
            "amplitude": [],
            "phase": [],
        }
        # Note: exact accessor depends on OrcaWave version; see
        # orcawave-to-orcaflex skill for the canonical extraction pattern

    return results