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 EXPLOSIONIBC_H_
00029 #define EXPLOSIONIBC_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 #include "SolidF_F.H"
00041
00043
00052 class ExplosionIBC:public PhysIBC
00053 {
00054 public:
00056
00058 ExplosionIBC();
00059
00061
00063 ExplosionIBC(const Real& a_gamma,
00064 const Real& a_ms,
00065 const RealVect& a_center,
00066 const Real& a_size,
00067 const RealVect& a_velocity,
00068 const Real& a_artvisc);
00069
00071
00073 ~ExplosionIBC() {};
00074
00076
00078 void setFortranCommon(const Real& a_gamma,
00079 const Real& a_ms,
00080 const RealVect& a_center,
00081 const Real& a_size,
00082 const RealVect& a_velocity,
00083 const Real& a_artvisc);
00084
00086
00091 void setFortranCommonSet();
00092
00094
00099 PhysIBC *new_physIBC();
00100
00102
00104 void fluxBC(FArrayBox& a_F,
00105 const FArrayBox& a_W,
00106 const FArrayBox& a_Wextrap,
00107 const int& a_dir,
00108 const Side::LoHiSide& a_side,
00109 const Real& a_time);
00110
00112
00117 void setBdrySlopes(FArrayBox& a_dW,
00118 const FArrayBox& a_W,
00119 const int& a_dir,
00120 const Real& a_time);
00121
00123
00125 void initialize(LevelData<FArrayBox>& a_U);
00126
00127 protected:
00128
00129 bool m_isFortranCommonSet;
00130 };
00131
00132 #endif