Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

p4zmeso.F90 @ 14697

Last change on this file since 14697 was 14075, checked in by cguiavarch, 4 years ago
UKMO/NEMO_4.0.4_mirror : Remove SVN keywords.
File size: 18.5 KB

Rev	Line
[3443]	1	MODULE p4zmeso
	2	!!======================================================================
	3	!! * MODULE p4zmeso *
	4	!! TOP : PISCES Compute the sources/sinks for mesozooplankton
	5	!!======================================================================
	6	!! History : 1.0 ! 2002 (O. Aumont) Original code
	7	!! 2.0 ! 2007-12 (C. Ethe, G. Madec) F90
	8	!! 3.4 ! 2011-06 (O. Aumont, C. Ethe) Quota model for iron
	9	!!----------------------------------------------------------------------
[9169]	10	!! p4z_meso : Compute the sources/sinks for mesozooplankton
	11	!! p4z_meso_init : Initialization of the parameters for mesozooplankton
[3443]	12	!!----------------------------------------------------------------------
[9169]	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
[3443]	19
	20	IMPLICIT NONE
	21	PRIVATE
	22
	23	PUBLIC p4z_meso ! called in p4zbio.F90
	24	PUBLIC p4z_meso_init ! called in trcsms_pisces.F90
	25
[4147]	26	REAL(wp), PUBLIC :: part2 !: part of calcite not dissolved in mesozoo guts
[10362]	27	REAL(wp), PUBLIC :: xpref2d !: mesozoo preference for diatoms
	28	REAL(wp), PUBLIC :: xpref2n !: mesozoo preference for nanophyto
	29	REAL(wp), PUBLIC :: xpref2z !: mesozoo preference for microzooplankton
	30	REAL(wp), PUBLIC :: xpref2c !: mesozoo preference for POC
[4147]	31	REAL(wp), PUBLIC :: xthresh2zoo !: zoo feeding threshold for mesozooplankton
	32	REAL(wp), PUBLIC :: xthresh2dia !: diatoms feeding threshold for mesozooplankton
	33	REAL(wp), PUBLIC :: xthresh2phy !: nanophyto feeding threshold for mesozooplankton
	34	REAL(wp), PUBLIC :: xthresh2poc !: poc feeding threshold for mesozooplankton
	35	REAL(wp), PUBLIC :: xthresh2 !: feeding threshold for mesozooplankton
	36	REAL(wp), PUBLIC :: resrat2 !: exsudation rate of mesozooplankton
	37	REAL(wp), PUBLIC :: mzrat2 !: microzooplankton mortality rate
	38	REAL(wp), PUBLIC :: grazrat2 !: maximal mesozoo grazing rate
	39	REAL(wp), PUBLIC :: xkgraz2 !: non assimilated fraction of P by mesozoo
	40	REAL(wp), PUBLIC :: unass2 !: Efficicency of mesozoo growth
	41	REAL(wp), PUBLIC :: sigma2 !: Fraction of mesozoo excretion as DOM
[10362]	42	REAL(wp), PUBLIC :: epsher2 !: growth efficiency
	43	REAL(wp), PUBLIC :: epsher2min !: minimum growth efficiency at high food for grazing 2
[4147]	44	REAL(wp), PUBLIC :: grazflux !: mesozoo flux feeding rate
[3443]	45
	46	!!----------------------------------------------------------------------
[10067]	47	!! NEMO/TOP 4.0 , NEMO Consortium (2018)
[14075]	48	!! $Id$
[10068]	49	!! Software governed by the CeCILL license (see ./LICENSE)
[3443]	50	!!----------------------------------------------------------------------
	51	CONTAINS
	52
[5385]	53	SUBROUTINE p4z_meso( kt, knt )
[3443]	54	!!---------------------------------------------------------------------
	55	!! * ROUTINE p4z_meso *
	56	!!
	57	!! ** Purpose : Compute the sources/sinks for mesozooplankton
	58	!!
	59	!! ** Method : - ???
	60	!!---------------------------------------------------------------------
[9169]	61	INTEGER, INTENT(in) :: kt, knt ! ocean time step and ???
	62	!
[3443]	63	INTEGER :: ji, jj, jk
	64	REAL(wp) :: zcompadi, zcompaph, zcompapoc, zcompaz, zcompam
