Merge pull request #14824 from mcgratta/master

FDS Source: Fix bug in PHI_LS interpolation

Author: mcgratta

Source/vege.f90

@@ -719,10 +719,10 @@ SUBROUTINE FILL_BOUNDARY_VALUES
       IIO_2 = IIO
       JJO_2 = JJO
       SELECT CASE(IOR)
-         CASE( 1) ; IIO_2 = IIO - 1
-         CASE(-1) ; IIO_2 = IIO + 1
-         CASE( 2) ; JJO_2 = JJO - 1
-         CASE(-2) ; JJO_2 = JJO + 1
+         CASE(-1) ; IIO_2 = IIO - 1
+         CASE( 1) ; IIO_2 = IIO + 1
+         CASE(-2) ; JJO_2 = JJO - 1
+         CASE( 2) ; JJO_2 = JJO + 1
       END SELECT
       IF (PREDICTOR) THEN
          PHI_LS_OTHER   = PHI_LS_OTHER   + OMESH(NOM)%PHI_LS(IIO  ,JJO  )