Merge pull request #492 from illorenzo7/fix_withcustom_for_heating

Fix with_custom for heating
geodynamics · Apr 3, 2024 · cd0d556 · cd0d556
2 parents a0c3e08 + 74228db
commit cd0d556
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 4 deletions.
diff --git a/src/Physics/PDE_Coefficients.F90 b/src/Physics/PDE_Coefficients.F90
@@ -1066,6 +1066,8 @@ Subroutine Augment_Reference()
             ref%heating(:) = ra_functions(:,6)/(ref%density*ref%temperature)*ra_constants(10)
             temp_functions(:,6) = ra_functions(:,6)
             temp_constants(10)  = ra_constants(10)
+            adjust_reference_heating = .false. ! don't distrust the user's choices. 
+            ! Q(r) = c_10*f_6 (whatever it is) WILL be implemented!
         Endif
 
         If (use_custom_constant(2) .and. use_custom_function(2)) Then

diff --git a/tests/custom_reference/augment/Benchmark_Reports/00000200 b/tests/custom_reference/augment/Benchmark_Reports/00000200
@@ -14,7 +14,7 @@
  ----------------------------------------------------------------------------
    Observable      |    Measured    | Suggested   | % Difference |  Std. Dev.
  ----------------------------------------------------------------------------
-   Kinetic Energy  :     29.150940     58.348000    -50.039522     10.837670
-   Temperature     :      0.321399      0.428120    -24.927858      0.020800
-   Vphi            :     -5.546394    -10.157100    -45.393917      1.983182
-   Drift Frequency :     13.356443      0.182400   7222.611493      0.000000
+   Kinetic Energy  :     29.100559     58.348000    -50.125867     10.770205
+   Temperature     :      0.469344      0.428120      9.629148      0.069833
+   Vphi            :     -5.387204    -10.157100    -46.961199      1.845794
+   Drift Frequency :     13.206574      0.182400   7140.446023      0.000000