Changeset 13167
- Timestamp:
- 2020-06-26T18:38:12+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/DYN/sshwzv.F90
r13164 r13167 307 307 ! 308 308 INTEGER :: ji, jj, jk ! dummy loop indices 309 REAL(wp) :: zCu, zcff, z1_e3t 309 REAL(wp) :: zCu, zcff, z1_e3t, zdt ! local scalars 310 310 REAL(wp) , PARAMETER :: Cu_min = 0.15_wp ! local parameters 311 311 REAL(wp) , PARAMETER :: Cu_max = 0.30_wp ! local parameters … … 324 324 ! 325 325 ! Calculate Courant numbers 326 zdt = 2._wp * rn_Dt ! 2*rn_Dt and not rDt (for restartability) 326 327 IF( ln_vvl_ztilde .OR. ln_vvl_layer ) THEN 327 328 DO_3D_00_00( 1, jpkm1 ) 328 329 z1_e3t = 1._wp / e3t(ji,jj,jk,Kmm) 329 Cu_adv(ji,jj,jk) = rDt * &330 Cu_adv(ji,jj,jk) = zdt * & 330 331 & ( ( MAX( ww(ji,jj,jk) , 0._wp ) - MIN( ww(ji,jj,jk+1) , 0._wp ) ) & 331 332 & + ( MAX( e2u(ji ,jj) * e3u(ji ,jj,jk,Kmm) & … … 344 345 DO_3D_00_00( 1, jpkm1 ) 345 346 z1_e3t = 1._wp / e3t(ji,jj,jk,Kmm) 346 Cu_adv(ji,jj,jk) = rDt *&347 Cu_adv(ji,jj,jk) = zdt * & 347 348 & ( ( MAX( ww(ji,jj,jk) , 0._wp ) - MIN( ww(ji,jj,jk+1) , 0._wp ) ) & 348 349 & + ( MAX( e2u(ji ,jj)*e3u(ji ,jj,jk,Kmm)*uu(ji ,jj,jk,Kmm), 0._wp ) - &
Note: See TracChangeset
for help on using the changeset viewer.