[FLASH-BUGS] Flash doesn't run stirturb test problem / Bug in visc_explicit.F90

[Achim Nonnenmacher]

Hello,

the "stirturb" test problem does not run in Flash 2.3. It aborts with the
error message (in flash.log) "Error: only cart, 2-d cyl, and 1-d sph
supported in diffusion", although cartesian coordinates are used.

To fix this bug, remove the "else" in visc_explicit.F90, line 152. Right
now, if it is cartesian, it aborts (see below). In therm_explicit.F90 the
code is correct.


Regards,

Achim Nonnenmacher



----------------------------------------------------------------------
visc_explicit.F90
Line 137 +

! check the geometry -- currently, only cartesian, 2-d cylindrical, and
! 1-d spherical geometry is supported

          if (.NOT. (igeomx == geom_cartesian .AND. &
                     igeomy == geom_cartesian .AND. &
                     igeomz == geom_cartesian) .AND.  &
              .NOT. (igeomx == geom_cylrad .AND. &
                     igeomy == geom_cartesian .AND. &
                     igeomz == geom_cylang .AND. &
                     ndim == 2) .AND. &
              .NOT. (igeomx == geom_sphrad .AND. &
                     igeomy == geom_sphtheta .AND. &
                     igeomz == geom_sphphi .AND. &
                     ndim == 1)) then
            print *, 'supported in explicit diffusion.'
! remove the following else
!          else
            call abort_flash( &
                   "Error: only cart, 2-d cyl, and 1-d sph supported in
diffusion")
           endif