Changeset 473 for trunk/NEMO/OPA_SRC/LDF

Timestamp:

2006-05-11T17:04:37+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_060: SM: IOM + 301 levels + CORE + begining of ctl_stop

File:

• trunk/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90 (modified) (4 diffs)

trunk/NEMO/OPA_SRC/LDF/ldfdyn_c3d.h90

@@ -190 +190 @@
       !! * Modules used
       USE ldftra_oce, ONLY : aht0
- 
+
       !! * Arguments
       LOGICAL, INTENT (in) :: ld_print   ! If true, output arrays on numout
@@ -237 +237 @@
       ENDIF
 111   CONTINUE
-      IF( iost /= 0 ) THEN
-         IF(lwp) THEN
-            WRITE(numout,*)
-            WRITE(numout,*) ' ===>>>> : bad opening file: ahmcoef,  we stop. verify the file '
-            WRITE(numout,*) ' =======   ===  '
-         ENDIF
-         nstop = nstop + 1
-      ENDIF
-
+      IF( iost /= 0 ) CALL ctl_stop( ' ',  &
+           &    ' ===>>>> : bad opening file: ahmcoef, verify the file ahmcoef', &
+           &    ' =======   ===  ' )
       REWIND inum
       READ(inum,9101) clexp, iim, ijm
@@ -375 +369 @@
       ! other level: re-increase the coef in the deep ocean
       
-      DO jk = 1, 21
+#if defined key_orca_lev10
+      DO jk = 1, 210
+         zcoef(jk) = 1.
+      END DO
+      DO jk= 211, 230
+         zcoef(jk) = 1. + 0.1 * FLOAT(jk-210)
+      END DO
+      DO jk= 231, 260
+         zcoef(jk) = 3. + 0.2 * FLOAT(jk-230)
+      END DO
+      DO jk= 261, 270
+         zcoef(jk) = 9. + 0.1 * FLOAT(jk-260)
+      END DO
+      DO jk= 271, jpk
+         zcoef(jk) = 10.
+      END DO
+      DO jk= 1, jpk
+         IF(lwp) WRITE(numout,*) 'k= ',jk, 'cof ', zcoef(jk)
+      END DO
+#else
+       DO jk = 1, 21
          zcoef(jk) = 1.
       END DO
@@ -386 +400 @@
          zcoef(jk) = 10.
       END DO
-      
+#endif
+
       DO jk = 2, jpk
          ahm1(:,:,jk) = MIN( zahm0(:,:), zcoef(jk) * ahm1(:,:,1) )