changed adjust_reference_heating to .false. if user specifies it usin…

…g with_custom
geodynamics · Nov 30, 2023 · e778002 · e778002
1 parent 4471ac2
commit e778002
Showing 1 changed file with 2 additions and 0 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