Changeset 620 for trunk/NEMO/OPA_SRC

Timestamp:

2007-02-22T08:57:39+01:00 (17 years ago)

Author:

opalod

Message:

nemo_v2_update_005 : CT : allow to run without horizontal diffusion in setting ln_traldf_lap & ln_traldf_bilap logicals to .false.

File:

• trunk/NEMO/OPA_SRC/TRA/traldf.F90 (modified) (3 diffs)

trunk/NEMO/OPA_SRC/TRA/traldf.F90

@@ -150 +150 @@
       IF( ln_traldf_lap   )   ioptio = ioptio + 1
       IF( ln_traldf_bilap )   ioptio = ioptio + 1
-      IF( ioptio /= 1 )   CALL ctl_stop( '          use ONE of the 2 lap/bilap operator type on tracer' )
+      IF( ioptio >  1 )   CALL ctl_stop( '          use ONE or NONE of the 2 lap/bilap operator type on tracer' )
+      IF( ioptio == 0 )   nldf = -2   ! No lateral diffusion
       ioptio = 0
       IF( ln_traldf_level )   ioptio = ioptio + 1
@@ -159 +160 @@
       ! defined the type of lateral diffusion from ln_traldf_... logicals
       ierr = 0
-      IF ( ln_traldf_lap ) THEN      ! laplacian operator
+      IF( ln_traldf_lap ) THEN       ! laplacian operator
          IF ( ln_zco ) THEN                ! z-coordinate
             IF ( ln_traldf_level )   nldf = 0      ! iso-level  (no rotation)
@@ -211 +212 @@
       IF(lwp) THEN
          WRITE(numout,*)
-         IF( nldf == -1 )   WRITE(numout,*) '              ESOPA test All scheme used'
-         IF( nldf ==  0 )   WRITE(numout,*) '              laplacian operator'
-         IF( nldf ==  1 )   WRITE(numout,*) '              Rotated laplacian operator'
-         IF( nldf ==  2 )   WRITE(numout,*) '              bilaplacian operator'
-         IF( nldf ==  3 )   WRITE(numout,*) '              Rotated bilaplacian'
+         IF( nldf == -2 )   WRITE(numout,*) '          NO lateral diffusion'
+         IF( nldf == -1 )   WRITE(numout,*) '          ESOPA test All scheme used'
+         IF( nldf ==  0 )   WRITE(numout,*) '          laplacian operator'
+         IF( nldf ==  1 )   WRITE(numout,*) '          Rotated laplacian operator'
+         IF( nldf ==  2 )   WRITE(numout,*) '          bilaplacian operator'
+         IF( nldf ==  3 )   WRITE(numout,*) '          Rotated bilaplacian'
       ENDIF