Changeset 15387

Timestamp:

2021-10-15T19:56:32+02:00 (3 years ago)

Author:

amoulin

Message:

use of charnock coefficient from wave model in sbcblk_algo_ecmwf -ticket #2686

File:

• NEMO/trunk/src/OCE/SBC/sbcblk_algo_ecmwf.F90 (modified) (3 diffs)

NEMO/trunk/src/OCE/SBC/sbcblk_algo_ecmwf.F90

@@ -31 +31 @@
    USE sbc_phy         ! Catalog of functions for physical/meteorological parameters in the marine boundary layer
    USE sbcblk_skin_ecmwf ! cool-skin/warm layer scheme !LB
+   USE sbcwave, ONLY : charn
+   USE sbc_oce, ONLY : ln_charn ! wave module
 
    IMPLICIT NONE
@@ -229 +231 @@
       u_star = 0.035_wp*Ubzu*ztmp1/ztmp0       ! (u* = 0.035*Un10)
 
-      z0     = charn0_ecmwf*u_star*u_star/grav + 0.11_wp*znu_a/u_star
+      IF (ln_charn) THEN ! Charnock value if wave coupling
+         z0 = charn*u_star*u_star/grav + 0.11_wp*znu_a/u_star
+      ELSE
+         z0 = charn0_ecmwf*u_star*u_star/grav + 0.11_wp*znu_a/u_star
+      ENDIF
       z0     = MIN( MAX(ABS(z0), 1.E-9) , 1._wp )                      ! (prevents FPE from stupid values from masked region later on)
 
@@ -296 +302 @@
          ztmp2  = u_star*u_star
          ztmp1  = znu_a/u_star
+         IF (ln_charn) THEN ! Charnock value if wave coupling
+            z0  = MIN( ABS( alpha_M*ztmp1 + charn*ztmp2/grav ) , 0.001_wp)
+         ELSE
+            z0  = MIN( ABS( alpha_M*ztmp1 + charn0_ecmwf*ztmp2/grav ) , 0.001_wp)
+         ENDIF
          z0     = MIN( ABS( alpha_M*ztmp1 + charn0_ecmwf*ztmp2/grav ) , 0.001_wp)
          z0t    = MIN( ABS( alpha_H*ztmp1                           ) , 0.001_wp)   ! eq.3.26, Chap.3, p.34, IFS doc - Cy31r1