DynamicDeepSeaGravityWaves.h

 
#pragma once
 
#include "DeepSeaGravityWaves.h"
 
//#include <vector>
//using std::vector;
//using std::copy;
 
 
struct WaveInputInfo_Spectrum {
    string waveEnergySpecter;
    double LongCrestedPower;
    bool isLongCrested;
    double meanWavePeriod;
    double significantWaveHeight;
    double nDimOmegaLowLim;
    double nDimOmegaUppLim;
    double mainWaveDirectionRad;
    int randomSeed;
};
 
 
class DynamicDeepSeaGravityWaves : public DeepSeaGravityWaves
{
public:
 
    DynamicDeepSeaGravityWaves(string simObjectName, ISimObjectCreator* creator);
    //~DynamicDeepSeaGravityWaves(void);
 
    //void OdeFcn(const double dT, const double* const adX, double* const adXDot, const bool bIsMajorTimeStep) { std::cout << "working"; };
 
    
    //Pre ode-function pass for handling discontinuous/discrete states - in this case, there are no states but environmental parameters.
    void PreOdeFcn(const double T, const double *const X, IStateUpdater* updater);
 
#ifdef FH_VISUALIZATION
    void RenderUpdate(const double T, const double* const X);
#endif
 
 
private:
    
    // Updates m_WaveTheory according to new parameters (replaces m_WaveTheory)
    void updateWaveSpectrum();
 
 
    // vectors keeping "old" values for current modelling
    vector<double> m_current_velocity;      
    vector<double> m_current_direction;     
    vector<double> m_depth_layers;          
 
    // wave parameters that can be updated in "real time" if m_updateWaveField is TRUE
    WaveInputInfo_Spectrum m_waveInfo_Spectrum;
 
    double  m_cutoff; 
    int     m_NumWavesWithoutPadding;
    double  m_updateStartTime;
    bool    m_printUpdatedParameters;
    bool    m_waveUpdateIsNeeded;
 
    // input ports for updating Current parameters
    ISignalPort*    m_In_current_velocity;
    ISignalPort*    m_In_current_direction;
 
    // input ports for updating wave spectrum (if spectrum is defined)
    ISignalPort*    m_In_significantWaveHeight;
    ISignalPort*    m_In_mainWaveDirectionRad;
    ISignalPort*    m_In_meanWavePeriod;
 
    // input ports for updating wave components (if components are defined)
    ISignalPort*    m_In_waveFrequency;
    ISignalPort*    m_In_waveAmplitude;
    ISignalPort*    m_In_waveDirection;
    ISignalPort*    m_In_phaseAngle;
 
#ifdef FH_VISUALIZATION
    // updates visualization if new wave parameters are set 
    void updateWaveRenderParams(); 
#endif
 
};