- Timestamp:
- 2021-07-07T14:42:35+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/DYN/sshwzv.F90
r15058 r15102 359 359 IF(lwp) WRITE(numout,*) 'wAimp : Courant number-based partitioning of now vertical velocity ' 360 360 IF(lwp) WRITE(numout,*) '~~~~~ ' 361 wi(:,:,:) = 0._wp362 361 ENDIF 363 362 ! … … 365 364 zdt = 2._wp * rn_Dt ! 2*rn_Dt and not rDt (for restartability) 366 365 IF( ln_vvl_ztilde .OR. ln_vvl_layer ) THEN 367 DO_3D( nn_hls-1, nn_hls -1, nn_hls-1, nn_hls-1, 1, jpkm1 )366 DO_3D( nn_hls-1, nn_hls, nn_hls-1, nn_hls, 1, jpkm1 ) 368 367 z1_e3t = 1._wp / e3t(ji,jj,jk,Kmm) 369 368 Cu_adv(ji,jj,jk) = zdt * & … … 382 381 END_3D 383 382 ELSE 384 DO_3D( nn_hls-1, nn_hls -1, nn_hls-1, nn_hls-1, 1, jpkm1 )383 DO_3D( nn_hls-1, nn_hls, nn_hls-1, nn_hls, 1, jpkm1 ) 385 384 z1_e3t = 1._wp / e3t(ji,jj,jk,Kmm) 386 385 Cu_adv(ji,jj,jk) = zdt * & … … 395 394 END_3D 396 395 ENDIF 397 IF (nn_hls==1) CALL lbc_lnk( 'sshwzv', Cu_adv, 'T', 1.0_wp )398 !399 396 CALL iom_put("Courant",Cu_adv) 400 397 ! 401 398 IF( MAXVAL( Cu_adv(:,:,:) ) > Cu_min ) THEN ! Quick check if any breaches anywhere 402 DO_3DS( nn_hls , nn_hls, nn_hls, nn_hls, jpkm1, 2, -1 )! or scan Courant criterion and partition ! w where necessary399 DO_3DS( nn_hls-1, nn_hls, nn_hls-1, nn_hls, jpkm1, 2, -1 ) ! or scan Courant criterion and partition ! w where necessary 403 400 ! 404 401 zCu = MAX( Cu_adv(ji,jj,jk) , Cu_adv(ji,jj,jk-1) ) … … 428 425 ELSE 429 426 ! Fully explicit everywhere 430 Cu_adv(:,:,:) = 0._wp 427 Cu_adv(:,:,:) = 0._wp ! Reuse array to output coefficient below and in stp_ctl 431 428 wi (:,:,:) = 0._wp 432 429 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.