add comments and correct initialization for XMU

umr-lops · Apr 2, 2024 · 4c81f1c · 4c81f1c
1 parent 5635032
commit 4c81f1c
Showing 1 changed file with 55 additions and 4 deletions.
diff --git a/model/src/w3iogomd.F90 b/model/src/w3iogomd.F90
@@ -4848,13 +4848,35 @@ SUBROUTINE SECONDHH(NKHF,FAC0,FAC1,FAC2,FAC3)
 
 
      CONTAINS
+
+!-----------------------------------------------------------------------
+
      REAL(KIND=4) FUNCTION VMIN_D(XI,XJ,XK,XIJ,XIK,XJK,XOI,XOJ,XOK)
 
-!***  *VMIN_D*  DETERMINES THE NONLINEAR TRANSFER COEFFICIENT FOR THREE
-!               WAVE INTERACTIONS OF DEEP WATER WAVES.
-
 !     PETER JANSSEN
 
+!     PURPOSE.
+!     --------
+
+!              GIVES NONLINEAR TRANSFER COEFFICIENT FOR THREE
+!              WAVE INTERACTIONS OF DEEP-WATER WAVES IN THE
+!              IDEAL CASE OF NO CURRENT. (CF.ZAKHAROV)
+
+!     INTERFACE.
+!     ----------
+!              *VMIN_D(XI,XJ,XK)*
+!                      *XI*  - WAVE NUMBER
+!                      *XJ*  - WAVE NUMBER
+!                      *XK*  - WAVE NUMBER
+!     METHOD.
+!     -------
+!              NONE
+
+!     EXTERNALS.
+!     ----------
+!              NONE.
+
+
 !***  1. DETERMINE NONLINEAR TRANSFER.
 !     --------------------------------
       IMPLICIT NONE
@@ -4875,8 +4897,36 @@ REAL(KIND=4) FUNCTION VMIN_D(XI,XJ,XK,XIJ,XIK,XJK,XOI,XOJ,XOK)
 
       END FUNCTION VMIN_D      
 
+!-----------------------------------------------------------------------
+
       REAL(KIND=4) FUNCTION VPLUS_D(XI,XJ,XK,XIJ,XIK,XJK,XOI,XOJ,XOK)
 
+!***  *VPLUS_D*  DETERMINES THE NONLINEAR TRANSFER COEFFICIENT FOR THREE
+!                WAVE INTERACTIONS OF DEEP-WATER WAVES.
+
+!     PETER JANSSEN
+
+!     PURPOSE.
+!     --------
+
+!              GIVES NONLINEAR TRANSFER COEFFICIENT FOR THREE
+!              WAVE INTERACTIONS OF GRAVITY-CAPILLARY WAVES IN THE
+!              IDEAL CASE OF NO CURRENT. (CF.ZAKHAROV)
+
+!     INTERFACE.
+!     ----------
+!              *VPLUS_D(XI,XJ,XK)*
+!                        *XI*  - WAVE NUMBER
+!                        *XJ*  - WAVE NUMBER
+!                        *XK*  - WAVE NUMBER
+!     METHOD.
+!     -------
+!              NONE
+
+!     EXTERNALS.
+!     ----------
+!              NONE.
+
 
 
 !***  1. DETERMINE NONLINEAR TRANSFER.
@@ -5003,12 +5053,13 @@ SUBROUTINE SKEWNESS(A)
 !     -----------------------------------------------
 
     MSTART = 1
-    XMU(:,:,:) = 0.0
+
 
 #ifdef W3_OMPG
         !$OMP PARALLEL DO PRIVATE(JSEA)
 #endif
     DO JSEA=1, NSEAL
+      XMU(:,:,:) = 0.0
       DO K=1,NTH
         DO M=1,NK
           CONX = TPIINV / SIG(M) * CG(M,JSEA)