Changeset 12537 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p5zmicro.F90
- Timestamp:
- 2020-03-11T16:02:54+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p5zmicro.F90
r12349 r12537 97 97 ENDIF 98 98 ! 99 ! Use of excess carbon for metabolism 99 100 zmetexcess = 0.0 100 101 IF ( bmetexc ) zmetexcess = 1.0 … … 109 110 zproport = min(1.0, exp(-1.1 * MAX(0., ( sized(ji,jj,jk) - 1.8 ))**0.8 )) 110 111 zproport2 = sizen(ji,jj,jk)**(-0.54) 111 ! Michaelis-Menten mortality rates of microzooplankton 112 ! ----------------------------------------------------- 112 ! linear mortality of mesozooplankton 113 ! A michaelis menten modulation term is used to avoid extinction of 114 ! microzooplankton at very low food concentrations. Mortality is 115 ! enhanced in low O2 waters 116 ! ----------------------------------------------------------------- 113 117 zrespz = resrat * zfact * ( trb(ji,jj,jk,jpzoo) / ( xkmort + trb(ji,jj,jk,jpzoo) ) & 114 118 & + 3. * nitrfac(ji,jj,jk) ) 115 119 116 ! Zooplankton mortality. A square function has been selected with 117 ! no real reason except that it seems to be more stable and may mimic predation. 118 ! ------------------------------------------------------------------------------ 120 ! Zooplankton quadratic mortality. A square function has been selected with 121 ! to mimic predation and disease (density dependent mortality). It also tends 122 ! to stabilise the model 123 ! ------------------------------------------------------------------------- 119 124 ztortz = mzrat * 1.e6 * zfact * trb(ji,jj,jk,jpzoo) * (1. - nitrfac(ji,jj,jk)) 120 125 121 126 ! Computation of the abundance of the preys 122 127 ! A threshold can be specified in the namelist 123 ! -------------------------------------------- 128 ! Nanophyto and diatoms have a specific treatment with 129 ! teir preference decreasing with size. 130 ! -------------------------------------------------------- 124 131 zcompadi = zproport * MAX( ( trb(ji,jj,jk,jpdia) - xthreshdia ), 0.e0 ) 125 132 zcompaph = zproport2 * MAX( ( trb(ji,jj,jk,jpphy) - xthreshphy ), 0.e0 ) … … 128 135 zcompapoc = MAX( ( trb(ji,jj,jk,jppoc) - xthreshpoc ), 0.e0 ) 129 136 130 ! Microzooplankton grazing 131 ! ------------------------ 137 ! Microzooplankton grazing 138 ! The total amount of food is the sum of all preys accessible to mesozooplankton 139 ! multiplied by their food preference 140 ! A threshold can be specified in the namelist (xthresh). However, when food 141 ! concentration is close to this threshold, it is decreased to avoid the 142 ! accumulation of food in the mesozoopelagic domain 143 ! ------------------------------------------------------------------------------- 132 144 zfood = xprefn * zcompaph + xprefc * zcompapoc + xprefd * zcompadi & 133 145 & + xprefz * zcompaz + xprefp * zcompapi … … 136 148 zgraze = grazrat * xstep * tgfunc2(ji,jj,jk) * trb(ji,jj,jk,jpzoo) * (1. - nitrfac(ji,jj,jk)) 137 149 138 ! An active switching parameterization is used here. 139 ! We don't use the KTW parameterization proposed by 140 ! Vallina et al. because it tends to produce to steady biomass 141 ! composition and the variance of Chl is too low as it grazes 142 ! too strongly on winning organisms. Thus, instead of a square 143 ! a 1.5 power value is used which decreases the pressure on the 144 ! most abundant species 145 ! ------------------------------------------------------------ 150 ! An active switching parameterization is used here. 151 ! We don't use the KTW parameterization proposed by 152 ! Vallina et al. because it tends to produce too steady biomass 153 ! composition and the variance of Chl is too low as it grazes 154 ! too strongly on winning organisms. We use a generalized 155 ! switching parameterization proposed by Morozov and 156 ! Petrovskii (2013) 157 ! ------------------------------------------------------------ 158 ! The width of the selection window is increased when preys 159 ! have low abundance, .i.e. zooplankton become less specific 160 ! to avoid starvation. 161 ! ---------------------------------------------------------- 146 162 zsigma = 1.0 - zdenom**2/(0.05**2+zdenom**2) 147 163 zsigma = 0.5 + 1.0*zsigma … … 189 205 zgrazing(ji,jj,jk) = zgraztotc 190 206 191 ! 192 ! 207 ! Stoichiometruc ratios of the food ingested by zooplanton 208 ! -------------------------------------------------------- 193 209 zgrasratf = (zgraztotf + rtrn) / ( zgraztotc + rtrn ) 194 210 zgrasratn = (zgraztotn + rtrn) / ( zgraztotc + rtrn ) 195 211 zgrasratp = (zgraztotp + rtrn) / ( zgraztotc + rtrn ) 196 212 197 ! Growth efficiency is made a function of the quality 198 ! and the quantity of the preys 199 ! --------------------------------------------------- 213 ! Mesozooplankton efficiency. 214 ! We adopt a formulation proposed by Mitra et al. (2007) 215 ! The gross growth efficiency is controled by the most limiting nutrient. 216 ! Growth is also further decreased when the food quality is poor. This is currently 217 ! hard coded : it can be decreased by up to 50% (zepsherq) 218 ! GGE can also be decreased when food quantity is high, zepsherf (Montagnes and 219 ! Fulton, 2012) 220 ! ----------------------------------------------------------------------------------- 200 221 zepshert = MIN( 1., zgrasratn/ no3rat3, zgrasratp/ po4rat3, zgrasratf / ferat3) 201 222 zbeta = MAX( 0., (epsher - epshermin) ) … … 204 225 zepsherv = zepsherf * zepshert * zepsherq 205 226 206 ! 207 ! 208 ! 227 ! Respiration of microzooplankton 228 ! Excess carbon in the food is used preferentially 229 ! ------------------------------------------------ 209 230 zexcess = zgraztotc * zepsherf * (1.0 - zepshert) * zmetexcess 210 231 zbasresb = MAX(0., zrespz - zexcess) … … 212 233 zrespirc = srespir * zepsherv * zgraztotc + zbasresb 213 234 214 ! 215 ! 216 ! 235 ! When excess carbon is used, the other elements in excess 236 ! are also used proportionally to their abundance 237 ! -------------------------------------------------------- 217 238 zexcess = ( zgrasratn/ no3rat3 - zepshert ) / ( 1.0 - zepshert + rtrn) 218 239 zbasresn = zbasresi * zexcess * zgrasratn … … 222 243 zbasresf = zbasresi * zexcess * zgrasratf 223 244 224 ! 225 ! 245 ! Voiding of the excessive elements as DOM 246 ! ---------------------------------------- 226 247 zgradoct = (1. - unassc - zepsherv) * zgraztotc - zbasresi 227 248 zgradont = (1. - unassn) * zgraztotn - zepsherv * no3rat3 * zgraztotc - zbasresn … … 229 250 zgrareft = (1. - unassc) * zgraztotf - zepsherv * ferat3 * zgraztotc - zbasresf 230 251 231 ! 232 ! 233 ! 234 ! 252 ! Since only semilabile DOM is represented in PISCES 253 ! part of DOM is in fact labile and is then released 254 ! as dissolved inorganic compounds (ssigma) 255 ! -------------------------------------------------- 235 256 zgradoc = zgradoct * ssigma 236 257 zgradon = zgradont * ssigma … … 241 262 zgraref = zgrareft 242 263 243 ! 244 ! 264 ! Defecation as a result of non assimilated products 265 ! -------------------------------------------------- 245 266 zgrapoc = zgraztotc * unassc 246 267 zgrapon = zgraztotn * unassn … … 248 269 zgrapof = zgraztotf * unassc 249 270 250 ! 251 ! 271 ! Addition of respiration to the release of inorganic nutrients 272 ! ------------------------------------------------------------- 252 273 zgrarem = zgrarem + zbasresi + zrespirc 253 274 zgraren = zgraren + zbasresn + zrespirc * no3rat3 … … 255 276 zgraref = zgraref + zbasresf + zrespirc * ferat3 256 277 257 ! 258 ! 278 ! Update of the TRA arrays 279 ! ------------------------ 259 280 tra(ji,jj,jk,jppo4) = tra(ji,jj,jk,jppo4) + zgrarep 260 281 tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + zgraren … … 296 317 tra(ji,jj,jk,jpsfe) = tra(ji,jj,jk,jpsfe) + ferat3 * ztortz + zgrapof - zgrazpof 297 318 ! 298 ! calcite production 319 ! Calcite production 320 ! Calcite remineralization due to zooplankton activity 321 ! part of the ingested calcite is dissolving in the acidic gut 322 ! ------------------------------------------------------------- 299 323 zprcaca = xfracal(ji,jj,jk) * zgraznc 300 324 prodcal(ji,jj,jk) = prodcal(ji,jj,jk) + zprcaca ! prodcal=prodcal(nanophy)+prodcal(microzoo)+prodcal(mesozoo) … … 345 369 !! ** Purpose : Initialization of microzooplankton parameters 346 370 !! 347 !! ** Method : Read the namp iszoo namelist and check the parameters371 !! ** Method : Read the namp5zzoo namelist and check the parameters 348 372 !! called at the first timestep (nittrc000) 349 373 !! 350 !! ** input : Namelist namp iszoo374 !! ** input : Namelist namp5zzoo 351 375 !! 352 376 !!---------------------------------------------------------------------- … … 359 383 !!---------------------------------------------------------------------- 360 384 ! 361 REWIND( numnatp_ref ) ! Namelist namp iszoo in reference namelist : Pisces microzooplankton385 REWIND( numnatp_ref ) ! Namelist namp5zzoo in reference namelist : Pisces microzooplankton 362 386 READ ( numnatp_ref, namp5zzoo, IOSTAT = ios, ERR = 901) 363 387 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namp5zzoo in reference namelist' ) 364 388 ! 365 REWIND( numnatp_cfg ) ! Namelist namp iszoo in configuration namelist : Pisces microzooplankton389 REWIND( numnatp_cfg ) ! Namelist namp5zzoo in configuration namelist : Pisces microzooplankton 366 390 READ ( numnatp_cfg, namp5zzoo, IOSTAT = ios, ERR = 902 ) 367 391 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'namp5zzoo in configuration namelist' )
Note: See TracChangeset
for help on using the changeset viewer.