ww3_gse: Rewrite code for diagonal term. Clean and konsolidate code. …

…After mergin with ww3_difraction the subroutine for the differentation will be removed at the given place.
erdc · Aug 25, 2023 · e3661de · e3661de
1 parent 12d8a8e
commit e3661de
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/model/src/w3profsmd_pdlib.F90 b/model/src/w3profsmd_pdlib.F90
@@ -7016,8 +7016,8 @@ SUBROUTINE BLOCK_SOLVER_DIFFUSION(DTG)
           END DO
           CALL PDLIB_exchange1DREAL(PHI_V)
           DO JSEA =1, NSEAL
-            DIFFTOT = PHI_V(JSEA) + DV2DXY(JSEA) * DIFFVEC(3,JSEA)
-            VA(ISP,JSEA) = MAX(0.,VA(ISP,JSEA) - DT_DIFF * DIFFTOT / PDLIB_SI(JSEA) * DFAC)
+            DIFFTOT = PHI_V(JSEA) + 2 * DV2DXY(JSEA) * DIFFVEC(3,JSEA)
+            VA(ISP,JSEA) = MAX(0.,VA(ISP,JSEA) - DT_DIFF * PHI_V(JSEA) / PDLIB_SI(JSEA) * DFAC) + 2 * DV2DXY(JSEA) * DIFFVEC(3,JSEA) * DFAC
             !IF (ABS(PHI_V(JSEA) .gt. 0.d0)) write(1040+myrank,*) JSEA, DT_DIFF, PHI_V(JSEA), DV2DXY(JSEA) * DIFFVEC(3,JSEA)
           END DO       
         END DO