Fix for couplings = 'none' in TERAPI-SH-S0

PHOTOX · Nov 11, 2024 · 9e4d031 · 9e4d031
1 parent 98e235c
commit 9e4d031
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/surfacehop.F90 b/src/surfacehop.F90
@@ -719,7 +719,7 @@ subroutine surfacehop(x, y, z, vx, vy, vz, vx_old, vy_old, vz_old, dt, eclas)
          pop0 = get_elpop(ist)
 
          ! INTERPOLATION
-         if (inac == 0) then
+         if ((inac == 0) .or. (inac == 2)) then
             fr = real(itp, DP) / real(substep, DP)
             call interpolate(vx, vy, vz, vx_old, vy_old, vz_old, vx_newint, vy_newint, vz_newint, &
                              nacx_newint, nacy_newint, nacz_newint, en_array_newint, &