var for obliquity

RevathiJambunathan · May 18, 2022 · ff9997b · ff9997b
1 parent 47e1f2f
commit ff9997b
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 0 deletions.
diff --git a/Source/Particles/PulsarParameters.H b/Source/Particles/PulsarParameters.H
@@ -529,6 +529,10 @@ public:
      ** Upper bound specified must be less than 1, and default is 0.1 (10%)
      */
     static amrex::Real m_ubound_reldiff_sigma0;
+    /** Oblique angle between magnetic field axis and pulsar rotation axis.
+     ** (Read in degrees and converted to radians)
+     */
+    static amrex::Real m_Chi;
 private:
 };
 

diff --git a/Source/Particles/PulsarParameters.cpp b/Source/Particles/PulsarParameters.cpp
@@ -91,6 +91,7 @@ int Pulsar::m_print_injected_celldata;
 int Pulsar::m_print_celldata_starttime;
 amrex::Real Pulsar::m_lbound_ndens_magnetization = 1.e-16;
 amrex::Real Pulsar::m_ubound_reldiff_sigma0 = 0.1;
+amrex::Real Pulsar::m_Chi;
 
 
 Pulsar::Pulsar ()
@@ -250,6 +251,11 @@ Pulsar::ReadParameters () {
     if (m_sigma_tune_method == "relative_difference") {
         pp.query("upperBound_reldiff_sigma0", m_ubound_reldiff_sigma0);
     }
+    // Obliquity of the pulsar
+    pp.get("Chi", m_Chi);
+    // convert degrees to radians
+    m_Chi = m_Chi * MathConst::pi / 180._rt;
+    amrex::Print() << "Oblique angle between B-axis and Omega-axis " << m_Chi << " radians \n";
 }