Changeset 1176 for trunk/NEMO/TOP_SRC/LOBSTER/trcexp.F90

Timestamp:

2008-09-11T18:26:59+02:00 (16 years ago)

Author:

cetlod

Message:

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

File:

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

trunk/NEMO/TOP_SRC/LOBSTER/trcexp.F90

@@ -21 +21 @@
    USE trctrp_lec
    USE prtctl_trc      ! Print control for debbuging
+   USE trdmld_trc
+   USE trdmld_trc_oce
 
    IMPLICIT NONE
@@ -52 +54 @@
       INTEGER, INTENT( in ) ::   kt      ! ocean time-step index      
       !!
-      INTEGER  ::   ji, jj, jk
+      INTEGER  ::   ji, jj, jk, jl
       REAL(wp) ::   zgeolpoc, zfact
       INTEGER , DIMENSION(jpi,jpj) ::   ikbot
       REAL(wp), DIMENSION(jpi,jpj) ::   zwork
-
+      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::   ztrbio
       CHARACTER (len=25) :: charout
       !!---------------------------------------------------------------------
@@ -66 +68 @@
       ENDIF
 
-! VERTICAL DISTRIBUTION OF NEWLY PRODUCED BIOGENIC
-! POC IN THE WATER COLUMN
-! (PARTS OF NEWLY FORMED MATTER REMAINING IN THE DIFFERENT
-! LAYERS IS DETERMINED BY DMIN3 DEFINED IN common.passivetrc.*.h
-! ----------------------------------------------------------------------
+      ! VERTICAL DISTRIBUTION OF NEWLY PRODUCED BIOGENIC
+      ! POC IN THE WATER COLUMN
+      ! (PARTS OF NEWLY FORMED MATTER REMAINING IN THE DIFFERENT
+      ! LAYERS IS DETERMINED BY DMIN3 DEFINED IN sms_lobster.F90
+      ! ----------------------------------------------------------------------
+
+      IF( l_trdtrc )THEN
+         ALLOCATE( ztrbio(jpi,jpj,jpk) )
+!CDIR COLLAPSE
+         ztrbio(:,:,:) = tra(:,:,:,jpno3)
+      ENDIF
 
       DO jk = 1, jpkm1
@@ -166 +174 @@
       ENDIF
       !
+      IF( l_trdtrc ) THEN
+!CDIRR COLLAPSE
+         ztrbio(:,:,:) = tra(:,:,:,jpno3) - ztrbio(:,:,:)
+         jl = 17
+         CALL trd_mod_trc( ztrbio, jl, kt )   ! handle the trend
+      ENDIF
 
       IF(ln_ctl)   THEN  ! print mean trends (used for debugging)