- Timestamp:
- 2017-12-26T17:32:56+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_merge_2017/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zrem.F90
r9125 r9169 29 29 PUBLIC p4z_rem_alloc 30 30 31 !! * Shared module variables 32 REAL(wp), PUBLIC :: xremikc !: remineralisation rate of DOC 33 REAL(wp), PUBLIC :: xremikn !: remineralisation rate of DON 34 REAL(wp), PUBLIC :: xremikp !: remineralisation rate of DOP 35 REAL(wp), PUBLIC :: xremik !: remineralisation rate of POC 36 REAL(wp), PUBLIC :: nitrif !: NH4 nitrification rate 37 REAL(wp), PUBLIC :: xsirem !: remineralisation rate of POC 38 REAL(wp), PUBLIC :: xsiremlab !: fast remineralisation rate of POC 39 REAL(wp), PUBLIC :: xsilab !: fraction of labile biogenic silica 40 REAL(wp), PUBLIC :: feratb !: Fe/C quota in bacteria 41 REAL(wp), PUBLIC :: xkferb !: Half-saturation constant for bacteria Fe/C 42 43 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: denitr !: denitrification array 31 REAL(wp), PUBLIC :: xremikc !: remineralisation rate of DOC 32 REAL(wp), PUBLIC :: xremikn !: remineralisation rate of DON 33 REAL(wp), PUBLIC :: xremikp !: remineralisation rate of DOP 34 REAL(wp), PUBLIC :: xremik !: remineralisation rate of POC 35 REAL(wp), PUBLIC :: nitrif !: NH4 nitrification rate 36 REAL(wp), PUBLIC :: xsirem !: remineralisation rate of POC 37 REAL(wp), PUBLIC :: xsiremlab !: fast remineralisation rate of POC 38 REAL(wp), PUBLIC :: xsilab !: fraction of labile biogenic silica 39 REAL(wp), PUBLIC :: feratb !: Fe/C quota in bacteria 40 REAL(wp), PUBLIC :: xkferb !: Half-saturation constant for bacteria Fe/C 41 42 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: denitr !: denitrification array 44 43 45 44 !!---------------------------------------------------------------------- … … 303 302 INTEGER :: ios ! Local integer output status for namelist read 304 303 !!---------------------------------------------------------------------- 305 304 ! 305 IF(lwp) THEN 306 WRITE(numout,*) 307 WRITE(numout,*) 'p4z_rem_init : Initialization of remineralization parameters' 308 WRITE(numout,*) '~~~~~~~~~~~~' 309 ENDIF 310 ! 306 311 REWIND( numnatp_ref ) ! Namelist nampisrem in reference namelist : Pisces remineralization 307 312 READ ( numnatp_ref, nampisrem, IOSTAT = ios, ERR = 901) 308 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisrem in reference namelist', lwp ) 309 313 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'nampisrem in reference namelist', lwp ) 310 314 REWIND( numnatp_cfg ) ! Namelist nampisrem in configuration namelist : Pisces remineralization 311 315 READ ( numnatp_cfg, nampisrem, IOSTAT = ios, ERR = 902 ) 312 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nampisrem in configuration namelist', lwp )313 IF(lwm) WRITE 316 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'nampisrem in configuration namelist', lwp ) 317 IF(lwm) WRITE( numonp, nampisrem ) 314 318 315 319 IF(lwp) THEN ! control print 316 WRITE(numout,*) ' ' 317 WRITE(numout,*) ' Namelist parameters for remineralization, nampisrem' 318 WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 320 WRITE(numout,*) ' Namelist parameters for remineralization, nampisrem' 319 321 IF( ln_p4z ) THEN 320 WRITE(numout,*) ' remineralization rate of DOC xremik =', xremik322 WRITE(numout,*) ' remineralization rate of DOC xremik =', xremik 321 323 ELSE 322 WRITE(numout,*) ' remineralization rate of DOC xremikc =', xremikc323 WRITE(numout,*) ' remineralization rate of DON xremikn =', xremikn324 WRITE(numout,*) ' remineralization rate of DOP xremikp =', xremikp324 WRITE(numout,*) ' remineralization rate of DOC xremikc =', xremikc 325 WRITE(numout,*) ' remineralization rate of DON xremikn =', xremikn 326 WRITE(numout,*) ' remineralization rate of DOP xremikp =', xremikp 325 327 ENDIF 326 WRITE(numout,*) ' remineralization rate of Si xsirem =', xsirem327 WRITE(numout,*) ' fast remineralization rate of Si xsiremlab =', xsiremlab328 WRITE(numout,*) ' fraction of labile biogenic silica xsilab =', xsilab329 WRITE(numout,*) ' NH4 nitrification rate nitrif =', nitrif330 WRITE(numout,*) ' Bacterial Fe/C ratio feratb =', feratb331 WRITE(numout,*) ' Half-saturation constant for bact. Fe/C xkferb =', xkferb328 WRITE(numout,*) ' remineralization rate of Si xsirem =', xsirem 329 WRITE(numout,*) ' fast remineralization rate of Si xsiremlab =', xsiremlab 330 WRITE(numout,*) ' fraction of labile biogenic silica xsilab =', xsilab 331 WRITE(numout,*) ' NH4 nitrification rate nitrif =', nitrif 332 WRITE(numout,*) ' Bacterial Fe/C ratio feratb =', feratb 333 WRITE(numout,*) ' Half-saturation constant for bact. Fe/C xkferb =', xkferb 332 334 ENDIF 333 335 ! 334 denitr 336 denitr(:,:,:) = 0._wp 335 337 ! 336 338 END SUBROUTINE p4z_rem_init
Note: See TracChangeset
for help on using the changeset viewer.