Simple logic fix for time interpolation of boundary nodes at the end …

…of W3XYPFSNIMP. (#1005)
NOAA-EMC · May 5, 2023 · a09335e · a09335e
1 parent 7229a4c
commit a09335e
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/model/src/w3profsmd.F90 b/model/src/w3profsmd.F90
@@ -1266,7 +1266,7 @@ SUBROUTINE W3XYPFSNIMP ( ISP, C, LCALC, RD10, RD20, DT, AC)
       DO IBI=1, NBI
         IP    = MAPSF(ISBPI(IBI),1)
         AC(IP) = ( RD1*BBPI0(ISP,IBI) + RD2*BBPIN(ISP,IBI) )   &
-             *IOBPA(IP)*(1-IOBPD(ITH,IP)) / CG(IK,ISBPI(IBI)) * CLATS(ISBPI(IBI))
+             *IOBPA(IP)*IOBPD(ITH,IP) / CG(IK,ISBPI(IBI)) * CLATS(ISBPI(IBI))
       END DO
     END IF