TsunamiEvent2d.h

Go to the documentation of this file.

 
#ifndef TSUNAMI_LAB_SETUPS_TSUNAMI_EVENT_2D_H
#define TSUNAMI_LAB_SETUPS_TSUNAMI_EVENT_2D_H
 
#include "../io/NetCdf.h"
#include "Setup.h"
#include <string>
 
namespace tsunami_lab {
  namespace setups {
    class TsunamiEvent2d;
  }
}
 
class tsunami_lab::setups::TsunamiEvent2d: public Setup {
  private:
    static constexpr t_real DELTA = 20;
 
    io::NetCdf::Grid m_bathy;
 
    io::NetCdf::Grid m_displ;
 
    t_real m_offsetX = 0;
    t_real m_offsetY = 0;
 
    t_real displacement( t_real i_xFile, t_real i_yFile ) const;
 
  public:
    TsunamiEvent2d( std::string i_bathyPath,
                    std::string i_displPath );
 
    t_real getHeight( t_real i_x, t_real i_y ) const;
 
    t_real getMomentumX( t_real, t_real ) const;
 
    t_real getMomentumY( t_real, t_real ) const;
 
    t_real getBathymetry( t_real i_x, t_real i_y ) const;
 
    t_real domainLengthX() const {
      return m_bathy.x.empty() ? 0 : (m_bathy.x.back() - m_bathy.x.front());
    }
 
    t_real domainLengthY() const {
      return m_bathy.y.empty() ? 0 : (m_bathy.y.back() - m_bathy.y.front());
    }
};
 
#endif