- Timestamp:
- 2017-12-13T15:58:53+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_algo_coare.F90
r7646 r9019 44 44 USE lib_fortran ! to use key_nosignedzero 45 45 46 47 46 IMPLICIT NONE 48 47 PRIVATE … … 50 49 PUBLIC :: TURB_COARE ! called by sbcblk.F90 51 50 52 !! COARE own values for given constants: 53 REAL(wp), PARAMETER :: & 54 & zi0 = 600., & !: scale height of the atmospheric boundary layer...1 55 & Beta0 = 1.25, & !: gustiness parameter 56 & rctv0 = 0.608 !: constant to obtain virtual temperature... 57 58 51 ! !! COARE own values for given constants: 52 REAL(wp), PARAMETER :: zi0 = 600._wp ! scale height of the atmospheric boundary layer... 53 REAL(wp), PARAMETER :: Beta0 = 1.250_wp ! gustiness parameter 54 REAL(wp), PARAMETER :: rctv0 = 0.608_wp ! constant to obtain virtual temperature... 55 56 !!---------------------------------------------------------------------- 59 57 CONTAINS 60 58 61 59 SUBROUTINE turb_coare( zt, zu, sst, t_zt, ssq, q_zt, U_zu, & 62 & Cd, Ch, Ce, t_zu, q_zu, U_blk ) 60 & Cd, Ch, Ce, t_zu, q_zu, U_blk, & 61 & Cdn, Chn, Cen ) 63 62 !!---------------------------------------------------------------------- 64 63 !! *** ROUTINE turb_coare *** … … 106 105 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: q_zu ! spec. humidity adjusted at zu [kg/kg] 107 106 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: U_blk ! bulk wind at 10m [m/s] 107 REAL(wp), INTENT( out), DIMENSION(jpi,jpj) :: Cdn, Chn, Cen ! neutral transfer coefficients 108 108 ! 109 109 INTEGER :: j_itt … … 246 246 Ce = ztmp0*q_star/dq_zu 247 247 ! 248 ztmp1 = zu + z0 249 Cdn = vkarmn*vkarmn / (log(ztmp1/z0 )*log(ztmp1/z0 )) 250 Chn = vkarmn*vkarmn / (log(ztmp1/z0t)*log(ztmp1/z0t)) 251 Cen = Chn 252 ! 248 253 CALL wrk_dealloc( jpi,jpj, u_star, t_star, q_star, zeta_u, dt_zu, dq_zu ) 249 254 CALL wrk_dealloc( jpi,jpj, znu_a, z0, z0t, ztmp0, ztmp1, ztmp2 )
Note: See TracChangeset
for help on using the changeset viewer.