Changeset 911 for trunk/NEMO/OPA_SRC/DYN/dynspg_flt.F90

Timestamp:

2008-04-28T11:31:32+02:00 (16 years ago)

Author:

ctlod

Message:

Implementation of the BDY package, see ticket: #126

File:

• trunk/NEMO/OPA_SRC/DYN/dynspg_flt.F90 (modified) (4 diffs)

trunk/NEMO/OPA_SRC/DYN/dynspg_flt.F90

@@ -11 +11 @@
    !!            " "  !  05-11  (V. Garnier) Surface pressure gradient organization
    !!            " "  !  06-07  (S. Masson)  distributed restart using iom
+   !!            " "  !  05-01  (J.Chanut, A.Sellar) Calls to BDY routines. 
    !!----------------------------------------------------------------------
 #if defined key_dynspg_flt   ||   defined key_esopa  
@@ -36 +37 @@
    USE obcdyn          ! ocean open boundary condition (obc_dyn routines)
    USE obcvol          ! ocean open boundary condition (obc_vol routines)
+   USE bdy_oce         ! Unstructured open boundaries condition
+   USE bdydyn          ! Unstructured open boundaries condition (bdy_dyn routine) 
+   USE bdyvol          ! Unstructured open boundaries condition (bdy_vol routine)
    USE cla_dynspg      ! cross land advection
    USE in_out_manager  ! I/O manager
@@ -249 +253 @@
       CALL obc_vol( kt )      ! Correction of the barotropic componant velocity to control the volume of the system
 #endif
+#if defined key_bdy
+      ! Update velocities on unstructured boundary using the Flow Relaxation Scheme
+      CALL bdy_dyn_frs( kt )
+
+      IF (ln_bdy_vol) THEN
+      ! Correction of the barotropic component velocity to control the volume of the system
+        CALL bdy_vol( kt )
+      ENDIF
+#endif
 #if defined key_agrif
       CALL Agrif_dyn( kt )    ! Update velocities on each coarse/fine interfaces 
@@ -377 +390 @@
             spgu(ji,jj) = z2dt * ztdgu * obcumask(ji,jj)
             spgv(ji,jj) = z2dt * ztdgv * obcvmask(ji,jj)
+#elif defined key_bdy
+            ! caution : grad D = 0 along open boundaries
+            ! Remark: The filtering force could be reduced here in the FRS zone
+            !         by multiplying spgu/spgv by (1-alpha) ??  
+            spgu(ji,jj) = z2dt * ztdgu * bdyumask(ji,jj)
+            spgv(ji,jj) = z2dt * ztdgv * bdyvmask(ji,jj)           
 #else
             spgu(ji,jj) = z2dt * ztdgu