Changeset 13759 for NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/ZDF

Timestamp:

2020-11-09T18:04:19+01:00 (4 years ago)

Author:

emanuelaclementi

Message:

minor changes to include writing statements and update n. of max coupling fields with oasis - tickets #2155 #2339

File:

• NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/ZDF/zdftke.F90 (modified) (3 diffs)

NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/ZDF/zdftke.F90

@@ -408 +408 @@
 
       IF ( cpl_phioc .and. ln_phioc )  THEN
-         SELECT CASE (nn_bc_surf)  ! Dirichlet Boundary Condition using surface TKE flux from waves
-
-         CASE ( 0 )
+         SELECT CASE (nn_bc_surf)  ! Boundary Condition using surface TKE flux from waves
+
+         CASE ( 0 ) ! Dirichlet BC
 
             DO_2D( 0, 0, 0, 0 )    ! en(1)   = rn_ebb taum / rho0  (min value rn_emin0)
@@ -420 +420 @@
             END_2D
 
-         CASE ( 1 )
+         CASE ( 1 ) ! Neumann BC
             DO_2D( 0, 0, 0, 0 )
                IF ( phioc(ji,jj) < 0 )  phioc(ji,jj) = 0._wp
@@ -769 +769 @@
          WRITE(numout,*) '      Langmuir cells parametrization              ln_lc     = ', ln_lc
          WRITE(numout,*) '         coef to compute vertical velocity of LC     rn_lc  = ', rn_lc
+         IF ( cpl_phioc .and. ln_phioc )  THEN
+            SELECT CASE( nn_bc_surf)             ! Type of scaling under sea-ice
+            CASE( 0 )   ;   WRITE(numout,*) '  nn_bc_surf=0  ==>>>   DIRICHLET Surface Boundary Condition using surface TKE flux from waves'
+            CASE( 1 )   ;   WRITE(numout,*) '  nn_bc_surf=1  ==>>>   NEUMANN Surface Boundary Condition using surface TKE flux from waves'
+            END SELECT
+         ENDIF 
          WRITE(numout,*) '      test param. to add tke induced by wind      nn_etau   = ', nn_etau
          WRITE(numout,*) '          type of tke penetration profile            nn_htau   = ', nn_htau