Changeset 15533
 Timestamp:
 20211124T13:07:20+01:00 (6 months ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

NEMO/branches/2021/dev_r15388_updated_zdfiwm/src/OCE/ZDF/zdfiwm.F90
r15526 r15533 42 42 LOGICAL :: ln_tsdiff ! account for differential T/S wavedriven mixing (=T) or not (=F) 43 43 44 REAL(wp):: z1_6 = 1._wp / 6._wp45 REAL(wp):: znu = 1.4e6_wp ! molecular kinematic viscosity44 REAL(wp):: r1_6 = 1._wp / 6._wp 45 REAL(wp):: rnu = 1.4e6_wp ! molecular kinematic viscosity 46 46 47 47 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: ebot_iwm ! bottomintensified dissipation above abyssal hills (W/m2) … … 219 219 ! Calculate turbulence intensity parameter Reb 220 220 DO_3D( nn_hls1, nn_hls1, nn_hls1, nn_hls1, 2, jpkm1 ) 221 zReb(ji,jj,jk) = zemx_iwm(ji,jj,jk) / MAX( 1.e20_wp, znu * rn2(ji,jj,jk) )221 zReb(ji,jj,jk) = zemx_iwm(ji,jj,jk) / MAX( 1.e20_wp, rnu * rn2(ji,jj,jk) ) 222 222 END_3D 223 223 ! 224 224 ! Define internal waveinduced diffusivity 225 225 DO_3D( nn_hls1, nn_hls1, nn_hls1, nn_hls1, 2, jpkm1 ) 226 zav_wave(ji,jj,jk) = zReb(ji,jj,jk) * z1_6 * znu ! This corresponds to a constant mixing efficiency of 1/6226 zav_wave(ji,jj,jk) = zReb(ji,jj,jk) * r1_6 * rnu ! This corresponds to a constant mixing efficiency of 1/6 227 227 END_3D 228 228 ! … … 230 230 DO_3D( nn_hls1, nn_hls1, nn_hls1, nn_hls1, 2, jpkm1 ) ! energetic (Reb > 480) and buoyancycontrolled (Reb <10.224) regimes 231 231 IF( zReb(ji,jj,jk) > 480.00_wp ) THEN 232 zav_wave(ji,jj,jk) = 3.6515_wp * znu * SQRT( zReb(ji,jj,jk) )232 zav_wave(ji,jj,jk) = 3.6515_wp * rnu * SQRT( zReb(ji,jj,jk) ) 233 233 ELSEIF( zReb(ji,jj,jk) < 10.224_wp ) THEN 234 zav_wave(ji,jj,jk) = 0.052125_wp * znu * zReb(ji,jj,jk) * SQRT( zReb(ji,jj,jk) )234 zav_wave(ji,jj,jk) = 0.052125_wp * rnu * zReb(ji,jj,jk) * SQRT( zReb(ji,jj,jk) ) 235 235 ENDIF 236 236 END_3D … … 248 248 DO_3D( nn_hls1, nn_hls1, nn_hls1, nn_hls1, 2, jpkm1 ) ! Calculate S/T diffusivity ratio as a function of Reb (else it is set to 1) 249 249 zav_ratio(ji,jj,jk) = ( 0.505_wp + & 250 & 0.495_wp * TANH( 0.92_wp * ( LOG10( MAX( 1.e20, zReb(ji,jj,jk) * 5._wp * z1_6 ) )  0.60_wp ) ) &250 & 0.495_wp * TANH( 0.92_wp * ( LOG10( MAX( 1.e20, zReb(ji,jj,jk) * 5._wp * r1_6 ) )  0.60_wp ) ) & 251 251 & ) * wmask(ji,jj,jk) 252 252 END_3D … … 372 372 ! ensures that avt remains larger than its molecular value (=1.4e7). Therefore, avtb should 373 373 ! be set here to a very small value, and avmb to its (uniform) molecular value (=1.4e6). 374 avmb(:) = znu ! molecular value374 avmb(:) = rnu ! molecular value 375 375 avtb(:) = 1.e10_wp ! very small diffusive minimum (background avt is specified in zdf_iwm) 376 376 avtb_2d(:,:) = 1._wp ! uniform
Note: See TracChangeset
for help on using the changeset viewer.