Changes between Version 12 and Version 13 of 2009WP/2009Stream3/VVL

Timestamp:

2009-06-21T12:56:54+02:00 (15 years ago)

Author:

gm

Comment:

--

2009WP/2009Stream3/VVL

@@ -79 +79 @@
             time derivative of e3 should not be not included in the left hand side of the momentum equation [[BR]]
             solve dU/dt = ...  and not d( e3 U)/dt = ... ). [[BR]]
-            Change to be done in both dynnxt AND dynspg (_ts and flt)
+            Change to be done in both dynnxt AND dynspg (_ts ; _flt ; _exp)
 
             '''In dynnxt.F90'''   perform the thickness weighted time leap-frog and  Asselin filter only in ln_dynadv_vec=FALSE and lk_vvl=TRUE., in other word unweighted operation if  n_dynadv_vec=T and NOT lk_vvl. So replace the 2 two [[BR]]
@@ -101 +101 @@
 
     In addition : (1) the masking of ua, va can be moved inside the key_bdy case ; (2) all the comment have been updated. The resulting subroutine can be found in attachment (https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/2009WP/2009Stream3/VVL/dynnxt.F90).
+
+   '''In dynspg_flt.F90''', the time stepping of (ua,va) has to be modified as follows: [[BR]]
+
+{{{
+      IF( lk_vvl ) THEN          ! variable volume 
+
+         ! Evaluate the masked next velocity (effect of the additional force not included)
+         ! -------------------   (thickness weighted velocity, surface pressure gradient already included in dyn_hpg)
+         DO jk = 1, jpkm1
+            DO jj = 2, jpjm1
+               DO ji = fs_2, fs_jpim1   ! vector opt.
+                  ua(ji,jj,jk) = (        ub(ji,jj,jk) * fse3u_b(ji,jj,jk)      &
+                     &           + z2dt * ua(ji,jj,jk) * fse3u_n(ji,jj,jk)  )   &
+                     &         / fse3u_a(ji,jj,jk) * umask(ji,jj,jk)
+                  va(ji,jj,jk) = (        vb(ji,jj,jk) * fse3v_b(ji,jj,jk)      &
+                     &           + z2dt * va(ji,jj,jk) * fse3v_n(ji,jj,jk)  )   &
+                     &         / fse3v_a(ji,jj,jk) * vmask(ji,jj,jk)
+               END DO
+            END DO
+         END DO
+
+      ELSE                       ! fixed volume 
+
+         ! Surface pressure gradient (now)
+         DO jj = 2, jpjm1
+            DO ji = fs_2, fs_jpim1   ! vector opt.
+               spgu(ji,jj) = - grav * ( sshn(ji+1,jj) - sshn(ji,jj) ) / e1u(ji,jj)
+               spgv(ji,jj) = - grav * ( sshn(ji,jj+1) - sshn(ji,jj) ) / e2v(ji,jj)
+            END DO 
+         END DO 
+         !
+         ! add the surface pressure trend to the general trend and
+         ! evaluate the masked next velocity (effect of the additional force not included)
+         DO jk = 1, jpkm1
+            DO jj = 2, jpjm1
+               DO ji = fs_2, fs_jpim1   ! vector opt.
+                  ua(ji,jj,jk) = (  ub(ji,jj,jk) + z2dt * ( ua(ji,jj,jk) + spgu(ji,jj) )  ) * umask(ji,jj,jk)
+                  va(ji,jj,jk) = (  vb(ji,jj,jk) + z2dt * ( va(ji,jj,jk) + spgv(ji,jj) )  ) * vmask(ji,jj,jk)
+               END DO
+            END DO
+         END DO
+         !
+      ENDIF
+}}}
+   becomes
+
+{{{
+      ! Next velocity :   Leap-frog time stepping  (effect of the additional force not included)
+      ! -------------
+      IF( lk_vvl ) THEN          !* variable volume : surface pressure gradient already included in dyn_hpg
+         !
+         IF( ln_dynadv_vec ) THEN       ! vector invariant advection form: Leap-Frog applied on velocity
+            DO jk = 1, jpkm1
+               ua(:,:,jk) = ( ub(:,:,jk) + z2dt * ua(:,:,jk) ) * umask(:,:,jk)
+               va(:,:,jk) = ( vb(:,:,jk) + z2dt * va(:,:,jk) ) * vmask(:,:,jk)
+            END DO
+         ELSE                           ! flux advection form: Leap-Frog applied on thickness weighted velocity
+            DO jk = 1, jpkm1
+               ua(:,:,jk) = (          ub(:,:,jk) * fse3u_b(:,:,jk)     &
+                  &           + z2dt * ua(:,:,jk) * fse3u_n(:,:,jk) )   &
+                  &         / fse3u_a(:,:,jk) * umask(:,:,jk)
+               va(:,:,jk) = (          vb(:,:,jk) * fse3v_b(:,:,jk)     &
+                  &           + z2dt * va(:,:,jk) * fse3v_n(:,:,jk) )   &
+                  &         / fse3v_a(:,:,jk) * vmask(:,:,jk)            
+            END DO
+         ENDIF
+         !
+      ELSE                       !* fixed volume : add the surface pressure gradient trend
+         !
+         DO jj = 2, jpjm1               ! now surface pressure gradient
+            DO ji = fs_2, fs_jpim1   ! vector opt.
+               spgu(ji,jj) = - grav * ( sshn(ji+1,jj) - sshn(ji,jj) ) / e1u(ji,jj)
+               spgv(ji,jj) = - grav * ( sshn(ji,jj+1) - sshn(ji,jj) ) / e2v(ji,jj)
+            END DO 
+         END DO 
+         DO jk = 1, jpkm1               ! Leap-Frog applied on velocity
+            DO jj = 2, jpjm1
+               DO ji = fs_2, fs_jpim1   ! vector opt.
+                  ua(ji,jj,jk) = (  ub(ji,jj,jk) + z2dt * ( ua(ji,jj,jk) + spgu(ji,jj) )  ) * umask(ji,jj,jk)
+                  va(ji,jj,jk) = (  vb(ji,jj,jk) + z2dt * ( va(ji,jj,jk) + spgv(ji,jj) )  ) * vmask(ji,jj,jk)
+               END DO
+            END DO
+         END DO
+         !
+      ENDIF
+}}}
+
+    In addition, all the comment have been updated. The resulting subroutine can be found in attachment (https://forge.ipsl.jussieu.fr/nemo/attachment/wiki/2009WP/2009Stream3/VVL/dynspg_flt.F90).
+
+   '''In dynspg_exp.F90'''   the whole module have to be modified to solve a velocity equation, not a thickness weighted velocity equation. Then in the specific case of ln_dynadv_vec=FALSE and lk_vvl=TRUE, an term associated with the time derivative of the ssh must be added to the equation solved. [[BR]]
 
    '''In dynspg_ts.F90'''   the whole module have to be modified to solve a velocity equation, not a thickness weighted velocity equation. Then in the specific case of ln_dynadv_vec=FALSE and lk_vvl=TRUE, an term associated with the time derivative of the ssh must be added to the equation solved. [[BR]]