Changeset 416

Timestamp:

2006-03-21T08:52:26+01:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_039 : CT : add some tests to catch user's mistakes

File:

• trunk/NEMO/OPA_SRC/OBC/obcini.F90 (modified) (1 diff)

trunk/NEMO/OPA_SRC/OBC/obcini.F90

@@ -197 +197 @@
       ! ...  a small arbitrary value, (to avoid division by zero further on). 
       ! ...  rdpxxx is not used anymore.
-      IF( lp_obc_east .AND. ( rdpein == 0 .AND. rdpeob == 0 ) )  THEN
-         lfbceast = .TRUE.
-         rdpein = 1e-3
-         rdpeob = 1e-3
-         inumfbc = inumfbc+1
-      END IF
-      IF( lp_obc_west .AND. ( rdpwin == 0 .AND. rdpwob == 0 ) )  THEN
-         lfbcwest = .TRUE.
-         rdpwin = 1e-3
-         rdpwob = 1e-3
-         inumfbc = inumfbc+1
-      END IF
-      IF( lp_obc_north .AND. ( rdpnin == 0 .AND. rdpnob == 0 ) )  THEN
-         lfbcnorth = .TRUE.
-         rdpnin = 1e-3
-         rdpnob = 1e-3
-         inumfbc = inumfbc+1
-      END IF
-      IF( lp_obc_south .AND. ( rdpsin == 0 .AND. rdpsob == 0 ) )  THEN
-         lfbcsouth = .TRUE.
-         rdpsin = 1e-3
-         rdpsob = 1e-3
-         inumfbc = inumfbc+1
+      IF( lp_obc_east )  THEN
+         IF( (rdpein+rdpeob) == 0 )  THEN
+            lfbceast = .TRUE.
+            rdpein = 1e-3
+            rdpeob = 1e-3
+            inumfbc = inumfbc+1
+         ELSEIF ( (rdpein*rdpeob) == 0 )  THEN
+            IF(lwp) THEN
+               WRITE(numout,cform_err)
+               WRITE(numout,*) 'obc_init : rdpein & rdpeob must be both zero or non zero'
+               nstop = nstop + 1
+            ENDIF
+         END IF
+      END IF
+      IF( lp_obc_west )  THEN
+         IF( (rdpwin + rdpwob) == 0 )  THEN
+            lfbcwest = .TRUE.
+            rdpwin = 1e-3
+            rdpwob = 1e-3
+            inumfbc = inumfbc+1
+         ELSEIF ( (rdpwin*rdpwob) == 0 )  THEN
+            IF(lwp) THEN
+               WRITE(numout,cform_err)
+               WRITE(numout,*) 'obc_init : rdpwin & rdpwob must be both zero or non zero'
+               nstop = nstop + 1
+            ENDIF
+         END IF
+      END IF
+      IF( lp_obc_north )  THEN
+         IF( (rdpnin + rdpnob) == 0 )  THEN
+            lfbcnorth = .TRUE.
+            rdpnin = 1e-3
+            rdpnob = 1e-3
+            inumfbc = inumfbc+1
+         ELSEIF ( (rdpnin*rdpnob) == 0 )  THEN
+            IF(lwp) THEN
+               WRITE(numout,cform_err)
+               WRITE(numout,*) 'obc_init : rdpnin & rdpnob must be both zero or non zero'
+               nstop = nstop + 1
+            ENDIF
+         END IF
+      END IF
+      IF( lp_obc_south )  THEN
+         IF( (rdpsin + rdpsob) == 0 )  THEN
+            lfbcsouth = .TRUE.
+            rdpsin = 1e-3
+            rdpsob = 1e-3
+            inumfbc = inumfbc+1
+         ELSEIF ( (rdpsin*rdpsob) == 0 )  THEN
+            IF(lwp) THEN
+               WRITE(numout,cform_err)
+               WRITE(numout,*) 'obc_init : rdpsin & rdpsob must be both zero or non zero'
+               nstop = nstop + 1
+            ENDIF
+         END IF
       END IF