- Timestamp:
- 2015-12-16T10:25:22+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_merge_2015/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zexp.F90
r5836 r6060 42 42 43 43 !! * Substitutions 44 # include "domzgr_substitute.h90"45 44 # include "vectopt_loop_substitute.h90" 46 45 !!---------------------------------------------------------------------- … … 49 48 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 50 49 !!---------------------------------------------------------------------- 51 52 50 CONTAINS 53 51 … … 95 93 DO jj = 2, jpjm1 96 94 DO ji = fs_2, fs_jpim1 97 ze3t = 1. / fse3t(ji,jj,jk)95 ze3t = 1. / e3t_n(ji,jj,jk) 98 96 tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + ze3t * dmin3(ji,jj,jk) * xksi(ji,jj) 99 97 END DO … … 110 108 DO ji = fs_2, fs_jpim1 111 109 ikt = mbkt(ji,jj) 112 tra(ji,jj,ikt,jpno3) = tra(ji,jj,ikt,jpno3) + sedlam * sedpocn(ji,jj) / fse3t(ji,jj,ikt)110 tra(ji,jj,ikt,jpno3) = tra(ji,jj,ikt,jpno3) + sedlam * sedpocn(ji,jj) / e3t_n(ji,jj,ikt) 113 111 ! Deposition of organic matter in the sediment 114 112 zwork = vsed * trn(ji,jj,ikt,jpdet) … … 121 119 DO jj = 2, jpjm1 122 120 DO ji = fs_2, fs_jpim1 123 tra(ji,jj,1,jpno3) = tra(ji,jj,1,jpno3) + zgeolpoc * cmask(ji,jj) / areacot / fse3t(ji,jj,1)121 tra(ji,jj,1,jpno3) = tra(ji,jj,1,jpno3) + zgeolpoc * cmask(ji,jj) / areacot / e3t_n(ji,jj,1) 124 122 END DO 125 123 END DO … … 212 210 DO jj = 1, jpj 213 211 DO ji = 1, jpi 214 zfluo = ( fsdepw(ji,jj,jk ) / fsdepw(ji,jj,jpkb) )**xhr215 zfluu = ( fsdepw(ji,jj,jk+1) / fsdepw(ji,jj,jpkb) )**xhr212 zfluo = ( gdepw_n(ji,jj,jk ) / gdepw_n(ji,jj,jpkb) )**xhr 213 zfluu = ( gdepw_n(ji,jj,jk+1) / gdepw_n(ji,jj,jpkb) )**xhr 216 214 IF( zfluo.GT.1. ) zfluo = 1._wp 217 215 zdm0(ji,jj,jk) = zfluo - zfluu
Note: See TracChangeset
for help on using the changeset viewer.