Changeset 1177 for trunk/NEMO/TOP_SRC/trclec.F90

Timestamp:

2008-09-11T18:27:52+02:00 (16 years ago)

Author:

cetlod

Message:

update modules to take into account new trends organization, see ticket:248

File:

• trunk/NEMO/TOP_SRC/trclec.F90 (modified) (4 diffs)

trunk/NEMO/TOP_SRC/trclec.F90

@@ -27 +27 @@
    USE trclsm_my_trc     ! MY_TRC SMS namelist
    USE in_out_manager    ! I/O manager
+   USE trdmld_trc_oce
 
    IMPLICIT NONE
@@ -69 +70 @@
       !!
       NAMELIST/namtoptrc/ ndttrc, nwritetrc, lrsttr, nrsttr, tracer
-#if defined key_trc_diatrd
-      NAMELIST/namtoptrd/ luttrd, nwritetrd                                             ! dynamical trends
+#if defined key_trdmld_trc || defined key_trc_diatrd
+      NAMELIST/namtoptrd/ ntrd_trc, nctls_trc, ucf_trc, &
+                         ln_trdmld_trc_restart, ln_trdmld_trc_instant, &
+                         luttrd
 #endif
 
@@ -131 +134 @@
       ENDIF
 
-#if defined key_trc_diatrd
+#if defined key_trdmld_trc || defined key_trc_diatrd
 
       ! Namelist natrtd (transport trends)
       ! ----------------------------------------------
-      nwritetrd = 10                 ! default values (no dynamical trend recording)
+      ntrd_trc  = 20
+      nctls_trc =  6
+      ucf_trc   =  1.
+      ln_trdmld_trc_instant = .TRUE.
+      ln_trdmld_trc_restart =.FALSE.
       DO jn = 1, jptra
          luttrd(jn) = .FALSE.
       END DO
 
-      REWIND( numnat )               ! read natrtd
+      REWIND( numnat )               !  namelist namtoptrd : passive tracer trends diagnostic
       READ  ( numnat, namtoptrd )
 
-      nkeep=0
-      ikeep(:)=0
+#if defined key_trc_diatrd
+      nkeep    = 0
+      ikeep(:) = 0
       DO jn = 1, jptra
          IF( luttrd(jn) ) THEN 
-             nkeep    = nkeep + 1
-             ikeep(jn)=nkeep
+             nkeep     = nkeep + 1
+             ikeep(jn) = nkeep
          END IF 
       END DO
@@ -156 +164 @@
       ENDIF 
 
-      IF(lwp) THEN                   ! control print
+      nwritetrd = ntrd_trc
+#endif
+
+     IF(lwp) THEN
          WRITE(numout,*)
-         WRITE(numout,*) ' Namelist : namtoptrd'
-         WRITE(numout,*) '    frequency of outputs for dynamical trends nwritetrd = ', nwritetrd
+         WRITE(numout,*) ' trd_mld_trc_init : read namelist namtoptrd                        '
+         WRITE(numout,*) ' ~~~~~~~~~~~~~~~~                                               '
+         WRITE(numout,*) '   * frequency of trends diagnostics   ntrd_trc               = ', ntrd_trc
+         WRITE(numout,*) '   * control surface type              nctls_trc              = ', nctls_trc
+         WRITE(numout,*) '   * restart for ML diagnostics        ln_trdmld_trc_restart  = ', ln_trdmld_trc_restart
+         WRITE(numout,*) '   * flag to diagnose trends of                                 '
+         WRITE(numout,*) '     instantantaneous or mean ML T/S   ln_trdmld_trc_instant  = ', ln_trdmld_trc_instant
+         WRITE(numout,*) '   * unit conversion factor            ucf_trc                = ', ucf_trc
+#if defined key_trdmld_trc
+         DO jn = 1, jptra
+            WRITE(numout,*) '    keep dynamical trends for tracer number :', jn, luttrd(jn)
+         END DO
+#else
          DO jn = 1, jptra
             WRITE(numout,*) '    keep dynamical trends for tracer number :', jn, luttrd(jn), ikeep(jn)
          END DO
+
+         WRITE(numout,*)
          WRITE(numout,*) '    total = ', nkeep, ' tracers dyn trends saved'
          WRITE(numout,*) '    size of trtrd = ', jpi*jpj*jpk*nkeep*jpdiatrc
-      ENDIF
+#endif
+      ENDIF
+
 #endif