Changeset 261 for trunk/NEMO/TOP_SRC/TRP/trcstp.F90

Timestamp:

2005-09-08T15:18:32+02:00 (19 years ago)

Author:

opalod

Message:

nemo_v1_update_05:RB+OA:Update and rewritting of (part of) the TOP component

File:

• trunk/NEMO/TOP_SRC/TRP/trcstp.F90 (modified) (5 diffs)

trunk/NEMO/TOP_SRC/TRP/trcstp.F90

@@ -12 +12 @@
    USE trc              ! ocean passive tracers variables 
    USE trctrp           ! passive tracers transport
+   USE trcsms           ! passive tracers sources and sinks
 
    IMPLICIT NONE
@@ -19 +20 @@
    PUBLIC trc_stp           ! called by step
    !!----------------------------------------------------------------------
-   !!   TOP 1.0 , LOCEAN-IPSL (2005) 
-   !! $Header$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
+   !!   OPA 9.0 , LODYC-IPSL (2003)
    !!----------------------------------------------------------------------
 
@@ -41 +40 @@
       !! * Arguments
       INTEGER, INTENT( in ) ::  kt  ! ocean time-step index
+      INTEGER :: jn
+      REAL(wp) :: ztra
 
 
@@ -49 +50 @@
 
 
-      CALL trcsms( kt )
+      CALL trc_sms( kt )
+
+       DO jn = 1, jptra
+         IF(ln_ctl) THEN         ! print mean trends (used for debugging)
+            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )
+            WRITE(numout,*) ' trc/sms  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)
+            tra_ctl(jn) = ztra
+         ENDIF
+      END DO
+
  
 
@@ -55 +65 @@
 
       CALL trc_trp( kt )
+
+       DO jn = 1, jptra
+         IF(ln_ctl) THEN         ! print mean trends (used for debugging)
+            ztra = SUM( tra(2:nictle,2:njctle,1:jpkm1,jn) * tmask(2:nictle,2:njctle,1:jpkm1) )
+            WRITE(numout,*) ' trc/trp  - ',ctrcnm(jn),' : ', ztra-tra_ctl(jn)
+            tra_ctl(jn) = ztra
+         ENDIF
+      END DO
+