Changeset 1353 for trunk/NEMO/TOP_SRC/TRP/trcldf_iso_zps.F90

Timestamp:

2009-03-31T10:15:02+02:00 (15 years ago)

Author:

cetlod

Message:

correction of calendar in TOP outputs, see ticket:368

File:

• trunk/NEMO/TOP_SRC/TRP/trcldf_iso_zps.F90 (modified) (3 diffs)

trunk/NEMO/TOP_SRC/TRP/trcldf_iso_zps.F90

@@ -81 +81 @@
       REAL(wp) ::   &
          zabe1, zabe2, zcof1, zcof2,   &  ! temporary scalars
-         zmsku, zmskv, zbtr, ztra,     &
-         ztagu, ztagv
+         zmsku, zmskv, zbtr, ztra
 
       REAL(wp), DIMENSION(jpi,jpj) ::   &
@@ -93 +92 @@
       REAL(wp), DIMENSION(jpi,jpj) ::   &
          zftug, zftvg                     ! temporary workspace
-      REAL(wp) ::   z_hdivn_x, z_hdivn_y
       REAL(wp) ::   &
-         zuwk, zvwk,                   &
-         zuwk1, zvwk1,                 &
-         zcg1,zcg2
+         z_hdivn_x, z_hdivn_y, zcg1, zcg2,  &
+         zuwk, zvwk, zuwk1, zvwk1
 #endif
       CHARACTER (len=22) :: charout
@@ -222 +219 @@
             DO jj = 2 , jpjm1
                DO ji = fs_2, fs_jpim1   ! vector opt.
-                  zbtr= 1. / ( e1t(ji,jj)*e2t(ji,jj)*fse3t(ji,jj,jk) )
+                  zbtr = 1. / ( e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) )
                   ztra = zbtr * ( zftu(ji,jj,jk) - zftu(ji-1,jj,jk) + zftv(ji,jj,jk) - zftv(ji,jj-1,jk)  )
-                  tra (ji,jj,jk,jn) = tra (ji,jj,jk,jn) + ztra
+                  tra(ji,jj,jk,jn) = tra (ji,jj,jk,jn) + ztra
+               END DO
+            END DO
+
 #if defined key_trc_diatrd
-                  IF (luttrd(jn)) trtrd (ji,jj,jk,ikeep(jn),4) = zbtr * ( zftu(ji,jj,jk) - zftu(ji-1,  jj,jk) )
-                  IF (luttrd(jn)) trtrd (ji,jj,jk,ikeep(jn),5) = zbtr * ( zftv(ji,jj,jk) - zftv(ji  ,jj-1,jk) )
-#endif
-               END DO
-            END DO
-#if defined key_trc_diatrd
-#   if defined key_trcldf_eiv
             DO jj = 2 , jpjm1
                DO ji = fs_2, fs_jpim1   ! vector opt.
-                  zbtr= 1. / ( e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) )
-                  ztagu = ( zftug(ji,jj) - zftug(ji-1,jj  ) ) * zbtr
-                  ztagv = ( zftvg(ji,jj) - zftvg(ji  ,jj-1) ) * zbtr
-                  IF (luttrd(jn)) trtrd (ji,jj,jk,ikeep(jn),4) = trtrd(ji,jj,jk,ikeep(jn),4) - ztagu
-                  IF (luttrd(jn)) trtrd (ji,jj,jk,ikeep(jn),5) = trtrd(ji,jj,jk,ikeep(jn),5) - ztagv
-               END DO
-            END DO
-#   endif
+                  zbtr = 1. / ( e1t(ji,jj) * e2t(ji,jj) * fse3t(ji,jj,jk) )
+                  IF( luttrd(jn) ) THEN 
+                     trtrd (ji,jj,jk,ikeep(jn),4) = zbtr * ( zftu(ji,jj,jk) - zftu(ji-1,  jj,jk) )
+                     trtrd (ji,jj,jk,ikeep(jn),5) = zbtr * ( zftv(ji,jj,jk) - zftv(ji  ,jj-1,jk) )
+                  ENDIF
+#   if defined key_trcldf_eiv
+                  IF( luttrd(jn) ) THEN 
+                     trtrd (ji,jj,jk,ikeep(jn),4) = trtrd(ji,jj,jk,ikeep(jn),4)    &
+                     &                            - zbtr * ( zftug(ji,jj) - zftug(ji-1,jj  ) ) 
+                     trtrd (ji,jj,jk,ikeep(jn),5) = trtrd(ji,jj,jk,ikeep(jn),5)    &
+                     &                            - zbtr * ( zftvg(ji,jj) - zftvg(ji  ,jj-1) )
+                  ENDIF
+#   endif
+               END DO
+            END DO
 #endif