[4529]	65	REAL(wp) :: zgraze2 , zdenom, zdenom2
[10362]	66	REAL(wp) :: zfact , zfood, zfoodlim, zproport, zbeta
[7646]	67	REAL(wp) :: zmortzgoc, zfrac, zfracfe, zratio, zratio2, zfracal, zgrazcal
[12837]	68	REAL(wp) :: zepsherf, zepshert, zepsherv, zepsherq
	69	REAL(wp) :: zgrarsig, zgraztotc, zgraztotn, zgraztotf
[10362]	70	REAL(wp) :: zgrarem2, zgrafer2, zgrapoc2, zprcaca, zmortz, zgrasrat, zgrasratn
	71	REAL(wp) :: zrespz, ztortz, zgrazd, zgrazz, zgrazpof
[3443]	72	REAL(wp) :: zgrazn, zgrazpoc, zgraznf, zgrazf
[4148]	73	REAL(wp) :: zgrazfffp, zgrazfffg, zgrazffep, zgrazffeg
[12232]	74	REAL(wp), DIMENSION(jpi,jpj,jpk) :: zgrazing2, zfezoo2, zz2ligprod
[4148]	75	CHARACTER (len=25) :: charout
[3443]	76	!!---------------------------------------------------------------------
	77	!
[9124]	78	IF( ln_timing ) CALL timing_start('p4z_meso')
[3443]	79	!
	80	DO jk = 1, jpkm1
	81	DO jj = 1, jpj
	82	DO ji = 1, jpi
[5385]	83	zcompam = MAX( ( trb(ji,jj,jk,jpmes) - 1.e-9 ), 0.e0 )
[7646]	84	zfact = xstep * tgfunc2(ji,jj,jk) * zcompam
[3443]	85
	86	! Respiration rates of both zooplankton
	87	! -------------------------------------
[10362]	88	zrespz = resrat2 * zfact * ( trb(ji,jj,jk,jpmes) / ( xkmort + trb(ji,jj,jk,jpmes) ) &
	89	& + 3. * nitrfac(ji,jj,jk) )
[3443]	90
	91	! Zooplankton mortality. A square function has been selected with
	92	! no real reason except that it seems to be more stable and may mimic predation
	93	! ---------------------------------------------------------------
[10362]	94	ztortz = mzrat2 * 1.e6 * zfact * trb(ji,jj,jk,jpmes) * (1. - nitrfac(ji,jj,jk) )
[3443]	95	!
[5385]	96	zcompadi = MAX( ( trb(ji,jj,jk,jpdia) - xthresh2dia ), 0.e0 )
	97	zcompaz = MAX( ( trb(ji,jj,jk,jpzoo) - xthresh2zoo ), 0.e0 )
[10362]	98	zcompapoc = MAX( ( trb(ji,jj,jk,jppoc) - xthresh2poc ), 0.e0 )
[4529]	99	! Size effect of nanophytoplankton on grazing : the smaller it is, the less prone
	100	! it is to predation by mesozooplankton
	101	! -------------------------------------------------------------------------------
[5385]	102	zcompaph = MAX( ( trb(ji,jj,jk,jpphy) - xthresh2phy ), 0.e0 ) &
[4529]	103	& * MIN(1., MAX( 0., ( quotan(ji,jj,jk) - 0.2) / 0.3 ) )
[3443]	104
[10362]	105	! Mesozooplankton grazing
	106	! ------------------------
	107	zfood = xpref2d * zcompadi + xpref2z * zcompaz + xpref2n * zcompaph + xpref2c * zcompapoc
[4148]	108	zfoodlim = MAX( 0., zfood - MIN( 0.5 * zfood, xthresh2 ) )
[3443]	109	zdenom = zfoodlim / ( xkgraz2 + zfoodlim )
	110	zdenom2 = zdenom / ( zfood + rtrn )
[10222]	111	zgraze2 = grazrat2 * xstep * tgfunc2(ji,jj,jk) * trb(ji,jj,jk,jpmes) * (1. - nitrfac(ji,jj,jk))
[3443]	112
[10362]	113	zgrazd = zgraze2 * xpref2d * zcompadi * zdenom2
	114	zgrazz = zgraze2 * xpref2z * zcompaz * zdenom2
	115	zgrazn = zgraze2 * xpref2n * zcompaph * zdenom2
	116	zgrazpoc = zgraze2 * xpref2c * zcompapoc * zdenom2
