- Timestamp:
- 2019-11-22T15:29:17+01:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11943_MERGE_2019/src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src
- Property svn:mergeinfo deleted
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/PISCES/P2Z/p2zexp.F90
r10425 r11949 46 46 CONTAINS 47 47 48 SUBROUTINE p2z_exp( kt )48 SUBROUTINE p2z_exp( kt, Kmm, Krhs ) 49 49 !!--------------------------------------------------------------------- 50 50 !! *** ROUTINE p2z_exp *** … … 60 60 !!--------------------------------------------------------------------- 61 61 !! 62 INTEGER, INTENT( in ) :: kt ! ocean time-step index 62 INTEGER, INTENT( in ) :: kt ! ocean time-step index 63 INTEGER, INTENT( in ) :: Kmm, Krhs ! time level indices 63 64 !! 64 65 INTEGER :: ji, jj, jk, jl, ikt … … 70 71 IF( ln_timing ) CALL timing_start('p2z_exp') 71 72 ! 72 IF( kt == nittrc000 ) CALL p2z_exp_init 73 IF( kt == nittrc000 ) CALL p2z_exp_init( Kmm ) 73 74 74 75 zsedpoca(:,:) = 0. … … 83 84 DO jj = 2, jpjm1 84 85 DO ji = fs_2, fs_jpim1 85 ze3t = 1. / e3t _n(ji,jj,jk)86 tr a(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + ze3t * dmin3(ji,jj,jk) * xksi(ji,jj)86 ze3t = 1. / e3t(ji,jj,jk,Kmm) 87 tr(ji,jj,jk,jpno3,Krhs) = tr(ji,jj,jk,jpno3,Krhs) + ze3t * dmin3(ji,jj,jk) * xksi(ji,jj) 87 88 END DO 88 89 END DO … … 98 99 DO ji = fs_2, fs_jpim1 99 100 ikt = mbkt(ji,jj) 100 tr a(ji,jj,ikt,jpno3) = tra(ji,jj,ikt,jpno3) + sedlam * sedpocn(ji,jj) / e3t_n(ji,jj,ikt)101 tr(ji,jj,ikt,jpno3,Krhs) = tr(ji,jj,ikt,jpno3,Krhs) + sedlam * sedpocn(ji,jj) / e3t(ji,jj,ikt,Kmm) 101 102 ! Deposition of organic matter in the sediment 102 zwork = vsed * tr n(ji,jj,ikt,jpdet)103 zwork = vsed * tr(ji,jj,ikt,jpdet,Kmm) 103 104 zsedpoca(ji,jj) = ( zwork + dminl(ji,jj) * xksi(ji,jj) & 104 105 & - sedlam * sedpocn(ji,jj) - sedlostpoc * sedpocn(ji,jj) ) * rdt … … 109 110 DO jj = 2, jpjm1 110 111 DO ji = fs_2, fs_jpim1 111 tr a(ji,jj,1,jpno3) = tra(ji,jj,1,jpno3) + zgeolpoc * cmask(ji,jj) / areacot / e3t_n(ji,jj,1)112 tr(ji,jj,1,jpno3,Krhs) = tr(ji,jj,1,jpno3,Krhs) + zgeolpoc * cmask(ji,jj) / areacot / e3t(ji,jj,1,Kmm) 112 113 END DO 113 114 END DO … … 149 150 WRITE(charout, FMT="('exp')") 150 151 CALL prt_ctl_trc_info(charout) 151 CALL prt_ctl_trc(tab4d=tr a, mask=tmask, clinfo=ctrcnm)152 CALL prt_ctl_trc(tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm) 152 153 ENDIF 153 154 ! … … 157 158 158 159 159 SUBROUTINE p2z_exp_init 160 SUBROUTINE p2z_exp_init( Kmm ) 160 161 !!---------------------------------------------------------------------- 161 162 !! *** ROUTINE p4z_exp_init *** 162 163 !! ** purpose : specific initialisation for export 163 164 !!---------------------------------------------------------------------- 165 INTEGER, INTENT(in) :: Kmm ! time level index 164 166 INTEGER :: ji, jj, jk 165 167 REAL(wp) :: zmaskt, zfluo, zfluu … … 184 186 DO jj = 1, jpj 185 187 DO ji = 1, jpi 186 zfluo = ( gdepw _n(ji,jj,jk ) / gdepw_n(ji,jj,jpkb) )**xhr187 zfluu = ( gdepw _n(ji,jj,jk+1) / gdepw_n(ji,jj,jpkb) )**xhr188 zfluo = ( gdepw(ji,jj,jk ,Kmm) / gdepw(ji,jj,jpkb,Kmm) )**xhr 189 zfluu = ( gdepw(ji,jj,jk+1,Kmm) / gdepw(ji,jj,jpkb,Kmm) )**xhr 188 190 IF( zfluo.GT.1. ) zfluo = 1._wp 189 191 zdm0(ji,jj,jk) = zfluo - zfluu
Note: See TracChangeset
for help on using the changeset viewer.