Merge branch 'photoionization' of https://github.com/quokka-astro/quokka into photoionization

Author: James471

src/problems/StromgrenSphere/testStromgrenSphere.cpp

@@ -120,10 +120,11 @@ void RadSystem<StromgrenSphere>::SetRadEnergySource(array_t &radEnergy, const am
 // Numerically integrate dR/dt = (Q - 4*pi*R^3*alpha_B*n_HI0^2/3) / (Q/c + 4*pi*R^2*n_HI0)
 // Integrates forward by `dt_target` starting from `R0` using RK4 with step-doubling.
 // Aborts the run if convergence is not reached within allowed iterations.
-static amrex::Real integrate_radius(amrex::Real dt_target, amrex::Real Q, amrex::Real alpha_B, amrex::Real n_HI0,
-								   amrex::Real c_light, amrex::Real R0, amrex::Real r_s_est)
+static amrex::Real integrate_radius(amrex::Real dt_target, amrex::Real Q, amrex::Real alpha_B, amrex::Real n_HI0, amrex::Real c_light, amrex::Real R0,
+				    amrex::Real r_s_est)
 {
-	if (dt_target <= 0.0_rt) return R0;
+	if (dt_target <= 0.0_rt)
+		return R0;
 
 	auto rhs = [&](amrex::Real R) -> amrex::Real {
 		const amrex::Real num = Q - (4.0_rt * M_PI * R * R * R * alpha_B * n_HI0 * n_HI0) / 3.0_rt;
@@ -146,7 +147,8 @@ static amrex::Real integrate_radius(amrex::Real dt_target, amrex::Real Q, amrex:
 			const amrex::Real k3 = rhs(R + 0.5_rt * dt * k2);
 			const amrex::Real k4 = rhs(R + dt * k3);
 			R += (dt / 6.0_rt) * (k1 + 2.0_rt * k2 + 2.0_rt * k3 + k4);
-			if (R < 0.0_rt) R = 0.0_rt;
+			if (R < 0.0_rt)
+				R = 0.0_rt;
 		}
 
 		if (iter > 0 && std::abs(R - R_prev) < tol) {
@@ -160,7 +162,6 @@ static amrex::Real integrate_radius(amrex::Real dt_target, amrex::Real Q, amrex:
 	return R_prev; // unreachable
 }
 
-
 template <> struct SimulationData<StromgrenSphere> {
 	amrex::Real small_temp;
 	amrex::Real small_dens;
@@ -432,8 +433,7 @@ template <> void QuokkaSimulation<StromgrenSphere>::computeAfterTimestep()
 	if (amrex::ParallelDescriptor::IOProcessor()) {
 		const amrex::Real n_HI0 = userData_.primary_species_2;
 		const amrex::Real alpha_B = 2.6e-13;
-		const amrex::Real r_s =
-			std::pow((3.0_rt * userData_.Q) / (4.0_rt * M_PI * alpha_B * n_HI0 * n_HI0), 1.0_rt / 3.0_rt);
+		const amrex::Real r_s = std::pow((3.0_rt * userData_.Q) / (4.0_rt * M_PI * alpha_B * n_HI0 * n_HI0), 1.0_rt / 3.0_rt);
 
 		amrex::Real dt = tNew_[lev] - userData_.r_analytical_last_t;
 		if (dt < 0.0_rt) {
@@ -457,7 +457,7 @@ template <> void QuokkaSimulation<StromgrenSphere>::computeAfterTimestep()
 	userData_.r_analytical_vec_.push_back(r_analytical);
 
 	if (amrex::ParallelDescriptor::IOProcessor()) {
-			userData_.output_file_ << tNew_[lev] << ',' << r16 << ',' << r50 << ',' << r84 << ',' << r_analytical << '\n';
+		userData_.output_file_ << tNew_[lev] << ',' << r16 << ',' << r50 << ',' << r84 << ',' << r_analytical << '\n';
 	}
 }