weno/applications/euler/3d/Wedge-Dimotakis/src/Problem.h

// -*- C++ -*-

#ifndef AMROC_PROBLEM_H
#define AMROC_PROBLEM_H

#define DIM  3
#define NEQUATIONS 10  // Euler equations for gases in 3D (5 fields), 
                       // +2 fields for passive scalars
                       // +1 feild for temperature
                       // +1 feild to output dcflag and 
                       // +1 fieldsgske
#define NVARS 7        // Fixup used only for vector of state and 2 scalars
  
#include "WENOProblem.h"
//#include "WENOStdProblem.h"

#define OWN_GFMAMRSOLVER
#include "WENOStdGFMProblem.h"

#define f_init_commongm FORTRAN_NAME(comblgm, COMBLGM)

extern "C" {
  void f_init_commongm(const INTEGER& dim, INTEGER* shape, DOUBLE* geom, 
                       INTEGER& maxlev, INTEGER* refine);
}

class SolverSpecific : 
  public AMRGFMSolver<VectorType,FixupType,FlagType,DIM> {
  typedef VectorType::InternalDataType DataType;
  typedef AMRGFMSolver<VectorType,FixupType,FlagType,DIM> base;
  typedef WENOF77GFMFileOutput<VectorType,FixupType,FlagType,DIM> output_type;
public:

  SolverSpecific(IntegratorSpecific& integ, 
                 base::initial_condition_type& init,
                 base::boundary_conditions_type& bc) : base(integ, init, bc) {
    SetLevelTransfer(new F77LevelTransfer<VectorType,DIM>(f_prolong, f_restrict));
#ifdef f_flgout
    SetFileOutput(new WENOF77GFMFileOutput<VectorType,FixupType,FlagType,DIM>(*this,f_flgout,f_bounds)); 
#else   
    SetFileOutput(new GFMFileOutput<VectorType,FixupType,FlagType,DIM>(*this)); 
#endif
    SetFixup(new FixupSpecific());
    SetFlagging(new FlaggingSpecific(*this)); 
    AddGFM(new GhostFluidMethod<VectorType,DIM>(
              new F77GFMBoundary<VectorType,DIM>(f_ibndrfl,f_itrans),
              new F77GFMLevelSet<DataType,DIM>(f_lset)));
  }  
 
  ~SolverSpecific() {
    DeleteGFM(_GFM[0]);
    delete _Flagging;
    delete _Fixup;
    delete _FileOutput;
    delete _LevelTransfer;
  }

  virtual void SetupData() {
    int dim = DIM;
    base::SetupData();
    f_init_commongm(dim,shape,geom,MaxLev,_RefineFactor);
  }


};