Changeset 232 for trunk/NEMO/OPA_SRC/lbclnk.F90

Timestamp:

2005-03-22T11:12:44+01:00 (19 years ago)

Author:

opalod

Message:

CT : BUGFIX163 : - introduction of a sign change for 'F' and 'G' points on 3D and 2D fields for the North fold boundary condition

File:

• trunk/NEMO/OPA_SRC/lbclnk.F90 (modified) (5 diffs)

trunk/NEMO/OPA_SRC/lbclnk.F90

@@ -1 +1 @@
 MODULE lbclnk
    !!======================================================================
-   !!                       ***  MODULE  zdftke  ***
+   !!                       ***  MODULE  lbclnk  ***
    !! Ocean        : lateral boundary conditions
    !!=====================================================================
@@ -156 +156 @@
                   DO ji = 1, jpi-1
                      iju = jpi-ji+1
-                     pt3d(ji,jpj-1,jk) = pt3d(iju,jpj-2,jk)
-                     pt3d(ji,jpj  ,jk) = pt3d(iju,jpj-3,jk)
+                     pt3d(ji,jpj-1,jk) = psgn * pt3d(iju,jpj-2,jk)
+                     pt3d(ji,jpj  ,jk) = psgn * pt3d(iju,jpj-3,jk)
                   END DO
             END SELECT
@@ -192 +192 @@
                   DO ji = 1, jpi-1
                      iju = jpi-ji
-                     pt3d(ji,jpj  ,jk) = pt3d(iju,jpj-2,jk)
+                     pt3d(ji,jpj  ,jk) = psgn * pt3d(iju,jpj-2,jk)
                   END DO
                   DO ji = jpi/2+1, jpi-1
                      iju = jpi-ji
-                     pt3d(ji,jpjm1,jk) = pt3d(iju,jpjm1,jk)
+                     pt3d(ji,jpjm1,jk) = psgn * pt3d(iju,jpjm1,jk)
                   END DO
             END SELECT
@@ -328 +328 @@
             DO ji = 1, jpi-1
                iju = jpi - ji + 1
-               pt2d(ji,jpj-1) = pt2d(iju,jpj-2)
-               pt2d(ji,jpj  ) = pt2d(iju,jpj-3)
+               pt2d(ji,jpj-1) = psgn * pt2d(iju,jpj-2)
+               pt2d(ji,jpj  ) = psgn * pt2d(iju,jpj-3)
             END DO
 
@@ -377 +377 @@
             DO ji = 1, jpi-1
                iju = jpi-ji
-               pt2d(ji,jpj  ) = pt2d(iju,jpj-2)
+               pt2d(ji,jpj  ) = psgn * pt2d(iju,jpj-2)
             END DO
             DO ji = jpi/2+1, jpi-1
                iju = jpi-ji
-               pt2d(ji,jpjm1) = pt2d(iju,jpjm1)
+               pt2d(ji,jpjm1) = psgn * pt2d(iju,jpjm1)
             END DO
 
          CASE ( 'I' )                                  ! ice U-V point
-            pt2d( 2 ,jpj) = 0.e0           !!bug  ???
-            DO ji = 1 , jpim1            !!bug rob= 2,jpim1
-               ijt = jpi - ji            !!bug rob= ijt=jpi-ji+2   ???
+            pt2d( 2 ,jpj) = 0.e0
+            DO ji = 2 , jpim1
+               ijt = jpi - ji + 2
                pt2d(ji,jpj)= 0.5 * ( pt2d(ji,jpjm1) + psgn * pt2d(ijt,jpjm1) )
             END DO