[3443]	117
[5385]	118	zgraznf = zgrazn * trb(ji,jj,jk,jpnfe) / ( trb(ji,jj,jk,jpphy) + rtrn)
	119	zgrazf = zgrazd * trb(ji,jj,jk,jpdfe) / ( trb(ji,jj,jk,jpdia) + rtrn)
	120	zgrazpof = zgrazpoc * trb(ji,jj,jk,jpsfe) / ( trb(ji,jj,jk,jppoc) + rtrn)
[3443]	121
	122	! Mesozooplankton flux feeding on GOC
	123	! ----------------------------------
[7646]	124	zgrazffeg = grazflux * xstep * wsbio4(ji,jj,jk) &
[8533]	125	& * tgfunc2(ji,jj,jk) * trb(ji,jj,jk,jpgoc) * trb(ji,jj,jk,jpmes) &
	126	& * (1. - nitrfac(ji,jj,jk))
[5385]	127	zgrazfffg = zgrazffeg * trb(ji,jj,jk,jpbfe) / (trb(ji,jj,jk,jpgoc) + rtrn)
[7646]	128	zgrazffep = grazflux * xstep * wsbio3(ji,jj,jk) &
[8533]	129	& * tgfunc2(ji,jj,jk) * trb(ji,jj,jk,jppoc) * trb(ji,jj,jk,jpmes) &
	130	& * (1. - nitrfac(ji,jj,jk))
[5385]	131	zgrazfffp = zgrazffep * trb(ji,jj,jk,jpsfe) / (trb(ji,jj,jk,jppoc) + rtrn)
[10362]	132	!
	133	zgraztotc = zgrazd + zgrazz + zgrazn + zgrazpoc + zgrazffep + zgrazffeg
	134	! Compute the proportion of filter feeders
	135	zproport = (zgrazffep + zgrazffeg)/(rtrn + zgraztotc)
	136	! Compute fractionation of aggregates. It is assumed that
	137	! diatoms based aggregates are more prone to fractionation
	138	! since they are more porous (marine snow instead of fecal pellets)
	139	zratio = trb(ji,jj,jk,jpgsi) / ( trb(ji,jj,jk,jpgoc) + rtrn )
	140	zratio2 = zratio * zratio
	141	zfrac = zproport * grazflux * xstep * wsbio4(ji,jj,jk) &
[5385]	142	& * trb(ji,jj,jk,jpgoc) * trb(ji,jj,jk,jpmes) &
[4800]	143	& * ( 0.2 + 3.8 * zratio2 / ( 1.**2 + zratio2 ) )
[10362]	144	zfracfe = zfrac * trb(ji,jj,jk,jpbfe) / (trb(ji,jj,jk,jpgoc) + rtrn)
[4148]	145
[10362]	146	zgrazffep = zproport * zgrazffep
	147	zgrazffeg = zproport * zgrazffeg
	148	zgrazfffp = zproport * zgrazfffp
	149	zgrazfffg = zproport * zgrazfffg
	150	zgraztotc = zgrazd + zgrazz + zgrazn + zgrazpoc + zgrazffep + zgrazffeg
	151	zgraztotn = zgrazd * quotad(ji,jj,jk) + zgrazz + zgrazn * quotan(ji,jj,jk) &
	152	& + zgrazpoc + zgrazffep + zgrazffeg
	153	zgraztotf = zgrazf + zgraznf + zgrazz * ferat3 + zgrazpof + zgrazfffp + zgrazfffg
[4148]	154
[10362]	155	! Total grazing ( grazing by microzoo is already computed in p4zmicro )
[12232]	156	zgrazing2(ji,jj,jk) = zgraztotc
[3446]	157
[12837]	158	! Mesozooplankton efficiency.
	159	! We adopt a formulation proposed by Mitra et al. (2007)
	160	! The gross growth efficiency is controled by the most limiting nutrient.
	161	! Growth is also further decreased when the food quality is poor. This is currently
	162	! hard coded : it can be decreased by up to 50% (zepsherq)
	163	! GGE can also be decreased when food quantity is high, zepsherf (Montagnes and
	164	! Fulton, 2012)
	165	! -----------------------------------------------------------------------------------
