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

Timestamp:

2004-04-22T15:50:27+02:00 (20 years ago)

Author:

opalod

Message:

CT : UPDATE057 : # General syntax, alignement, comments corrections

# l_ctl alone replace the set (l_ctl .AND. lwp)
# Add of diagnostics which are activated when using l_ctl logical

File:

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

trunk/NEMO/OPA_SRC/lbclnk.F90

@@ -72 +72 @@
          !             !   =-1 , the sign is changed if north fold boundary
          !             !   = 1 , no sign change
+         !             !   = 0 , no sign change and > 0 required (use the inner
+         !             !         row/column if closed boundary)
 
       !! * Local declarations
@@ -89 +91 @@
 
          CASE ( 1 , 4 , 6 )                    ! * cyclic east-west
-
             pt3d( 1 ,:,jk) = pt3d(jpim1,:,jk)          ! all points
             pt3d(jpi,:,jk) = pt3d(  2  ,:,jk)
 
          CASE DEFAULT                          ! * closed
-
             SELECT CASE ( cd_type )
-
             CASE ( 'T' , 'U' , 'V' , 'W' )             ! T-, U-, V-, W-points
                pt3d( 1 ,:,jk) = 0.e0
                pt3d(jpi,:,jk) = 0.e0
-
             CASE ( 'F' )                               ! F-point
                pt3d(jpi,:,jk) = 0.e0
-
             END SELECT
 
@@ -115 +112 @@
 
             SELECT CASE ( cd_type )
-
             CASE ( 'T' , 'U' , 'W' )                   ! T-, U-, W-points
                pt3d(:, 1 ,jk) = pt3d(:,3,jk)
                pt3d(:,jpj,jk) = 0.e0
-
             CASE ( 'V' , 'F' )                         ! V-, F-points
                pt3d(:, 1 ,jk) = psgn * pt3d(:,2,jk)
                pt3d(:,jpj,jk) = 0.e0
-
             END SELECT
 
-         CASE ( 3 , 4 )                            ! *  North fold  T-point pivot
+         CASE ( 3 , 4 )                        ! *  North fold  T-point pivot
 
             pt3d( 1 ,jpj,jk) = 0.e0
@@ -132 +126 @@
 
             SELECT CASE ( cd_type )
-
             CASE ( 'T' , 'W' )                         ! T-, W-point
                DO ji = 2, jpi
@@ -143 +136 @@
                   pt3d(ji,jpjm1,jk) = psgn * pt3d(ijt,jpjm1,jk)
                END DO
-
             CASE ( 'U' )                               ! U-point
                DO ji = 1, jpi-1
@@ -154 +146 @@
                   pt3d(ji,jpjm1,jk) = psgn * pt3d(iju,jpjm1,jk)
                END DO
-
             CASE ( 'V' )                               ! V-point
                   DO ji = 2, jpi
@@ -162 +153 @@
                      pt3d(ji,jpj  ,jk) = psgn * pt3d(ijt,jpj-3,jk)
                   END DO
-
             CASE ( 'F' )                               ! F-point
                   DO ji = 1, jpi-1
@@ -169 +159 @@
                      pt3d(ji,jpj  ,jk) = pt3d(iju,jpj-3,jk)
                   END DO
-
             END SELECT
 
@@ -178 +167 @@
 
             SELECT CASE ( cd_type )
-
             CASE ( 'T' , 'W' )                         ! T-, W-point
                DO ji = 1, jpi
@@ -185 +173 @@
                   pt3d(ji,jpj,jk) = psgn * pt3d(ijt,jpj-1,jk)
                END DO
-
             CASE ( 'U' )                               ! U-point
                   DO ji = 1, jpi-1
@@ -192 +179 @@
                      pt3d(ji,jpj,jk) = psgn * pt3d(iju,jpj-1,jk)
                   END DO
-
             CASE ( 'V' )                               ! V-point
                   DO ji = 1, jpi
@@ -203 +189 @@
                      pt3d(ji,jpjm1,jk) = psgn * pt3d(ijt,jpjm1,jk)
                   END DO
-
             CASE ( 'F' )                               ! F-point
                   DO ji = 1, jpi-1
@@ -215 +200 @@
             END SELECT
 
-         CASE DEFAULT                           ! *  closed
+         CASE DEFAULT                          ! *  closed
 
             SELECT CASE ( cd_type )
-
             CASE ( 'T' , 'U' , 'V' , 'W' )             ! T-, U-, V-, W-points
                pt3d(:, 1 ,jk) = 0.e0
                pt3d(:,jpj,jk) = 0.e0
-
             CASE ( 'F' )                               ! F-point
                pt3d(:,jpj,jk) = 0.e0
-
             END SELECT
 
@@ -273 +255 @@
 
       CASE ( 1 , 4 , 6 )                       ! * cyclic east-west
-
          pt2d( 1 ,:) = pt2d(jpim1,:)
          pt2d(jpi,:) = pt2d(  2  ,:)
 
       CASE DEFAULT                             ! * closed 
-
          SELECT CASE ( cd_type )
-
          CASE ( 'T' , 'U' , 'V' , 'W' )                ! T-, U-, V-, W-points
             pt2d( 1 ,:) = 0.e0
             pt2d(jpi,:) = 0.e0
-
-         CASE ( 'F' , 'I' )                            ! F-point, ice U-V point
+         CASE ( 'F' )                                  ! F-point, ice U-V point
             pt2d(jpi,:) = 0.e0 
-
+         CASE ( 'I' )                                  ! F-point, ice U-V point
+            pt2d( 1 ,:) = 0.e0 
+            pt2d(jpi,:) = 0.e0 
          END SELECT
 
@@ -299 +279 @@
 
          SELECT CASE ( cd_type )
-
          CASE ( 'T' , 'U' , 'W' )                      ! T-, U-, W-points
             pt2d(:, 1 ) = pt2d(:,3)
             pt2d(:,jpj) = 0.e0
-
-         CASE ( 'V' , 'F' )                            ! V-, F-points, ice U-V point
+         CASE ( 'V' , 'F' , 'I' )                      ! V-, F-points, ice U-V point
             pt2d(:, 1 ) = psgn * pt2d(:,2)
             pt2d(:,jpj) = 0.e0
-
          END SELECT
 
@@ -419 +396 @@
 
          SELECT CASE ( cd_type )
-
          CASE ( 'T' , 'U' , 'V' , 'W' )                ! T-, U-, V-, W-points
             pt2d(:, 1 ) = 0.e0
             pt2d(:,jpj) = 0.e0
-
          CASE ( 'F' )                                  ! F-point
             pt2d(:,jpj) = 0.e0
-
          CASE ( 'I' )                                  ! ice U-V point
             pt2d(:, 1 ) = 0.e0
             pt2d(:,jpj) = 0.e0
-
          END SELECT