Changeset 3637

Timestamp:

2012-11-23T12:40:53+01:00 (11 years ago)

Author:

cetlod

Message:

bugfix in PISCES/p4zsink.F90 : replace trb with ztrb in p4z_sink2 subroutine, see ticket #1011

File:

• trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zsink.F90 (modified) (4 diffs)

trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zsink.F90

@@ -615 +615 @@
       INTEGER  ::   ji, jj, jk, jn
       REAL(wp) ::   zigma,zew,zign, zflx, zstep
-      REAL(wp), POINTER, DIMENSION(:,:,:) :: ztraz, zakz, zwsink2 
+      REAL(wp), POINTER, DIMENSION(:,:,:) :: ztraz, zakz, zwsink2, ztrb 
       !!---------------------------------------------------------------------
       !
@@ -621 +621 @@
       !
       ! Allocate temporary workspace
-      CALL wrk_alloc( jpi, jpj, jpk, ztraz, zakz, zwsink2 )
+      CALL wrk_alloc( jpi, jpj, jpk, ztraz, zakz, zwsink2, ztrb )
 
       zstep = rfact2 / 2.
@@ -627 +627 @@
       ztraz(:,:,:) = 0.e0
       zakz (:,:,:) = 0.e0
+      ztrb (:,:,:) = trn(:,:,:,jp_tra)
 
       DO jk = 1, jpkm1
@@ -696 +697 @@
             DO ji = 1, jpi
                zflx = ( psinkflx(ji,jj,jk) - psinkflx(ji,jj,jk+1) ) / fse3t(ji,jj,jk)
-               trb(ji,jj,jk,jp_tra) = trb(ji,jj,jk,jp_tra) + 2. * zflx
+               ztrb(ji,jj,jk,jp_tra) = ztrb(ji,jj,jk,jp_tra) + 2. * zflx
             END DO
          END DO
       END DO
 
-      trn     (:,:,:,jp_tra) = trb(:,:,:,jp_tra)
+      trn     (:,:,:,jp_tra) = ztrb(:,:,:,jp_tra)
       psinkflx(:,:,:)        = 2. * psinkflx(:,:,:)
       !
-      CALL wrk_dealloc( jpi, jpj, jpk, ztraz, zakz, zwsink2 )
+      CALL wrk_dealloc( jpi, jpj, jpk, ztraz, zakz, zwsink2, ztrb )
       !
       IF( nn_timing == 1 )  CALL timing_stop('p4z_sink2')