Changeset 198

Timestamp:

2004-12-22T16:48:32+01:00 (19 years ago)

Author:

opalod

Message:

CT : BUGFIX139 : set the temporary work array zavt and zavs to 3D instead of 2D to be able to save all the 3D array avt(,:,:) and avs(:,:,:)

File:

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

trunk/NEMO/OPA_SRC/TRA/trazdf_iso.F90

@@ -109 +109 @@
       !!   8.5  !  02-08  (G. Madec)  F90: Free form and module
       !!---------------------------------------------------------------------
+      !! * Modules used
+      USE oce                ,   &
+# if defined key_zdfddm
+         zavs => va,  &  ! use va as workspace
+# endif
+         zavt => ua      ! use ua as workspace
+
       !! * Arguments
       INTEGER, INTENT( in ) ::   kt           ! ocean time-step index
@@ -134 +141 @@
       REAL(wp), DIMENSION(jpi,jpk) ::   &
          ztfw, zdit, zdjt, zdj1t,   &
-         zsfw, zdis, zdjs, zdj1s,   &
-         zavt
+         zsfw, zdis, zdjs, zdj1s
 #if defined key_traldf_eiv   ||   defined key_esopa
       REAL(wp), DIMENSION(jpi,jpk) ::   &
@@ -143 +149 @@
          zuwk, zvwk,                    &
          zuwki, zvwki
-#endif
-#if defined key_zdfddm   ||   defined key_esopa
-      REAL(wp), DIMENSION(jpi,jpk) ::   &
-         zavs 
 #endif
       !!---------------------------------------------------------------------
@@ -270 +272 @@
 
                ! save avt in zavt to recover avt for mixed layer depth diag.
-               zavt(ji,jk) = avt(ji,jj,jk)
+               zavt(ji,jj,jk) = avt(ji,jj,jk)
                ! add isopycnal vertical coeff. to avt
                avt(ji,jj,jk) = avt(ji,jj,jk) + zavi
@@ -276 +278 @@
 #if defined key_zdfddm
                ! save avs in zavs to recover avs in output files
-               zavs(ji,jk) = fsavs(ji,jj,jk)
+               zavs(ji,jj,jk) = fsavs(ji,jj,jk)
                ! add isopycnal vertical coeff. to avs
                fsavs(ji,jj,jk) = fsavs(ji,jj,jk) + zavi
@@ -522 +524 @@
          DO jk = 2, jpkm1
             DO ji = 2, jpim1
-               avt(ji,jj,jk) = zavt(ji,jk)
+               avt(ji,jj,jk) = zavt(ji,jj,jk)
 #if defined key_zdfddm
-               fsavs(ji,jj,jk) = zavs(ji,jk)
+               fsavs(ji,jj,jk) = zavs(ji,jj,jk)
 #endif
             END DO