Changeset 64

Timestamp:

2004-04-22T14:02:25+02:00 (20 years ago)

Author:

opalod

Message:

CT : BUGFIX038 : Indices correction which concerns all ORCA configurations with NO closed Seas

File:

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

trunk/NEMO/OPA_SRC/DOM/closea.F90

@@ -77 +77 @@
       
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*)' dom_clo : closed seas '
-      IF(lwp) WRITE(numout,*)' ~~~~~~~'
+      IF(lwp) WRITE(numout,*)'dom_clo : closed seas '
+      IF(lwp) WRITE(numout,*)'~~~~~~~'
 
       ! initial values
@@ -153 +153 @@
          ncsj1(jc)   = mj0( ncsj1(jc) )
 
-         ncsi2(jc)   = mi1( ncsi1(jc) )
-         ncsj2(jc)   = mj1( ncsj1(jc) )
+         ncsi2(jc)   = mi1( ncsi2(jc) )   
+         ncsj2(jc)   = mj1( ncsj2(jc) )  
       END DO
          
@@ -191 +191 @@
       IF( kt == nit000 ) THEN 
          IF(lwp) WRITE(numout,*)
-         IF(lwp) WRITE(numout,*)' flx_clo : closed seas '
-         IF(lwp) WRITE(numout,*)' ~~~~~~~'
+         IF(lwp) WRITE(numout,*)'flx_clo : closed seas '
+         IF(lwp) WRITE(numout,*)'~~~~~~~'
 
          ! Total surface of ocean
@@ -233 +233 @@
       ! 2 - Computation
       ! ---------------
-      zemp = 0.0e0
+      zemp = 0.e0
 
       DO jc = 1, jpncs
@@ -258 +258 @@
          ELSEIF( ncstt(jc) == 1 ) THEN 
             ! Excess water in open sea, at outflow location, excess evap shared
-            IF ( zemp(jc) <= 0.0e0 ) THEN 
+            IF ( zemp(jc) <= 0.e0 ) THEN 
                 DO jn = 1, ncsnr(jc)
                   ji = mi0(ncsir(jc,jn))