Changeset 1368 for branches/dev_004_VVL/NEMO/OPA_SRC/DYN/dynspg_flt.F90

Timestamp:

2009-04-01T14:02:17+02:00 (15 years ago)

Author:

rblod

Message:

dev_004_VVL: change time stepping for ssh, see ticket #394

File:

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

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

@@ -427 +427 @@
       END DO
 
-      ! Sea surface elevation time stepping
-      ! -----------------------------------
-      IF( lk_vvl ) THEN   ! after free surface elevation
-         zssha(:,:) = ssha(:,:)
-      ELSE
-         zssha(:,:) = sshb(:,:) + z2dt * ( wn(:,:,1) - emp(:,:) * zraur ) * tmask(:,:,1)
-      ENDIF
-#if defined key_obc
-# if defined key_agrif
-      IF ( Agrif_Root() ) THEN
-# endif
-         zssha(:,:)=zssha(:,:)*obctmsk(:,:)
-         CALL lbc_lnk(zssha,'T',1.)  ! absolutly compulsory !! (jmm)
-# if defined key_agrif
-      ENDIF
-# endif
-#endif
-
-      IF( neuler == 0 .AND. kt == nit000 ) THEN      ! Euler (forward) time stepping, no time filter
-         ! swap of arrays
-         sshb(:,:) = sshn (:,:)
-         sshn(:,:) = zssha(:,:)
-      ELSE                                           ! Leap-frog time stepping and time filter
-         ! time filter and array swap
-         sshb(:,:) = atfp * ( sshb(:,:) + zssha(:,:) ) + atfp1 * sshn(:,:)
-         sshn(:,:) = zssha(:,:)
-      ENDIF
-
       ! write filtered free surface arrays in restart file
       ! --------------------------------------------------