Don't test using ZHI_COMP unless you know you're doing inhomo_reion

AMReX-Astro · Mar 8, 2018 · 4bac257 · 4bac257
1 parent b98ac0e
commit 4bac257
Showing 1 changed file with 3 additions and 4 deletions.
diff --git a/Source/Src_3d/compute_temp_3d.f90 b/Source/Src_3d/compute_temp_3d.f90
@@ -69,10 +69,9 @@ subroutine fort_compute_temp(lo,hi, &
 
                    eint = state(i,j,k,UEINT) * rhoInv
 
-                   if ((inhomogeneous_on) .and. (z .gt. diag_eos(i,j,k,ZHI_COMP))) then
-                       JH = 0
-                   else
-                       JH = 1
+                   JH = 1
+                   if (inhomogeneous_on) then
+                       if (z .gt. diag_eos(i,j,k,ZHI_COMP)) JH = 0
                    endif
 
                    call nyx_eos_T_given_Re(JH, JHe, diag_eos(i,j,k,TEMP_COMP), diag_eos(i,j,k,NE_COMP), &