- Timestamp:
- 2015-10-31T08:40:45+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/ZDF/zdfbfr.F90
r5332 r5845 56 56 !! * Substitutions 57 57 # include "vectopt_loop_substitute.h90" 58 # include "domzgr_substitute.h90"59 58 !!---------------------------------------------------------------------- 60 59 !! NEMO/OPA 4.0 , NEMO Consortium (2011) … … 117 116 ikbt = mbkt(ji,jj) 118 117 !! JC: possible WAD implementation should modify line below if layers vanish 119 ztmp = tmask(ji,jj,ikbt) * ( vkarmn / LOG( 0.5_wp * fse3t_n(ji,jj,ikbt) / rn_bfrz0 ))**2._wp118 ztmp = tmask(ji,jj,ikbt) * ( vkarmn / LOG( 0.5_wp * e3t_n(ji,jj,ikbt) / rn_bfrz0 ))**2._wp 120 119 zbfrt(ji,jj) = MAX(bfrcoef2d(ji,jj), ztmp) 121 120 zbfrt(ji,jj) = MIN(zbfrt(ji,jj), rn_bfri2_max) … … 128 127 ikbt = mikt(ji,jj) 129 128 ! JC: possible WAD implementation should modify line below if layers vanish 130 ztmp = (1-tmask(ji,jj,1)) * ( vkarmn / LOG( 0.5_wp * fse3t_n(ji,jj,ikbt) / rn_bfrz0 ))**2._wp129 ztmp = (1-tmask(ji,jj,1)) * ( vkarmn / LOG( 0.5_wp * e3t_n(ji,jj,ikbt) / rn_bfrz0 ))**2._wp 131 130 ztfrt(ji,jj) = MAX(tfrcoef2d(ji,jj), ztmp) 132 131 ztfrt(ji,jj) = MIN(ztfrt(ji,jj), rn_tfri2_max) … … 375 374 DO ji = 1, jpi 376 375 ikbt = mbkt(ji,jj) 377 ztmp = tmask(ji,jj,ikbt) * ( vkarmn / LOG( 0.5_wp * fse3t_n(ji,jj,ikbt) / rn_bfrz0 ))**2._wp376 ztmp = tmask(ji,jj,ikbt) * ( vkarmn / LOG( 0.5_wp * e3t_n(ji,jj,ikbt) / rn_bfrz0 ))**2._wp 378 377 bfrcoef2d(ji,jj) = MAX(bfrcoef2d(ji,jj), ztmp) 379 378 bfrcoef2d(ji,jj) = MIN(bfrcoef2d(ji,jj), rn_bfri2_max) … … 384 383 DO ji = 1, jpi 385 384 ikbt = mikt(ji,jj) 386 ztmp = tmask(ji,jj,ikbt) * ( vkarmn / LOG( 0.5_wp * fse3t_n(ji,jj,ikbt) / rn_tfrz0 ))**2._wp385 ztmp = tmask(ji,jj,ikbt) * ( vkarmn / LOG( 0.5_wp * e3t_n(ji,jj,ikbt) / rn_tfrz0 ))**2._wp 387 386 tfrcoef2d(ji,jj) = MAX(tfrcoef2d(ji,jj), ztmp) 388 387 tfrcoef2d(ji,jj) = MIN(tfrcoef2d(ji,jj), rn_tfri2_max) … … 424 423 ikbu = mbku(ji,jj) ! deepest ocean level at u- and v-points 425 424 ikbv = mbkv(ji,jj) 426 zfru = 0.5 * fse3u(ji,jj,ikbu) / rdt427 zfrv = 0.5 * fse3v(ji,jj,ikbv) / rdt425 zfru = 0.5 * e3u_n(ji,jj,ikbu) / rdt 426 zfrv = 0.5 * e3v_n(ji,jj,ikbv) / rdt 428 427 IF( ABS( bfrcoef2d(ji,jj) ) > zfru ) THEN 429 428 IF( ln_ctl ) THEN … … 446 445 ikbu = miku(ji,jj) ! 1st wet ocean level at u- and v-points 447 446 ikbv = mikv(ji,jj) 448 zfru = 0.5 * fse3u(ji,jj,ikbu) / rdt449 zfrv = 0.5 * fse3v(ji,jj,ikbv) / rdt447 zfru = 0.5 * e3u_n(ji,jj,ikbu) / rdt 448 zfrv = 0.5 * e3v_n(ji,jj,ikbv) / rdt 450 449 IF( ABS( tfrcoef2d(ji,jj) ) > zfru ) THEN 451 450 IF( ln_ctl ) THEN
Note: See TracChangeset
for help on using the changeset viewer.