Changeset 367 for trunk/NEMO/OPA_SRC/DYN/dynnxt.F90

Timestamp:

2005-12-28T10:25:10+01:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_035 : CT : OBCs adapted to the new surface pressure gradient algorithms

File:

• trunk/NEMO/OPA_SRC/DYN/dynnxt.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/DYN/dynnxt.F90

@@ -12 +12 @@
    USE dom_oce         ! ocean space and time domain
    USE in_out_manager  ! I/O manager
+   USE obc_oce         ! ocean open boundary conditions
    USE obcdyn          ! open boundary condition for momentum (obc_dyn routine)
+   USE obcdyn_bt       ! 2D open boundary condition for momentum (obc_dyn_bt routine)
+   USE obcvol          ! ocean open boundary condition (obc_vol routines)
+   USE dynspg_oce      ! type of surface pressure gradient
    USE lbclnk          ! lateral boundary condition (or mpp link)
    USE prtctl          ! Print control
@@ -105 +109 @@
       ! Update (ua,va) along open boundaries (only in the rigid-lid case)
       CALL obc_dyn( kt )
+
+      IF ( lk_dynspg_exp .OR. lk_dynspg_ts ) THEN
+         !Flather boundary condition :
+         !        - Update sea surface height on each open boundary
+         !                 sshn (= after ssh) for explicit case
+         !                 sshn_b (= after ssha_b) for time-splitting case
+         !        - Correct the barotropic velocities
+         CALL obc_dyn_bt( kt )
+
+         !Boundary conditions on sshn ( after ssh)
+         CALL lbc_lnk( sshn, 'T', 1. )
+
+         IF(ln_ctl) THEN         ! print sum trends (used for debugging)
+            CALL prt_ctl(tab2d_1=sshn, clinfo1=' ssh      : ', mask1=tmask)
+         ENDIF
+
+         IF ( ln_vol_cst ) CALL obc_vol( kt )
+
+      ENDIF
+
       !                                                ! ===============
       DO jk = 1, jpkm1                                 ! Horizontal slab