Changeset 6861 for branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/DOM

Timestamp:

2016-08-09T17:55:10+02:00 (8 years ago)

Author:

clevy

Message:

Add north south cyclic boundaries option, see ticket 1761

Location:

branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/DOM

Files:

• domcfg.F90 (modified) (1 diff)
• domzgr.F90 (modified) (4 diffs)

branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/DOM/domcfg.F90

@@ -56 +56 @@
          IF( jperio == 5 )   WRITE(numout,*) '      jperio= 5, north fold with F-point pivot'
          IF( jperio == 6 )   WRITE(numout,*) '      jperio= 6, cyclic east-west and north fold with F-point pivot'
-      ENDIF
-      !
-      IF( jperio <  0 .OR. jperio > 6 )   CALL ctl_stop( 'jperio is out of range' )
+         IF( jperio == 7 )   WRITE(numout,*) '      jperio= 7, cyclic east-west and north-south'
+      ENDIF
+      !
+      IF( jperio <  0 .OR. jperio > 7 )   CALL ctl_stop( 'jperio is out of range' )
       !
       CALL dom_glo                   ! global domain versus zoom and/or local domain

branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/DOM/domzgr.F90

@@ -466 +466 @@
             idta( 1    , :    ) = 0                 ;      zdta( 1    , :    ) =  0._wp
             idta(jpidta, :    ) = 0                 ;      zdta(jpidta, :    ) =  0._wp
+         ELSEIF( jperio == 7 ) THEN
+!           Nothing to do here
          ELSE
             ih = 0                                  ;      zh = 0._wp
@@ -730 +732 @@
          IF( lk_mpp ) THEN
             IF( nbondi == -1 .OR. nbondi == 2 ) THEN
-               IF( jperio /= 1 )   mbathy(1,:) = 0
+               IF( jperio /= 1 .AND. jperio /= 7 )   mbathy(1,:) = 0
             ENDIF
             IF( nbondi == 1 .OR. nbondi == 2 ) THEN
-               IF( jperio /= 1 )   mbathy(nlci,:) = 0
+               IF( jperio /= 1 .AND. jperio /= 7 )   mbathy(nlci,:) = 0
             ENDIF
          ELSE
@@ -748 +750 @@
          mbathy( 1 ,:) = mbathy(jpim1,:)
          mbathy(jpi,:) = mbathy(  2  ,:)
+         IF (jperio == 7) THEN
+            IF(lwp) WRITE(numout,*)' north south boundary conditions on mbathy: jperio = ', jperio
+            mbathy( : ,1) = mbathy(:, jpjm1)
+            mbathy(:, jpj)= mbathy(:,2)
+         ENDIF
       ELSEIF( nperio == 2 ) THEN
          IF(lwp) WRITE(numout,*) '   equatorial boundary conditions on mbathy: nperio = ', nperio
@@ -779 +786 @@
       ENDIF
 
-      IF( lwp .AND. nprint == 1 ) THEN      ! control print
+!!      IF( lwp .AND. nprint == 1 ) THEN      ! control print
+      IF( lwp ) THEN
          WRITE(numout,*)
          WRITE(numout,*) ' bathymetric field :   number of non-zero T-levels '