- 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/P4Z/p4zprod.F90
r5836 r6060 54 54 REAL(wp) :: texcret2 !: 1 - excret2 55 55 56 !! * Substitutions57 # include "domzgr_substitute.h90"58 56 !!---------------------------------------------------------------------- 59 57 !! NEMO/TOP 3.3 , NEMO Consortium (2010) … … 277 275 DO jj = 1, jpj 278 276 DO ji = 1, jpi 279 IF( fsdepw(ji,jj,jk+1) <= hmld(ji,jj) ) THEN277 IF( gdepw_n(ji,jj,jk+1) <= hmld(ji,jj) ) THEN 280 278 zprbio(ji,jj,jk) = zprbio(ji,jj,jk) * zmixnano(ji,jj) 281 279 zprdia(ji,jj,jk) = zprdia(ji,jj,jk) * zmixdiat(ji,jj) … … 321 319 DO jj = 1, jpj 322 320 DO ji = 1, jpi 323 IF( fsdepw(ji,jj,jk+1) <= hmld(ji,jj) ) THEN321 IF( gdepw_n(ji,jj,jk+1) <= hmld(ji,jj) ) THEN 324 322 zprnch(ji,jj,jk) = zprnch(ji,jj,jk) * zmixnano(ji,jj) 325 323 zprdch(ji,jj,jk) = zprdch(ji,jj,jk) * zmixdiat(ji,jj) … … 462 460 zw2d(:,:) = 0. 463 461 DO jk = 1, jpkm1 464 zw2d(:,:) = zw2d(:,:) + zprorca (:,:,jk) * fse3t(:,:,jk) * zfact * tmask(:,:,jk) ! vert. integrated primary produc. by nano462 zw2d(:,:) = zw2d(:,:) + zprorca (:,:,jk) * e3t_n(:,:,jk) * zfact * tmask(:,:,jk) ! vert. integrated primary produc. by nano 465 463 ENDDO 466 464 CALL iom_put( "INTPPPHY" , zw2d ) … … 468 466 zw2d(:,:) = 0. 469 467 DO jk = 1, jpkm1 470 zw2d(:,:) = zw2d(:,:) + zprorcad(:,:,jk) * fse3t(:,:,jk) * zfact * tmask(:,:,jk) ! vert. integrated primary produc. by diatom468 zw2d(:,:) = zw2d(:,:) + zprorcad(:,:,jk) * e3t_n(:,:,jk) * zfact * tmask(:,:,jk) ! vert. integrated primary produc. by diatom 471 469 ENDDO 472 470 CALL iom_put( "INTPPPHY2" , zw2d ) … … 475 473 zw2d(:,:) = 0. 476 474 DO jk = 1, jpkm1 477 zw2d(:,:) = zw2d(:,:) + ( zprorca(:,:,jk) + zprorcad(:,:,jk) ) * fse3t(:,:,jk) * zfact * tmask(:,:,jk) ! vert. integrated pp475 zw2d(:,:) = zw2d(:,:) + ( zprorca(:,:,jk) + zprorcad(:,:,jk) ) * e3t_n(:,:,jk) * zfact * tmask(:,:,jk) ! vert. integrated pp 478 476 ENDDO 479 477 CALL iom_put( "INTPP" , zw2d ) … … 482 480 zw2d(:,:) = 0. 483 481 DO jk = 1, jpkm1 484 zw2d(:,:) = zw2d(:,:) + ( zpronew(:,:,jk) + zpronewd(:,:,jk) ) * fse3t(:,:,jk) * zfact * tmask(:,:,jk) ! vert. integrated new prod482 zw2d(:,:) = zw2d(:,:) + ( zpronew(:,:,jk) + zpronewd(:,:,jk) ) * e3t_n(:,:,jk) * zfact * tmask(:,:,jk) ! vert. integrated new prod 485 483 ENDDO 486 484 CALL iom_put( "INTPNEW" , zw2d ) … … 489 487 zw2d(:,:) = 0. 490 488 DO jk = 1, jpkm1 491 zw2d(:,:) = zw2d(:,:) + ( zprofen(:,:,jk) + zprofed(:,:,jk) ) * fse3t(:,:,jk) * zfact * tmask(:,:,jk) ! vert integr. bfe prod489 zw2d(:,:) = zw2d(:,:) + ( zprofen(:,:,jk) + zprofed(:,:,jk) ) * e3t_n(:,:,jk) * zfact * tmask(:,:,jk) ! vert integr. bfe prod 492 490 ENDDO 493 491 CALL iom_put( "INTPBFE" , zw2d ) … … 496 494 zw2d(:,:) = 0. 497 495 DO jk = 1, jpkm1 498 zw2d(:,:) = zw2d(:,:) + zprorcad(:,:,jk) * zysopt(:,:,jk) * fse3t(:,:,jk) * zfact * tmask(:,:,jk) ! vert integr. bsi prod496 zw2d(:,:) = zw2d(:,:) + zprorcad(:,:,jk) * zysopt(:,:,jk) * e3t_n(:,:,jk) * zfact * tmask(:,:,jk) ! vert integr. bsi prod 499 497 ENDDO 500 498 CALL iom_put( "INTPBSI" , zw2d )
Note: See TracChangeset
for help on using the changeset viewer.