darwin3 |
|
|
|
|||
File indexing completed on 2025-12-02 20:42:53 UTC
view on githubraw file Latest commit a092808e on 2025-12-02 20:09:45 UTC8fbfd1f382 Oliv*0001 #ifdef ALLOW_DARWIN 0002 0003 CBOP 0004 C !ROUTINE: DARWIN_PARAMS.h 0005 C !INTERFACE: 0006 C #include DARWIN_PARAMS.h 0007 0008 C !DESCRIPTION: 0009 C Contains run-time parameters for the darwin package 0010 C 0011 C Requires: DARWIN_SIZE.h 0012 0013 _RL DARWIN_UNINIT_RL 0014 PARAMETER(DARWIN_UNINIT_RL=-999999999 _d 0) 0015 5910c7e5d2 Oliv*0016 C-- COMMON/darwin_forcing_params_l/ darwin parameters related to forcing 0017 C darwin_chlInitBalanced :: Initialize Chlorophyll to a balanced value following Geider 6ed0cf4850 Oliv*0018 C darwin_haveSurfPAR :: whether PAR forcing fields/constant were given 0019 C darwin_haveVentHe3 :: whether hydrothermal vent forcing was given 5910c7e5d2 Oliv*0020 C darwin_useSEAICE :: whether to use ice area from seaice pkg 0021 C darwin_useQsw :: whether to use model shortwave radiation 0022 C darwin_useEXFwind :: whether to use wind speed from exf package 8fbfd1f382 Oliv*0023 COMMON/darwin_forcing_params_l/ 0024 & darwin_chlInitBalanced, 0025 & darwin_haveSurfPAR, 6ed0cf4850 Oliv*0026 & darwin_haveVentHe3, 8fbfd1f382 Oliv*0027 & darwin_useSEAICE, 1342b05049 Oliv*0028 & darwin_useQsw, 8fbfd1f382 Oliv*0029 & darwin_useEXFwind 0030 LOGICAL darwin_chlInitBalanced 0031 LOGICAL darwin_haveSurfPAR 6ed0cf4850 Oliv*0032 LOGICAL darwin_haveVentHe3 8fbfd1f382 Oliv*0033 LOGICAL darwin_useSEAICE 1342b05049 Oliv*0034 LOGICAL darwin_useQsw 8fbfd1f382 Oliv*0035 LOGICAL darwin_useEXFwind 0036 5910c7e5d2 Oliv*0037 C-- COMMON/darwin_forcing_params_i/ darwin parameters related to forcing 0038 C darwin_chlIter0 :: Iteration number when to initialize Chlorophyll 8fbfd1f382 Oliv*0039 COMMON/darwin_forcing_params_i/ 0040 & darwin_chlIter0 0041 INTEGER darwin_chlIter0 0042 5910c7e5d2 Oliv*0043 C-- COMMON /DARWIN_CONSTANTS_r/ 0044 C rad2deg :: 8fbfd1f382 Oliv*0045 COMMON /DARWIN_CONSTANTS_r/ 0046 & rad2deg 0047 _RL rad2deg 0048 0049 #ifdef DARWIN_ALLOW_CARBON 5910c7e5d2 Oliv*0050 C-- COMMON /CARBON_CONSTANTS_r/ Coefficients for DIC chemistry 0051 C Pa2Atm :: Convert pressure in Pascal to atm 0052 C ptr2mol :: convert ptracers (in mmol/m3) to mol/m3 0053 C- 0054 C sca1 :: Schmidt no. coefficient for CO2 0055 C sca2 :: Schmidt no. coefficient for CO2 0056 C sca3 :: Schmidt no. coefficient for CO2 0057 C sca4 :: Schmidt no. coefficient for CO2 0058 C- 0059 C sox1 :: [] Schmidt no. coefficient for O2 [Keeling et al, GBC, 12, 141, (1998)] 0060 C sox2 :: [] Schmidt no. coefficient for O2 [Keeling et al, GBC, 12, 141, (1998)] 0061 C sox3 :: [] Schmidt no. coefficient for O2 [Keeling et al, GBC, 12, 141, (1998)] 0062 C sox4 :: [] Schmidt no. coefficient for O2 [Keeling et al, GBC, 12, 141, (1998)] 0063 C- 0064 C oA0 :: Coefficient for determining saturation O2 0065 C oA1 :: Coefficient for determining saturation O2 0066 C oA2 :: Coefficient for determining saturation O2 0067 C oA3 :: Coefficient for determining saturation O2 0068 C oA4 :: Coefficient for determining saturation O2 0069 C oA5 :: Coefficient for determining saturation O2 0070 C oB0 :: Coefficient for determining saturation O2 0071 C oB1 :: Coefficient for determining saturation O2 0072 C oB2 :: Coefficient for determining saturation O2 0073 C oB3 :: Coefficient for determining saturation O2 0074 C oC0 :: Coefficient for determining saturation O2 8fbfd1f382 Oliv*0075 COMMON /CARBON_CONSTANTS_r/ 0076 & Pa2Atm, 0077 & ptr2mol, 0078 & sca1, 0079 & sca2, 0080 & sca3, 0081 & sca4, 0082 & sox1, 0083 & sox2, 0084 & sox3, 0085 & sox4, 0086 & oA0, 0087 & oA1, 0088 & oA2, 0089 & oA3, 0090 & oA4, 0091 & oA5, 0092 & oB0, 0093 & oB1, 0094 & oB2, 0095 & oB3, 0096 & oC0 0097 _RL Pa2Atm 0098 _RL ptr2mol 0099 _RL sca1 0100 _RL sca2 0101 _RL sca3 0102 _RL sca4 0103 _RL sox1 0104 _RL sox2 0105 _RL sox3 0106 _RL sox4 0107 _RL oA0 0108 _RL oA1 0109 _RL oA2 0110 _RL oA3 0111 _RL oA4 0112 _RL oA5 0113 _RL oB0 0114 _RL oB1 0115 _RL oB2 0116 _RL oB3 0117 _RL oC0 6aa5674e10 Oliv*0118 8a10163480 Oliv*0119 # ifdef DARWIN_ALLOW_RADIv1 6aa5674e10 Oliv*0120 C-- COMMON /DARWIN_SED_CONSTANTS_R/ Coefficients for sediment model 0121 C sed_a1 :: coefficient for sediment model 0122 C sed_b1 :: coefficient for sediment model 0123 C sed_c1 :: coefficient for sediment model 0124 C sed_d1 :: coefficient for sediment model 0125 C sed_a2 :: coefficient for sediment model 0126 C sed_b2 :: coefficient for sediment model 0127 C sed_c2 :: coefficient for sediment model 0128 C sed_d2 :: coefficient for sediment model 0129 C sed_a3 :: coefficient for sediment model 0130 C sed_b3 :: coefficient for sediment model 0131 C sed_c3 :: coefficient for sediment model 0132 C sed_d3 :: coefficient for sediment model 0133 C sed_a4 :: coefficient for sediment model 0134 C sed_b4 :: coefficient for sediment model 0135 C sed_c4 :: coefficient for sediment model 0136 C sed_d4 :: coefficient for sediment model 0137 C sed_a5 :: coefficient for sediment model 0138 C sed_b5 :: coefficient for sediment model 0139 C sed_c5 :: coefficient for sediment model 0140 C sed_d5 :: coefficient for sediment model 0141 C sed_c :: coefficient for sediment model 0142 COMMON /DARWIN_SED_CONSTANTS_R/ 0143 & sed_a1, 0144 & sed_b1, 0145 & sed_c1, 0146 & sed_d1, 0147 & sed_a2, 0148 & sed_b2, 0149 & sed_c2, 0150 & sed_d2, 0151 & sed_a3, 0152 & sed_b3, 0153 & sed_c3, 0154 & sed_d3, 0155 & sed_a4, 0156 & sed_b4, 0157 & sed_c4, 0158 & sed_d4, 0159 & sed_a5, 0160 & sed_b5, 0161 & sed_c5, 0162 & sed_d5, 0163 & sed_c 0164 _RL sed_a1 0165 _RL sed_b1 0166 _RL sed_c1 0167 _RL sed_d1 0168 _RL sed_a2 0169 _RL sed_b2 0170 _RL sed_c2 0171 _RL sed_d2 0172 _RL sed_a3 0173 _RL sed_b3 0174 _RL sed_c3 0175 _RL sed_d3 0176 _RL sed_a4 0177 _RL sed_b4 0178 _RL sed_c4 0179 _RL sed_d4 0180 _RL sed_a5 0181 _RL sed_b5 0182 _RL sed_c5 0183 _RL sed_d5 0184 _RL sed_c 0185 # endif 8a10163480 Oliv*0186 0187 # ifdef DARWIN_ALLOW_RADIv2 0188 C-- COMMON /DARWIN_SED_CONSTANTS_r/ Coefficients for sediment model 0189 C sed_depth_threshold :: depth threshold for deep/coast equations 5938ff335d Oliv*0190 C sed_globala1 :: coefficient for sediment model 0191 C sed_globalb1 :: coefficient for sediment model 0192 C sed_globalc1 :: coefficient for sediment model 0193 C sed_globald1 :: coefficient for sediment model 0194 C sed_globale1 :: coefficient for sediment model 0195 C sed_globala2 :: coefficient for sediment model 0196 C sed_globalb2 :: coefficient for sediment model 0197 C sed_globalc2 :: coefficient for sediment model 0198 C sed_globald2 :: coefficient for sediment model 0199 C sed_globale2 :: coefficient for sediment model 0200 C sed_globala3 :: coefficient for sediment model 0201 C sed_globalb3 :: coefficient for sediment model 0202 C sed_globalc3 :: coefficient for sediment model 0203 C sed_globald3 :: coefficient for sediment model 0204 C sed_globale3 :: coefficient for sediment model 0205 C sed_globala4 :: coefficient for sediment model 0206 C sed_globalb4 :: coefficient for sediment model 0207 C sed_globalc4 :: coefficient for sediment model 0208 C sed_globald4 :: coefficient for sediment model 0209 C sed_globale4 :: coefficient for sediment model 0210 C sed_globalf4 :: coefficient for sediment model 0211 C sed_globalg4 :: coefficient for sediment model 0212 C sed_globala5 :: coefficient for sediment model 0213 C sed_globalb5 :: coefficient for sediment model 0214 C sed_globalc5 :: coefficient for sediment model 0215 C sed_globald5 :: coefficient for sediment model 0216 C sed_globale5 :: coefficient for sediment model 0217 C sed_globala6 :: coefficient for sediment model 0218 C sed_globalb6 :: coefficient for sediment model 0219 C sed_globalc6 :: coefficient for sediment model 0220 C sed_globald6 :: coefficient for sediment model 0221 C sed_globale6 :: coefficient for sediment model 8a10163480 Oliv*0222 0223 COMMON /DARWIN_SED_CONSTANTS_r/ 0224 & sed_globala1, 0225 & sed_globalb1, 0226 & sed_globalc1, 0227 & sed_globald1, 0228 & sed_globale1, 0229 & sed_globala2, 0230 & sed_globalb2, 0231 & sed_globalc2, 0232 & sed_globald2, 0233 & sed_globale2, 0234 & sed_globala3, 0235 & sed_globalb3, 0236 & sed_globalc3, 0237 & sed_globald3, 0238 & sed_globale3, 0239 & sed_globala4, 0240 & sed_globalb4, 0241 & sed_globalc4, 0242 & sed_globald4, 0243 & sed_globale4, 5938ff335d Oliv*0244 & sed_globalf4, 0245 & sed_globalg4, 8a10163480 Oliv*0246 & sed_globala5, 0247 & sed_globalb5, 0248 & sed_globalc5, 0249 & sed_globald5, 0250 & sed_globale5, 0251 & sed_globala6, 0252 & sed_globalb6, 0253 & sed_globalc6, 0254 & sed_globald6, 0255 & sed_globale6 0256 _RL sed_globala1 0257 _RL sed_globalb1 0258 _RL sed_globalc1 0259 _RL sed_globald1 0260 _RL sed_globale1 0261 _RL sed_globala2 0262 _RL sed_globalb2 0263 _RL sed_globalc2 0264 _RL sed_globald2 0265 _RL sed_globale2 0266 _RL sed_globala3 0267 _RL sed_globalb3 0268 _RL sed_globalc3 0269 _RL sed_globald3 0270 _RL sed_globale3 0271 _RL sed_globala4 0272 _RL sed_globalb4 0273 _RL sed_globalc4 0274 _RL sed_globald4 0275 _RL sed_globale4 5938ff335d Oliv*0276 _RL sed_globalf4 0277 _RL sed_globalg4 8a10163480 Oliv*0278 _RL sed_globala5 0279 _RL sed_globalb5 0280 _RL sed_globalc5 0281 _RL sed_globald5 0282 _RL sed_globale5 0283 _RL sed_globala6 0284 _RL sed_globalb6 0285 _RL sed_globalc6 0286 _RL sed_globald6 0287 _RL sed_globale6 0288 # endif 8fbfd1f382 Oliv*0289 #endif 0290 5910c7e5d2 Oliv*0291 C COMMON /DARWIN_PARAMS_c/ General parameters (same for all plankton) 0292 C darwin_pickupSuff :: pickup suffix for darwin; set to ' ' to disable reading at PTRACERS_Iter0 f17484eeb1 Oliv*0293 COMMON /DARWIN_PARAMS_c/ darwin_pickupSuff 0294 CHARACTER*10 darwin_pickupSuff 5910c7e5d2 Oliv*0295 C darwin_strict_check :: stop instead of issuing warnings 0296 C darwin_linFSConserve :: correct non-conservation due to linear free surface (globally) 0297 C darwin_read_phos :: initial conditions for plankton biomass are in mmol P/m3 0298 C-- COMMON /DARWIN_PARAMS_l/ General parameters (same for all plankton) 8fbfd1f382 Oliv*0299 COMMON /DARWIN_PARAMS_l/ 0300 & darwin_strict_check, 0301 & darwin_linFSConserve, 1342b05049 Oliv*0302 & darwin_read_phos 8fbfd1f382 Oliv*0303 LOGICAL darwin_strict_check 0304 LOGICAL darwin_linFSConserve 0305 LOGICAL darwin_read_phos 5910c7e5d2 Oliv*0306 0307 C-- COMMON /DARWIN_PARAMS_i/ General parameters (same for all plankton) 0308 C darwin_seed :: seed for random number generator (for DARWIN_RANDOM_TRAITS) 0309 C iDEBUG :: index in x dimension for debug prints 0310 C jDEBUG :: index in y dimension for debug prints 0311 C kDEBUG :: index in z dimension for debug prints ba0b6d5d33 Oliv*0312 C DARWIN_disscSelect :: Switch for PIC dissolution rate formulation 0313 C :: 0: fixed dissolution rate Kdissc 0314 C :: 1: Keir 1980 power law b5b44fea51 Oliv*0315 C :: 2: Naviaux et al. 2019 8fbfd1f382 Oliv*0316 COMMON /DARWIN_PARAMS_i/ 0317 & darwin_seed, ba0b6d5d33 Oliv*0318 & darwin_disscSelect, 8fbfd1f382 Oliv*0319 & iDEBUG, 0320 & jDEBUG, 0321 & kDEBUG 0322 INTEGER darwin_seed ba0b6d5d33 Oliv*0323 INTEGER darwin_disscSelect 8fbfd1f382 Oliv*0324 INTEGER iDEBUG 0325 INTEGER jDEBUG 0326 INTEGER kDEBUG 5910c7e5d2 Oliv*0327 0328 C-- COMMON /DARWIN_PARAMS_r/ General parameters (same for all plankton) 0329 C katten_w :: [1/m] atten coefficient water 0330 C katten_chl :: [m2/mg Chl] atten coefficient chl 0331 C 0332 C parfrac :: [] fraction Qsw that is PAR 0333 C parconv :: [uEin/s/W] conversion from W/m2 to uEin/m2/s 0334 C tempnorm :: [] set temperature function (was 1.0) 0335 C TempAeArr :: [K] slope for pseudo-Arrhenius (TEMP_VERSION 2) 0336 C TemprefArr :: [K] reference temp for pseudo-Arrhenius (TEMP_VERSION 2) 0337 C TempCoeffArr :: [] pre-factor for pseudo-Arrhenius (TEMP_VERSION 2) a092808e6b shlo*0338 C TempAeArrMacromol :: [K] slope for pseudo-Arrhenius for macromolecular (TEMP_VERSION 2) 09b16dd232 Oliv*0339 C reminTempAe :: [1/K] temperature coefficient for remineralization (TEMP_VERSION 4) 0340 C mortTempAe :: [1/K] temperature coefficient for linear mortality (TEMP_VERSION 4) 0341 C mort2TempAe :: [1/K] temperature coefficient for quadr. mortality (TEMP_VERSION 4) 0342 C uptakeTempAe :: [1/K] temperature coefficient for uptake (TEMP_VERSION 4) 5910c7e5d2 Oliv*0343 C 0344 C- Iron parameters bbcb16c4fd Oliv*0345 C alpfe :: [] solubility of Fe dust 66ade11a9e Oliv*0346 C scav :: [1/s] fixed iron scavenging rate (#undef DARWIN_PART_SCAV) bbcb16c4fd Oliv*0347 C ligand_tot :: [mol/m3] total ligand concentration 0348 C ligand_stab :: [m3/mol] ligand stability rate ratio 0349 C freefemax :: [mol/m3] max concentration of free iron 66ade11a9e Oliv*0350 C scav_tau :: [1] factor to go from Th scavenging rate to iron 0351 C scav_inter :: [L^e mg^-e s^-1] intercept of scavenging power law (e=scav_exp) 0352 C scav_exp :: [1] exponent of scavenging power law 0353 C scav_POC_wgt :: [g / mmol C] weight POC contributes to POM for scavenging 28ac947a3e Oliv*0354 C scav_PSi_wgt :: [g / mmol Si] weight PSi contributes to POM for scavenging 66ade11a9e Oliv*0355 C scav_PIC_wgt :: [g / mmol C] weight PIC contributes to POM for scavenging 28ac947a3e Oliv*0356 C scav_degrPOM :: [g/m3] concentration of non-labile POM for scavenging 66ade11a9e Oliv*0357 C scav_rat :: [1] factor Th to iron for DARWIN_PART_SCAV_POP 0358 C scav_R_POPPOC :: [mmol P / g C] POP:POC ratio for DARWIN_PART_SCAV_POP bbcb16c4fd Oliv*0359 C depthfesed :: [m] depth above which to add sediment source (was -1000) 0360 C fesedflux :: [mmol Fe /m2/s] fixed iron flux from sediment 0361 C fesedflux_pcm :: [mmol Fe / mmol C] iron input per POC sinking into bottom for DARWIN_IRON_SED_SOURCE_VARIABLE 0362 C fesedflux_min :: [mmol Fe /s] min iron input rate subtracted from fesedflux_pcm*wc_sink*POC 0363 C R_CP_fesed :: [mmol C / mmol P] POC:POP conversion for DARWIN_IRON_SED_SOURCE_POP 2c93eb88ef Oliv*0364 C depthFeVent :: [m] depth below which iron from hydrothermal vents is added 0365 C solFeVent :: [] solubility of iron from hydrothermal vents 0366 C R_FeHe3_vent :: [mmol Fe / mmol He3] Fe:He3 ratio for hydrothermal vents 5910c7e5d2 Oliv*0367 C bbcb16c4fd Oliv*0368 C Knita :: [1/s] ammonia oxidation rate 0369 C Knitb :: [1/s] nitrite oxidation rate 0370 C PAR_oxi :: [uEin/m2/s] critical light level after which oxidation starts 5910c7e5d2 Oliv*0371 C bbcb16c4fd Oliv*0372 C Kdoc :: [1/s] DOC remineralization rate 0373 C Kdop :: [1/s] DON remineralization rate 0374 C Kdon :: [1/s] DOP remineralization rate 0375 C KdoFe :: [1/s] DOFe remineralization rate 0376 C KPOC :: [1/s] POC remineralization rate 0377 C KPON :: [1/s] PON remineralization rate 0378 C KPOP :: [1/s] POP remineralization rate 0379 C KPOFe :: [1/s] POFe remineralization rate 0380 C KPOSi :: [1/s] POSi remineralization rate a092808e6b shlo*0381 C ksatO2remin :: [mmol O2/m3] half-saturation conc. of O2 for remineralization 5910c7e5d2 Oliv*0382 C bbcb16c4fd Oliv*0383 C wC_sink :: [m/s] sinking velocity for POC 0384 C wN_sink :: [m/s] sinking velocity for PON 0385 C wP_sink :: [m/s] sinking velocity for POP 0386 C wFe_sink :: [m/s] sinking velocity for POFe 0387 C wSi_sink :: [m/s] sinking velocity for POSi 0388 C wPIC_sink :: [m/s] sinking velocity for PIC ba0b6d5d33 Oliv*0389 C darwin_KeirCoeff :: [1/s] Keir PIC dissolution rate coefficient 0390 C darwin_KeirExp :: [1] Keir PIC dissolution rate exponent bbcb16c4fd Oliv*0391 C Kdissc :: [1/s] dissolution rate for PIC 5910c7e5d2 Oliv*0392 C 0393 C- Carbon chemistry parameters bbcb16c4fd Oliv*0394 C R_OP :: [mmol O2 / mmol P] O:P ratio for respiration and consumption 0395 C R_OC :: [mmol O2 / mmol C] NOT USED 0396 C m3perkg :: [m3/kg] constant for converting per kg to per m^3 0397 C surfSaltMinInit :: [ppt] minimum salt for carbon solver at initialization 0398 C surfSaltMaxInit :: [ppt] maximum salt for carbon solver at initialization 0399 C surfTempMinInit :: [degrees C] minimum temp for carbon solver at initialization 0400 C surfTempMaxInit :: [degrees C] maximum temp for carbon solver at initialization 0401 C surfDICMinInit :: [mmol C m^-3] minimum DIC for carbon solver at initialization 0402 C surfDICMaxInit :: [mmol C m^-3] maximum DIC for carbon solver at initialization 0403 C surfALKMinInit :: [meq m^-3] minimum alkalinity for carbon solver at initialization 0404 C surfALKMaxInit :: [meq m^-3] maximum alkalinity for carbon solver at initialization 0405 C surfPO4MinInit :: [mmol P m^-3] minimum PO4 for carbon solver at initialization 0406 C surfPO4MaxInit :: [mmol P m^-3] maximum PO4 for carbon solver at initialization 0407 C surfSiMinInit :: [mmol Si m^-3] minimum SiO2 for carbon solver at initialization 0408 C surfSiMaxInit :: [mmol Si m^-3] maximum SiO2 for carbon solver at initialization 0409 C surfSaltMin :: [ppt] minimum salt for carbon solver during run 0410 C surfSaltMax :: [ppt] maximum salt for carbon solver during run 0411 C surfTempMin :: [degrees C] minimum temp for carbon solver during run 0412 C surfTempMax :: [degrees C] maximum temp for carbon solver during run 0413 C surfDICMin :: [mmol C m^-3] minimum DIC for carbon solver during run 0414 C surfDICMax :: [mmol C m^-3] maximum DIC for carbon solver during run 0415 C surfALKMin :: [meq m^-3] minimum alkalinity for carbon solver during run 0416 C surfALKMax :: [meq m^-3] maximum alkalinity for carbon solver during run 0417 C surfPO4Min :: [mmol P m^-3] minimum PO4 for carbon solver during run 0418 C surfPO4Max :: [mmol P m^-3] maximum PO4 for carbon solver during run 0419 C surfSiMin :: [mmol Si m^-3] minimum SiO2 for carbon solver during run 0420 C surfSiMax :: [mmol Si m^-3] maximum SiO2 for carbon solver during run 5910c7e5d2 Oliv*0421 C 0422 C diaz_ini_fac :: reduce tracer concentrations by this factor on initialization 0423 C 0424 C- Denitrification 0425 C O2crit :: [mmol O2 m-3] critical oxygen for O2/NO3 remineralization 0426 C denit_NP :: [mmol N / mmol P] ratio of n to p in denitrification process 0427 C denit_NO3 :: [mmol N / mmol P] ratio of NO3 uptake to phos remineralization in denitrification 0428 C NO3crit :: [mmol N m-3] critical nitrate below which no denit (or remin) happens 0429 C 0430 C- These should probably be traits 0431 C PARmin :: [uEin/m2/s] minimum light for photosynthesis; for non-Geider: 1.0 b3165fc3b9 Oliv*0432 C aphy_chl_ave :: [m2/mg Chl] Chl-specific absorption coefficient 5910c7e5d2 Oliv*0433 C chl2nmax :: [mg Chl / mmol N] max Chl:N ratio for Chl synthesis following Moore 2002 0434 C synthcost :: [mmol C / mmol N] cost of biosynthesis 0435 C inhib_graz :: [(mmol C m-3)-1] inverse decay scale for grazing inhibition 0436 C inhib_graz_exp :: [] exponent for grazing inhibition (0 to turn off inhibition) 0437 C hillnumGraz :: [] exponent for limiting quota uptake in grazing 0438 C hollexp :: [] grazing exponential 1= "Holling 2", 2= "Holling 3" 0439 C phygrazmin :: [mmol C m-3] minimum total prey conc for grazing to occur 0440 C 0441 C- Bacteria 0442 C pmaxDIN :: [1/s] max DIN uptake rate for denitrifying bacteria 0443 C pcoefO2 :: [m3/mmol O2/s] max O2-specific O2 uptake rate for aerobic bacteria 0444 C ksatDIN :: [mmol N m-3] half-saturation conc of dissolved inorganic nitrogen 0445 C alpha_hydrol :: [] increase in POM needed due to hydrolysis 0446 C yod :: [] organic matter yield of aerobic bacteria 0447 C yoe :: [] energy yield of aerobic bacteria 0448 C ynd :: [] organic matter yield of denitrifying bacteria 0449 C yne :: [] energy yield of denitrifying bacteria 0450 C fnh4 :: [] not implemented (for ammonia-oxidizing bacteria) 0451 C ynh4 :: [] not implemented (for ammonia-oxidizing bacteria) 0452 C yonh4 :: [] not implemented (for ammonia-oxidizing bacteria) 0453 C fno2 :: [] not implemented (for nitrite-oxidizing bacteria) 0454 C yno2 :: [] not implemented (for nitrite-oxidizing bacteria) 0455 C yono2 :: [] not implemented (for nitrite-oxidizing bacteria) 0456 C 0457 C- To be implemented 0458 C depthdenit :: [m] not implemented (depth for denitrification relaxation to start) 8fbfd1f382 Oliv*0459 COMMON /DARWIN_PARAMS_r/ 0460 & katten_w, 0461 & katten_chl, 0462 & parfrac, 0463 & parconv, 0464 & tempnorm, 0465 & TempAeArr, 0466 & TemprefArr, 0467 & TempCoeffArr, a092808e6b shlo*0468 & TempAeArrMacromol, 09b16dd232 Oliv*0469 & reminTempAe, 0470 & mortTempAe, 0471 & mort2TempAe, 0472 & uptakeTempAe, 8fbfd1f382 Oliv*0473 & alpfe, 0474 & ligand_tot, 0475 & ligand_stab, 0476 & freefemax, 66ade11a9e Oliv*0477 #ifdef DARWIN_PART_SCAV_POP 8fbfd1f382 Oliv*0478 & scav_rat, 0479 & scav_inter, 0480 & scav_exp, 0481 & scav_R_POPPOC, 66ade11a9e Oliv*0482 #elif defined(DARWIN_PART_SCAV) 0483 & scav_tau, 0484 & scav_inter, 0485 & scav_exp, 0486 & scav_POC_wgt, 28ac947a3e Oliv*0487 & scav_PSi_wgt, 66ade11a9e Oliv*0488 & scav_PIC_wgt, 28ac947a3e Oliv*0489 & scav_degrPOM, 66ade11a9e Oliv*0490 #else 0491 & scav, 0492 #endif 8fbfd1f382 Oliv*0493 & depthfesed, 0494 & fesedflux, 0495 & fesedflux_pcm, ad5342a8d6 Oliv*0496 & fesedflux_min, 8fbfd1f382 Oliv*0497 & R_CP_fesed, 2c93eb88ef Oliv*0498 & depthFeVent, 0499 & solFeVent, 0500 & R_FeHe3_vent, 8fbfd1f382 Oliv*0501 & Knita, 0502 & Knitb, 0503 & PAR_oxi, 0504 & Kdoc, 0505 & Kdop, 0506 & Kdon, 0507 & KdoFe, 0508 & KPOC, 0509 & KPON, 5910c7e5d2 Oliv*0510 & KPOP, 8fbfd1f382 Oliv*0511 & KPOFe, 0512 & KPOSi, a092808e6b shlo*0513 & ksatO2remin, 8fbfd1f382 Oliv*0514 & wC_sink, 0515 & wN_sink, 5910c7e5d2 Oliv*0516 & wP_sink, 8fbfd1f382 Oliv*0517 & wFe_sink, 0518 & wSi_sink, 0519 & wPIC_sink, 0520 & Kdissc, ba0b6d5d33 Oliv*0521 & DARWIN_KeirCoeff, 0522 & DARWIN_KeirExp, 8fbfd1f382 Oliv*0523 #ifdef DARWIN_ALLOW_CARBON 0524 & R_OP, 0525 & R_OC, 0526 & m3perkg, 0527 & surfSaltMinInit, 0528 & surfSaltMaxInit, 0529 & surfTempMinInit, 0530 & surfTempMaxInit, 0531 & surfDICMinInit, 0532 & surfDICMaxInit, 0533 & surfALKMinInit, 0534 & surfALKMaxInit, 0535 & surfPO4MinInit, 0536 & surfPO4MaxInit, 0537 & surfSiMinInit, 0538 & surfSiMaxInit, 0539 & surfSaltMin, 0540 & surfSaltMax, 0541 & surfTempMin, 0542 & surfTempMax, 0543 & surfDICMin, 0544 & surfDICMax, 0545 & surfALKMin, 0546 & surfALKMax, 0547 & surfPO4Min, 0548 & surfPO4Max, 0549 & surfSiMin, 0550 & surfSiMax, 0551 #endif 0552 & diaz_ini_fac, 0553 & O2crit, 0554 & denit_NP, 0555 & denit_NO3, 0556 & NO3crit, 0557 & PARmin, b3165fc3b9 Oliv*0558 & aphy_chl_ave, 8fbfd1f382 Oliv*0559 & chl2nmax, 0560 & synthcost, 0561 & inhib_graz, 0562 & inhib_graz_exp, 0563 & hillnumGraz, 0564 & hollexp, 0565 & phygrazmin, 0566 & pcoefO2, 0567 & pmaxDIN, 0568 & ksatDIN, 0569 & alpha_hydrol, 0570 & yod, 0571 & yoe, 0572 & ynd, 0573 & yne, 5910c7e5d2 Oliv*0574 C & fnh4, 0575 C & ynh4, 0576 C & yonh4, 0577 C & fno2, 0578 C & yno2, 0579 C & yono2, 8fbfd1f382 Oliv*0580 & depthdenit 0581 _RL katten_w 0582 _RL katten_chl 0583 _RL parfrac 0584 _RL parconv 0585 _RL tempnorm 0586 _RL TempAeArr 0587 _RL TemprefArr 0588 _RL TempCoeffArr a092808e6b shlo*0589 _RL TempAeArrMacromol 09b16dd232 Oliv*0590 _RL reminTempAe 0591 _RL mortTempAe 0592 _RL mort2TempAe 0593 _RL uptakeTempAe 8fbfd1f382 Oliv*0594 _RL alpfe 0595 _RL ligand_tot 0596 _RL ligand_stab 0597 _RL freefemax 66ade11a9e Oliv*0598 #ifdef DARWIN_PART_SCAV_POP 8fbfd1f382 Oliv*0599 _RL scav_rat 0600 _RL scav_inter 0601 _RL scav_exp 0602 _RL scav_R_POPPOC 66ade11a9e Oliv*0603 #elif defined(DARWIN_PART_SCAV) 0604 _RL scav_tau 0605 _RL scav_inter 0606 _RL scav_exp 0607 _RL scav_POC_wgt 28ac947a3e Oliv*0608 _RL scav_PSi_wgt 66ade11a9e Oliv*0609 _RL scav_PIC_wgt 28ac947a3e Oliv*0610 _RL scav_degrPOM 66ade11a9e Oliv*0611 #else 0612 _RL scav 0613 #endif 8fbfd1f382 Oliv*0614 _RL depthfesed 0615 _RL fesedflux 0616 _RL fesedflux_pcm ad5342a8d6 Oliv*0617 _RL fesedflux_min 8fbfd1f382 Oliv*0618 _RL R_CP_fesed 2c93eb88ef Oliv*0619 _RL depthFeVent 0620 _RL solFeVent 0621 _RL R_FeHe3_vent 8fbfd1f382 Oliv*0622 _RL Knita 0623 _RL Knitb 0624 _RL PAR_oxi 0625 _RL Kdoc 0626 _RL Kdop 0627 _RL Kdon 0628 _RL KdoFe 0629 _RL KPOC 0630 _RL KPON 5910c7e5d2 Oliv*0631 _RL KPOP 8fbfd1f382 Oliv*0632 _RL KPOFe 0633 _RL KPOSi a092808e6b shlo*0634 _RL ksatO2remin 8fbfd1f382 Oliv*0635 _RL wC_sink 0636 _RL wN_sink 5910c7e5d2 Oliv*0637 _RL wP_sink 8fbfd1f382 Oliv*0638 _RL wFe_sink 0639 _RL wSi_sink 0640 _RL wPIC_sink 0641 _RL Kdissc ba0b6d5d33 Oliv*0642 _RL DARWIN_KeirCoeff 0643 _RL DARWIN_KeirExp 8fbfd1f382 Oliv*0644 #ifdef DARWIN_ALLOW_CARBON 0645 _RL R_OP 0646 _RL R_OC 0647 _RL m3perkg 0648 _RL surfSaltMinInit 0649 _RL surfSaltMaxInit 0650 _RL surfTempMinInit 0651 _RL surfTempMaxInit 0652 _RL surfDICMinInit 0653 _RL surfDICMaxInit 0654 _RL surfALKMinInit 0655 _RL surfALKMaxInit 0656 _RL surfPO4MinInit 0657 _RL surfPO4MaxInit 0658 _RL surfSiMinInit 0659 _RL surfSiMaxInit 0660 _RL surfSaltMin 0661 _RL surfSaltMax 0662 _RL surfTempMin 0663 _RL surfTempMax 0664 _RL surfDICMin 0665 _RL surfDICMax 0666 _RL surfALKMin 0667 _RL surfALKMax 0668 _RL surfPO4Min 0669 _RL surfPO4Max 0670 _RL surfSiMin 0671 _RL surfSiMax 0672 #endif 0673 _RL diaz_ini_fac 0674 _RL O2crit 0675 _RL denit_NP 0676 _RL denit_NO3 0677 _RL NO3crit 0678 _RL PARmin b3165fc3b9 Oliv*0679 _RL aphy_chl_ave 8fbfd1f382 Oliv*0680 _RL chl2nmax 0681 _RL synthcost 0682 _RL inhib_graz 0683 _RL inhib_graz_exp 0684 _RL hillnumGraz 0685 _RL hollexp 0686 _RL phygrazmin 0687 _RL pcoefO2 0688 _RL pmaxDIN 0689 _RL ksatDIN 0690 _RL alpha_hydrol 0691 _RL yod 0692 _RL yoe 0693 _RL ynd 0694 _RL yne 5910c7e5d2 Oliv*0695 C _RL fnh4 0696 C _RL ynh4 0697 C _RL yonh4 0698 C _RL fno2 0699 C _RL yno2 0700 C _RL yono2 8fbfd1f382 Oliv*0701 _RL depthdenit 0702 f0a72e2151 Oliv*0703 #ifdef DARWIN_SOLVESAPHE 0704 C If using Solvesaphe routines (Munhoven, 2013) then in addition, 0705 C selectBTconst :: estimates borate concentration from salinity: 0706 C =1 :: use default formulation of Uppström (1974)(same as S/R CARBON_COEFFS) 0707 C =2 :: use new formulation from Lee et al (2010) 0708 C 0709 C selectFTconst :: estimates fluoride concentration from salinity: 0710 C =1 :: use default formulation of Riley (1965) (same as S/R CARBON_COEFFS) 0711 C =2 :: use new formulation from Culkin (1965) 0712 C 0713 C selectHFconst :: sets the first dissociation constant for hydrogen fluoride: 0714 C =1 :: use default Dickson and Riley (1979) (same as S/R CARBON_COEFFS) 0715 C =2 :: use new formulation of Perez and Fraga (1987) 0716 C 0717 C selectK1K2const :: sets the 1rst & 2nd dissociation constants of carbonic acid: 0718 C =1 :: use default formulation of Millero (1995) with data 0719 C from Mehrbach et al. (1973) (same as S/R CARBON_COEFFS) 0720 C =2 :: use formulation of Roy et al. (1993) 0721 C =3 :: use "combination" formulation of Millero (1995) 0722 C =4 :: use formulation of Luecker et al. (2000) 0723 C =5 :: use formulation of Millero (2010, Mar. Fresh Wat. Res.) 0724 C =6 :: use formulation of Waters, Millero, Woosley (2014, Mar. Chem.) 0725 C selectPHsolver :: sets the pH solver to use: 0726 C =1 :: use the GENERAL solver ; =2 :: use SEC solver ; 0727 C =3 :: use FAST solver routine. 0728 0729 COMMON /DARWIN_SOLVESAPHE_I/ 0730 & at_maxniter, 0731 & selectBTconst,selectFTconst, 0732 & selectHFconst,selectK1K2const, 0733 & selectPHsolver 0734 0735 INTEGER at_maxniter 0736 INTEGER selectBTconst 0737 INTEGER selectFTconst 0738 INTEGER selectHFconst 0739 INTEGER selectK1K2const 0740 INTEGER selectPHsolver 0741 #endif /* DARWIN_SOLVESAPHE */ 0742 8fbfd1f382 Oliv*0743 #ifdef DARWIN_ALLOW_CDOM 5910c7e5d2 Oliv*0744 C-- COMMON /DARWIN_CDOM_PARAMS_r/ 0745 C fracCDOM :: [] fraction of remineralized POP contributing to CDOM 0746 C CDOMdegrd :: [1/s] CDOM degradation rate 0747 C CDOMbleach :: [1/s] CDOM bleaching rate 0748 C PARCDOM :: [uEin/m2/s] PAR where CDOM bleaching becomes maximal 1c72adde5c Oliv*0749 C R_NC_CDOM :: [mmol N / mmol C] CDOM N:C ratio 0750 C R_FeC_CDOM :: [mmol Fe / mmol C] CDOM Fe:C ratio 0751 C R_PC_CDOM :: [mmol P / mmol C] CDOM P:C ratio 5910c7e5d2 Oliv*0752 C R_NP_CDOM :: [mmol N / mmol P] CDOM N:P ratio 0753 C R_FeP_CDOM :: [mmol Fe / mmol P] CDOM Fe:P ratio 0754 C R_CP_CDOM :: [mmol C / mmol P] CDOM C:P ratio 1c72adde5c Oliv*0755 C CDOMcoeff :: [m2 / mmol C] C-specific absorption coefficient of CDOM 0756 C (with #define DARWIN_CDOM_UNITS_CARBON) 5910c7e5d2 Oliv*0757 C CDOMcoeff :: [m2 / mmol P] P-specific absorption coefficient of CDOM 1c72adde5c Oliv*0758 C (with #undef DARWIN_CDOM_UNITS_CARBON) 8fbfd1f382 Oliv*0759 COMMON /DARWIN_CDOM_PARAMS_r/ 0760 & fracCDOM, 0761 & CDOMdegrd, 0762 & CDOMbleach, 0763 & PARCDOM, 1c72adde5c Oliv*0764 # ifdef DARWIN_CDOM_UNITS_CARBON 0765 & R_NC_CDOM, 0766 & R_FeC_CDOM, 0767 & R_PC_CDOM, 0768 # else 8fbfd1f382 Oliv*0769 & R_NP_CDOM, 0770 & R_FeP_CDOM, 0771 & R_CP_CDOM, 1c72adde5c Oliv*0772 # endif 8fbfd1f382 Oliv*0773 & CDOMcoeff 0774 _RL fracCDOM 0775 _RL CDOMdegrd 0776 _RL CDOMbleach 0777 _RL PARCDOM 1c72adde5c Oliv*0778 # ifdef DARWIN_CDOM_UNITS_CARBON 0779 _RL R_NC_CDOM 0780 _RL R_FeC_CDOM 0781 _RL R_PC_CDOM 0782 # else 8fbfd1f382 Oliv*0783 _RL R_NP_CDOM 0784 _RL R_FeP_CDOM 0785 _RL R_CP_CDOM 1c72adde5c Oliv*0786 # endif 8fbfd1f382 Oliv*0787 _RL CDOMcoeff 0788 #endif 0789 5910c7e5d2 Oliv*0790 C-- COMMON /DARWIN_DEPENDENT_PARAMS_i/ 2c93eb88ef Oliv*0791 C laCDOM :: index of reference waveband for CDOM absorption spectrum 0792 C kMinFeSed :: minimum level index for iron sedimentation 0793 C kMaxFeSed :: maximum level index for iron sedimentation 0794 C kMinFeVent :: minimum level index for hydrothermal vents 8fbfd1f382 Oliv*0795 COMMON /DARWIN_DEPENDENT_PARAMS_i/ 0796 & darwin_dependent_i_dummy, 0797 #ifdef ALLOW_RADTRANS 0798 #ifdef DARWIN_ALLOW_CDOM 0799 #else 0800 & laCDOM, 0801 #endif 0802 #endif 0803 & kMinFeSed, 2c93eb88ef Oliv*0804 & kMaxFeSed, 0805 & kMinFeVent 8fbfd1f382 Oliv*0806 INTEGER darwin_dependent_i_dummy 0807 #ifdef ALLOW_RADTRANS 0808 #ifdef DARWIN_ALLOW_CDOM 0809 #else 0810 INTEGER laCDOM 0811 #endif 0812 #endif 0813 INTEGER kMinFeSed 0814 INTEGER kMaxFeSed 2c93eb88ef Oliv*0815 INTEGER kMinFeVent 8fbfd1f382 Oliv*0816 0817 0818 #endif /* ALLOW_DARWIN */ 0819
| [ Source navigation ] | [ Diff markup ] | [ Identifier search ] | [ general search ] |
|
This page was automatically generated from https://github.com/darwinproject/darwin3 by the 2.3.7-MITgcm-0.1 LXR engine. The LXR team |
|