- Timestamp:
- 2017-07-20T15:39:18+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zrem.F90
r6618 r8356 44 44 REAL(wp), PUBLIC :: xsiremlab !: fast remineralisation rate of POC 45 45 REAL(wp), PUBLIC :: xsilab !: fraction of labile biogenic silica 46 REAL(wp), PUBLIC :: oxymin !: halk saturation constant for anoxia47 48 46 49 47 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: denitr !: denitrification array … … 111 109 zdepprod(ji,jj,jk) = zdepmin**0.273 112 110 ENDIF 113 END DO114 END DO115 END DO116 117 DO jk = 1, jpkm1118 DO jj = 1, jpj119 DO ji = 1, jpi120 ! denitrification factor computed from O2 levels121 nitrfac(ji,jj,jk) = MAX( 0.e0, 0.4 * ( 6.e-6 - trb(ji,jj,jk,jpoxy) ) &122 & / ( oxymin + trb(ji,jj,jk,jpoxy) ) )123 nitrfac(ji,jj,jk) = MIN( 1., nitrfac(ji,jj,jk) )124 111 END DO 125 112 END DO … … 357 344 !! 358 345 !!---------------------------------------------------------------------- 359 NAMELIST/nampisrem/ xremik, xremip, nitrif, xsirem, xsiremlab, xsilab, & 360 & oxymin 346 NAMELIST/nampisrem/ xremik, xremip, nitrif, xsirem, xsiremlab, xsilab 361 347 INTEGER :: ios ! Local integer output status for namelist read 362 348 … … 380 366 WRITE(numout,*) ' fraction of labile biogenic silica xsilab =', xsilab 381 367 WRITE(numout,*) ' NH4 nitrification rate nitrif =', nitrif 382 WRITE(numout,*) ' halk saturation constant for anoxia oxymin =', oxymin383 368 ENDIF 384 369 ! 385 nitrfac (:,:,:) = 0._wp386 370 denitr (:,:,:) = 0._wp 387 371 denitnh4(:,:,:) = 0._wp
Note: See TracChangeset
for help on using the changeset viewer.