- Timestamp:
- 2009-04-06T19:48:42+02:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_004_VVL/NEMO/OPA_SRC/DYN/wzvmod.F90
r1384 r1386 59 59 !!---------------------------------------------------------------------- 60 60 INTEGER, INTENT(in) :: kt 61 62 ! Empty routine 63 64 WRITE(*,*) 'wzv : you should not be here : error ?' 65 61 !! 62 INTEGER :: jk ! dummy loop indices 63 REAL(wp) :: z2dt ! temporary scalar 64 !!---------------------------------------------------------------------- 65 66 IF( kt == nit000 ) THEN 67 IF(lwp) WRITE(numout,*) 68 IF(lwp) WRITE(numout,*) 'wzv : vertical velocity from continuity eq.' 69 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 70 71 ! bottom boundary condition: w=0 (set once for all) 72 wn(:,:,jpk) = 0.e0 73 ENDIF 74 75 IF( lk_vvl ) THEN ! Variable volume 76 ! 77 z2dt = 2. * rdt ! time step: leap-frog 78 IF( neuler == 0 .AND. kt == nit000 ) z2dt = rdt ! time step: Euler if restart from rest 79 80 81 ! Vertical velocity computed from bottom 82 ! -------------------------------------- 83 DO jk = jpkm1, 1, -1 84 wn(:,:,jk) = wn(:,:,jk+1) - fse3t(:,:,jk) * hdivn(:,:,jk) & 85 & - ( fse3t_a(:,:,jk) - fse3t_b(:,:,jk) ) / z2dt 86 END DO 87 88 ELSE ! Fixed volume 89 90 ! Vertical velocity computed from bottom 91 ! -------------------------------------- 92 DO jk = jpkm1, 1, -1 93 wn(:,:,jk) = wn(:,:,jk+1) - fse3t_n(:,:,jk) * hdivn(:,:,jk) 94 END DO 95 96 ENDIF 97 98 IF(ln_ctl) CALL prt_ctl(tab3d_1=wn, clinfo1=' w**2 - : ', mask1=wn) 99 ! 66 100 END SUBROUTINE wzv 101 67 102 68 103 SUBROUTINE ssh_wzv( kt )
Note: See TracChangeset
for help on using the changeset viewer.