Back to home page

darwin3

 
 

    


File indexing completed on 2025-12-02 20:42:54 UTC

view on githubraw file Latest commit a092808e on 2025-12-02 20:09:45 UTC
8fbfd1f382 Oliv*0001 #ifdef ALLOW_DARWIN
                0002 
                0003 CBOP
                0004 C     !ROUTINE: DARWIN_TRAITS.h
                0005 C     !INTERFACE:
                0006 C #include DARWIN_TRAITS.h
                0007 
                0008 C     !DESCRIPTION:
                0009 C Contains run-time parameters for the darwin package
                0010 C the parameters in this file are traits
                0011 C
                0012 C Requires: DARWIN_SIZE.h
                0013 
5910c7e5d2 Oliv*0014 C--   COMMON /DARWIN_TRAITS_i/ Per-plankton traits (generated, maybe overwritten by data.traits)
                0015 C     isPhoto   :: 1: does photosynthesis, 0: not
                0016 C     bactType  :: 1: particle associated, 2: free living bacteria, 0: not bacteria
                0017 C     isAerobic :: 1: is aerobic bacteria (also set bactType), 0: not
                0018 C     isDenit   :: 1: is dentrifying bacteria (also set (bactType), 0: not
                0019 C     hasSi     :: 1: uses silica (Diatom), 0: not
                0020 C     hasPIC    :: 1: calcifying, 0: set R_PICPOC to zero
                0021 C     diazo     :: 1: use molecular instead of mineral nitrogen, 0: not
                0022 C     useNH4    :: 1: can use ammonia, 0: not
                0023 C     useNO2    :: 1: can use nitrite, 0: not
                0024 C     useNO3    :: 1: can use nitrate, 0: not
                0025 C     combNO    :: 1: combined nitrite/nitrate limitation, 0: not
                0026 C     isPrey    :: 1: can be grazed, 0: not
                0027 C     isPred    :: 1: can graze, 0: not
                0028 C     tempMort  :: 1: mortality is temperature dependent, 0: turn dependence off
                0029 C     tempMort2 :: 1: quadratic mortality is temperature dependent, 0: turn dependence off
                0030 C     tempGraz  :: 1: grazing is temperature dependent, 0: turn dependence off
8fbfd1f382 Oliv*0031       COMMON /DARWIN_TRAITS_i/
                0032      &    isPhoto,
                0033      &    bactType,
                0034      &    isAerobic,
                0035      &    isDenit,
                0036      &    hasSi,
                0037      &    hasPIC,
                0038      &    diazo,
                0039      &    useNH4,
                0040      &    useNO2,
                0041      &    useNO3,
                0042      &    combNO,
                0043      &    isPrey,
                0044      &    isPred,
                0045      &    tempMort,
14aae40637 Oliv*0046      &    tempMort2,
                0047      &    tempGraz
8fbfd1f382 Oliv*0048       INTEGER isPhoto(nplank)
                0049       INTEGER bactType(nplank)
                0050       INTEGER isAerobic(nplank)
                0051       INTEGER isDenit(nplank)
                0052       INTEGER hasSi(nplank)
                0053       INTEGER hasPIC(nplank)
                0054       INTEGER diazo(nplank)
                0055       INTEGER useNH4(nplank)
                0056       INTEGER useNO2(nplank)
                0057       INTEGER useNO3(nplank)
                0058       INTEGER combNO(nplank)
                0059       INTEGER isPrey(nplank)
                0060       INTEGER isPred(nplank)
                0061       INTEGER tempMort(nplank)
                0062       INTEGER tempMort2(nplank)
14aae40637 Oliv*0063       INTEGER tempGraz(nplank)
5910c7e5d2 Oliv*0064 
                0065 C--   COMMON /DARWIN_TRAITS_r/ Per-plankton traits (generated, maybe overwritten by data.traits)
                0066 C     Xmin               :: [mmol C m^-3]              minimum abundance for mortality, respiration and exudation
                0067 C     amminhib           :: [(mmol N m^-3)^-1]         coefficient for NH4 inhibition of NO uptake
                0068 C     acclimtimescl      :: [s^-1]                     rate of chlorophyll acclimation
                0069 C
                0070 C     mort               :: [s^-1]                     linear mortality rate
                0071 C     mort2              :: [(mmol C m^-3)^-1 s^-1]    quadratic mortality coefficient
                0072 C     ExportFracMort     :: []                         fraction of linear mortality to POM
                0073 C     ExportFracMort2    :: []                         fraction of quadratic mortality to POM
                0074 C     ExportFracExude    :: []                         fraction of exudation to POM
c7b6c66d45 Oliv*0075 C     FracExudeC         :: []                         fraction of excess C exuded
5910c7e5d2 Oliv*0076 C
                0077 C- temperature dependence:
                0078 C     phytoTempCoeff     :: []                         see :numref:`pkg_darwin_temperature_params`
                0079 C     phytoTempExp1      :: [exp(1/degrees C)]         see :numref:`pkg_darwin_temperature_params`
09b16dd232 Oliv*0080 C     phytoTempAe        :: [1/degrees C]              see :numref:`pkg_darwin_temperature_params`
5910c7e5d2 Oliv*0081 C     phytoTempExp2      :: []                         see :numref:`pkg_darwin_temperature_params`
                0082 C     phytoTempOptimum   :: [degrees C]                see :numref:`pkg_darwin_temperature_params`
                0083 C     phytoDecayPower    :: []                         see :numref:`pkg_darwin_temperature_params`
09b16dd232 Oliv*0084 C     hetTempAe          :: [1/degrees C]              see :numref:`pkg_darwin_temperature_params`
                0085 C     hetTempExp2        :: []                         see :numref:`pkg_darwin_temperature_params`
                0086 C     hetTempOptimum     :: [degrees C]                see :numref:`pkg_darwin_temperature_params`
                0087 C     hetDecayPower      :: []                         see :numref:`pkg_darwin_temperature_params`
                0088 C     grazTempAe         :: [1/degrees C]              see :numref:`pkg_darwin_temperature_params`
                0089 C     grazTempExp2       :: []                         see :numref:`pkg_darwin_temperature_params`
                0090 C     grazTempOptimum    :: [degrees C]                see :numref:`pkg_darwin_temperature_params`
                0091 C     grazDecayPower     :: []                         see :numref:`pkg_darwin_temperature_params`
5910c7e5d2 Oliv*0092 C
                0093 C     R_NC               :: [mmol N (mmol C)^-1]       nitrogen-carbon ratio (not used with DARWIN_ALLOW_NQUOTA)
                0094 C     R_PC               :: [mmol P (mmol C)^-1]       phosphorus-carbon ratio (not used with DARWIN_ALLOW_PQUOTA)
                0095 C     R_SiC              :: [mmol Si (mmol C)^-1]      silica-carbon ratio (not used with DARWIN_ALLOW_SIQUOTA)
                0096 C     R_FeC              :: [mmol Fe (mmol C)^-1]      iron-carbon ratio (not used with DARWIN_ALLOW_FEQUOTA)
                0097 C     R_ChlC             :: [mg Chl (mmol C)^-1]       chlorophyll-carbon ratio (not used with DARWIN_ALLOW_CHLQUOTA)
                0098 C     R_PICPOC           :: [mmol PIC (mmol POC)^-1]   inorganic-organic carbon ratio
                0099 C
                0100 C     biosink            :: [m s^-1]                   sinking velocity (positive downwards)
                0101 C     bioswim            :: [m s^-1]                   upward swimming velocity (positive upwards)
5e7acb36b1 daat*0102 C     bioswimDVM         :: [m s^-1]                   swimming speed for diel vertically migrating plankton (pos. upwards)
5910c7e5d2 Oliv*0103 C
                0104 C     respRate           :: [s^-1]                     respiration rate
                0105 C     PCmax              :: [s^-1]                     maximum carbon-specific growth rate
                0106 C
                0107 C     Qnmax              :: [mmol N (mmol C)^-1]       maximum nitrogen quota (only with DARWIN_ALLOW_NQUOTA)
                0108 C     Qnmin              :: [mmol N (mmol C)^-1]       minimum nitrogen quota (only with DARWIN_ALLOW_NQUOTA)
                0109 C     Qpmax              :: [mmol P (mmol C)^-1]       maximum phosphorus quota (only with DARWIN_ALLOW_PQUOTA)
                0110 C     Qpmin              :: [mmol P (mmol C)^-1]       minimum phosphorus quota (only with DARWIN_ALLOW_PQUOTA)
                0111 C     Qsimax             :: [mmol Si (mmol C)^-1]      maximum silica quota (only with DARWIN_ALLOW_SIQUOTA)
                0112 C     Qsimin             :: [mmol Si (mmol C)^-1]      minimum silica quota (only with DARWIN_ALLOW_SIQUOTA)
                0113 C     Qfemax             :: [mmol Fe (mmol C)^-1]      maximum iron quota (only with DARWIN_ALLOW_FEQUOTA)
                0114 C     Qfemin             :: [mmol Fe (mmol C)^-1]      minimum iron quota (only with DARWIN_ALLOW_FEQUOTA)
                0115 C
                0116 C     VmaxNH4            :: [mmol N (mmol C)^-1 s^-1]  maximum ammonia uptake rate (only with DARWIN_ALLOW_NQUOTA)
                0117 C     VmaxNO2            :: [mmol N (mmol C)^-1 s^-1]  maximum nitrite uptake rate (only with DARWIN_ALLOW_NQUOTA)
                0118 C     VmaxNO3            :: [mmol N (mmol C)^-1 s^-1]  maximum nitrate uptake rate (only with DARWIN_ALLOW_NQUOTA)
                0119 C     VmaxN              :: [mmol N (mmol C)^-1 s^-1]  maximum nitrogen uptake rate for diazotrophs (only with DARWIN_ALLOW_NQUOTA)
                0120 C                        ::                            has to be >= vmaxNO3 + vmaxNO2 + vmaxNH4
                0121 C     VmaxPO4            :: [mmol P (mmol C)^-1 s^-1]  maximum phosphate uptake rate (only with DARWIN_ALLOW_PQUOTA)
                0122 C     VmaxSiO2           :: [mmol Si (mmol C)^-1 s^-1] maximum silica uptake rate (only with DARWIN_ALLOW_SIQUOTA)
                0123 C     VmaxFeT            :: [mmol Fe (mmol C)^-1 s^-1] maximum iron uptake rate (only with DARWIN_ALLOW_FEQUOTA)
                0124 C
                0125 C     ksatNH4            :: [mmol N m^-3]              half-saturation conc. for ammonia uptake/limitation
                0126 C     ksatNO2            :: [mmol N m^-3]              half-saturation conc. for nitrite uptake/limitation
                0127 C     ksatNO3            :: [mmol N m^-3]              half-saturation conc. for nitrate uptake/limitation
                0128 C     ksatPO4            :: [mmol P m^-3]              half-saturation conc. for phosphate uptake/limitation
                0129 C     ksatSiO2           :: [mmol Si m^-3]             half-saturation conc. for silica uptake/limitation
                0130 C     ksatFeT            :: [mmol Fe m^-3]             half-saturation conc. for iron uptake/limitation
                0131 C
                0132 C     kexcc              :: [s^-1]                 exudation rate for carbon
                0133 C     kexcn              :: [s^-1]                 exudation rate for nitrogen
                0134 C     kexcp              :: [s^-1]                 exudation rate for phosphorus
                0135 C     kexcsi             :: [s^-1]                 exudation rate for silica
                0136 C     kexcfe             :: [s^-1]                 exudation rate for iron
                0137 C
                0138 C- Geider
                0139 C     inhibGeider        :: []                     photo-inhibition coefficient for Geider growth
                0140 C- not Geider
                0141 C     ksatPAR            :: [(uEin m^-2 s^-1)^-1]  saturation coefficient for PAR (w/o GEIDER)
                0142 C     kinhPAR            :: [(uEin m^-2 s^-1)^-1]  inhibition coefficient for PAR (w/o GEIDER)
                0143 C- always defined
                0144 C     mQyield            :: [mmol C (uEin)^-1]     maximum quantum yield
                0145 C     chl2cmax           :: [mg Chl (mmol C)^-1]   maximum Chlorophyll-carbon ratio
                0146 C
                0147 C- Grazing
                0148 C     grazemax           :: [s^-1]          maximum grazing rate
                0149 C     kgrazesat          :: [mmol C m^-3]   grazing half-saturation concentration
                0150 C     palat              :: []              palatability matrix
                0151 C     asseff             :: []              assimilation efficiency matrix
                0152 C     ExportFracPreyPred :: []              fraction of unassimilated prey becoming particulate organic matter
                0153 C
                0154 C- Bacteria
                0155 C     yield              :: []              bacterial growth yield for all organic matter
                0156 C     yieldO2            :: []              bacterial growth yield for oxygen
                0157 C     yieldNO3           :: []              bacterial growth yield for nitrate
                0158 C     ksatPON            :: [mmol N m^-3]   half-saturation of PON for bacterial growth
                0159 C     ksatPOC            :: [mmol C m^-3]   half-saturation of POC for bacterial growth
                0160 C     ksatPOP            :: [mmol P m^-3]   half-saturation of POP for bacterial growth
                0161 C     ksatPOFe           :: [mmol Fe m^-3]  half-saturation of POFe for bacterial growth
                0162 C     ksatDON            :: [mmol N m^-3]   half-saturation of DON for bacterial growth
                0163 C     ksatDOC            :: [mmol C m^-3]   half-saturation of DOC for bacterial growth
                0164 C     ksatDOP            :: [mmol P m^-3]   half-saturation of DOP for bacterial growth
                0165 C     ksatDOFe           :: [mmol Fe m^-3]  half-saturation of DOFe for bacterial growth
5e7acb36b1 daat*0166 C- PAR preference
                0167 C     PARpref            :: [uEin m^-2 s^-1] preferred PAR isolume for determining DVM
                0168 C- Light-mediated mortality for DVM
                0169 C     mortmaxDVM         :: [s^-1]           maximum mortality rate related to PAR for DVM
                0170 C     ksatDVM            :: [mmol C m^-3]    half saturation for DVM mortality
                0171 C     ksatPARDVM         :: [uEin m^-2 s^-1] half sat for light limitation for DVM
                0172 C     fracPARmort        :: []               fraction of mortality from light-dependent mortality
                0173 C     ExportFracDVM      :: []               fraction of light-dep mortality from DVM to POM
a092808e6b shlo*0174 C
                0175 C- MACROMOLECULAR_GROWTH parameters
                0176 C     Y_CP_Plip          :: [molC molP^-1]   C/P molar ratio of thylacoid membrane
                0177 C     Y_CN_protein       :: [molC molN^-1]   C/N molar ratio in protein
                0178 C     Y_NC_chl           :: [molN molC^-1]   N/C molar ratio in chlorophyll
                0179 C     Y_CN_cyano         :: [molC molN^-1]   C/N molar ratio of cyanophycin
                0180 C     Y_PN_nucacid       :: [molP molN^-1]   P/N molar ratio of RNA
                0181 C     Y_CN_DNA           :: [molC molN^-1]   C/N molar ratio of DNA
                0182 C     Y_CN_RNA           :: [molC molN^-1]   C/N molar ratio of RNA
                0183 C     Y_THY_P            :: [(molP)/(molC in chl)^-1] the stoichiometric ratio for cell phosphorus in thylakoid membrane to chlorophyll
                0184 C     Y_FeN_photo        :: [molFe mol N^-1] Fe/N ratio in photosystem iron
                0185 C     ECo2Prod           :: [dimensionless]  CO2 production ratio
                0186 C     maintConsum        :: [s^-1]           maintenance carbohydrate consumption
                0187 C     VI_max             :: [molC (molC in Chl)^-1 s^-1] carbon fixing rate
                0188 C     A_I                :: [umol^-1 m2 s]   coefficient characterizing the absorption cross section
                0189 C     QC_other           :: [molC molC^-1]   essential carbon (lipid membrane, etc.)
                0190 C     QC_pro_other       :: [molC molC^-1]
                0191 C     QP_other           :: [molP molC^-1]   constant part of phosphorus
                0192 C     QP_RNA_min         :: [molP molC^-1]
                0193 C     QC_DNA             :: [molC molC^-1]   constant part of DNA in carbon
                0194 C     QN_pro_other       :: [molN molC^-1]
                0195 C     QN_RNA_min         :: [molN molC^-1]   constant part of RNA in nitrogen
                0196 C     QC_RNA_min         :: [molN molC^-1]   constant part of RNA in carbon
                0197 C     QN_DNA             :: [molN molC^-1]   DNA in nitrogen
                0198 C     QP_DNA             :: [molP molC^-1]   DNA in phosphorous
                0199 C     QN_sto_max         :: [molN molC^-1]   maximum nitrogen storage
                0200 C     Qp_max             :: [molP molC^-1]   maximum phosphorus quota
                0201 C     Qfe_max            :: [molFe molC^-1]  maximum iron quota
                0202 C     A_pho              :: [(molC)/(molC in chl)^-1] A constant of proportionalty
                0203 C     A_bio              :: [molC molC^-1 s] constant for variable part of biosynthesis protein
                0204 C     AP_RNA             :: [molP molC^-1 s] constant for Variable part of RNA
                0205 C     AN_RNA             :: [molN molN^-1 s] constant for Variable part of RNA
                0206 C     A_thy              :: [molC (molC in chl)^-1]
                0207 C     Sf                 :: [unitless]       enhancement of photosynthesis due to size
                0208 C     QC_const           :: [molC molC^-1]   constant portion of the cell
                0209 C     VI_min             :: [molC (molC in Chl)^-1 s^-1] minimum photosynthesis rate
                0210 C     QC_chlMax          :: [molC molC^-1]   maximum chlorophyll concentration at minimum light
                0211 C     QnNoChl            :: [molN molC^-1]   minimum QN at zero growth rate
                0212 C     QpNoChl            :: [molP molC^-1]   minimum QP at zero growth rate
                0213 C     QfeNoChl           :: [molFe molC^-1]  minimum QFe at zero growth rate
8fbfd1f382 Oliv*0214       COMMON /DARWIN_TRAITS_r/
                0215      &    Xmin,
                0216      &    amminhib,
                0217      &    acclimtimescl,
                0218      &    mort,
                0219      &    mort2,
                0220      &    ExportFracMort,
                0221      &    ExportFracMort2,
                0222      &    ExportFracExude,
c7b6c66d45 Oliv*0223      &    FracExudeC,
8fbfd1f382 Oliv*0224      &    phytoTempCoeff,
                0225      &    phytoTempExp1,
09b16dd232 Oliv*0226      &    phytoTempAe,
8fbfd1f382 Oliv*0227      &    phytoTempExp2,
                0228      &    phytoTempOptimum,
                0229      &    phytoDecayPower,
09b16dd232 Oliv*0230      &    hetTempAe,
                0231      &    hetTempExp2,
                0232      &    hetTempOptimum,
                0233      &    hetDecayPower,
                0234      &    grazTempAe,
                0235      &    grazTempExp2,
                0236      &    grazTempOptimum,
                0237      &    grazDecayPower,
8fbfd1f382 Oliv*0238      &    R_NC,
                0239      &    R_PC,
                0240      &    R_SiC,
                0241      &    R_FeC,
                0242      &    R_ChlC,
                0243      &    R_PICPOC,
                0244      &    biosink,
                0245      &    bioswim,
5e7acb36b1 daat*0246      &    bioswimDVM,
8fbfd1f382 Oliv*0247      &    respRate,
                0248      &    PCmax,
                0249      &    Qnmax,
                0250      &    Qnmin,
                0251      &    Qpmax,
                0252      &    Qpmin,
                0253      &    Qsimax,
                0254      &    Qsimin,
                0255      &    Qfemax,
                0256      &    Qfemin,
                0257      &    VmaxNH4,
                0258      &    VmaxNO2,
                0259      &    VmaxNO3,
                0260      &    VmaxN,
                0261      &    VmaxPO4,
                0262      &    VmaxSiO2,
                0263      &    VmaxFeT,
                0264      &    ksatNH4,
                0265      &    ksatNO2,
                0266      &    ksatNO3,
                0267      &    ksatPO4,
                0268      &    ksatSiO2,
                0269      &    ksatFeT,
a092808e6b shlo*0270      &    hillnumDIN,
                0271      &    hillnumPO4,
                0272      &    hillnumFeT,
                0273      &    hillnumSiO2,
8fbfd1f382 Oliv*0274      &    kexcc,
                0275      &    kexcn,
                0276      &    kexcp,
                0277      &    kexcsi,
                0278      &    kexcfe,
                0279 #ifdef DARWIN_ALLOW_GEIDER
                0280      &    inhibGeider,
                0281 #else
                0282      &    ksatPAR,
                0283      &    kinhPAR,
                0284 #endif
                0285      &    mQyield,
                0286      &    chl2cmax,
                0287      &    grazemax,
                0288      &    kgrazesat,
                0289      &    palat,
                0290      &    asseff,
                0291      &    ExportFracPreyPred,
                0292      &    yield,
                0293      &    yieldO2,
                0294      &    yieldNO3,
                0295      &    ksatPON,
                0296      &    ksatPOC,
                0297      &    ksatPOP,
                0298      &    ksatPOFe,
                0299      &    ksatDON,
                0300      &    ksatDOC,
                0301      &    ksatDOP,
5e7acb36b1 daat*0302      &    ksatDOFe,
                0303      &    PARpref,
                0304      &    mortmaxDVM,
                0305      &    ksatDVM,
                0306      &    ksatPARDVM,
                0307      &    fracPARmort,
                0308      &    ExportFracDVM
a092808e6b shlo*0309 #ifdef DARWIN_MACROMOLECULAR_GROWTH
                0310      &   ,Y_CN_protein,
                0311      &    Y_CP_Plip,
                0312      &    Y_NC_chl,
                0313      &    Y_CN_cyano,
                0314      &    Y_PN_nucacid,
                0315      &    Y_CN_DNA,
                0316      &    Y_CN_RNA,
                0317      &    Y_THY_P,
                0318      &    Y_FeN_photo,
                0319      &    ECo2Prod,
                0320      &    maintConsum,
                0321      &    VI_max,
                0322      &    A_I,
                0323      &    QC_other,
                0324      &    QC_pro_other,
                0325      &    QP_other,
                0326      &    QP_RNA_min,
                0327      &    QC_DNA,
                0328      &    QN_pro_other,
                0329      &    QN_RNA_min,
                0330      &    QC_RNA_min,
                0331      &    QN_DNA,
                0332      &    QP_DNA,
                0333      &    QN_sto_max,
                0334      &    Qp_max,
                0335      &    Qfe_max,
                0336      &    A_pho,
                0337      &    A_bio,
                0338      &    AP_RNA,
                0339      &    AN_RNA,
                0340      &    A_thy,
                0341      &    Sf,
                0342      &    VI_min,
                0343      &    QC_chlMax,
                0344      &    QnNoChl,
                0345      &    QpNoChl,
                0346      &    QfeNoChl,
                0347      &    QC_const
                0348 #endif
8fbfd1f382 Oliv*0349       _RL Xmin(nplank)
                0350       _RL amminhib(nplank)
                0351       _RL acclimtimescl(nplank)
                0352       _RL mort(nplank)
                0353       _RL mort2(nplank)
                0354       _RL ExportFracMort(nplank)
                0355       _RL ExportFracMort2(nplank)
                0356       _RL ExportFracExude(nplank)
c7b6c66d45 Oliv*0357       _RL FracExudeC(nplank)
8fbfd1f382 Oliv*0358       _RL phytoTempCoeff(nplank)
                0359       _RL phytoTempExp1(nplank)
09b16dd232 Oliv*0360       _RL phytoTempAe(nplank)
8fbfd1f382 Oliv*0361       _RL phytoTempExp2(nplank)
                0362       _RL phytoTempOptimum(nplank)
                0363       _RL phytoDecayPower(nplank)
09b16dd232 Oliv*0364       _RL hetTempAe(nplank)
                0365       _RL hetTempExp2(nplank)
                0366       _RL hetTempOptimum(nplank)
                0367       _RL hetDecayPower(nplank)
                0368       _RL grazTempAe(nplank)
                0369       _RL grazTempExp2(nplank)
                0370       _RL grazTempOptimum(nplank)
                0371       _RL grazDecayPower(nplank)
8fbfd1f382 Oliv*0372       _RL R_NC(nplank)
                0373       _RL R_PC(nplank)
                0374       _RL R_SiC(nplank)
                0375       _RL R_FeC(nplank)
                0376       _RL R_ChlC(nplank)
                0377       _RL R_PICPOC(nplank)
                0378       _RL biosink(nplank)
                0379       _RL bioswim(nplank)
5e7acb36b1 daat*0380       _RL bioswimDVM(nplank)
8fbfd1f382 Oliv*0381       _RL respRate(nplank)
                0382       _RL PCmax(nplank)
                0383       _RL Qnmax(nplank)
                0384       _RL Qnmin(nplank)
                0385       _RL Qpmax(nplank)
                0386       _RL Qpmin(nplank)
                0387       _RL Qsimax(nplank)
                0388       _RL Qsimin(nplank)
                0389       _RL Qfemax(nplank)
                0390       _RL Qfemin(nplank)
                0391       _RL VmaxNH4(nplank)
                0392       _RL VmaxNO2(nplank)
                0393       _RL VmaxNO3(nplank)
                0394       _RL VmaxN(nplank)
                0395       _RL VmaxPO4(nplank)
                0396       _RL VmaxSiO2(nplank)
                0397       _RL VmaxFeT(nplank)
                0398       _RL ksatNH4(nplank)
                0399       _RL ksatNO2(nplank)
                0400       _RL ksatNO3(nplank)
                0401       _RL ksatPO4(nplank)
                0402       _RL ksatSiO2(nplank)
                0403       _RL ksatFeT(nplank)
a092808e6b shlo*0404       _RL hillnumDIN(nplank)
                0405       _RL hillnumPO4(nplank)
                0406       _RL hillnumFeT(nplank)
                0407       _RL hillnumSiO2(nplank)
8fbfd1f382 Oliv*0408       _RL kexcc(nplank)
                0409       _RL kexcn(nplank)
                0410       _RL kexcp(nplank)
                0411       _RL kexcsi(nplank)
                0412       _RL kexcfe(nplank)
                0413 #ifdef DARWIN_ALLOW_GEIDER
                0414       _RL inhibGeider(nplank)
                0415 #else
                0416       _RL ksatPAR(nplank)
                0417       _RL kinhPAR(nplank)
                0418 #endif
                0419       _RL mQyield(nplank)
                0420       _RL chl2cmax(nplank)
                0421       _RL grazemax(nplank)
                0422       _RL kgrazesat(nplank)
                0423       _RL palat(nplank,nplank)
                0424       _RL asseff(nplank,nplank)
                0425       _RL ExportFracPreyPred(nplank,nplank)
                0426       _RL yield(nplank)
                0427       _RL yieldO2(nplank)
                0428       _RL yieldNO3(nplank)
                0429       _RL ksatPON(nplank)
                0430       _RL ksatPOC(nplank)
                0431       _RL ksatPOP(nplank)
                0432       _RL ksatPOFe(nplank)
                0433       _RL ksatDON(nplank)
                0434       _RL ksatDOC(nplank)
                0435       _RL ksatDOP(nplank)
                0436       _RL ksatDOFe(nplank)
5e7acb36b1 daat*0437       _RL PARpref(nplank)
                0438       _RL mortmaxDVM(nplank)
                0439       _RL ksatDVM(nplank)
                0440       _RL ksatPARDVM(nplank)
                0441       _RL fracPARmort(nplank)
                0442       _RL ExportFracDVM(nplank)
a092808e6b shlo*0443 #ifdef DARWIN_MACROMOLECULAR_GROWTH
                0444       _RL Y_CN_protein(nplank)
                0445       _RL Y_CP_Plip(nplank)
                0446       _RL Y_NC_chl(nplank)
                0447       _RL Y_CN_cyano(nplank)
                0448       _RL Y_PN_nucacid(nplank)
                0449       _RL Y_CN_DNA(nplank)
                0450       _RL Y_CN_RNA(nplank)
                0451       _RL Y_THY_P(nplank)
                0452       _RL Y_FeN_photo(nplank)
                0453       _RL ECo2Prod(nplank)
                0454       _RL maintConsum(nplank)
                0455       _RL VI_max(nplank)
                0456       _RL A_I(nplank)
                0457       _RL QC_other(nplank)
                0458       _RL QC_pro_other(nplank)
                0459       _RL QP_other(nplank)
                0460       _RL QP_RNA_min(nplank)
                0461       _RL QC_DNA(nplank)
                0462       _RL QN_pro_other(nplank)
                0463       _RL QN_RNA_min(nplank)
                0464       _RL QC_RNA_min(nplank)
                0465       _RL QN_DNA(nplank)
                0466       _RL QP_DNA(nplank)
                0467       _RL QN_sto_max(nplank)
                0468       _RL Qp_max(nplank)
                0469       _RL Qfe_max(nplank)
                0470       _RL A_pho(nplank)
                0471       _RL A_bio(nplank)
                0472       _RL AP_RNA(nplank)
                0473       _RL AN_RNA(nplank)
                0474       _RL A_thy(nplank)
                0475       _RL Sf(nplank)
                0476       _RL VI_min(nplank)
                0477       _RL QC_chlMax(nplank)
                0478       _RL QnNoChl(nplank)
                0479       _RL QpNoChl(nplank)
                0480       _RL QfeNoChl(nplank)
                0481       _RL QC_const(nplank)
                0482 #endif
8fbfd1f382 Oliv*0483 
5910c7e5d2 Oliv*0484 C--   COMMON /DARWIN_DEPENDENT_TRAITS_i/ Dependent and constant (not read-in) parameters
                0485 C     group  :: which group this type belongs to
                0486 C     igroup :: index within group
                0487 C
                0488 C- Radtrans only:
                0489 C     aptype :: optical type (for absorption/scattering spectra)
8fbfd1f382 Oliv*0490       COMMON /DARWIN_DEPENDENT_TRAITS_i/
                0491 #ifdef ALLOW_RADTRANS
                0492      &    aptype,
                0493 #endif
                0494      &    group,
                0495      &    igroup
                0496 #ifdef ALLOW_RADTRANS
                0497       INTEGER aptype(nplank)
                0498 #endif
                0499       INTEGER group(nplank)
                0500       INTEGER igroup(nplank)
5910c7e5d2 Oliv*0501 
                0502 C--   COMMON /DARWIN_DEPENDENT_TRAITS_r/ Dependent and constant (not read-in) parameters
                0503 C     normI            :: []                    normalization factor for non-Geider light curve
                0504 C     biovol           :: [um^3]                volume
                0505 C     qcarbon          :: [mmol C/cell]         cellular carbon content
                0506 C     biovol_bygroup   :: [um^3]                volume of types in each group
                0507 C     chl2cmin         :: [mg Chl (mmol C)^-1]  minimum Chl:C ratio (function of chl2cmax and alpha_mean)
                0508 C     alpha_mean       :: [mmol C s-1 (uEin m^-2 s^-1)^-1 (mg Chl)^-1]  mean initial slope of light curve (over wavebands)
8fbfd1f382 Oliv*0509       COMMON /DARWIN_DEPENDENT_TRAITS_r/
                0510 #ifndef DARWIN_ALLOW_GEIDER
                0511      &    normI,
                0512 #endif
                0513      &    biovol,
                0514      &    qcarbon,
                0515      &    biovol_bygroup,
                0516      &    alpha_mean,
7e97ec2956 Oliv*0517      &    chl2cmin
8fbfd1f382 Oliv*0518 #ifndef DARWIN_ALLOW_GEIDER
                0519       _RL normI(nplank)
                0520 #endif
                0521       _RL biovol(nplank)
                0522       _RL qcarbon(nplank)
                0523       _RL biovol_bygroup(nplank,ngroup)
                0524       _RL alpha_mean(nplank)
                0525       _RL chl2cmin(nplank)
                0526 
                0527 
                0528 #endif /* ALLOW_DARWIN */
                0529