Changeset 1256

Timestamp:

2009-01-13T11:32:36+01:00 (15 years ago)

Author:

cetlod

Message:

Add a namelist parameter to use or not the trcrad routine, see ticket:300

Location:

trunk/NEMO/TOP_SRC/TRP

Files:

• trctrp.F90 (modified) (2 diffs)
• trctrp_lec.F90 (modified) (3 diffs)

trunk/NEMO/TOP_SRC/TRP/trctrp.F90

@@ -13 +13 @@
    !!----------------------------------------------------------------------
    USE oce_trc         ! ocean dynamics and active tracers variables
-   USE trp_trc             ! ocean passive tracers variables 
+   USE trp_trc         ! ocean passive tracers variables 
    USE trctrp_lec      ! passive tracers transport parameters
-   USE prtctl_trc     ! Print control for debbuging
+   USE prtctl_trc      ! Print control for debbuging
 
    USE trcbbl          ! bottom boundary layer               (trc_bbl routine)
@@ -112 +112 @@
 
                                CALL trc_nxt( kt )            ! tracer fields at next time step
- 
-                               CALL trc_rad( kt )            ! Correct artificial negative concentrations
+     
+      IF( ln_trcrad )          CALL trc_rad( kt )            ! Correct artificial negative concentrations
       !                                                      ! especially useful when isopycnal mixing is used
       !                                                      

trunk/NEMO/TOP_SRC/TRP/trctrp_lec.F90

@@ -38 +38 @@
    LOGICAL , PUBLIC ::   ln_trcldf_hor   = .FALSE.    !: horizontal (geopotential) direction
    LOGICAL , PUBLIC ::   ln_trcldf_iso   = .TRUE.     !: iso-neutral direction
+
+   !                                                 !!: ** Treatment of Negative concentrations ( nam_trcrad )
+   LOGICAL , PUBLIC ::   ln_trcrad       = .TRUE.     !: flag to artificially correct negative concentrations
 
    !                                                 !!: flag of the lateral diff. scheme used
@@ -96 +99 @@
          &                 ahtrc0, ahtrb0, aeivtr0, trcrat
       NAMELIST/namtopzdf/ ln_trczdf_exp, n_trczdf_exp
+      NAMELIST/namtoprad/ ln_trcrad
 #if defined key_trcdmp
       NAMELIST/namtopdmp/ ndmptr, ndmpftr, nmldmptr, sdmptr, bdmptr, hdmptr
@@ -169 +173 @@
       ENDIF
 
+      !
+      REWIND( numnat )                 !   Read Namelist namtoprad
+      READ  ( numnat, namtoprad )
+
+      IF(lwp) THEN                     !   ! Control print
+         WRITE(numout,*)
+         WRITE(numout,*) '   Namelist namtoprad : treatment of negative concentrations'
+         WRITE(numout,*) '      correct artificially negative concen. or not ln_trcrad = ', ln_trcrad
+      ENDIF
+
+
 # if defined key_trcdmp
       !                                ! passive tracres damping term