Changeset 216 for trunk/NEMO/OPA_SRC/TRA/traldf_bilapg.F90

Timestamp:

2005-03-17T15:02:38+01:00 (19 years ago)

Author:

opalod

Message:

CT : UPDATE151 : New trends organization

File:

• trunk/NEMO/OPA_SRC/TRA/traldf_bilapg.F90 (modified) (5 diffs)

trunk/NEMO/OPA_SRC/TRA/traldf_bilapg.F90

@@ -16 +16 @@
    USE dom_oce         ! ocean space and time domain variables
    USE ldftra_oce      ! ocean active tracers: lateral physics
-   USE trdtra_oce      ! ocean active tracers: trend variables
+   USE trdmod          ! ocean active tracers trends 
+   USE trdmod_oce      ! ocean variables trends
    USE in_out_manager  ! I/O manager
    USE ldfslp          ! iso-neutral slopes available
@@ -61 +62 @@
       !! ** Action : - Update (ta,sa) arrays with the before geopotential 
       !!               biharmonic mixing trend.
-      !!             - Save the trends  in (ttrd,strd) ('key_diatrends')
+      !!             - Save the trends  in (ztdta,ztdsa) ('key_trdtra')
       !!
       !! History :
       !!   8.0  !  97-07  (G. Madec)  Original code
       !!   8.5  !  02-08  (G. Madec)  F90: Free form and module
+      !!   9.0  !  04-08  (C. Talandier) New trends organization
       !!----------------------------------------------------------------------
       !! * Modules used
@@ -76 +78 @@
       !! * Local declarations
       INTEGER ::   ji, jj, jk                 ! dummy loop indices
+      REAL(wp) ::   zta, zsa                  ! temporary scalars 
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   &
          wk3, wk4                ! work array used for rotated biharmonic
@@ -116 +119 @@
                ta(ji,jj,jk) = ta(ji,jj,jk) + wk3(ji,jj,jk)
                sa(ji,jj,jk) = sa(ji,jj,jk) + wk4(ji,jj,jk)
-#if defined key_trdtra || defined key_trdmld
-               ! save the horizontal diffusive trends
-               ttrd(ji,jj,jk,3) = wk3(ji,jj,jk)
-               strd(ji,jj,jk,3) = wk4(ji,jj,jk)
-#endif
             END DO
          END DO
@@ -126 +124 @@
       END DO                                           !   End of slab
       !                                                ! ===============
+
+      ! save the trends for diagnostic
+      ! save the horizontal diffusive trends
+      IF( l_trdtra )   THEN
+
+         CALL trd_mod(wk3, wk4, jpttdldf, 'TRA', kt)
+      ENDIF
+
+      IF(l_ctl) THEN         ! print mean trends (used for debugging)
+         zta = SUM( ta(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) )
+         zsa = SUM( sa(2:nictl,2:njctl,1:jpkm1) * tmask(2:nictl,2:njctl,1:jpkm1) )
+         WRITE(numout,*) ' ldf  - Ta: ', zta-t_ctl, ' Sa: ', zsa-s_ctl
+         t_ctl = zta   ;   s_ctl = zsa
+      ENDIF
+
    END SUBROUTINE tra_ldf_bilapg