FDS Source: make level set flux limiter flags consistent with other scalar transport

Author: ericvmueller

Source/vege.f90

@@ -197,9 +197,9 @@ SUBROUTINE INITIALIZE_LEVEL_SET_FIRESPREAD_2(NM,MODE)
 LSET_TAN2    = .FALSE. ! Flag for ROS proportional to Tan(slope)^2
 
 ! Flux limiters
-! LIMITER_LS=1 MINMOD
+! LIMITER_LS=1 First order upwinding
 ! LIMITER_LS=2 SUPERBEE
-! LIMITER_LS=3 First order upwinding
+! LIMITER_LS=3 MINMOD
 
 LIMITER_LS = I_FLUX_LIMITER
 IF (LIMITER_LS > 3) LIMITER_LS = 1
@@ -1042,7 +1042,7 @@ END SUBROUTINE LEVEL_SET_ADVECT_FLUX
 !>
 !> \param SR_XY If positive, indicates that flow is from left to right
 !> \param Z Scalar quantity
-!> \param LIMITER Flux limiter (1) MINMOD, (2) SUPERBEE, (3) first-order upwinding (monotone)
+!> \param LIMITER Flux limiter (1) first-order upwinding (monotone), (2) SUPERBEE, (3) MINMOD
 !> \details
 !                           face
 !    |     o     |     o     |     o     |     o     |
@@ -1081,11 +1081,11 @@ REAL(EB) FUNCTION SCALAR_FACE_VALUE_LS(SR_XY,Z,LIMITER)
 ENDIF
 
 IF (LIMITER==1) THEN
-   B = MAX(0._EB,MIN(1._EB,R))
+   B = 0._EB
 ELSEIF (LIMITER==2) THEN
    B = MAX(0._EB,MIN(2._EB*R,1._EB),MIN(R,2._EB))
 ELSEIF (LIMITER==3) THEN
-   B = 0._EB
+   B = MAX(0._EB,MIN(1._EB,R))
 ENDIF
 
 SCALAR_FACE_VALUE_LS = ZUP + 0.5_EB * B * ( ZDWN - ZUP )