Changeset 732 for tags/nemo_v2_3/NEMO/TOP_SRC/trcdit.F90

Timestamp:

2007-11-19T16:19:15+01:00 (17 years ago)

Author:

cetlod

Message:

Update modules for passive tracers transport trends computation, see ticket:13

File:

• tags/nemo_v2_3/NEMO/TOP_SRC/trcdit.F90 (modified) (5 diffs)

tags/nemo_v2_3/NEMO/TOP_SRC/trcdit.F90

@@ -2 +2 @@
    !!----------------------------------------------------------------------
    !!  TOP 1.0,  LOCEAN-IPSL (2005) 
-   !! $Header$ 
+   !! $Header: /home/opalod/NEMOCVSROOT/NEMO/TOP_SRC/trcdit.F90,v 1.9 2007/10/12 09:22:19 opalod Exp $ 
    !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
    !!----------------------------------------------------------------------
@@ -418 +418 @@
                   IF (jl.eq.7) THEN
 ! short and long title for x gent velocity for tracer
-                      WRITE (cltra,'("XGV",16a)') ctrcnm(jn)
+                      WRITE (cltra,'("XGV_",16a)') ctrcnm(jn)
                       WRITE (cltral,'("X gent velocity trend for ",53a)')  &
                       &      ctrcnl(jn)(1:53)
@@ -436 +436 @@
 # endif
 # if defined key_trcdmp
-                  IF (jl.eq.jpdiatrc)  THEN
+                  IF (jl.eq.jpdiatrc-1)  THEN
 ! last trends for tracer damping : short and long title
                       WRITE (cltra,'("TDM_",16a)') ctrcnm(jn)
@@ -443 +443 @@
                   END IF
 # endif
+                  IF (jl.eq.jpdiatrc)  THEN
+! last trends for tracer damping : short and long title
+                      WRITE (cltra,'("SBC_",16a)') ctrcnm(jn)
+                      WRITE (cltral,'("Surface boundary flux ",58a)')  &
+                      &      ctrcnl(jn)(1:58)
+                  END IF
+
                   call flush(numout)
                   cltrau=ctrcun(jn) ! UNIT for tracer /trends
@@ -518 +525 @@
 # endif
 # if defined key_trcdmp
-                  IF (jl.eq.jpdiatrc) THEN
+                  IF (jl.eq.jpdiatrc-1) THEN
 ! short for x gent velocity for tracer
                       WRITE (cltra,'("TDM_",16a)') ctrcnm(jn)
                   END IF
 # endif
+                  IF (jl.eq.jpdiatrc) THEN
+! short for surface boundary conditions for tracer
+                      WRITE (cltra,'("SBC_",a)') ctrcnm(jn)
+                  END IF
 
                   CALL histwrite(nit6(jn), cltra, it, trtrd(:,:,:,ikeep(jn),jl)  &