Changeset 377 for trunk/NEMO/TOP_SRC

Timestamp:

2006-01-09T17:03:02+01:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_bugfix018:RB+CE: save kz in a correct way for passive tracers

File:

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

trunk/NEMO/TOP_SRC/TRP/trczdf_iso.F90

@@ -161 +161 @@
       ENDIF
 
-      ! 0. Local constant initialization
+      ! 0.0  Local constant initialization
       ! --------------------------------
       IF( ln_trcadv_cen2 .OR. ln_trcadv_tvd ) THEN
@@ -173 +173 @@
          rdttrc(:) =  rdttra(:) * FLOAT(ndttrc)      
       ENDIF
+
+
+
+      ! 0.1 Save avs in zavs to recover avs in output files
+      !---------------------------------------------------
+      zavs(:,:,:) = fstravs(:,:,:)
+
 
 
@@ -262 +269 @@
                      &                     +wslpj(ji,jj,jk)*wslpj(ji,jj,jk) )
 
-                  ! save avs in zavs to recover avs in output files
-                  zavs(ji,jj,jk) = fstravs(ji,jj,jk)
                   ! add isopycnal vertical coeff. to avs
                   fstravs(ji,jj,jk) = fstravs(ji,jj,jk) + zavi
@@ -423 +428 @@
                END DO
             END DO
-
-
-            ! III. recover the avt (avs) resulting from vertical physics only
-            ! ===============================================================
-
-            DO jk = 2, jpkm1
-               DO ji = 2, jpim1
-                  fstravs(ji,jj,jk) = zavs(ji,jj,jk)
-               END DO
-            END DO
-
             !                                             ! ===============
          END DO                                           !   End of slab
@@ -440 +434 @@
       END DO
 
+
+
+      ! III. recover the avt (avs) resulting from vertical physics only
+      !---------------------------------------------------------------
+      fstravs(:,:,:) = zavs(:,:,:)
+
+      
       IF(ln_ctl)   THEN  ! print mean trends (used for debugging)
          WRITE(charout, FMT="('zdf - 2')")