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

Timestamp:

2004-04-22T14:38:42+02:00 (20 years ago)

Author:

opalod

Message:

CT : BUGFIX048 : Bug correction of Poleward Transport diagnostics

File:

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

trunk/NEMO/OPA_SRC/TRA/traldf_bilapg.F90

@@ -15 +15 @@
    USE oce             ! ocean dynamics and tracers variables
    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 ldftra_oce      ! ocean active tracers: lateral physics
+   USE trdtra_oce      ! ocean active tracers: trend variables
    USE in_out_manager  ! I/O manager
    USE ldfslp          ! iso-neutral slopes available
    USE lbclnk          ! ocean lateral boundary condition (or mpp link)
+   USE ptr             ! poleward transport diagnostics 
 
    IMPLICIT NONE
@@ -87 +88 @@
       ! and multiply by aht (output in (wk1,wk2) )
 
-      CALL ldfght ( tb, sb, wk1, wk2, 1 )
+      CALL ldfght ( kt, tb, sb, wk1, wk2, 1 )
 
 
@@ -98 +99 @@
       ! (output in (wk3,wk4) )
 
-      CALL ldfght ( wk1, wk2, wk3, wk4, 2 )
+      CALL ldfght ( kt, wk1, wk2, wk3, wk4, 2 )
 
 
@@ -124 +125 @@
 
 
-   SUBROUTINE ldfght ( pt, ps, plt, pls, kaht )
+   SUBROUTINE ldfght ( kt, pt, ps, plt, pls, kaht )
       !!----------------------------------------------------------------------
       !!                  ***  ROUTINE ldfght  ***
@@ -169 +170 @@
       !!----------------------------------------------------------------------
       !! * Arguments
+      INTEGER, INTENT( in ) ::   kt           ! ocean time-step index
       REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT( in  ) ::   &
          pt, ps           ! tracer fields (before t and s for 1st call
@@ -271 +273 @@
       !                                                ! ===============
 
+      IF( kaht == 2 .AND. MOD( kt, nf_ptr ) == 0 ) THEN
 #if defined key_diaptr
-      IF( kaht == 2 .AND. MOD( kt, nf_ptr ) == 0 ) THEN
          ! "zonal" mean diffusive heat and salt transport
-         pht_ldf(:,:) = prt_vj( zftv(:,:,:) )
-         pst_ldf(:,:) = prt_vj( zfsv(:,:,:) )
+!         pht_ldf(:) = prt_vj( zftv(:,:) )
+!         pst_ldf(:) = prt_vj( zfsv(:,:) )
+        write(numout,cform_err)
+        nstop = nstop + 1
+#endif
       ENDIF
-#endif