Changeset 3637


Ignore:
Timestamp:
2012-11-23T12:40:53+01:00 (8 years ago)
Author:
cetlod
Message:

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zsink.F90

    r3558 r3637  
    615615      INTEGER  ::   ji, jj, jk, jn 
    616616      REAL(wp) ::   zigma,zew,zign, zflx, zstep 
    617       REAL(wp), POINTER, DIMENSION(:,:,:) :: ztraz, zakz, zwsink2  
     617      REAL(wp), POINTER, DIMENSION(:,:,:) :: ztraz, zakz, zwsink2, ztrb  
    618618      !!--------------------------------------------------------------------- 
    619619      ! 
     
    621621      ! 
    622622      ! Allocate temporary workspace 
    623       CALL wrk_alloc( jpi, jpj, jpk, ztraz, zakz, zwsink2 ) 
     623      CALL wrk_alloc( jpi, jpj, jpk, ztraz, zakz, zwsink2, ztrb ) 
    624624 
    625625      zstep = rfact2 / 2. 
     
    627627      ztraz(:,:,:) = 0.e0 
    628628      zakz (:,:,:) = 0.e0 
     629      ztrb (:,:,:) = trn(:,:,:,jp_tra) 
    629630 
    630631      DO jk = 1, jpkm1 
     
    696697            DO ji = 1, jpi 
    697698               zflx = ( psinkflx(ji,jj,jk) - psinkflx(ji,jj,jk+1) ) / fse3t(ji,jj,jk) 
    698                trb(ji,jj,jk,jp_tra) = trb(ji,jj,jk,jp_tra) + 2. * zflx 
     699               ztrb(ji,jj,jk,jp_tra) = ztrb(ji,jj,jk,jp_tra) + 2. * zflx 
    699700            END DO 
    700701         END DO 
    701702      END DO 
    702703 
    703       trn     (:,:,:,jp_tra) = trb(:,:,:,jp_tra) 
     704      trn     (:,:,:,jp_tra) = ztrb(:,:,:,jp_tra) 
    704705      psinkflx(:,:,:)        = 2. * psinkflx(:,:,:) 
    705706      ! 
    706       CALL wrk_dealloc( jpi, jpj, jpk, ztraz, zakz, zwsink2 ) 
     707      CALL wrk_dealloc( jpi, jpj, jpk, ztraz, zakz, zwsink2, ztrb ) 
    707708      ! 
    708709      IF( nn_timing == 1 )  CALL timing_stop('p4z_sink2') 
Note: See TracChangeset for help on using the changeset viewer.