00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028 #ifndef WAVEMHDIBC_H_
00029 #define WAVEMHDIBC_H_
00030
00031 #include <iostream>
00032
00033 #include "LevelData.H"
00034 #include "FArrayBox.H"
00035 #include "Vector.H"
00036 #include "RealVect.H"
00037 #include "AMRIO.H"
00038
00039 #include "PhysIBC.H"
00040
00042
00052 class WaveMHDIBC:public PhysIBC
00053 {
00054 public:
00056
00058 WaveMHDIBC();
00059
00061
00063 WaveMHDIBC(const Real& a_gamma,
00064 const int& a_pdir,
00065 const int& a_kratio,
00066 const int& a_waveNumber,
00067 const Real& a_pertAmplitude);
00068
00070
00072 ~WaveMHDIBC() {};
00073
00075
00077 void setParameters(const Real& a_gamma,
00078 const int& a_pdir,
00079 const int& a_kratio,
00080 const int& a_waveNumber,
00081 const Real& a_pertAmplitude
00082 );
00083
00085
00090 void setParameterSet();
00091
00093
00098 PhysIBC *new_physIBC();
00099
00101
00103 void fluxBC(FArrayBox& a_F,
00104 const FArrayBox& a_W,
00105 const FArrayBox& a_Wextrap,
00106 const int& a_dir,
00107 const Side::LoHiSide& a_side,
00108 const Real& a_time);
00109
00111
00116 void setBdrySlopes(FArrayBox& a_dW,
00117 const FArrayBox& a_W,
00118 const int& a_dir,
00119 const Real& a_time);
00120
00122
00124 void initialize(LevelData<FArrayBox>& a_U);
00125
00126 protected:
00127
00128 bool m_isParameterSet;
00129
00130
00131
00132
00133
00134 };
00135
00136 #endif