Changeset 2968
- Timestamp:
- 2011-10-20T13:35:13+02:00 (13 years ago)
- Location:
- branches/2011/dev_r2787_PISCES_improvment/NEMOGCM/NEMO/TOP_SRC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2787_PISCES_improvment/NEMOGCM/NEMO/TOP_SRC/PISCES/p4zprod.F90
r2963 r2968 90 90 REAL(wp) :: zratio, zmax, zsilim, ztn, zadap 91 91 REAL(wp) :: zlim, zsilfac2, zsiborn, zprod, zetot2, zproreg, zproreg2 92 REAL(wp) :: zmxltst, zmxlday 92 REAL(wp) :: zmxltst, zmxlday, zmaxday 93 93 REAL(wp) :: zpislopen , zpislope2n 94 94 REAL(wp) :: zrum, zcodel, zargu, zval … … 295 295 zprbio(ji,jj,jk) = zprbio(ji,jj,jk) * zmixnano(ji,jj) 296 296 zprdia(ji,jj,jk) = zprdia(ji,jj,jk) * zmixdiat(ji,jj) 297 zprnch(ji,jj,jk) = zprnch(ji,jj,jk) * zmixnano(ji,jj)298 zprdch(ji,jj,jk) = zprdch(ji,jj,jk) * zmixdiat(ji,jj)299 297 ENDIF 300 298 END DO … … 302 300 END DO 303 301 304 IF( ln_newprod ) THEN 305 zprochln(:,:,:) = chlcmin * 12. * zprorca (:,:,:) 306 zprochld(:,:,:) = chlcmin * 12. * zprorcad(:,:,:) 307 ELSE 308 zprochln(:,:,:) = 0._wp 309 zprochld(:,:,:) = 0._wp 310 ENDIF 311 302 ! Computation of the various production terms 312 303 !CDIR NOVERRCHK 313 304 DO jk = 1, jpkm1 … … 316 307 !CDIR NOVERRCHK 317 308 DO ji = 1, jpi 318 319 309 IF( etot(ji,jj,jk) > 1.E-3 ) THEN 320 ! Computation of the various production terms for nanophyto. 321 zetot2 = enano(ji,jj,jk) * zstrn(ji,jj) 310 ! production terms for nanophyto. 322 311 zprorca(ji,jj,jk) = zprbio(ji,jj,jk) * xlimphy(ji,jj,jk) * trn(ji,jj,jk,jpphy) * rfact2 323 312 zpronew(ji,jj,jk) = zprorca(ji,jj,jk) * xnanono3(ji,jj,jk) / ( xnanono3(ji,jj,jk) + xnanonh4(ji,jj,jk) + rtrn ) 324 !325 zprod = rday * zprorca(ji,jj,jk) * zprnch(ji,jj,jk) * xlimphy(ji,jj,jk)326 zprochln(ji,jj,jk) = zprochln(ji,jj,jk) + chlcnm * 12. * zprod / ( zpislopead(ji,jj,jk) * zetot2 +rtrn)327 313 ! 328 314 zratio = trn(ji,jj,jk,jpnfe) / ( trn(ji,jj,jk,jpphy) + rtrn ) … … 333 319 & * trn(ji,jj,jk,jpfer) / ( trn(ji,jj,jk,jpfer) + concnfe(ji,jj,jk) ) & 334 320 & * zmax * trn(ji,jj,jk,jpphy) * rfact2 335 ENDIF 336 END DO 337 END DO 338 END DO 339 340 !CDIR NOVERRCHK 341 DO jk = 1, jpkm1 342 !CDIR NOVERRCHK 343 DO jj = 1, jpj 344 !CDIR NOVERRCHK 345 DO ji = 1, jpi 346 IF( etot(ji,jj,jk) > 1.E-3 ) THEN 347 ! Computation of the various production terms for diatoms 348 zetot2 = ediat(ji,jj,jk) * zstrn(ji,jj) 321 ! production terms for diatomees 349 322 zprorcad(ji,jj,jk) = zprdia(ji,jj,jk) * xlimdia(ji,jj,jk) * trn(ji,jj,jk,jpdia) * rfact2 350 323 zpronewd(ji,jj,jk) = zprorcad(ji,jj,jk) * xdiatno3(ji,jj,jk) / ( xdiatno3(ji,jj,jk) + xdiatnh4(ji,jj,jk) + rtrn ) 351 !352 zprod = rday * zprorcad(ji,jj,jk) * zprdch(ji,jj,jk) * xlimdia(ji,jj,jk)353 zprochld(ji,jj,jk) = zprochld(ji,jj,jk) + chlcdm * 12. * zprod / ( zpislopead2(ji,jj,jk) * zetot2 +rtrn )354 324 ! 355 325 zratio = trn(ji,jj,jk,jpdfe) / ( trn(ji,jj,jk,jpdia) + rtrn ) … … 364 334 END DO 365 335 END DO 336 337 IF( ln_newprod ) THEN 338 !CDIR NOVERRCHK 339 DO jk = 1, jpkm1 340 !CDIR NOVERRCHK 341 DO jj = 1, jpj 342 !CDIR NOVERRCHK 343 DO ji = 1, jpi 344 IF( fsdepw(ji,jj,jk+1) <= hmld(ji,jj) ) THEN 345 zprnch(ji,jj,jk) = zprnch(ji,jj,jk) * zmixnano(ji,jj) 346 zprdch(ji,jj,jk) = zprdch(ji,jj,jk) * zmixdiat(ji,jj) 347 ENDIF 348 IF( etot(ji,jj,jk) > 1.E-3 ) THEN 349 ! production terms for nanophyto. ( chlorophyll ) 350 znanotot = enano(ji,jj,jk) * zstrn(ji,jj) 351 zprod = rday * zprorca(ji,jj,jk) * zprnch(ji,jj,jk) * xlimphy(ji,jj,jk) 352 zprochln(ji,jj,jk) = chlcmin * 12. * zprorca (ji,jj,jk) 353 zprochln(ji,jj,jk) = zprochln(ji,jj,jk) + chlcnm * 12. * zprod / ( zpislopead(ji,jj,jk) * znanotot +rtrn) 354 ! production terms for diatomees ( chlorophyll ) 355 zdiattot = ediat(ji,jj,jk) * zstrn(ji,jj) 356 zprod = rday * zprorcad(ji,jj,jk) * zprdch(ji,jj,jk) * xlimdia(ji,jj,jk) 357 zprochld(ji,jj,jk) = chlcmin * 12. * zprorcad(ji,jj,jk) 358 zprochld(ji,jj,jk) = zprochld(ji,jj,jk) + chlcdm * 12. * zprod / ( zpislopead2(ji,jj,jk) * zdiattot +rtrn ) 359 ENDIF 360 END DO 361 END DO 362 END DO 363 ELSE 364 !CDIR NOVERRCHK 365 DO jk = 1, jpkm1 366 !CDIR NOVERRCHK 367 DO jj = 1, jpj 368 !CDIR NOVERRCHK 369 DO ji = 1, jpi 370 IF( etot(ji,jj,jk) > 1.E-3 ) THEN 371 ! production terms for nanophyto. ( chlorophyll ) 372 znanotot = enano(ji,jj,jk) * zstrn(ji,jj) 373 zprod = rday * zprorca(ji,jj,jk) * zprnch(ji,jj,jk) * trn(ji,jj,jk,jpphy) * xlimphy(ji,jj,jk) 374 zprochln(ji,jj,jk) = chlcnm * 144. * zprod / ( zpislopead(ji,jj,jk) * trn(ji,jj,jk,jpnch) * znanotot +rtrn) 375 ! production terms for diatomees ( chlorophyll ) 376 zdiattot = ediat(ji,jj,jk) * zstrn(ji,jj) 377 zprod = rday * zprorcad(ji,jj,jk) * zprdch(ji,jj,jk) * trn(ji,jj,jk,jpdia) * xlimdia(ji,jj,jk) 378 zprochld(ji,jj,jk) = chlcdm * 144. * zprod / ( zpislopead2(ji,jj,jk) * trn(ji,jj,jk,jpdch) * zdiattot +rtrn ) 379 ENDIF 380 END DO 381 END DO 382 END DO 383 ENDIF 366 384 367 385 ! Update the arrays TRA which contain the biological sources and sinks -
branches/2011/dev_r2787_PISCES_improvment/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r2963 r2968 230 230 ENDIF 231 231 ! Control of date 232 IF( nit000 - NINT( zkt ) /= 1.AND. nn_rsttr /= 0 ) &232 IF( nit000 - NINT( zkt ) /= nn_dttrc .AND. nn_rsttr /= 0 ) & 233 233 & CALL ctl_stop( ' ===>>>> : problem with nit000 for the restart', & 234 234 & ' verify the restart file or rerun with nn_rsttr = 0 (namelist)' )
Note: See TracChangeset
for help on using the changeset viewer.