Changeset 10368 for NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/PISCES/P4Z/p4zopt.F90
- Timestamp:
- 2018-12-03T12:45:01+01:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2018/dev_r10164_HPC09_ESIWACE_PREP_MERGE/src/TOP/PISCES/P4Z/p4zopt.F90
r10069 r10368 109 109 DO jk = 1, nksrp 110 110 etot_ndcy(:,:,jk) = ze1(:,:,jk) + ze2(:,:,jk) + ze3(:,:,jk) 111 enano (:,:,jk) = 2.1 * ze1(:,:,jk) + 0.42 * ze2(:,:,jk) + 0.4* ze3(:,:,jk)112 ediat (:,:,jk) = 1.6 * ze1(:,:,jk) + 0.69 * ze2(:,:,jk) + 0.7* ze3(:,:,jk)111 enano (:,:,jk) = 1.85 * ze1(:,:,jk) + 0.69 * ze2(:,:,jk) + 0.46 * ze3(:,:,jk) 112 ediat (:,:,jk) = 1.62 * ze1(:,:,jk) + 0.74 * ze2(:,:,jk) + 0.63 * ze3(:,:,jk) 113 113 END DO 114 114 IF( ln_p5z ) THEN 115 115 DO jk = 1, nksrp 116 epico (:,:,jk) = 2.1 * ze1(:,:,jk) + 0.42* ze2(:,:,jk) + 0.4 * ze3(:,:,jk)116 epico (:,:,jk) = 1.94 * ze1(:,:,jk) + 0.66 * ze2(:,:,jk) + 0.4 * ze3(:,:,jk) 117 117 END DO 118 118 ENDIF … … 134 134 DO jk = 1, nksrp 135 135 etot (:,:,jk) = ze1(:,:,jk) + ze2(:,:,jk) + ze3(:,:,jk) 136 enano(:,:,jk) = 2.1 * ze1(:,:,jk) + 0.42 * ze2(:,:,jk) + 0.4* ze3(:,:,jk)137 ediat(:,:,jk) = 1.6 * ze1(:,:,jk) + 0.69 * ze2(:,:,jk) + 0.7* ze3(:,:,jk)136 enano(:,:,jk) = 1.85 * ze1(:,:,jk) + 0.69 * ze2(:,:,jk) + 0.46 * ze3(:,:,jk) 137 ediat(:,:,jk) = 1.62 * ze1(:,:,jk) + 0.74 * ze2(:,:,jk) + 0.63 * ze3(:,:,jk) 138 138 END DO 139 139 IF( ln_p5z ) THEN 140 140 DO jk = 1, nksrp 141 epico(:,:,jk) = 2.1 * ze1(:,:,jk) + 0.42* ze2(:,:,jk) + 0.4 * ze3(:,:,jk)141 epico(:,:,jk) = 1.94 * ze1(:,:,jk) + 0.66 * ze2(:,:,jk) + 0.4 * ze3(:,:,jk) 142 142 END DO 143 143 ENDIF … … 182 182 zetmp1 (:,:) = 0.e0 183 183 zetmp2 (:,:) = 0.e0 184 zetmp3 (:,:) = 0.e0185 zetmp4 (:,:) = 0.e0186 184 187 185 DO jk = 1, nksrp … … 191 189 zetmp1 (ji,jj) = zetmp1 (ji,jj) + etot (ji,jj,jk) * e3t_n(ji,jj,jk) ! remineralisation 192 190 zetmp2 (ji,jj) = zetmp2 (ji,jj) + etot_ndcy(ji,jj,jk) * e3t_n(ji,jj,jk) ! production 193 zetmp3 (ji,jj) = zetmp3 (ji,jj) + enano (ji,jj,jk) * e3t_n(ji,jj,jk) ! production194 zetmp4 (ji,jj) = zetmp4 (ji,jj) + ediat (ji,jj,jk) * e3t_n(ji,jj,jk) ! production195 191 zdepmoy(ji,jj) = zdepmoy(ji,jj) + e3t_n(ji,jj,jk) 196 192 ENDIF … … 209 205 emoy (ji,jj,jk) = zetmp1(ji,jj) * z1_dep 210 206 zpar (ji,jj,jk) = zetmp2(ji,jj) * z1_dep 211 enano(ji,jj,jk) = zetmp3(ji,jj) * z1_dep 212 ediat(ji,jj,jk) = zetmp4(ji,jj) * z1_dep 207 ENDIF 208 END DO 209 END DO 210 END DO 211 ! 212 zdepmoy(:,:) = 0.e0 213 zetmp3 (:,:) = 0.e0 214 zetmp4 (:,:) = 0.e0 215 ! 216 DO jk = 1, nksrp 217 DO jj = 1, jpj 218 DO ji = 1, jpi 219 IF( gdepw_n(ji,jj,jk+1) <= MIN(hmld(ji,jj), heup_01(ji,jj)) ) THEN 220 zetmp3 (ji,jj) = zetmp3 (ji,jj) + enano (ji,jj,jk) * e3t_n(ji,jj,jk) ! production 221 zetmp4 (ji,jj) = zetmp4 (ji,jj) + ediat (ji,jj,jk) * e3t_n(ji,jj,jk) ! production 222 zdepmoy(ji,jj) = zdepmoy(ji,jj) + e3t_n(ji,jj,jk) 223 ENDIF 224 END DO 225 END DO 226 END DO 227 enanom(:,:,:) = enano(:,:,:) 228 ediatm(:,:,:) = ediat(:,:,:) 229 ! 230 DO jk = 1, nksrp 231 DO jj = 1, jpj 232 DO ji = 1, jpi 233 IF( gdepw_n(ji,jj,jk+1) <= hmld(ji,jj) ) THEN 234 z1_dep = 1. / ( zdepmoy(ji,jj) + rtrn ) 235 enanom(ji,jj,jk) = zetmp3(ji,jj) * z1_dep 236 ediatm(ji,jj,jk) = zetmp4(ji,jj) * z1_dep 213 237 ENDIF 214 238 END DO … … 221 245 DO jj = 1, jpj 222 246 DO ji = 1, jpi 223 IF( gdepw_n(ji,jj,jk+1) <= hmld(ji,jj) ) THEN 247 IF( gdepw_n(ji,jj,jk+1) <= MIN(hmld(ji,jj), heup_01(ji,jj)) ) THEN 248 zetmp5(ji,jj) = zetmp5 (ji,jj) + epico(ji,jj,jk) * e3t_n(ji,jj,jk) ! production 249 ENDIF 250 END DO 251 END DO 252 END DO 253 ! 254 epicom(:,:,:) = epico(:,:,:) 255 ! 256 DO jk = 1, nksrp 257 DO jj = 1, jpj 258 DO ji = 1, jpi 259 IF( gdepw_n(ji,jj,jk+1) <= hmld(ji,jj) ) THEN 224 260 z1_dep = 1. / ( zdepmoy(ji,jj) + rtrn ) 225 zetmp5(ji,jj) = zetmp5 (ji,jj) + epico(ji,jj,jk) * e3t_n(ji,jj,jk) ! production 226 epico(ji,jj,jk) = zetmp5(ji,jj) * z1_dep 261 epicom(ji,jj,jk) = zetmp5(ji,jj) * z1_dep 227 262 ENDIF 228 263 END DO
Note: See TracChangeset
for help on using the changeset viewer.