[10362]	166	zgrasrat = ( zgraztotf + rtrn )/ ( zgraztotc + rtrn )
	167	zgrasratn = ( zgraztotn + rtrn )/ ( zgraztotc + rtrn )
[4529]	168	zepshert = MIN( 1., zgrasratn, zgrasrat / ferat3)
[10362]	169	zbeta = MAX(0., (epsher2 - epsher2min) )
	170	zepsherf = epsher2min + zbeta / ( 1.0 + 0.04E6 * 12. * zfood * zbeta )
[12837]	171	zepsherq = 0.5 + (1.0 - 0.5) * zepshert * ( 1.0 + 1.0 ) / ( zepshert + 1.0 )
	172	zepsherv = zepsherf * zepshert * zepsherq
[3443]	173
[10362]	174	zgrarem2 = zgraztotc * ( 1. - zepsherv - unass2 ) &
	175	& + ( 1. - epsher2 - unass2 ) / ( 1. - epsher2 ) * ztortz
	176	zgrafer2 = zgraztotc * MAX( 0. , ( 1. - unass2 ) * zgrasrat - ferat3 * zepsherv ) &
	177	& + ferat3 * ( ( 1. - epsher2 - unass2 ) /( 1. - epsher2 ) * ztortz )
	178	zgrapoc2 = zgraztotc * unass2
	179
[3443]	180	! Update the arrays TRA which contain the biological sources and sinks
	181	zgrarsig = zgrarem2 * sigma2
	182	tra(ji,jj,jk,jppo4) = tra(ji,jj,jk,jppo4) + zgrarsig
	183	tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + zgrarsig
	184	tra(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) + zgrarem2 - zgrarsig
[7646]	185	!
[10362]	186	IF( ln_ligand ) THEN
	187	tra(ji,jj,jk,jplgw) = tra(ji,jj,jk,jplgw) + (zgrarem2 - zgrarsig) * ldocz
	188	zz2ligprod(ji,jj,jk) = (zgrarem2 - zgrarsig) * ldocz
	189	ENDIF
[7646]	190	!
[3443]	191	tra(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy) - o2ut * zgrarsig
	192	tra(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) + zgrafer2
[10362]	193	zfezoo2(ji,jj,jk) = zgrafer2
[3443]	194	tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) + zgrarsig
	195	tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) + rno3 * zgrarsig
[4148]	196
[10362]	197	zmortz = ztortz + zrespz
	198	zmortzgoc = unass2 / ( 1. - epsher2 ) * ztortz + zrespz
	199	tra(ji,jj,jk,jpmes) = tra(ji,jj,jk,jpmes) - zmortz + zepsherv * zgraztotc
[3443]	200	tra(ji,jj,jk,jpdia) = tra(ji,jj,jk,jpdia) - zgrazd
	201	tra(ji,jj,jk,jpzoo) = tra(ji,jj,jk,jpzoo) - zgrazz
	202	tra(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) - zgrazn
[5385]	203	tra(ji,jj,jk,jpnch) = tra(ji,jj,jk,jpnch) - zgrazn * trb(ji,jj,jk,jpnch) / ( trb(ji,jj,jk,jpphy) + rtrn )
	204	tra(ji,jj,jk,jpdch) = tra(ji,jj,jk,jpdch) - zgrazd * trb(ji,jj,jk,jpdch) / ( trb(ji,jj,jk,jpdia) + rtrn )
	205	tra(ji,jj,jk,jpdsi) = tra(ji,jj,jk,jpdsi) - zgrazd * trb(ji,jj,jk,jpdsi) / ( trb(ji,jj,jk,jpdia) + rtrn )
	206	tra(ji,jj,jk,jpgsi) = tra(ji,jj,jk,jpgsi) + zgrazd * trb(ji,jj,jk,jpdsi) / ( trb(ji,jj,jk,jpdia) + rtrn )
