make Gconst a runtime parameter with a default value

quokka-astro · Aug 4, 2023 · c8ac6fd · c8ac6fd
1 parent b990f17
commit c8ac6fd
Show file tree

Hide file tree

Showing 3 changed files with 11 additions and 1 deletion.
diff --git a/src/RadhydroSimulation.hpp b/src/RadhydroSimulation.hpp
@@ -123,7 +123,7 @@ template <typename problem_t> class RadhydroSimulation : public AMRSimulation<pr
 
 	amrex::Long radiationCellUpdates_ = 0; // total number of radiation cell-updates
 
-	amrex::Real Gconst_ = 1.0; // gravitational constant G (code units)
+	amrex::Real Gconst_ = C::Gconst; // gravitational constant G
 
 	// member functions
 	explicit RadhydroSimulation(amrex::Vector<amrex::BCRec> &BCs_cc, amrex::Vector<amrex::BCRec> &BCs_fc) : AMRSimulation<problem_t>(BCs_cc, BCs_fc)
@@ -316,6 +316,12 @@ template <typename problem_t> void RadhydroSimulation<problem_t>::readParmParse(
 		hpp.query("artificial_viscosity_coefficient", artificialViscosityK_);
 	}
 
+	// set gravity runtime parameter
+	{
+		amrex::ParmParse hpp("gravity");
+		hpp.query("Gconst", Gconst_);
+	}
+
 	// set cooling runtime parameters
 	{
 		amrex::ParmParse hpp("cooling");

diff --git a/tests/StarCluster.in b/tests/StarCluster.in
@@ -25,6 +25,8 @@ cfl = 0.2
 max_timesteps = 300
 stop_time = 0.5             # t_ff = 0.55
 
+gravity.Gconst = 1.0    #gravitational constant
+
 do_reflux = 1
 do_subcycle = 0
 

diff --git a/tests/StarCluster_AMR.in b/tests/StarCluster_AMR.in
@@ -25,6 +25,8 @@ cfl = 0.2
 max_timesteps = 10000
 stop_time = 0.5             # t_ff = 0.55
 
+gravity.Gconst = 1.0    #gravitational constant
+
 do_reflux = 1
 do_subcycle = 0