Changeset 4328 for branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_flt.F90

Timestamp:

2013-12-06T11:25:13+01:00 (10 years ago)

Author:

davestorkey

Message:

Remove OBC module at NEMO 3.6. See ticket #1189.

File:

• branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_flt.F90 (modified) (4 diffs)

branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DYN/dynspg_flt.F90

@@ -25 +25 @@
    USE zdf_oce         ! ocean vertical physics
    USE sbc_oce         ! surface boundary condition: ocean
-   USE obc_oce         ! Lateral open boundary condition
    USE bdy_oce         ! Lateral open boundary condition
    USE sol_oce         ! ocean elliptic solver
@@ -34 +33 @@
    USE solpcg          ! preconditionned conjugate gradient solver
    USE solsor          ! Successive Over-relaxation solver
-   USE obcdyn          ! ocean open boundary condition on dynamics
-   USE obcvol          ! ocean open boundary condition (obc_vol routine)
    USE bdydyn          ! ocean open boundary condition on dynamics
    USE bdyvol          ! ocean open boundary condition (bdy_vol routine)
@@ -184 +181 @@
       ENDIF
 
-#if defined key_obc
-      IF( lk_obc ) CALL obc_dyn( kt )   ! Update velocities on each open boundary with the radiation algorithm
-      IF( lk_obc ) CALL obc_vol( kt )   ! Correction of the barotropic componant velocity to control the volume of the system
-#endif
 #if defined key_bdy
       IF( lk_bdy ) CALL bdy_dyn( kt )   ! Update velocities on each open boundary
@@ -304 +297 @@
             ztdgv = z2dtg * (gcx(ji  ,jj+1) - gcx(ji,jj) ) / e2v(ji,jj)
             ! multiplied by z2dt
-#if defined key_obc
-            IF(lk_obc) THEN
-            ! 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 * obcumask(ji,jj)
-               spgv(ji,jj) = z2dt * ztdgv * obcvmask(ji,jj)
-            ELSE
-               spgu(ji,jj) = z2dt * ztdgu
-               spgv(ji,jj) = z2dt * ztdgv
-            ENDIF
-#elif defined key_bdy
+#if defined key_bdy
             IF(lk_bdy) THEN
             ! caution : grad D = 0 along open boundaries