Changeset 12489 for NEMO/trunk/tests/BENCH
- Timestamp:
- 2020-02-28T16:55:11+01:00 (4 years ago)
- Location:
- NEMO/trunk/tests/BENCH
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/tests/BENCH/EXPREF/namelist_cfg_orca025_like
r11536 r12489 36 36 &namdom ! time and space domain 37 37 !----------------------------------------------------------------------- 38 rn_ rdt = 900. ! time step for the dynamics (and tracer if nn_acc=0)38 rn_Dt = 900. ! time step for the dynamics (and tracer if nn_acc=0) 39 39 rn_atfp = 0.05 ! asselin time filter parameter 40 40 / … … 189 189 ln_dynspg_ts = .true. ! split-explicit free surface 190 190 ln_bt_auto = .false. ! Number of sub-step defined from: 191 nn_ baro = 30 ! =F : the number of sub-step in rn_rdt seconds191 nn_e = 30 ! =F : the number of sub-step in rn_Dt seconds 192 192 / 193 193 !----------------------------------------------------------------------- -
NEMO/trunk/tests/BENCH/EXPREF/namelist_cfg_orca12_like
r11536 r12489 36 36 &namdom ! time and space domain 37 37 !----------------------------------------------------------------------- 38 rn_ rdt = 300. ! time step for the dynamics (and tracer if nn_acc=0)38 rn_Dt = 300. ! time step for the dynamics (and tracer if nn_acc=0) 39 39 rn_atfp = 0.05 ! asselin time filter parameter 40 40 / … … 188 188 ln_dynspg_ts = .true. ! split-explicit free surface 189 189 ln_bt_auto = .false. ! Number of sub-step defined from: 190 nn_ baro = 30 ! =F : the number of sub-step in rn_rdt seconds190 nn_e = 30 ! =F : the number of sub-step in rn_Dt seconds 191 191 / 192 192 !----------------------------------------------------------------------- -
NEMO/trunk/tests/BENCH/EXPREF/namelist_cfg_orca1_like
r11536 r12489 36 36 &namdom ! time and space domain 37 37 !----------------------------------------------------------------------- 38 rn_ rdt = 3600. ! time step for the dynamics (and tracer if nn_acc=0)38 rn_Dt = 3600. ! time step for the dynamics (and tracer if nn_acc=0) 39 39 rn_atfp = 0.05 ! asselin time filter parameter 40 40 / … … 188 188 ln_dynspg_ts = .true. ! split-explicit free surface 189 189 ln_bt_auto = .false. ! Number of sub-step defined from: 190 nn_ baro = 30 ! =F : the number of sub-step in rn_rdt seconds190 nn_e = 30 ! =F : the number of sub-step in rn_Dt seconds 191 191 / 192 192 !----------------------------------------------------------------------- -
NEMO/trunk/tests/BENCH/MY_SRC/zdfiwm.F90
r12487 r12489 87 87 !! This is divided into three components: 88 88 !! 1. Bottom-intensified low-mode dissipation at critical slopes 89 !! zemx_iwm(z) = ( ecri_iwm / r au0 ) * EXP( -(H-z)/hcri_iwm )89 !! zemx_iwm(z) = ( ecri_iwm / rho0 ) * EXP( -(H-z)/hcri_iwm ) 90 90 !! / ( 1. - EXP( - H/hcri_iwm ) ) * hcri_iwm 91 91 !! where hcri_iwm is the characteristic length scale of the bottom 92 92 !! intensification, ecri_iwm a map of available power, and H the ocean depth. 93 93 !! 2. Pycnocline-intensified low-mode dissipation 94 !! zemx_iwm(z) = ( epyc_iwm / r au0 ) * ( sqrt(rn2(z))^nn_zpyc )94 !! zemx_iwm(z) = ( epyc_iwm / rho0 ) * ( sqrt(rn2(z))^nn_zpyc ) 95 95 !! / SUM( sqrt(rn2(z))^nn_zpyc * e3w(z) ) 96 96 !! where epyc_iwm is a map of available power, and nn_zpyc … … 98 98 !! energy dissipation. 99 99 !! 3. WKB-height dependent high mode dissipation 100 !! zemx_iwm(z) = ( ebot_iwm / r au0 ) * rn2(z) * EXP(-z_wkb(z)/hbot_iwm)100 !! zemx_iwm(z) = ( ebot_iwm / rho0 ) * rn2(z) * EXP(-z_wkb(z)/hbot_iwm) 101 101 !! / SUM( rn2(z) * EXP(-z_wkb(z)/hbot_iwm) * e3w(z) ) 102 102 !! where hbot_iwm is the characteristic length scale of the WKB bottom … … 151 151 DO_2D_11_11 152 152 zhdep(ji,jj) = gdepw_0(ji,jj,mbkt(ji,jj)+1) ! depth of the ocean 153 zfact(ji,jj) = r au0 * ( 1._wp - EXP( -zhdep(ji,jj) / hcri_iwm(ji,jj) ) )153 zfact(ji,jj) = rho0 * ( 1._wp - EXP( -zhdep(ji,jj) / hcri_iwm(ji,jj) ) ) 154 154 IF( zfact(ji,jj) /= 0._wp ) zfact(ji,jj) = ecri_iwm(ji,jj) / zfact(ji,jj) 155 155 END_2D … … 181 181 ! 182 182 DO_2D_11_11 183 IF( zfact(ji,jj) /= 0 ) zfact(ji,jj) = epyc_iwm(ji,jj) / ( r au0 * zfact(ji,jj) )183 IF( zfact(ji,jj) /= 0 ) zfact(ji,jj) = epyc_iwm(ji,jj) / ( rho0 * zfact(ji,jj) ) 184 184 END_2D 185 185 ! … … 196 196 ! 197 197 DO_2D_11_11 198 IF( zfact(ji,jj) /= 0 ) zfact(ji,jj) = epyc_iwm(ji,jj) / ( r au0 * zfact(ji,jj) )198 IF( zfact(ji,jj) /= 0 ) zfact(ji,jj) = epyc_iwm(ji,jj) / ( rho0 * zfact(ji,jj) ) 199 199 END_2D 200 200 ! … … 243 243 ! 244 244 DO_2D_11_11 245 IF( zfact(ji,jj) /= 0 ) zfact(ji,jj) = ebot_iwm(ji,jj) / ( r au0 * zfact(ji,jj) )245 IF( zfact(ji,jj) /= 0 ) zfact(ji,jj) = ebot_iwm(ji,jj) / ( rho0 * zfact(ji,jj) ) 246 246 END_2D 247 247 ! … … 255 255 ! Calculate molecular kinematic viscosity 256 256 znu_t(:,:,:) = 1.e-4_wp * ( 17.91_wp - 0.53810_wp * ts(:,:,:,jp_tem,Kmm) + 0.00694_wp * ts(:,:,:,jp_tem,Kmm) * ts(:,:,:,jp_tem,Kmm) & 257 & + 0.02305_wp * ts(:,:,:,jp_sal,Kmm) ) * tmask(:,:,:) * r1_r au0257 & + 0.02305_wp * ts(:,:,:,jp_sal,Kmm) ) * tmask(:,:,:) * r1_rho0 258 258 DO jk = 2, jpkm1 259 259 znu_w(:,:,jk) = 0.5_wp * ( znu_t(:,:,jk-1) + znu_t(:,:,jk) ) * wmask(:,:,jk) … … 293 293 END_3D 294 294 CALL mpp_sum( 'zdfiwm', zztmp ) 295 zztmp = r au0 * zztmp ! Global integral of rauo * Kz * N^2 = power contributing to mixing295 zztmp = rho0 * zztmp ! Global integral of rauo * Kz * N^2 = power contributing to mixing 296 296 ! 297 297 IF(lwp) THEN … … 337 337 !* output useful diagnostics: Kz*N^2 , 338 338 !!gm Kz*N2 should take into account the ratio avs/avt if it is used.... (see diaar5) 339 ! vertical integral of r au0 * Kz * N^2 , energy density (zemx_iwm)339 ! vertical integral of rho0 * Kz * N^2 , energy density (zemx_iwm) 340 340 IF( iom_use("bflx_iwm") .OR. iom_use("pcmap_iwm") ) THEN 341 341 ALLOCATE( z2d(jpi,jpj) , z3d(jpi,jpj,jpk) ) … … 345 345 z2d(:,:) = z2d(:,:) + e3w(:,:,jk,Kmm) * z3d(:,:,jk) * wmask(:,:,jk) 346 346 END DO 347 z2d(:,:) = r au0 * z2d(:,:)347 z2d(:,:) = rho0 * z2d(:,:) 348 348 CALL iom_put( "bflx_iwm", z3d ) 349 349 CALL iom_put( "pcmap_iwm", z2d )
Note: See TracChangeset
for help on using the changeset viewer.