New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 15548 for NEMO/branches/2021/ticket2632_r14588_theta_sbcblk/src/OCE/SBC/sbcblk_algo_ecmwf.F90 – NEMO

Ignore:
Timestamp:
2021-11-28T18:59:49+01:00 (3 years ago)
Author:
gsamson
Message:

update branch to the head of the trunk (r15547); ticket #2632

Location:
NEMO/branches/2021/ticket2632_r14588_theta_sbcblk
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • NEMO/branches/2021/ticket2632_r14588_theta_sbcblk

    • Property svn:externals

      •  

        old new  
        99 
        1010# SETTE 
        11 ^/utils/CI/sette@14244        sette 
         11^/utils/CI/sette@HEAD        sette 
         12 

  • NEMO/branches/2021/ticket2632_r14588_theta_sbcblk/src/OCE/SBC/sbcblk_algo_ecmwf.F90

    r14592 r15548  
    3131   USE sbc_phy         ! Catalog of functions for physical/meteorological parameters in the marine boundary layer 
    3232   USE sbcblk_skin_ecmwf ! cool-skin/warm layer scheme !LB 
     33   USE sbcwave, ONLY : charn 
     34   USE sbc_oce, ONLY : ln_charn ! wave module 
    3335 
    3436   IMPLICIT NONE 
    … …  
    230232      u_star = 0.035_wp*Ubzu*ztmp1/ztmp0       ! (u* = 0.035*Un10) 
    231233 
    232       z0     = charn0_ecmwf*u_star*u_star/grav + 0.11_wp*znu_a/u_star 
     234      IF (ln_charn) THEN ! Charnock value if wave coupling 
     235         z0 = charn*u_star*u_star/grav + 0.11_wp*znu_a/u_star 
     236      ELSE 
     237         z0 = charn0_ecmwf*u_star*u_star/grav + 0.11_wp*znu_a/u_star 
     238      ENDIF 
    233239      z0     = MIN( MAX(ABS(z0), 1.E-9) , 1._wp )                      ! (prevents FPE from stupid values from masked region later on) 
    234240 
    … …  
    297303         ztmp2  = u_star*u_star 
    298304         ztmp1  = znu_a/u_star 
     305         IF (ln_charn) THEN ! Charnock value if wave coupling 
     306            z0  = MIN( ABS( alpha_M*ztmp1 + charn*ztmp2/grav ) , 0.001_wp) 
     307         ELSE 
     308            z0  = MIN( ABS( alpha_M*ztmp1 + charn0_ecmwf*ztmp2/grav ) , 0.001_wp) 
     309         ENDIF 
    299310         z0     = MIN( ABS( alpha_M*ztmp1 + charn0_ecmwf*ztmp2/grav ) , 0.001_wp) 
    300311         z0t    = MIN( ABS( alpha_H*ztmp1                           ) , 0.001_wp)   ! eq.3.26, Chap.3, p.34, IFS doc - Cy31r1 
Note: See TracChangeset for help on using the changeset viewer.