Changeset 3637
- Timestamp:
- 2012-11-23T12:40:53+01:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zsink.F90
r3558 r3637 615 615 INTEGER :: ji, jj, jk, jn 616 616 REAL(wp) :: zigma,zew,zign, zflx, zstep 617 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztraz, zakz, zwsink2 617 REAL(wp), POINTER, DIMENSION(:,:,:) :: ztraz, zakz, zwsink2, ztrb 618 618 !!--------------------------------------------------------------------- 619 619 ! … … 621 621 ! 622 622 ! Allocate temporary workspace 623 CALL wrk_alloc( jpi, jpj, jpk, ztraz, zakz, zwsink2 )623 CALL wrk_alloc( jpi, jpj, jpk, ztraz, zakz, zwsink2, ztrb ) 624 624 625 625 zstep = rfact2 / 2. … … 627 627 ztraz(:,:,:) = 0.e0 628 628 zakz (:,:,:) = 0.e0 629 ztrb (:,:,:) = trn(:,:,:,jp_tra) 629 630 630 631 DO jk = 1, jpkm1 … … 696 697 DO ji = 1, jpi 697 698 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. * zflx699 ztrb(ji,jj,jk,jp_tra) = ztrb(ji,jj,jk,jp_tra) + 2. * zflx 699 700 END DO 700 701 END DO 701 702 END DO 702 703 703 trn (:,:,:,jp_tra) = trb(:,:,:,jp_tra)704 trn (:,:,:,jp_tra) = ztrb(:,:,:,jp_tra) 704 705 psinkflx(:,:,:) = 2. * psinkflx(:,:,:) 705 706 ! 706 CALL wrk_dealloc( jpi, jpj, jpk, ztraz, zakz, zwsink2 )707 CALL wrk_dealloc( jpi, jpj, jpk, ztraz, zakz, zwsink2, ztrb ) 707 708 ! 708 709 IF( nn_timing == 1 ) CALL timing_stop('p4z_sink2')
Note: See TracChangeset
for help on using the changeset viewer.