Changeset 5951 for branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90
- Timestamp:
- 2015-11-30T12:48:01+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4650_UKMO14.4_OBS_GENERAL_VINTERP/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zsms.F90
r5950 r5951 38 38 39 39 REAL(wp) :: alkbudget, no3budget, silbudget, ferbudget, po4budget 40 REAL(wp) :: xfact1, xfact2 40 REAL(wp) :: xfact1, xfact2, xfact3 41 41 INTEGER :: numco2, numnut, numnit !: logical unit for co2 budget 42 42 … … 45 45 46 46 47 !! * Substitutions48 # include "top_substitute.h90"49 47 !!---------------------------------------------------------------------- 50 48 !! NEMO/TOP 3.3 , NEMO Consortium (2010) … … 474 472 !!--------------------------------------------------------------------- 475 473 ! 476 INTEGER , INTENT( in ) :: kt ! ocean time-step index 477 REAL(wp) :: zfact 478 REAL(wp) :: zrdenittot, zsdenittot, znitrpottot 474 INTEGER, INTENT( in ) :: kt ! ocean time-step index 475 REAL(wp) :: zrdenittot, zsdenittot, znitrpottot 479 476 CHARACTER(LEN=100) :: cltxt 480 477 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zvol … … 492 489 xfact1 = rfact2r * 12. / 1.e15 * ryyss ! conversion molC/kt --> PgC/yr 493 490 xfact2 = 1.e+3 * rno3 * 14. / 1.e12 * ryyss ! conversion molC/l/s ----> TgN/m3/yr 491 xfact3 = 1.e+3 * rfact2r * rno3 ! conversion molC/l/kt ----> molN/m3/s 494 492 cltxt='time-step Alkalinity Nitrate Phosphorus Silicate Iron' 495 493 IF( lwp ) WRITE(numnut,*) TRIM(cltxt) … … 574 572 IF( iom_use( "Sdenit" ) .OR. ( ln_check_mass .AND. kt == nitend ) ) THEN 575 573 zsdenittot = glob_sum ( sdenit(:,:) * e1e2t(:,:) ) 576 CALL iom_put( "Sdenit", sdenit(:,:) * zfact* tmask(:,:,1) ) ! Nitrate reduction in the sediments574 CALL iom_put( "Sdenit", sdenit(:,:) * xfact3 * tmask(:,:,1) ) ! Nitrate reduction in the sediments 577 575 ENDIF 578 576
Note: See TracChangeset
for help on using the changeset viewer.