DeepSeaGravityWaves.h

 
#pragma once
 
#include <CEnvironment.h>
#include <vector>
#include <CNCField.h>   
using std::vector;
using std::copy;
 
namespace NetCDFsinmod
{
    class Parser;
    struct Grid;
}
 
class LinearWaveTheory;
 
 
#ifdef FH_VISUALIZATION
    #include <CFhCamera.h>
    #include "sfh/ogre/C3DLine.h"
#endif
 
class DeepSeaGravityWaves : public CEnvironment
{
public:
    DeepSeaGravityWaves(string simObjectName, ISimObjectCreator* creator);
    ~DeepSeaGravityWaves(void);
 
   virtual void FinalSetup(const double T, const double *const X, ISimObjectCreator *const creator); 
 
    double  GetSurfaceElevation     (const double dT, const double* const dX, const double adPos[2]);
    void    GetSurfaceUnitNormal    (const double dT, const double* const dX, const double adPos[2], double *adUnitNormalOut);
    void    GetParticleVelocity     (const double dT, const double* const dX, const double adPos[3], double *adVelOut);
    
    void    GetParticleAcceleration (const double dT, const double* const dX, const double adPos[3], double *adAccOut);
    double  GetPressure             (const double dT, const double* const dX, const double adPos[3]);
    double  GetTotalPressure        (const double dT, const double* const dX, const double adPos[3]);
    double  GetSeadepth             (const double adPos[2]);
    double  GetSeabedSigma          (const double adPos[3]);
    double  GetSeabedTau            (const double adPos[3]);
    double  GetSeabedRho            (const double adPos[3]);
    void    RegisterCurrentField    (ICurrentField* currentField);
    void GetWaveFrequencies (std::vector<double> & vec);
    void GetWaveAmplitudes  (std::vector<double> & vec);
    void GetWavePhaseAngles (std::vector<double> & vec);
    void GetWaveDirections  (std::vector<double> & vec);
    void GetWaveNumbers     (std::vector<double> & vec);
    void GetCurrentVelocity (const double dT, const double* const dX, const double adPos[3], double *adVelOut);
    void GetCurrentProducerVelocity (const double dT, const double* const dX, const double adPos[3], double *adVelOut);
 
    double  GetDensity(const double T, const double *const X, const double pos[3]);
 
 
    EnvironmentQuery PointEnvironmentQuery(const double T, const double* const X, const double Position[3]);
    void PointEnvironmentQuery2(const double T, const double* const X, const double Position[3], EnvironmentQuery &q);
 
    #ifdef FH_VISUALIZATION
        void    RenderInit(Ogre::Root* const ogreRoot, ISimObjectCreator* const creator);
        void    RenderUpdate(const double T, const double* const X);
    #endif
 
protected:
    LinearWaveTheory*   m_WaveTheory;   
    string m_waveTheoryName;
    string m_spectrum;
    int m_MaxNumWaves;
    int m_NumWaves;
    int m_NumWaves4Multiplier;          
 
    vector<double> m_adWaveFrequency;       
    vector<double> m_adWaveAmplitude;       
    vector<double> m_adWaveDirection;       
    vector<double> m_adPhaseAngle;          
    vector<double> m_adWaveNumber;          
 
    NetCDFsinmod::Grid *m_grid;
    double* m_depthMap;
    ICurrentField* m_globalCurrentField;
    vector<ICurrentField*> m_registeredCurrentFields;
    double m_length_x;
    double m_length_y;
   
   // ==================== Use NetCDFField (PC)
   bool        m_useNetCDFField;    
   CNCField*    m_ncField;           
   // ==================== Use dummy/artificial velocity field for tanks.
   bool        m_useArtificialTankVelField;  
   double      m_UxyArtField;          
   double      m_UzArtField;           
   double      m_tankRadius;           
   
 
    // ==================== Wake effect (BiaoSu)
    bool             m_wakeEffect;
    int           m_wakeObjNum;
    ISignalPort**  m_wakeObjPort;
    ISignalPort*   m_wakeDirPort;
   ISignalPort*   m_wakeRatioPort;
    double        m_wakeObjRadius[3];
    double        m_wakeObjDepth[3];
   double         m_wakeObjPower;
 
    double         GetWakeRatio(const double T, const double *const X, const double adPos[3]);   
    // ==================== Wake effect (BiaoSu)
 
   // ==================== Avoid crashes due to initial condition (BiaoSu)
   double m_RampT;
   // ==================== Avoid crashes due to initial condition (BiaoSu)
 
    #ifdef FH_VISUALIZATION
        string m_sMaterial;
        Ogre::SceneNode*        m_pWorldSceneNode;      
        Ogre::SceneNode*        m_RenderOceanNode;      
        Ogre::ManualObject* m_pOceanRenderObj;      
        Ogre::Light*            m_pLight;                  
        C3DLine*                   m_lines;                   
        Ogre::ManualObject**    m_depthPolygons;
        Ogre::SceneNode**       m_depthPolygonsNodes;
        CFhCamera*              m_FhCamera;
        Ogre::Camera*           m_OgreCamera;
 
        bool                       m_BoolVisualization;    
    #endif
private:
 
 
 
};