fix: Cray ftn continuation and AMD flang case-opt array mismatch

m_time_steppers.fpp: remove invalid trailing & on call to
s_compute_moment_of_inertia (Cray ftn rejects blank continuation lines).

m_re_visc.fpp: add #:include 'case.fpp'. Without it, MFC_CASE_OPTIMIZATION
is undefined (False) during Fypp preprocessing, so AMD case-opt builds take
the USING_AMD dimension(3) path while callers that include case.fpp take the
dimension(num_fluids) path — causing 'actual argument has fewer elements'
errors when num_fluids != 3.

Author: sbryngelson

src/simulation/m_re_visc.fpp

@@ -2,6 +2,7 @@
 !! @file m_re_visc.f90
 !! @brief Contains module m_re_visc
 
+#:include 'case.fpp'
 #:include 'macros.fpp'
 
 !> @brief The module contains routines that compute viscosity-related quantities for both Newtonian and non-Newtonian fluids.

src/simulation/m_time_steppers.fpp

@@ -762,9 +762,8 @@ contains
                     ! update the angular velocity with the torque value
                     patch_ib(i)%angular_vel = (patch_ib(i)%angular_vel*patch_ib(i)%moment) + (rk_coef(s, &
                              & 3)*dt*patch_ib(i)%torque/rk_coef(s, 4))  ! add the torque to the angular momentum
-                    call s_compute_moment_of_inertia(i, &
-                                                     & patch_ib(i)%angular_vel) &
-                                                     &  ! update the moment of inertia to be based on the direction of the angular momentum
+                    ! update the moment of inertia based on the direction of the angular momentum
+                    call s_compute_moment_of_inertia(i, patch_ib(i)%angular_vel)
                     patch_ib(i)%angular_vel = patch_ib(i)%angular_vel/patch_ib(i) &
                              & %moment  ! convert back to angular velocity with the new moment of inertia
                 end if