File indexing completed on 2025-12-21 17:50:42 UTC
view on githubraw file Latest commit 1f645d6f on 2025-11-26 17:01:07 UTC
8fbfd1f382 Oliv*0001 #include "DARWIN_OPTIONS.h"
0002
0003 SUBROUTINE DARWIN_READ_TRAITPARAMS(iUnit,oUnit,myThid)
0004 IMPLICIT NONE
0005 #include "SIZE.h"
0006 #include "EEPARAMS.h"
0007 #include "PARAMS.h"
0008 #include "DARWIN_SIZE.h"
0009 #include "DARWIN_PARAMS.h"
0010 #include "DARWIN_TRAITPARAMS.h"
0011
0012 CHARACTER*(MAX_LEN_MBUF) msgBuf
0013 INTEGER iUnit,oUnit,myThid
0014
0015 #ifdef ALLOW_DARWIN
0016
0017 NAMELIST/DARWIN_RANDOM_PARAMS/
54fc1c2e36 Oliv*0018 & phymin,
8fbfd1f382 Oliv*0019 & Smallgrow,
0020 & Biggrow,
0021 & Smallgrowrange,
0022 & Biggrowrange,
0023 & diaz_growfac,
0024 & cocco_growfac,
0025 & diatom_growfac,
0026 & Smallmort,
0027 & Bigmort,
0028 & Smallmortrange,
0029 & Bigmortrange,
0030 & Smallexport,
0031 & Bigexport,
0032 & tempcoeff1,
0033 & tempcoeff2_small,
0034 & tempcoeff2_big,
0035 & tempcoeff3,
0036 & tempmax,
0037 & temprange,
0038 & tempdecay,
0039 & val_R_NC,
0040 & val_R_NC_diaz,
0041 & val_R_PC,
0042 & val_R_SiC_diatom,
0043 & val_R_FeC,
0044 & val_R_FeC_diaz,
0045 & val_R_PICPOC,
0046 & val_R_ChlC,
0047 & val_R_NC_zoo,
0048 & val_R_PC_zoo,
0049 & val_R_SiC_zoo,
0050 & val_R_FeC_zoo,
0051 & val_R_PICPOC_zoo,
0052 & val_R_ChlC_zoo,
0053 & SmallSink,
0054 & BigSink,
0055 & SmallPsat,
0056 & BigPsat,
0057 & ProcPsat,
0058 & UniDzPsat,
0059 & CoccoPsat,
0060 & SmallPsatrange,
0061 & BigPsatrange,
0062 & ProcPsatrange,
0063 & UniDzPsatrange,
0064 & CoccoPsatrange,
0065 & ksatNH4fac,
0066 & ksatNO2fac,
0067 & val_amminhib,
0068 & val_ksatsio2,
0069 & smallksatpar,
0070 & smallksatparstd,
0071 & smallkinhpar,
0072 & smallkinhparstd,
0073 & Bigksatpar,
0074 & Bigksatparstd,
0075 & Bigkinhpar,
0076 & Bigkinhparstd,
0077 & LLProkinhpar,
0078 & Coccokinhpar,
0079 & inhibcoef_geid_val,
0080 & smallmQyield,
0081 & smallmQyieldrange,
0082 & BigmQyield,
0083 & BigmQyieldrange,
0084 & smallchl2cmax,
0085 & smallchl2cmaxrange,
0086 & Bigchl2cmax,
0087 & Bigchl2cmaxrange,
0088 & val_acclimtimescl,
0089 & oldTwoGrazers,
0090 & GrazeFast,
0091 & GrazeSlow,
0092 & ZooexfacSmall,
0093 & ZooexfacBig,
0094 & ZoomortSmall,
0095 & ZoomortBig,
0096 & ZoomortSmall2,
0097 & ZoomortBig2,
0098 & ExGrazfracbig,
0099 & ExGrazfracsmall,
0100 & palathi,
0101 & palatlo,
0102 & diatomgraz,
0103 & coccograz,
0104 & olargegraz,
0105 & GrazeEfflow,
0106 & GrazeEffmod,
0107 & GrazeEffhi,
0108 & GrazeRate,
0109 & ExGrazfrac,
0110 & val_palat,
0111 & val_ass_eff,
0112 & kgrazesat_val,
0113 & Zoomort,
0114 & Zoomort2,
0115 & Zooexfac,
0116 & ZooDM
0117
0118 NAMELIST/DARWIN_TRAIT_PARAMS/
0119 & darwin_sort_biovol,
5910c7e5d2 Oliv*0120 & darwin_effective_ksat,
8fbfd1f382 Oliv*0121 & darwin_select_kn_allom,
0122 & logvolbase,
0123 & logvolinc,
0124 & biovol0,
0125 & biovolfac,
0126 & logvol0ind,
0127 & grp_logvolind,
0128 & grp_biovol,
0129 & grp_names,
0130 & grp_nplank,
0131 & grp_photo,
0132 & grp_bacttype,
0133 & grp_aerobic,
0134 & grp_denit,
0135 & grp_pred,
0136 & grp_prey,
0137 & grp_hasSi,
0138 & grp_hasPIC,
0139 & grp_diazo,
0140 & grp_useNH4,
0141 & grp_useNO2,
0142 & grp_useNO3,
0143 & grp_combNO,
0144 & grp_aptype,
0145 & grp_tempMort,
0146 & grp_tempMort2,
14aae40637 Oliv*0147 & grp_tempGraz,
8fbfd1f382 Oliv*0148 & a_Xmin,
0149 & a_R_NC,
0150 & a_R_PC,
0151 & a_R_SiC,
0152 & a_R_FeC,
0153 & a_R_ChlC,
0154 & a_R_PICPOC,
0155 & a_ExportFracMort,
0156 & a_ExportFracMort2,
0157 & a_ExportFracExude,
c7b6c66d45 Oliv*0158 & a_FracExudeC,
8fbfd1f382 Oliv*0159 & a_mort,
0160 & a_mort2,
0161 & a_phytoTempCoeff,
0162 & a_phytoTempExp1,
09b16dd232 Oliv*0163 & a_phytoTempAe,
0164 & a_phytoTempExp2,
8fbfd1f382 Oliv*0165 & a_phytoTempOptimum,
0166 & a_phytoDecayPower,
09b16dd232 Oliv*0167 & a_hetTempAe,
0168 & a_hetTempExp2,
0169 & a_hetTempOptimum,
0170 & a_hetDecayPower,
0171 & a_grazTempAe,
0172 & a_grazTempExp2,
0173 & a_grazTempOptimum,
0174 & a_grazDecayPower,
a092808e6b shlo*0175 #ifdef DARWIN_MACROMOLECULAR_GROWTH
0176 & a_Y_CP_Plip,
0177 & a_Y_CN_protein,
0178 & a_Y_NC_chl,
0179 & a_Y_CN_cyano,
0180 & a_Y_PN_nucacid,
0181 & a_Y_CN_DNA,
0182 & a_Y_CN_RNA,
0183 & a_Y_THY_P,
0184 & a_Y_FeN_photo,
0185 & a_ECo2Prod,
0186 & a_maintConsum,
0187 & a_VI_max,
0188 & b_VI_max,
0189 & a_A_I,
0190 & a_QC_other,
0191 & a_QC_pro_other,
0192 & a_QP_other,
0193 & a_QP_RNA_min,
0194 & a_QC_DNA,
0195 & a_QN_sto_max,
0196 & b_QN_sto_max,
0197 & a_Qp_max,
0198 & b_Qp_max,
0199 & a_Qfe_max,
0200 & b_Qfe_max,
0201 & a_A_pho,
0202 & a_A_bio,
0203 & a_AP_RNA,
0204 #endif
8fbfd1f382 Oliv*0205 #ifdef DARWIN_ALLOW_GEIDER
0206 & a_mQyield,
0207 & a_chl2cmax,
0208 & a_inhibGeider,
0209 #else
0210 & a_ksatPAR,
0211 & a_kinhPAR,
0212 #endif
0213 & a_amminhib,
0214 & a_acclimtimescl,
0215 & a_acclimtimescl_denom,
367178819a Oliv*0216 & a_ksatPON,
0217 & a_ksatDON,
a092808e6b shlo*0218 & a_hillnumDIN,
0219 & a_hillnumPO4,
0220 & a_hillnumFeT,
0221 & a_hillnumSiO2,
8fbfd1f382 Oliv*0222 & a_grazemax,
0223 & a_grazemax_denom,
0224 & b_grazemax,
0225 & a_kgrazesat,
0226 & b_kgrazesat,
0227 & a_biosink,
0228 & a_biosink_denom,
0229 & b_biosink,
0230 & a_bioswim,
0231 & a_bioswim_denom,
0232 & b_bioswim,
5e7acb36b1 daat*0233 & a_bioswimDVM,
0234 & b_bioswimDVM,
1d634e507f Oliv*0235 & a_ppSig,
0236 & a_ppOpt,
0237 & b_ppOpt,
faa67d1773 Oliv*0238 & palat_min,
8fbfd1f382 Oliv*0239 & a_PCmax,
0240 & a_PCmax_denom,
0241 & b_PCmax,
0242 & a_qcarbon,
0243 & b_qcarbon,
0244 & a_respRate_c,
0245 & a_respRate_c_denom,
0246 & b_respRate_c,
0247 & a_kexcC,
0248 & b_kexcC,
0249 & a_vmaxNO3,
0250 & a_vmaxNO3_denom,
0251 & b_vmaxNO3,
0252 & a_ksatNO3,
0253 & b_ksatNO3,
0254 & a_Qnmin,
0255 & b_Qnmin,
0256 & a_Qnmax,
0257 & b_Qnmax,
0258 & a_kexcN,
0259 & b_kexcN,
0260 & a_vmaxNO2,
0261 & a_vmaxNO2_denom,
0262 & b_vmaxNO2,
0263 & a_ksatNO2,
0264 & b_ksatNO2,
0265 & a_ksatNO2fac,
0266 & a_vmaxNH4,
0267 & a_vmaxNH4_denom,
0268 & b_vmaxNH4,
0269 & a_ksatNH4,
0270 & b_ksatNH4,
0271 & a_ksatNH4fac,
0272 & a_vmaxN,
0273 & a_vmaxN_denom,
0274 & b_vmaxN,
0275 & a_vmaxPO4,
0276 & a_vmaxPO4_denom,
0277 & b_vmaxPO4,
0278 & a_ksatPO4,
0279 & b_ksatPO4,
0280 & a_Qpmin,
0281 & b_Qpmin,
0282 & a_Qpmax,
0283 & b_Qpmax,
0284 & a_kexcP,
0285 & b_kexcP,
0286 & a_vmaxSiO2,
0287 & a_vmaxSiO2_denom,
0288 & b_vmaxSiO2,
0289 & a_ksatSiO2,
0290 & b_ksatSiO2,
0291 & a_Qsimin,
0292 & b_Qsimin,
0293 & a_Qsimax,
0294 & b_Qsimax,
0295 & a_kexcSi,
0296 & b_kexcSi,
0297 & a_vmaxFeT,
0298 & a_vmaxFeT_denom,
0299 & b_vmaxFeT,
0300 & a_ksatFeT,
0301 & b_ksatFeT,
0302 & a_Qfemin,
0303 & b_Qfemin,
0304 & a_Qfemax,
0305 & b_Qfemax,
0306 & a_kexcFe,
0307 & b_kexcFe,
0308 & grp_ExportFracPreyPred,
5e7acb36b1 daat*0309 & grp_ass_eff,
0310 & a_PARpref,
0311 & b_PARpref,
0312 & a_mortmaxDVM,
0313 & b_mortmaxDVM,
0314 & a_ksatDVM,
0315 & b_ksatDVM,
0316 & a_ksatPARDVM,
0317 & b_ksatPARDVM,
0318 & a_fracPARmort,
0319 & b_fracPARmort,
0320 & a_ExportFracDVM
8fbfd1f382 Oliv*0321
0322
0323
0324 _RL pday
0325
0326
0327
0328
0329
0330 pday = 86400.0 _d 0
0331
0332
0333
0334
0335
54fc1c2e36 Oliv*0336 phymin = 120 _d -20
8fbfd1f382 Oliv*0337 Smallgrow = .7 _d 0
0338 Biggrow = .4 _d 0
0339 Smallgrowrange = 0. _d 0
0340 Biggrowrange = 0. _d 0
0341 diaz_growfac = 2. _d 0
0342 cocco_growfac = 1.3 _d 0
0343 diatom_growfac = 0.95 _d 0
0344 Smallmort = 10. _d 0
0345 Bigmort = 10. _d 0
0346 Smallmortrange = 0. _d 0
0347 Bigmortrange = 0. _d 0
0348 Smallexport = 0.2 _d 0
0349 Bigexport = 0.5 _d 0
0350 tempcoeff1 = 1. _d 0/3. _d 0
0351 tempcoeff2_small = 0.001 _d 0
0352 tempcoeff2_big = 0.0003 _d 0
0353 tempcoeff3 = 1.04 _d 0
0354 tempmax = 30. _d 0
0355 temprange = 32. _d 0
0356 tempdecay = 4. _d 0
0357 val_R_NC = 16.0 _d 0/120.0 _d 0
0358 val_R_NC_diaz = 40.0 _d 0/120.0 _d 0
0359 val_R_PC = 1.0 _d 0/120.0 _d 0
0360 val_R_SiC_diatom = 16.0 _d 0/120.0 _d 0
0361 val_R_FeC = 1.0 _d -3/120.0 _d 0
0362 val_R_FeC_diaz = 30.0 _d 0 * val_R_FeC
0363 val_R_PICPOC = 0.8 _d 0
0364 val_R_ChlC = 16.0 _d 0/120 _d 0
0365 val_R_NC_zoo = 16.0 _d 0/120.0 _d 0
0366 val_R_PC_zoo = 1 _d 0/120.0 _d 0
0367 val_R_SiC_zoo = 0.0 _d 0
0368 val_R_FeC_zoo = 1.0 _d -3/120.0 _d 0
0369 val_R_PICPOC_zoo = 0.0 _d 0
0370 val_R_ChlC_zoo = 0 _d 0/120 _d 0
0371 SmallSink = 0.0 _d 0/pday
0372 BigSink = 0.5 _d 0/pday
0373 SmallPsat = 0.015 _d 0
0374 BigPsat = 0.035 _d 0
0375 ProcPsat = 0.01 _d 0
0376 UniDzPsat = 0.012 _d 0
0377 CoccoPsat = 0.035 _d 0
0378 SmallPsatrange = 0.02 _d 0
0379 BigPsatrange = 0.02 _d 0
0380 ProcPsatrange = 0.005 _d 0
0381 UniDzPsatrange = 0.02 _d 0
0382 CoccoPsatrange = 0.02 _d 0
0383 ksatNH4fac = .50 _d 0
0384 ksatNO2fac = 1.0 _d 0
0385 val_amminhib = 4.6 _d 0
0386 val_ksatsio2 = 1. _d 0
0387 smallksatpar = 0.12 _d -1
0388 smallksatparstd = 0.20 _d -1
0389 smallkinhpar = 6.0 _d -3
0390 smallkinhparstd = 0.10 _d -3
0391 Bigksatpar = 0.12 _d -1
0392 Bigksatparstd = 0.06 _d -1
0393 Bigkinhpar = 1.0 _d -3
0394 Bigkinhparstd = 0.05 _d -3
0395 LLProkinhpar = 6.0 _d -3
0396 Coccokinhpar = 0.5 _d -3
0397 inhibcoef_geid_val = 0 _d 0
0398 smallmQyield = 5. _d -5
0399 smallmQyieldrange = 1. _d -4
0400 BigmQyield = 3. _d -5
0401 BigmQyieldrange = 4. _d -5
0402 smallchl2cmax = 0.2 _d 0
0403 smallchl2cmaxrange = 0.3 _d 0
0404 Bigchl2cmax = 0.5 _d 0
0405 Bigchl2cmaxrange = 0.3 _d 0
0406 val_acclimtimescl = 1./(60. _d 0*60. _d 0*24. _d 0*20. _d 0)
0407 oldTwoGrazers = .FALSE.
0408 GrazeFast = 1.0 _d 0/(2.0 _d 0*pday)
0409 GrazeSlow = 1.0 _d 0/(7.0 _d 0*pday)
0410 ZooexfacSmall = 0.2 _d 0
0411 ZooexfacBig = 0.7 _d 0
0412 ZoomortSmall = 1.0 _d 0/(30.0 _d 0*pday)
0413 ZoomortBig = 1.0 _d 0/(30.0 _d 0*pday)
0414 ZoomortSmall2 = 0. _d 0
0415 ZoomortBig2 = 0. _d 0
0416 ExGrazfracbig = 0.8 _d 0
0417 ExGrazfracsmall = 0.8 _d 0
0418 palathi = 1.0 _d 0
0419 palatlo = 0.2 _d 0
0420 diatomgraz = 0.7 _d 0
0421 coccograz = 0.6 _d 0
0422 olargegraz = 1.0 _d 0
0423 GrazeEfflow = 0.2 _d 0
0424 GrazeEffmod = 0.5 _d 0
0425 GrazeEffhi = 0.7 _d 0
0426 GrazeRate = 1.0 _d 0/(2.0 _d 0*pday)
0427 ExGrazfrac = 0.8 _d 0
0428 val_palat = 0.0 _d 0
0429 val_ass_eff = 0.70 _d 0
0430 kgrazesat_val = 12 _d 0
0431 Zoomort = 1.0 _d 0/(30.0 _d 0*pday)
0432 Zoomort2 = 0. _d 0
0433 Zooexfac = 0.7 _d 0
0434 ZooDM = 100 _d 0
0435
0436 READ(UNIT=iUnit, NML=DARWIN_RANDOM_PARAMS)
0437
0438
0439
0440
0441
0442 darwin_sort_biovol = .FALSE.
5910c7e5d2 Oliv*0443 darwin_effective_ksat = .FALSE.
8fbfd1f382 Oliv*0444 darwin_select_kn_allom = 2
0445 logvolbase = 0.0 _d 0
0446 logvolinc = 0.0 _d 0
0447 biovol0 = 0.0 _d 0
0448 biovolfac = 1.0 _d 0
0449 logvol0ind = 0
0450 grp_logvolind = 0 _d 0
0451 grp_biovol = 0 _d 0
0452 grp_names = ''
0453 grp_nplank = 0
0454 grp_photo = 1
0455 grp_bacttype = 0
0456 grp_aerobic = 0
0457 grp_denit = 0
0458 grp_pred = 0
0459 grp_prey = 1
0460 grp_hasSi = 0
0461 grp_hasPIC = 0
0462 grp_diazo = 0
0463 grp_useNH4 = 1
0464 grp_useNO2 = 1
0465 grp_useNO3 = 1
0466 grp_combNO = 1
0467 grp_aptype = 0
0468 grp_tempMort = 1
0469 grp_tempMort2 = 1
14aae40637 Oliv*0470 grp_tempGraz = 1
5910c7e5d2 Oliv*0471
8fbfd1f382 Oliv*0472 a_Xmin = 0 _d 0
0473 a_R_NC = 16.0 _d 0/120.0 _d 0
0474 a_R_PC = 1 _d 0/120.0 _d 0
0475 a_R_SiC = 0 _d 0
0476 a_R_FeC = 1.0 _d -3/120.0 _d 0
0477 a_R_ChlC = 16 _d 0/120 _d 0
0478 a_R_PICPOC = 0.8 _d 0
0479 a_ExportFracMort = 0.5 _d 0
0480 a_ExportFracMort2 = 0.5 _d 0
0481 a_ExportFracExude = DARWIN_UNINIT_RL
c7b6c66d45 Oliv*0482 a_FracExudeC = 0.3 _d 0
8fbfd1f382 Oliv*0483 a_mort = 0.02 _d 0 / pday
0484 a_mort2 = 0.0 _d 0
0485 a_phytoTempCoeff = 1. _d 0/3. _d 0
0486 a_phytoTempExp1 = 1.04 _d 0
09b16dd232 Oliv*0487 a_phytoTempAe = 0.0438 _d 0
0488 a_phytoTempExp2 = 0.001 _d 0
8fbfd1f382 Oliv*0489 a_phytoTempOptimum = 2. _d 0
0490 a_phytoDecayPower = 4. _d 0
09b16dd232 Oliv*0491 a_hetTempAe = 0.0438 _d 0
0492 a_hetTempExp2 = 0.001 _d 0
0493 a_hetTempOptimum = 2. _d 0
0494 a_hetDecayPower = 4. _d 0
0495 a_grazTempAe = 0.0438 _d 0
0496 a_grazTempExp2 = 0.001 _d 0
0497 a_grazTempOptimum = 2. _d 0
0498 a_grazDecayPower = 4. _d 0
a092808e6b shlo*0499 #ifdef DARWIN_MACROMOLECULAR_GROWTH
0500 a_A_pho = 16.0 _d 0
0501 a_Y_THY_P = 0.028163 _d 0
0502 a_Y_CN_protein = 5.3 _d 0 /1.4 _d 0
0503 a_Y_FeN_photo = 0.00163 _d 0
0504 a_Y_NC_chl = 4.0 _d 0/55.0 _d 0
0505 a_Y_CN_cyano = 2.0 _d 0
0506 a_Y_PN_nucacid = 1. _d 0 /3.75 _d 0
0507 a_Y_CN_DNA = 9.75 _d 0 /3.75 _d 0
0508 a_Y_CN_RNA = 9.50 _d 0 /3.75 _d 0
0509 a_Y_CP_Plip = 40.0 _d 0
0510 a_ECo2Prod = 0.774 _d 0
0511 a_maintConsum = 0.393 _d 0/3600/24
0512 a_VI_max = 277 _d 0/3600/24
0513 b_VI_max = 0 _d 0
0514 a_A_I = 0.008633641 _d 0
0515 a_QC_other = 0.0182 _d 0
0516 a_QC_pro_other = 0.24 _d 0
0517 a_QP_other = 6.5344 _d -4
0518 a_QP_RNA_min = 2.23 _d -4
0519 a_QC_DNA = 9.41 _d -4
0520 a_QN_sto_max = 0.035 _d 0
0521 b_QN_sto_max = 0 _d 0
0522 a_Qp_max = 0.0052 _d 0
0523 b_Qp_max = 0 _d 0
0524 a_Qfe_max = 2.436 _d -4
0525 b_Qfe_max = 0 _d 0
0526 a_A_bio = 0.2711 _d 0*3600*24
0527 a_AP_RNA = 0.00423 _d 0*3600*24
0528 #endif
8fbfd1f382 Oliv*0529 #ifdef DARWIN_ALLOW_GEIDER
5910c7e5d2 Oliv*0530 a_mQyield = 75 _d -6
8fbfd1f382 Oliv*0531 a_chl2cmax = .3 _d 0
0532 a_inhibGeider = 0 _d 0
0533 #else
0534 a_ksatPAR = 0.012 _d 0
0535 a_kinhPAR = 6.0 _d -3
0536 #endif
0537 a_amminhib = 4.6 _d 0
0538 a_acclimtimescl = 1 _d 0/(20*pday)
0539 a_acclimtimescl_denom = 1 _d 0
367178819a Oliv*0540 a_ksatPON = 1 _d 0
0541 a_ksatDON = 1 _d 0
a092808e6b shlo*0542 a_hillnumDIN = 1 _d 0
0543 a_hillnumPO4 = 1 _d 0
0544 a_hillnumFeT = 1 _d 0
0545 a_hillnumSiO2 = 1 _d 0
8fbfd1f382 Oliv*0546 a_grazemax = 21.9 _d 0 / pday
0547 a_grazemax_denom = 1. _d 0
0548 b_grazemax = -0.16 _d 0
5910c7e5d2 Oliv*0549
8fbfd1f382 Oliv*0550 a_kgrazesat = 1.00 _d 0
0551 b_kgrazesat = 0.00 _d 0
0552 a_biosink = 0.28 _d -1 / pday
0553 a_biosink_denom = 1. _d 0
0554 b_biosink = 0.39 _d 0
0555 a_bioswim = 0.00 _d 0 / pday
0556 a_bioswim_denom = 1. _d 0
0557 b_bioswim = 0.18 _d 0
1f645d6fac Oliv*0558 a_bioswimDVM = 0.0 _d 0
0559 b_bioswimDVM = 0.0 _d 0
1d634e507f Oliv*0560 a_ppSig = 1.0 _d 0
0561 a_ppOpt = 1024. _d 0
0562 b_ppOpt = 0.00 _d 0
faa67d1773 Oliv*0563 palat_min = 0 _d 0
8fbfd1f382 Oliv*0564 a_PCmax = 1.00 _d 0 / pday
0565 a_PCmax_denom = 1. _d 0
0566 b_PCmax = -0.15 _d 0
0567 a_qcarbon = 1.80 _d -11
0568 b_qcarbon = 0.94 _d 0
5910c7e5d2 Oliv*0569
8fbfd1f382 Oliv*0570 a_respRate_c = 0.00 _d 0
0571 a_respRate_c_denom = 1. _d 0
0572 b_respRate_c = 0.93 _d 0
5910c7e5d2 Oliv*0573
8fbfd1f382 Oliv*0574 a_kexcC = 0.00 _d 0
0575 b_kexcC = -0.33 _d 0
0576 a_vmaxNO3 = 0.51 _d 0 / pday
0577 a_vmaxNO3_denom = 1. _d 0
0578 b_vmaxNO3 = -0.27 _d 0
0579 a_ksatNO3 = 0.17 _d 0
0580 b_ksatNO3 = 0.27 _d 0
0581 a_Qnmin = 0.07 _d 0
0582 b_Qnmin = -0.17 _d 0
0583 a_Qnmax = 0.25 _d 0
0584 b_Qnmax = -0.13 _d 0
5910c7e5d2 Oliv*0585
8fbfd1f382 Oliv*0586 a_kexcN = 0.00 _d 0
0587 b_kexcN = -0.33 _d 0
0588 a_vmaxNO2 = 0.51 _d 0 / pday
0589 a_vmaxNO2_denom = 1.0 _d 0
0590 b_vmaxNO2 = -0.27 _d 0
0591 a_ksatNO2 = 0.17 _d 0
0592 b_ksatNO2 = 0.27 _d 0
0593 a_ksatNO2fac = 1.0 _d 0
0594 a_vmaxNH4 = 0.26 _d 0 / pday
0595 a_vmaxNH4_denom = 1.0 _d 0
0596 b_vmaxNH4 = -0.27 _d 0
0597 a_ksatNH4 = 0.85 _d -1
0598 b_ksatNH4 = 0.27 _d 0
0599 a_ksatNH4fac = 0.5 _d 0
0600 a_vmaxN = 1.28 _d 0 / pday
0601 a_vmaxN_denom = 1. _d 0
0602 b_vmaxN = -0.27 _d 0
0603 a_vmaxPO4 = 0.77 _d -1 / pday
0604 a_vmaxPO4_denom = 1.0 _d 0
0605 b_vmaxPO4 = -0.27 _d 0
0606 a_ksatPO4 = 0.26 _d -1
0607 b_ksatPO4 = 0.27 _d 0
0608 a_Qpmin = 2.00 _d -3
0609 b_Qpmin = 0.00 _d 0
0610 a_Qpmax = 0.01 _d 0
0611 b_Qpmax = 0.00 _d 0
0612 a_kexcP = 0.24 _d -1 / pday
0613 b_kexcP = -0.33 _d 0
0614 a_vmaxSiO2 = 0.77 _d -1 / pday
0615 a_vmaxSiO2_denom = 1.0 _d 0
0616 b_vmaxSiO2 = -0.27 _d 0
0617 a_ksatSiO2 = 0.24 _d -1
0618 b_ksatSiO2 = 0.27 _d 0
5910c7e5d2 Oliv*0619
8fbfd1f382 Oliv*0620 a_Qsimin = 2.00 _d -3
5910c7e5d2 Oliv*0621
8fbfd1f382 Oliv*0622 b_Qsimin = 0.00 _d 0
5910c7e5d2 Oliv*0623
8fbfd1f382 Oliv*0624 a_Qsimax = 4.00 _d -3
5910c7e5d2 Oliv*0625
8fbfd1f382 Oliv*0626 b_Qsimax = 0.00 _d 0
0627 a_kexcSi = 0.00 _d 0 / pday
0628 b_kexcSi = 0.00 _d 0
5910c7e5d2 Oliv*0629
8fbfd1f382 Oliv*0630 a_vmaxFeT = 14.0 _d -6 / pday
0631 a_vmaxFeT_denom = 1.0 _d 0
0632 b_vmaxFeT = -0.27 _d 0
5910c7e5d2 Oliv*0633
8fbfd1f382 Oliv*0634 a_ksatFeT = 80.0 _d -6
0635 b_ksatFeT = 0.27 _d 0
5910c7e5d2 Oliv*0636
8fbfd1f382 Oliv*0637 a_Qfemin = 1.50 _d -6
0638 b_Qfemin = 0.00 _d 0
5910c7e5d2 Oliv*0639
8fbfd1f382 Oliv*0640 a_Qfemax = 80.0 _d -6
0641 b_Qfemax = 0.00 _d 0
0642 a_kexcFe = 0.00 _d 0 / pday
0643 b_kexcFe = 0.00 _d 0
0644 grp_ExportFracPreyPred = 0.5 _d 0
5910c7e5d2 Oliv*0645 grp_ass_eff = 0.7 _d 0
5e7acb36b1 daat*0646 a_PARpref = 0.0046 _d 0
0647 b_PARpref = 0.00 _d 0
0648 a_mortmaxDVM = 8.389 _d -5
0649
0650
0651 b_mortmaxDVM = -0.16 _d 0
0652 a_ksatDVM = 1.0 _d 0/a_R_NC
0653 b_ksatDVM = 0.0 _d 0
0654 a_ksatPARDVM = 0.46 _d 0
0655 b_ksatPARDVM = 0.0 _d 0
0656 a_fracPARmort = 0.9 _d 0
0657 b_fracPARmort = 0.0 _d 0
0658 a_ExportFracDVM = 0.5 _d 0
8fbfd1f382 Oliv*0659
0660 READ(UNIT=iUnit, NML=DARWIN_TRAIT_PARAMS)
0661
0662
0663 a_acclimtimescl = a_acclimtimescl/a_acclimtimescl_denom
0664 a_grazemax = a_grazemax/a_grazemax_denom
0665 a_biosink = a_biosink/a_biosink_denom
0666 a_bioswim = a_bioswim/a_bioswim_denom
0667 a_PCmax = a_PCmax/a_PCmax_denom
0668 a_respRate_c = a_respRate_c/a_respRate_c_denom
0669 a_vmaxNO3 = a_vmaxNO3/a_vmaxNO3_denom
0670 a_vmaxNO2 = a_vmaxNO2/a_vmaxNO2_denom
0671 a_vmaxNH4 = a_vmaxNH4/a_vmaxNH4_denom
0672 a_vmaxN = a_vmaxN/a_vmaxN_denom
0673 a_vmaxPO4 = a_vmaxPO4/a_vmaxPO4_denom
0674 a_vmaxSiO2 = a_vmaxSiO2/a_vmaxSiO2_denom
0675 a_vmaxFeT = a_vmaxFeT/a_vmaxFeT_denom
0676
0677
0678
0679
0680 IF ( oUnit .GE. 0 ) THEN
0681 WRITE(UNIT=oUnit, NML=DARWIN_RANDOM_PARAMS)
0682 WRITE(UNIT=oUnit, NML=DARWIN_TRAIT_PARAMS)
0683 ENDIF
0684
0685 #endif /* ALLOW_DARWIN */
0686
0687 RETURN
0688 END