[3443]	207	tra(ji,jj,jk,jpnfe) = tra(ji,jj,jk,jpnfe) - zgraznf
	208	tra(ji,jj,jk,jpdfe) = tra(ji,jj,jk,jpdfe) - zgrazf
	209
[10362]	210	tra(ji,jj,jk,jppoc) = tra(ji,jj,jk,jppoc) - zgrazpoc - zgrazffep + zfrac
	211	prodpoc(ji,jj,jk) = prodpoc(ji,jj,jk) + zfrac
	212	conspoc(ji,jj,jk) = conspoc(ji,jj,jk) - zgrazpoc - zgrazffep
	213	tra(ji,jj,jk,jpgoc) = tra(ji,jj,jk,jpgoc) + zmortzgoc - zgrazffeg + zgrapoc2 - zfrac
	214	prodgoc(ji,jj,jk) = prodgoc(ji,jj,jk) + zmortzgoc + zgrapoc2
	215	consgoc(ji,jj,jk) = consgoc(ji,jj,jk) - zgrazffeg - zfrac
	216	tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) - zgrazpof - zgrazfffp + zfracfe
	217	tra(ji,jj,jk,jpbfe) = tra(ji,jj,jk,jpbfe) + ferat3 * zmortzgoc - zgrazfffg &
[4529]	218	& + zgraztotf * unass2 - zfracfe
[10362]	219	zfracal = trb(ji,jj,jk,jpcal) / (trb(ji,jj,jk,jppoc) + trb(ji,jj,jk,jpgoc) + rtrn )
	220	zgrazcal = (zgrazffeg + zgrazpoc) * (1. - part2) * zfracal
	221	! calcite production
	222	zprcaca = xfracal(ji,jj,jk) * zgrazn
	223	prodcal(ji,jj,jk) = prodcal(ji,jj,jk) + zprcaca ! prodcal=prodcal(nanophy)+prodcal(microzoo)+prodcal(mesozoo)
	224	!
	225	zprcaca = part2 * zprcaca
	226	tra(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) + zgrazcal - zprcaca
	227	tra(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) - 2. * ( zgrazcal + zprcaca )
	228	tra(ji,jj,jk,jpcal) = tra(ji,jj,jk,jpcal) - zgrazcal + zprcaca
[3443]	229	END DO
	230	END DO
	231	END DO
	232	!
[5385]	233	IF( lk_iomput .AND. knt == nrdttrc ) THEN
[12276]	234	CALL iom_put( "PCAL" , prodcal(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) ) ! Calcite production
	235	IF( iom_use("GRAZ2") ) THEN ! Total grazing of phyto by zooplankton
	236	zgrazing2(:,:,jpk) = 0._wp ; CALL iom_put( "GRAZ2" , zgrazing2(:,:,:) * 1.e+3 * rfact2r * tmask(:,:,:) )
[4996]	237	ENDIF
[12276]	238	IF( iom_use("FEZOO2") ) THEN
	239	zfezoo2 (:,:,jpk) = 0._wp ; CALL iom_put( "FEZOO2", zfezoo2(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) )
	240	ENDIF
	241	IF( ln_ligand ) THEN
	242	zz2ligprod(:,:,jpk) = 0._wp ; CALL iom_put( "LPRODZ2", zz2ligprod(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) )
	243	ENDIF
[3443]	244	ENDIF
	245	!
	246	IF(ln_ctl) THEN ! print mean trends (used for debugging)
	247	WRITE(charout, FMT="('meso')")
	248	CALL prt_ctl_trc_info(charout)
	249	CALL prt_ctl_trc(tab4d=tra, mask=tmask, clinfo=ctrcnm)
	250	ENDIF
	251	!
[9124]	252	IF( ln_timing ) CALL timing_stop('p4z_meso')
[3443]	253	!
	254	END SUBROUTINE p4z_meso
	255
[9124]	256
[3443]	257	SUBROUTINE p4z_meso_init
	258	!!----------------------------------------------------------------------
	259	!! * ROUTINE p4z_meso_init *
	260	!!
	261	!! ** Purpose : Initialization of mesozooplankton parameters
	262	!!
	263	!! ** Method : Read the nampismes namelist and check the parameters
	264	!! called at the first timestep (nittrc000)
	265	!!
	266	!! ** input : Namelist nampismes
	267	!!----------------------------------------------------------------------
