From 0e58b0457c3f3c533e333a6b4aca1def06d43b26 Mon Sep 17 00:00:00 2001
From: Jooheon Yoo <jy884@cornell.edu>
Date: Tue, 12 Sep 2023 06:47:00 -0700
Subject: [PATCH] Add LimitLorentzFactor for ValenciaDivClean

---
 .../ValenciaDivClean/EvolveValenciaDivClean.hpp     | 13 +++++++++++++
 1 file changed, 13 insertions(+)

diff --git a/src/Evolution/Executables/GrMhd/ValenciaDivClean/EvolveValenciaDivClean.hpp b/src/Evolution/Executables/GrMhd/ValenciaDivClean/EvolveValenciaDivClean.hpp
index 18a2e9db4fd42..58b66184c7978 100644
--- a/src/Evolution/Executables/GrMhd/ValenciaDivClean/EvolveValenciaDivClean.hpp
+++ b/src/Evolution/Executables/GrMhd/ValenciaDivClean/EvolveValenciaDivClean.hpp
@@ -85,6 +85,7 @@
 #include "Evolution/Systems/GrMhd/ValenciaDivClean/Tags.hpp"
 #include "Evolution/VariableFixing/Actions.hpp"
 #include "Evolution/VariableFixing/FixToAtmosphere.hpp"
+#include "Evolution/VariableFixing/LimitLorentzFactor.hpp"
 #include "Evolution/VariableFixing/ParameterizedDeleptonization.hpp"
 #include "Evolution/VariableFixing/Tags.hpp"
 #include "IO/Observer/Actions/RegisterEvents.hpp"
@@ -423,6 +424,8 @@ struct EvolutionMetavars {
       parameterized_deleptonization,
       VariableFixing::Actions::FixVariables<
           VariableFixing::FixToAtmosphere<volume_dim>>,
+      VariableFixing::Actions::FixVariables<
+          VariableFixing::LimitLorentzFactor>,
       Actions::UpdateConservatives>>;
 
   using dg_subcell_step_actions = tmpl::flatten<tmpl::list<
@@ -448,6 +451,8 @@ struct EvolutionMetavars {
       parameterized_deleptonization,
       VariableFixing::Actions::FixVariables<
           VariableFixing::FixToAtmosphere<volume_dim>>,
+      VariableFixing::Actions::FixVariables<
+          VariableFixing::LimitLorentzFactor>,
       Actions::UpdateConservatives,
       Actions::Goto<evolution::dg::subcell::Actions::Labels::EndOfSolvers>,
 
@@ -487,6 +492,8 @@ struct EvolutionMetavars {
       parameterized_deleptonization,
       VariableFixing::Actions::FixVariables<
           VariableFixing::FixToAtmosphere<volume_dim>>,
+      VariableFixing::Actions::FixVariables<
+          VariableFixing::LimitLorentzFactor>,
       Actions::UpdateConservatives,
 
       Actions::Label<evolution::dg::subcell::Actions::Labels::EndOfSolvers>>>;
@@ -522,6 +529,8 @@ struct EvolutionMetavars {
               parameterized_deleptonization,
               VariableFixing::Actions::FixVariables<
                   VariableFixing::FixToAtmosphere<volume_dim>>,
+              VariableFixing::Actions::FixVariables<
+                  VariableFixing::LimitLorentzFactor>,
               Actions::UpdateConservatives,
               evolution::dg::subcell::Actions::SetAndCommunicateInitialRdmpData<
                   volume_dim,
@@ -539,12 +548,16 @@ struct EvolutionMetavars {
               parameterized_deleptonization,
               VariableFixing::Actions::FixVariables<
                   VariableFixing::FixToAtmosphere<volume_dim>>,
+              VariableFixing::Actions::FixVariables<
+                  VariableFixing::LimitLorentzFactor>,
               Actions::UpdateConservatives>,
           tmpl::list<evolution::Initialization::Actions::SetVariables<
                          domain::Tags::Coordinates<3, Frame::ElementLogical>>,
                      parameterized_deleptonization,
                      VariableFixing::Actions::FixVariables<
                          VariableFixing::FixToAtmosphere<volume_dim>>,
+                     VariableFixing::Actions::FixVariables<
+                         VariableFixing::LimitLorentzFactor>,
                      Actions::UpdateConservatives>>,
 
       Initialization::Actions::AddComputeTags<