- 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/p4zmeso.F90
r9125 r9169 8 8 !! 3.4 ! 2011-06 (O. Aumont, C. Ethe) Quota model for iron 9 9 !!---------------------------------------------------------------------- 10 !! p4z_meso : 11 !! p4z_meso_init : 12 !!---------------------------------------------------------------------- 13 USE oce_trc ! 14 USE trc ! 15 USE sms_pisces ! 16 USE p4zprod ! 17 USE prtctl_trc ! 18 USE iom ! 10 !! p4z_meso : Compute the sources/sinks for mesozooplankton 11 !! p4z_meso_init : Initialization of the parameters for mesozooplankton 12 !!---------------------------------------------------------------------- 13 USE oce_trc ! shared variables between ocean and passive tracers 14 USE trc ! passive tracers common variables 15 USE sms_pisces ! PISCES Source Minus Sink variables 16 USE p4zprod ! production 17 USE prtctl_trc ! print control for debugging 18 USE iom ! I/O manager 19 19 20 20 IMPLICIT NONE … … 24 24 PUBLIC p4z_meso_init ! called in trcsms_pisces.F90 25 25 26 !! * Shared module variables27 26 REAL(wp), PUBLIC :: part2 !: part of calcite not dissolved in mesozoo guts 28 27 REAL(wp), PUBLIC :: xprefc !: mesozoo preference for POC … … 49 48 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 50 49 !!---------------------------------------------------------------------- 51 52 50 CONTAINS 53 51 … … 60 58 !! ** Method : - ??? 61 59 !!--------------------------------------------------------------------- 62 INTEGER, INTENT(in) :: kt, knt ! ocean time step 60 INTEGER, INTENT(in) :: kt, knt ! ocean time step and ??? 61 ! 63 62 INTEGER :: ji, jj, jk 64 63 REAL(wp) :: zcompadi, zcompaph, zcompapoc, zcompaz, zcompam … … 73 72 CHARACTER (len=25) :: charout 74 73 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zgrazing 75 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zw3d 76 74 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: zw3d 77 75 !!--------------------------------------------------------------------- 78 76 ! … … 122 120 123 121 ! Mesozooplankton flux feeding on GOC 124 ! ----------------------------------125 122 ! ---------------------------------- 126 123 zgrazffeg = grazflux * xstep * wsbio4(ji,jj,jk) & … … 253 250 !! 254 251 !! ** input : Namelist nampismes 255 !!256 252 !!---------------------------------------------------------------------- 257 253 INTEGER :: ios ! Local integer … … 262 258 !!---------------------------------------------------------------------- 263 259 ! 260 IF(lwp) THEN 261 WRITE(numout,*) 262 WRITE(numout,*) 'p4z_meso_init : Initialization of mesozooplankton parameters' 263 WRITE(numout,*) '~~~~~~~~~~~~~' 264 ENDIF 265 ! 264 266 REWIND( numnatp_ref ) ! Namelist nampismes in reference namelist : Pisces mesozooplankton 265 267 READ ( numnatp_ref, namp4zmes, IOSTAT = ios, ERR = 901) 266 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp4zmes in reference namelist', lwp ) 267 ! 268 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp4zmes in reference namelist', lwp ) 268 269 REWIND( numnatp_cfg ) ! Namelist nampismes in configuration namelist : Pisces mesozooplankton 269 270 READ ( numnatp_cfg, namp4zmes, IOSTAT = ios, ERR = 902 ) 270 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'namp4zmes in configuration namelist', lwp )271 IF(lwm) WRITE 271 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'namp4zmes in configuration namelist', lwp ) 272 IF(lwm) WRITE( numonp, namp4zmes ) 272 273 ! 273 274 IF(lwp) THEN ! control print 274 WRITE(numout,*) ' ' 275 WRITE(numout,*) ' Namelist parameters for mesozooplankton, namp4zmes' 276 WRITE(numout,*) ' ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~' 277 WRITE(numout,*) ' part of calcite not dissolved in mesozoo guts part2 =', part2 278 WRITE(numout,*) ' mesozoo preference for phyto xprefc =', xprefc 279 WRITE(numout,*) ' mesozoo preference for POC xprefp =', xprefp 280 WRITE(numout,*) ' mesozoo preference for zoo xprefz =', xprefz 281 WRITE(numout,*) ' mesozoo preference for poc xprefpoc =', xprefpoc 282 WRITE(numout,*) ' microzoo feeding threshold for mesozoo xthresh2zoo =', xthresh2zoo 283 WRITE(numout,*) ' diatoms feeding threshold for mesozoo xthresh2dia =', xthresh2dia 284 WRITE(numout,*) ' nanophyto feeding threshold for mesozoo xthresh2phy =', xthresh2phy 285 WRITE(numout,*) ' poc feeding threshold for mesozoo xthresh2poc =', xthresh2poc 286 WRITE(numout,*) ' feeding threshold for mesozooplankton xthresh2 =', xthresh2 287 WRITE(numout,*) ' exsudation rate of mesozooplankton resrat2 =', resrat2 288 WRITE(numout,*) ' mesozooplankton mortality rate mzrat2 =', mzrat2 289 WRITE(numout,*) ' maximal mesozoo grazing rate grazrat2 =', grazrat2 290 WRITE(numout,*) ' mesozoo flux feeding rate grazflux =', grazflux 291 WRITE(numout,*) ' non assimilated fraction of P by mesozoo unass2 =', unass2 292 WRITE(numout,*) ' Efficicency of Mesozoo growth epsher2 =', epsher2 293 WRITE(numout,*) ' Fraction of mesozoo excretion as DOM sigma2 =', sigma2 294 WRITE(numout,*) ' half sturation constant for grazing 2 xkgraz2 =', xkgraz2 275 WRITE(numout,*) ' Namelist : namp4zmes' 276 WRITE(numout,*) ' part of calcite not dissolved in mesozoo guts part2 =', part2 277 WRITE(numout,*) ' mesozoo preference for phyto xprefc =', xprefc 278 WRITE(numout,*) ' mesozoo preference for POC xprefp =', xprefp 279 WRITE(numout,*) ' mesozoo preference for zoo xprefz =', xprefz 280 WRITE(numout,*) ' mesozoo preference for poc xprefpoc =', xprefpoc 281 WRITE(numout,*) ' microzoo feeding threshold for mesozoo xthresh2zoo =', xthresh2zoo 282 WRITE(numout,*) ' diatoms feeding threshold for mesozoo xthresh2dia =', xthresh2dia 283 WRITE(numout,*) ' nanophyto feeding threshold for mesozoo xthresh2phy =', xthresh2phy 284 WRITE(numout,*) ' poc feeding threshold for mesozoo xthresh2poc =', xthresh2poc 285 WRITE(numout,*) ' feeding threshold for mesozooplankton xthresh2 =', xthresh2 286 WRITE(numout,*) ' exsudation rate of mesozooplankton resrat2 =', resrat2 287 WRITE(numout,*) ' mesozooplankton mortality rate mzrat2 =', mzrat2 288 WRITE(numout,*) ' maximal mesozoo grazing rate grazrat2 =', grazrat2 289 WRITE(numout,*) ' mesozoo flux feeding rate grazflux =', grazflux 290 WRITE(numout,*) ' non assimilated fraction of P by mesozoo unass2 =', unass2 291 WRITE(numout,*) ' Efficicency of Mesozoo growth epsher2 =', epsher2 292 WRITE(numout,*) ' Fraction of mesozoo excretion as DOM sigma2 =', sigma2 293 WRITE(numout,*) ' half sturation constant for grazing 2 xkgraz2 =', xkgraz2 295 294 ENDIF 296 295 !
Note: See TracChangeset
for help on using the changeset viewer.