[9124]	268	INTEGER :: ios ! Local integer
	269	!
[10362]	270	NAMELIST/namp4zmes/ part2, grazrat2, resrat2, mzrat2, xpref2n, xpref2d, xpref2z, &
	271	& xpref2c, xthresh2dia, xthresh2phy, xthresh2zoo, xthresh2poc, &
	272	& xthresh2, xkgraz2, epsher2, epsher2min, sigma2, unass2, grazflux
[9124]	273	!!----------------------------------------------------------------------
	274	!
[9169]	275	IF(lwp) THEN
	276	WRITE(numout,*)
	277	WRITE(numout,*) 'p4z_meso_init : Initialization of mesozooplankton parameters'
	278	WRITE(numout,*) '~~~~~~~~~~~~~'
	279	ENDIF
	280	!
[12276]	281	REWIND( numnatp_ref )
[7646]	282	READ ( numnatp_ref, namp4zmes, IOSTAT = ios, ERR = 901)
[11536]	283	901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp4zmes in reference namelist' )
[12276]	284
	285	REWIND( numnatp_cfg )
[7646]	286	READ ( numnatp_cfg, namp4zmes, IOSTAT = ios, ERR = 902 )
[11536]	287	902 IF( ios > 0 ) CALL ctl_nam ( ios , 'namp4zmes in configuration namelist' )
[9169]	288	IF(lwm) WRITE( numonp, namp4zmes )
[9124]	289	!
[3443]	290	IF(lwp) THEN ! control print
[9169]	291	WRITE(numout,*) ' Namelist : namp4zmes'
	292	WRITE(numout,*) ' part of calcite not dissolved in mesozoo guts part2 =', part2
[10362]	293	WRITE(numout,*) ' mesozoo preference for phyto xpref2n =', xpref2n
	294	WRITE(numout,*) ' mesozoo preference for diatoms xpref2d =', xpref2d
	295	WRITE(numout,*) ' mesozoo preference for zoo xpref2z =', xpref2z
	296	WRITE(numout,*) ' mesozoo preference for poc xpref2c =', xpref2c
[9169]	297	WRITE(numout,*) ' microzoo feeding threshold for mesozoo xthresh2zoo =', xthresh2zoo
	298	WRITE(numout,*) ' diatoms feeding threshold for mesozoo xthresh2dia =', xthresh2dia
	299	WRITE(numout,*) ' nanophyto feeding threshold for mesozoo xthresh2phy =', xthresh2phy
	300	WRITE(numout,*) ' poc feeding threshold for mesozoo xthresh2poc =', xthresh2poc
	301	WRITE(numout,*) ' feeding threshold for mesozooplankton xthresh2 =', xthresh2
	302	WRITE(numout,*) ' exsudation rate of mesozooplankton resrat2 =', resrat2
	303	WRITE(numout,*) ' mesozooplankton mortality rate mzrat2 =', mzrat2
	304	WRITE(numout,*) ' maximal mesozoo grazing rate grazrat2 =', grazrat2
	305	WRITE(numout,*) ' mesozoo flux feeding rate grazflux =', grazflux
	306	WRITE(numout,*) ' non assimilated fraction of P by mesozoo unass2 =', unass2
[10362]	307	WRITE(numout,*) ' Efficiency of Mesozoo growth epsher2 =', epsher2
	308	WRITE(numout,*) ' Minimum Efficiency of Mesozoo growth epsher2min =', epsher2min
[9169]	309	WRITE(numout,*) ' Fraction of mesozoo excretion as DOM sigma2 =', sigma2
	310	WRITE(numout,*) ' half sturation constant for grazing 2 xkgraz2 =', xkgraz2
[3443]	311	ENDIF
[9124]	312	!
[3443]	313	END SUBROUTINE p4z_meso_init
	314
	315	!!======================================================================
[5656]	316	END MODULE p4zmeso

Note: See TracBrowser for help on using the repository browser.

New URL for NEMO forge! http://forge.nemo-ocean.eu

Context Navigation

source: NEMO/branches/UKMO/NEMO_4.0.4_mirror/src/TOP/PISCES/P4Z/p4zmeso.F90 @ 14697

Download in other formats: