Ignore:
Timestamp:
2017-10-03T17:09:40+02:00 (3 years ago)
Author:
clem
Message:

add option Lupkes2015 for ice-atm drag

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_algo_coare.F90

    r7646 r8585  
    5252   !! COARE own values for given constants: 
    5353   REAL(wp), PARAMETER :: & 
    54       &   zi0     = 600.,  &   !: scale height of the atmospheric boundary layer...1 
     54      &   zi0     = 600., &   !: scale height of the atmospheric boundary layer...1 
    5555      &  Beta0    = 1.25, &   !: gustiness parameter 
    5656      &  rctv0    = 0.608     !: constant to obtain virtual temperature... 
     
    6060 
    6161   SUBROUTINE turb_coare( zt, zu, sst, t_zt, ssq, q_zt, U_zu, & 
    62       &                   Cd, Ch, Ce, t_zu, q_zu, U_blk ) 
     62      &                   Cd, Ch, Ce, t_zu, q_zu, U_blk,      & 
     63      &                   Cdn, Chn, Cen                       ) 
     64 
    6365      !!---------------------------------------------------------------------- 
    6466      !!                      ***  ROUTINE  turb_coare  *** 
     
    106108      REAL(wp), INTENT(  out), DIMENSION(jpi,jpj) ::   q_zu     ! spec. humidity adjusted at zu           [kg/kg] 
    107109      REAL(wp), INTENT(  out), DIMENSION(jpi,jpj) ::   U_blk    ! bulk wind at 10m                          [m/s] 
     110      REAL(wp), INTENT(  out), DIMENSION(jpi,jpj) ::   Cdn, Chn, Cen ! neutral transfer coefficients 
    108111      ! 
    109112      INTEGER :: j_itt 
     
    246249      Ce   = ztmp0*q_star/dq_zu 
    247250      ! 
     251      ztmp1 = zu + z0 
     252      Cdn = vkarmn*vkarmn / (log(ztmp1/z0 )*log(ztmp1/z0 )) 
     253      Chn = vkarmn*vkarmn / (log(ztmp1/z0t)*log(ztmp1/z0t)) 
     254      Cen = Chn 
     255      ! 
    248256      CALL wrk_dealloc( jpi,jpj, u_star, t_star, q_star, zeta_u, dt_zu, dq_zu ) 
    249257      CALL wrk_dealloc( jpi,jpj, znu_a, z0, z0t, ztmp0, ztmp1, ztmp2 ) 
Note: See TracChangeset for help on using the changeset viewer.