Changeset 1386

Timestamp:

2009-04-06T19:48:42+02:00 (15 years ago)

Author:

rblod

Message:

Add wn computation for vvl case

File:

• branches/dev_004_VVL/NEMO/OPA_SRC/DYN/wzvmod.F90 (modified) (1 diff)

branches/dev_004_VVL/NEMO/OPA_SRC/DYN/wzvmod.F90

@@ -59 +59 @@
       !!----------------------------------------------------------------------
       INTEGER, INTENT(in) :: kt
-      
-      ! Empty routine
-
-      WRITE(*,*) 'wzv : you should not be here : error ?'  
-
+      !!
+      INTEGER  ::   jk     ! dummy loop indices
+      REAL(wp) ::   z2dt   ! temporary scalar
+      !!----------------------------------------------------------------------
+
+      IF( kt == nit000 ) THEN
+         IF(lwp) WRITE(numout,*)
+         IF(lwp) WRITE(numout,*) 'wzv     : vertical velocity from continuity eq.'
+         IF(lwp) WRITE(numout,*) '~~~~~~~ ' 
+
+         ! bottom boundary condition: w=0 (set once for all)
+         wn(:,:,jpk) = 0.e0
+      ENDIF
+
+      IF( lk_vvl ) THEN                ! Variable volume
+         !
+         z2dt = 2. * rdt                                       ! time step: leap-frog
+         IF( neuler == 0 .AND. kt == nit000 ) z2dt = rdt       ! time step: Euler if restart from rest
+
+
+         ! Vertical velocity computed from bottom
+         ! --------------------------------------
+         DO jk = jpkm1, 1, -1
+            wn(:,:,jk) = wn(:,:,jk+1) - fse3t(:,:,jk) * hdivn(:,:,jk) &
+              &        - ( fse3t_a(:,:,jk) - fse3t_b(:,:,jk) ) / z2dt
+         END DO
+
+      ELSE                             ! Fixed volume
+
+         ! Vertical velocity computed from bottom
+         ! --------------------------------------
+         DO jk = jpkm1, 1, -1
+            wn(:,:,jk) = wn(:,:,jk+1) - fse3t_n(:,:,jk) * hdivn(:,:,jk)
+         END DO
+
+      ENDIF
+
+      IF(ln_ctl)   CALL prt_ctl(tab3d_1=wn, clinfo1=' w**2 -   : ', mask1=wn)
+      !
    END SUBROUTINE wzv
+
 
    SUBROUTINE ssh_wzv( kt )