/import/home/rettenbs/src/SWE/src/blocks/SWE_WavePropagationBlock.hh

Go to the documentation of this file.

#ifndef SWEWAVEPROPAGATIONBLOCK_HH_
#define SWEWAVEPROPAGATIONBLOCK_HH_

#include "blocks/SWE_Block.hh"
#ifdef DYNAMIC_DISPLACEMENTS
#include "scenarios/Asagi.hpp"
#endif
#include "tools/help.hh"

#include <string>

//which wave propagation solver should be used
//  0: Hybrid
//  1: f-Wave
//  2: Approximate Augmented Riemann solver
//  3: Approximate Augmented Riemann solver which uses the underlying Fortran routines of GeoClaw directly.
#if WAVE_PROPAGATION_SOLVER==1
#include "solvers/FWave.hpp"
#elif WAVE_PROPAGATION_SOLVER==2
#include "solvers/AugRie.hpp"
#elif WAVE_PROPAGATION_SOLVER==3
#include "solvers/AugRieGeoClaw.hpp"
#elif WAVE_PROPAGATION_SOLVER==4
#include "solvers/FWaveVec.hpp"
#else
#include "solvers/Hybrid.hpp"
#endif

class SWE_WavePropagationBlock: public SWE_Block {
//private:
  //OpenMp: Every task defines it own solver -> SWE_WavePropagationBlock.cpp
#ifndef LOOP_OPENMP
    //specify the wave propagation solver
#if WAVE_PROPAGATION_SOLVER==1

    solver::FWave<float> wavePropagationSolver;
#elif WAVE_PROPAGATION_SOLVER==2

    solver::AugRie<float> wavePropagationSolver;
#elif WAVE_PROPAGATION_SOLVER==3

    solver::AugRieGeoClaw<double> wavePropagationSolver;
#elif WAVE_PROPAGATION_SOLVER==4

    solver::FWaveVec<float> wavePropagationSolver;
#else

    solver::Hybrid<float> wavePropagationSolver;
#endif
#endif

    Float2D hNetUpdatesLeft;
    Float2D hNetUpdatesRight;

    Float2D huNetUpdatesLeft;
    Float2D huNetUpdatesRight;

    #if WAVE_PROPAGATION_SOLVER==3

    Float2D hvNetUpdatesLeft;
    Float2D hvNetUpdatesRight;
    #endif

    Float2D hNetUpdatesBelow;
    Float2D hNetUpdatesAbove;

    Float2D hvNetUpdatesBelow;
    Float2D hvNetUpdatesAbove;

    #if WAVE_PROPAGATION_SOLVER==3

    Float2D huNetUpdatesBelow;
    Float2D huNetUpdatesAbove;
    #endif

  public:
    //constructor of a SWE_WavePropagationBlock.
    SWE_WavePropagationBlock(int l_nx, int l_ny,
                                        float l_dx, float l_dy);

    //executes a single timestep.
    virtual void simulateTimestep(float dt);

    //computes the net-updates for the block
    void computeNumericalFluxes();

    //update the cells
    void updateUnknowns(float dt);

    //runs the simulation until i_tEnd is reached.
    float simulate(float i_tStart, float i_tEnd);

    //updates the bathymetry with the current displacment values
#ifdef DYNAMIC_DISPLACEMENTS
    bool updateBathymetryWithDynamicDisplacement(scenarios::Asagi &i_asagiScenario, float time);
#endif

    virtual ~SWE_WavePropagationBlock() {}
};

#endif /* SWEWAVEPROPAGATIONBLOCK_HH_ */