- Timestamp:
- 2009-11-06T15:01:56+01:00 (14 years ago)
- Location:
- trunk/NEMO/OPA_SRC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/DYN/dynbfr.F90
r1708 r1719 75 75 ikbvm1 = MAX( ikbv-1, 1 ) 76 76 ! 77 ! Apply stability criteria 77 ! Apply stability criteria on absolute value : Min abs(bfr) => Max (bfr) 78 78 zbfru = MAX( bfrua(ji,jj), fse3u(ji,jj,ikbu)*zinv ) 79 79 zbfrv = MAX( bfrva(ji,jj), fse3v(ji,jj,ikbv)*zinv ) -
trunk/NEMO/OPA_SRC/ZDF/zdftke.F90
r1705 r1719 213 213 ! ! Bottom boundary condition on tke 214 214 ! !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 215 ! en(bot) = 0.5*sqrt(u_botfr^2+v_botfr^2) (min value rn_emin) 216 !CDIR NOVERRCHK 217 DO jj = 2, jpjm1 218 !CDIR NOVERRCHK 219 DO ji = fs_2, fs_jpim1 ! vector opt. 220 ikbu = MIN( mbathy(ji+1,jj), mbathy(ji,jj) ) 221 ikbv = MIN( mbathy(ji,jj+1), mbathy(ji,jj) ) 222 ikbum1 = MAX( ikbu-1, 1 ) 223 ikbvm1 = MAX( ikbv-1, 1 ) 224 ikbu = MIN( mbathy(ji,jj), mbathy(ji-1,jj) ) 225 ikbv = MIN( mbathy(ji,jj), mbathy(ji,jj-1) ) 226 ikbumm1 = MAX( ikbu-1, 1 ) 227 ikbvmm1 = MAX( ikbv-1, 1 ) 228 ikbt = MAX( mbathy(ji,jj), 1 ) 229 ztx2 = bfrua(ji-1,jj) * fse3u(ji-1,jj ,ikbumm1) + & 230 bfrua(ji ,jj) * fse3u(ji ,jj ,ikbum1 ) 231 zty2 = bfrva(ji,jj ) * fse3v(ji ,jj ,ikbvm1) + & 232 bfrva(ji,jj-1) * fse3v(ji ,jj-1,ikbvmm1 ) 233 zebot = 0.25_wp * SQRT( ztx2 * ztx2 + zty2 * zty2 ) 234 en (ji,jj,ikbt) = MAX( zebot, rn_emin ) * tmask(ji,jj,1) 235 END DO 236 END DO 215 ! 216 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 217 ! Tests to date have found the bottom boundary condition on tke to have very little effect. 218 ! The condition is coded here for completion but commented out until there is proof that the 219 ! computational cost is justified 220 !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 221 ! 222 ! en(bot) = (rn_ebb0/rau0)*0.5*sqrt(u_botfr^2+v_botfr^2) (min value rn_emin) 223 !CDIR NOVERRCHK 224 !! DO jj = 2, jpjm1 225 !CDIR NOVERRCHK 226 !! DO ji = fs_2, fs_jpim1 ! vector opt. 227 !! ikbu = MIN( mbathy(ji+1,jj), mbathy(ji,jj) ) 228 !! ikbv = MIN( mbathy(ji,jj+1), mbathy(ji,jj) ) 229 !! ikbum1 = MAX( ikbu-1, 1 ) 230 !! ikbvm1 = MAX( ikbv-1, 1 ) 231 !! ikbu = MIN( mbathy(ji,jj), mbathy(ji-1,jj) ) 232 !! ikbv = MIN( mbathy(ji,jj), mbathy(ji,jj-1) ) 233 !! ikbumm1 = MAX( ikbu-1, 1 ) 234 !! ikbvmm1 = MAX( ikbv-1, 1 ) 235 !! ikbt = MAX( mbathy(ji,jj), 1 ) 236 !! ztx2 = bfrua(ji-1,jj) * ub(ji-1,jj,ikbumm1) + & 237 !! bfrua(ji ,jj) * ub(ji ,jj ,ikbum1 ) 238 !! zty2 = bfrva(ji,jj ) * vb(ji ,jj ,ikbvm1) + & 239 !! bfrva(ji,jj-1) * vb(ji ,jj-1,ikbvmm1 ) 240 !! zebot = 0.001875_wp * SQRT( ztx2 * ztx2 + zty2 * zty2 ) ! where 0.001875 = (rn_ebb0/rau0) * 0.5 = 3.75*0.5/1000. 241 !! en (ji,jj,ikbt) = MAX( zebot, rn_emin ) * tmask(ji,jj,1) 242 !! END DO 243 !! END DO 237 244 ! 238 245 !
Note: See TracChangeset
for help on using the changeset viewer.