Changeset 1176 for trunk/NEMO/TOP_SRC/LOBSTER/trcsed.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/trcsed.F90 (modified) (4 diffs)

trunk/NEMO/TOP_SRC/LOBSTER/trcsed.F90

@@ -18 +18 @@
    USE sms_lobster
    USE lbclnk
+   USE trdmld_trc
+   USE trdmld_trc_oce
 
    USE prtctl_trc      ! Print control for debbuging
@@ -56 +58 @@
       INTEGER, INTENT( in ) ::   kt      ! ocean time-step index      
       !!
-      INTEGER  ::   ji, jj, jk
+      INTEGER  ::   ji, jj, jk, jl
       REAL(wp) ::   ztra
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zwork
-
+      REAL(wp), DIMENSION(:,:,:), ALLOCATABLE ::   ztrbio
       CHARACTER (len=25) :: charout
       !!---------------------------------------------------------------------
@@ -76 +78 @@
       zwork(:,:,1  ) = 0.e0      ! surface value set to zero
       zwork(:,:,jpk) = 0.e0      ! bottom value  set to zero
+
+      IF( l_trdtrc )THEN
+         ALLOCATE( ztrbio(jpi,jpj,jpk) )
+!CDIR COLLAPSE
+         ztrbio(:,:,:) = tra(:,:,:,jpdet)
+      ENDIF
 
       ! tracer flux at w-point: we use -vsed (downward flux)  with simplification : no e1*e2
@@ -108 +116 @@
       !
 
+      IF( l_trdtrc ) THEN
+!CDIRR COLLAPSE
+         ztrbio(:,:,:) = tra(:,:,:,jpdet) - ztrbio(:,:,:)
+         jl = 8
+         CALL trd_mod_trc( ztrbio, jl, kt )   ! handle the trend
+      ENDIF
+
       IF(ln_ctl)   THEN  ! print mean trends (used for debugging)
          WRITE(charout, FMT="('sed')")