![]() |
|
|||
File indexing completed on 2024-08-30 05:10:52 UTC
view on githubraw file Latest commit ae2be615 on 2024-08-29 19:00:27 UTC3c9c51750e Jean*0001 C *==========================================================* 1507ef23b7 Davi*0002 C | DIC_VARS.h cffdbccbc4 Jona*0003 C | o Abiotic Carbon Variables 3c9c51750e Jean*0004 C *==========================================================* 1507ef23b7 Davi*0005 cffdbccbc4 Jona*0006 C AtmospCO2 :: Atmospheric pCO2 (atm). 0007 C AtmosP :: Atmospheric Pressure loaded from file (atm). 2e3e8c330d Jona*0008 C pH :: surface ocean pH (acidity) for pCO2 calculations. cffdbccbc4 Jona*0009 C pCO2 :: surface ocean partial pressure of CO2 (atm). 0010 C FluxCO2 :: Air-sea flux of CO2 (mol/m2/s). 6acab690ae Jona*0011 C wind :: Wind speed loaded from file for air-sea 02609aa0cd Jona*0012 C flux calculations (m/s). a1d0e455fd Hann*0013 C fIce :: Fraction of sea ice cover loaded from file 3cf5faaa8f Jona*0014 C (or set by thice/seaice) 02609aa0cd Jona*0015 C for air-sea flux calculations. 6acab690ae Jona*0016 C Kwexch_Pre :: Common part of piston velocity used for cffdbccbc4 Jona*0017 C for air-sea CO2 and O2 flux calculations. 2e3e8c330d Jona*0018 C silicaSurf :: Surface ocean concentration of silicate for 0019 C pCO2 calculations. Read in from file (mol/m3). 0020 C zca :: Scale depth for CaCO3 remineralization [m] 0021 C useCalciteSaturation :: Dissolve calcium carbonate only below saturation 0022 C horizon (needs DIC_CALCITE_SAT to be defined) 0023 C calcOmegaCalciteFreq :: Frequency at which 3D calcite saturation state, 0024 C omegaC, is updated (s). 0025 C nIterCO3 :: Number of iterations of the Follows 3D pH solver to 0026 C calculate deep carbonate ion concenetration (no 0027 C effect when using the Munhoven/SolveSapHe solvers). a2c35952f2 Jona*0028 C calciteDissolRate :: Rate constant (%) for calcite dissolution 0029 C from Keir (1980) Geochem. Cosmochem. Acta. 0030 C calciteDissolExp :: Rate exponent for calcite dissolution 0031 C from Keir (1980) Geochem. Cosmochem. Acta. 2e3e8c330d Jona*0032 C WsinkPIC :: sinking speed (m/s) of particulate inorganic carbon for 0033 C calculation of calcite dissolution through the watercolumn 0034 C selectCalciteBottomRemin :: to either remineralize in bottom or top layer 0035 C if flux reaches bottom layer; =0 : bottom, =1 : top a2c35952f2 Jona*0036 C selectCalciteDissolution :: flag to control calcite dissolution rate method: 0037 C =0 : Constant dissolution rate; 0038 C =1 : Follows (default) ; 0039 C =2 : Keir (1980) Geochem. Cosmochem. Acta. ; 0040 C =3 : Naviaux et al. 2019, Marine Chemistry 51e381e9c9 Jean*0041 C pH_isLoaded(1) :: = T when surface pH is loaded from pickup file 0042 C pH_isLoaded(2) :: = T when 3-D pH is loaded from pickup file 2e3e8c330d Jona*0043 1507ef23b7 Davi*0044 COMMON /CARBON_NEEDS/ 0045 & AtmospCO2, AtmosP, pH, pCO2, FluxCO2, 2e3e8c330d Jona*0046 & wind, fIce, Kwexch_Pre, silicaSurf, a2c35952f2 Jona*0047 & calciteDissolRate, calciteDissolExp, 0048 & calcOmegaCalciteFreq, zca, 0049 & WsinkPIC, selectCalciteBottomRemin, 51e381e9c9 Jean*0050 & selectCalciteDissolution, nIterCO3, 0051 & useCalciteSaturation, pH_isLoaded 2e3e8c330d Jona*0052 1507ef23b7 Davi*0053 _RL AtmospCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0054 _RL AtmosP(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0055 _RL pCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) a1d0e455fd Hann*0056 _RL pH(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 1507ef23b7 Davi*0057 _RL FluxCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0058 _RL wind(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) a1d0e455fd Hann*0059 _RL fIce(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 1507ef23b7 Davi*0060 _RL Kwexch_Pre(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) a1d0e455fd Hann*0061 _RL silicaSurf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) a2c35952f2 Jona*0062 _RL calciteDissolRate(2), calciteDissolExp(2) 2e3e8c330d Jona*0063 _RL calcOmegaCalciteFreq a2c35952f2 Jona*0064 _RL zca 2e3e8c330d Jona*0065 _RL WsinkPIC 0066 INTEGER selectCalciteBottomRemin a2c35952f2 Jona*0067 INTEGER selectCalciteDissolution 2e3e8c330d Jona*0068 INTEGER nIterCO3 0069 LOGICAL useCalciteSaturation 51e381e9c9 Jean*0070 LOGICAL pH_isLoaded(2) 2e3e8c330d Jona*0071 0072 #ifdef DIC_CALCITE_SAT 0073 C silicaDeep :: 3D-field of silicate concentration for pH and 0074 C carbonate calculations. Read in from file (mol/m3). 0075 C omegaC :: Local saturation state with respect to calcite 0076 COMMON /DIC_CALCITE_SAT_FIELDS/ ae2be6150b Jona*0077 & pH3D, silicaDeep, omegaC 0078 _RL pH3D (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) a1d0e455fd Hann*0079 _RL silicaDeep(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) 2e3e8c330d Jona*0080 _RL omegaC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) 0081 #endif 1507ef23b7 Davi*0082 6acab690ae Jona*0083 C Store dissociation and carbon chemistry coefficients for 0084 C pCO2 solvers (see carbon_chem.F). 0085 C ak0 = [H2CO2]/pCO2 cffdbccbc4 Jona*0086 C [Weiss 1974] 0087 C ak1 = [H][HCO3]/[H2CO3] 0088 C [Millero p.664 (1995) using Mehrbach et al. data on seawater scale] 0089 C ak2 = [H][CO3]/[HCO3] 0090 C [Millero p.664 (1995) using Mehrbach et al. data on seawater scale] 0091 C akw = [H][OH] 0092 C [Millero p.670 (1995) using composite data] 0093 C akb = [H][BO2]/[HBO2] 0094 C [Millero p.669 (1995) using data from dickson (1990)] 0095 C aks = [H][SO4]/[HSO4] 0096 C [dickson (1990, J. chem. Thermodynamics 22, 113)] 0097 C akf = [H][F]/[HF] 0098 C [dickson and Riley (1979)] 0099 C ak1p = [H][H2PO4]/[H3PO4] 0100 C [DOE(1994) eq 7.2.20 with footnote using data from Millero (1974)] 0101 C ak2p = [H][HPO4]/[H2PO4] 0102 C [DOE(1994) eq 7.2.23 with footnote using data from Millero (1974)] 0103 C ak3p = [H][PO4]/[HPO4] 0104 C [DOE(1994) eq 7.2.26 with footnote using data from Millero (1974)] 0105 C aksi = [H][SiO(OH)3]/[Si(OH)4] 0106 C [Millero p.671 (1995) using data from Yao and Millero (1995)] 6acab690ae Jona*0107 C ft = estimated fluoride concentration cffdbccbc4 Jona*0108 C [Riley (1965)] 6acab690ae Jona*0109 C st = estimated sulphate concentration cffdbccbc4 Jona*0110 C [Morris & Riley (1966)] 6acab690ae Jona*0111 C bt = estimated borate concentration cffdbccbc4 Jona*0112 C [Uppstrom (1974)] 6acab690ae Jona*0113 C fugf :: correct for non-ideality in ocean cffdbccbc4 Jona*0114 C [Weiss (1974) Marine Chemistry] 6acab690ae Jona*0115 C ff :: used for water vapor and pressure correction cffdbccbc4 Jona*0116 C [Weiss & Price (1980, Mar. Chem., 8, 347-359; Eq 13 with table 6 values)] 6acab690ae Jona*0117 C Ksp_TP_Calc :: solubility product for calcite cffdbccbc4 Jona*0118 C [Following Mucci (1983) with pressure dependence from Ingle (1975)] 6acab690ae Jona*0119 1507ef23b7 Davi*0120 COMMON /CARBON_CHEM/ 0121 & ak0,ak1,ak2,akw,akb,aks,akf, 13adbda1b1 Jean*0122 & ak1p,ak2p,ak3p,aksi, fugf, 1507ef23b7 Davi*0123 & ff,ft,st,bt, Ksp_TP_Calc 6acab690ae Jona*0124 1507ef23b7 Davi*0125 _RL ak0(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0126 _RL ak1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0127 _RL ak2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0128 _RL akw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0129 _RL akb(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0130 _RL aks(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0131 _RL akf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0132 _RL ak1p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0133 _RL ak2p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0134 _RL ak3p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0135 _RL aksi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0136 _RL ff(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) d0092a57ac Step*0137 C Fugacity Factor added by Val Bennington Nov. 2010 0138 _RL fugf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 1507ef23b7 Davi*0139 _RL ft(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0140 _RL st(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0141 _RL bt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0142 _RL Ksp_TP_Calc(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0143 6acab690ae Jona*0144 #ifdef CARBONCHEM_SOLVESAPHE 0145 C If using Solvesaphe routines (Munhoven, 2013) then in adittion, 0146 C cat = total calcium concentration, Ref.: Culkin (1965) 0147 C akn = the dissociation constant of ammonium [H][NH3]/[NH4] 0148 C Ref.: Yao and Millero (1995) 0149 C akhs = the dissociation constant of hydrogen sulfide [H][HS]/[H2S] 0150 C Ref.: Millero et al. (1988) 0151 C (cited by Millero (1995) and Yao and Millero (1995)) 0152 C aphscale = pH scale conversion factor ; convert from the total to the free 0153 C scale for solvesaphe calculations; Ref.: Munhoven, 2013 0154 C Ksp_TP_Arag = solubility product for aragonite, Ref.: Mucci (1983) 0155 C---- 0156 C selectBTconst :: estimates borate concentration from salinity: 0157 C =1 :: use default formulation of Uppström (1974)(same as S/R CARBON_COEFFS) 0158 C =2 :: use new formulation from Lee et al (2010) 0159 C 0160 C selectFTconst :: estimates fluoride concentration from salinity: 0161 C =1 :: use default formulation of Riley (1965) (same as S/R CARBON_COEFFS) 0162 C =2 :: use new formulation from Culkin (1965) 0163 C 0164 C selectHFconst :: sets the first dissociation constant for hydrogen fluoride: 0165 C =1 :: use default Dickson and Riley (1979) (same as S/R CARBON_COEFFS) 0166 C =2 :: use new formulation of Perez and Fraga (1987) 0167 C 0168 C selectK1K2const :: sets the 1rst & 2nd dissociation constants of carbonic acid: 0169 C =1 :: use default formulation of Millero (1995) with data 0170 C from Mehrbach et al. (1973) (same as S/R CARBON_COEFFS) 0171 C =2 :: use formulation of Roy et al. (1993) 0172 C =3 :: use "combination" formulation of Millero (1995) 0173 C =4 :: use formulation of Luecker et al. (2000) 0174 C =5 :: use formulation of Millero (2010, Mar. Fresh Wat. Res.) 0175 C =6 :: use formulation of Waters, Millero, Woosley (2014, Mar. Chem.) 0176 C selectPHsolver :: sets the pH solver to use: 0177 C =1 :: use the GENERAL solver ; =2 :: use SEC solver ; 0178 C =3 :: use FAST solver routine. 0179 0180 COMMON /CARBONCHEM_SOLVESAPHE/ 0181 & cat, akn, akhs, aphscale, Ksp_TP_Arag, 0182 & at_maxniter, 0183 & selectBTconst,selectFTconst, 0184 & selectHFconst,selectK1K2const, 0185 & selectPHsolver 0186 0187 _RL cat(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0188 _RL akn(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0189 _RL akhs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0190 _RL aphscale(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0191 _RL Ksp_TP_Arag(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0192 0193 INTEGER at_maxniter 0194 INTEGER selectBTconst 0195 INTEGER selectFTconst 0196 INTEGER selectHFconst 0197 INTEGER selectK1K2const 0198 INTEGER selectPHsolver 0199 #endif /* CARBONCHEM_SOLVESAPHE */ 0200 cffdbccbc4 Jona*0201 C Store dissociation coefficients for O2 solubility 1507ef23b7 Davi*0202 COMMON /OXYGEN_CHEM/ 0203 & oA0,oA1,oA2,oA3,oA4,oA5, 0204 & oB0,oB1,oB2,oB3, 946d3aa393 Jean*0205 & oC0 1507ef23b7 Davi*0206 _RL oA0,oA1,oA2,oA3,oA4,oA5 0207 _RL oB0,oB1,oB2,oB3 0208 _RL oC0 0209 cffdbccbc4 Jona*0210 C Store surface mean tracer values for virtual FW fluxes 1507ef23b7 Davi*0211 C permil : is conversion factor for mol/m3 to mol/kg 0212 C assumes uniform (surface) density 0213 C Pa2Atm : for conversion of atmospheric pressure 946d3aa393 Jean*0214 C when coming from atmospheric model 1507ef23b7 Davi*0215 COMMON /GLOBAL_SURF_MEAN/ 0216 & gsm_alk,gsm_s,gsm_t,gsm_dic, 0217 & gsm_c14,permil,Pa2Atm 0218 _RL gsm_alk 0219 _RL gsm_s 0220 _RL gsm_t 0221 _RL gsm_DIC 0222 _RL gsm_C14 0223 _RL permil 0224 _RL Pa2Atm 0225 cffdbccbc4 Jona*0226 C Schmidt number coefficients for air sea gas exchange of 0227 C carbon and oxygen 946d3aa393 Jean*0228 COMMON /DIC_SCHMIDT_NO/ 1507ef23b7 Davi*0229 & sca1, sca2, sca3, sca4, 0230 & sox1, sox2, sox3, sox4 0231 _RL sca1 0232 _RL sca2 0233 _RL sca3 0234 _RL sca4 0235 _RL sox1 0236 _RL sox2 0237 _RL sox3 0238 _RL sox4 0239 254eed5bf2 Davi*0240 C-- COMMON /DIC_FILENAMES/ 946d3aa393 Jean*0241 C DIC_windFile :: file name of wind speeds 0242 C DIC_atmospFile :: file name of atmospheric pressure 0243 C DIC_iceFile :: file name of seaice fraction 0244 C DIC_ironFile :: file name of aeolian iron flux 0245 C DIC_silicaFile :: file name of surface silica 2e3e8c330d Jona*0246 C DIC_deepSilicaFile :: file name of 3D silica fields 10167de8af Step*0247 C DIC_parFile :: file name of photosynthetically available radiation (PAR) 13adbda1b1 Jean*0248 C DIC_chlaFile :: file name of chlorophyll climatology 946d3aa393 Jean*0249 C DIC_forcingPeriod :: periodic forcing parameter specific for dic (seconds) 0250 C DIC_forcingCycle :: periodic forcing parameter specific for dic (seconds) 6acab690ae Jona*0251 C dic_int* :: Handle the atmospheric boundary condition for pCO2 cffdbccbc4 Jona*0252 C dic_int1: 0253 C 0=use default pCO2 (278.d-6) 0254 C 1=use constant value of dic_pCO2 read in from data.dic 0255 C 2=read in from file named co2atmos.dat 6acab690ae Jona*0256 C 3=use well mixed atmospheric box (can use dic_pCO2 as initial value) cffdbccbc4 Jona*0257 C dic_int2 :: number pCO2 entries to read from file 0258 C dic_int3 :: start timestep 0259 C dic_int4 :: timestep between file entries 0260 C dic_pCO2 :: atmospheric pCO2 to be read from data.dic 254eed5bf2 Davi*0261 COMMON /DIC_FILENAMES/ a1d0e455fd Hann*0262 & DIC_windFile, DIC_atmospFile, DIC_silicaFile, 2e3e8c330d Jona*0263 & DIC_deepSilicaFile, 0264 & DIC_iceFile, DIC_parFile, a1d0e455fd Hann*0265 & DIC_chlaFile, DIC_ironFile, 946d3aa393 Jean*0266 & DIC_forcingPeriod, DIC_forcingCycle, 0267 & dic_pCO2, dic_int1, dic_int2, dic_int3, dic_int4 0268 0269 CHARACTER*(MAX_LEN_FNAM) DIC_windFile 0270 CHARACTER*(MAX_LEN_FNAM) DIC_atmospFile 0271 CHARACTER*(MAX_LEN_FNAM) DIC_silicaFile a1d0e455fd Hann*0272 CHARACTER*(MAX_LEN_FNAM) DIC_deepSilicaFile 0273 CHARACTER*(MAX_LEN_FNAM) DIC_iceFile 10167de8af Step*0274 CHARACTER*(MAX_LEN_FNAM) DIC_parFile 13adbda1b1 Jean*0275 CHARACTER*(MAX_LEN_FNAM) DIC_chlaFile a1d0e455fd Hann*0276 CHARACTER*(MAX_LEN_FNAM) DIC_ironFile 946d3aa393 Jean*0277 _RL DIC_forcingPeriod 0278 _RL DIC_forcingCycle cffdbccbc4 Jona*0279 _RL dic_pCO2 254eed5bf2 Davi*0280 INTEGER dic_int1 0281 INTEGER dic_int2 0282 INTEGER dic_int3 0283 INTEGER dic_int4 0284 1507ef23b7 Davi*0285 #ifdef DIC_BIOTIC 3c9c51750e Jean*0286 C *==========================================================* 1507ef23b7 Davi*0287 C | o Biological Carbon Variables 3c9c51750e Jean*0288 C *==========================================================* 1507ef23b7 Davi*0289 6acab690ae Jona*0290 C Values for biotic biogeochemistry 0291 C (many set in data.dic, defaults to values in dic_readparms.F) 0292 C par :: photosynthetically active radiation (light available 0293 C for phytoplankton growth) [W/m2] 0294 C alpha :: maximum rate of biological activity [mol P/m3/s] 0295 C alphaUniform :: read in alphaUniform to fill in 2d array alpha 0296 C rain_ratio :: inorganic/organic particulate carbon rain ratio (PIC/POC) 0297 C rainRatioUniform :: read in rainRatioUniform to fill in 2d array rain_ratio 0298 C InputFe :: aeolian deposition of TOTAL IRON in dust [mol/m2/s] 0299 C CHL :: chlorophyll climatology data for self-shading effect [mg/m3] 0300 C Kpo4, KFE, lit0 :: half saturation constants for phosphate [mol P/m3], 0301 C iron [mol Fe/m3] and light [W/m2] 0302 C DOPfraction :: fraction of new production going to DOP 0303 C zcrit :: Minimum Depth (m) over which biological activity is computed 0304 C nlev :: level index just below -zcrit 0305 C KRemin :: remineralization power law coeffient 0306 C KDOPremin :: DOP remineralization rate [1/s] 0307 C R_op, R_cp :: stochiometric ratios of nutrients 0308 C R_NP, R_FeP :: stochiometric ratios of nutrients 0309 C (assumption of stoichometry of plankton and 0310 C particulate and dissolved organic matter) 0311 C O2crit :: critical oxygen level [mol/m3] 0312 C alpfe :: solubility of aeolian fe [fraction] 0313 C KScav :: iron scavenging rate [1/s] 0314 C ligand_stab :: ligand-free iron stability constant [m3/mol] 0315 C ligand_tot :: uniform, invariant total free ligand conc [mol/m3] 0316 C freefemax :: max soluble free iron [mol/m3] 0317 C fesedflux_pcm :: ratio of sediment iron to sinking organic matter 0318 C FeIntSec :: Sediment Fe flux, intersect value in: 0319 C Fe_flux = fesedflux_pcm*pflux + FeIntSec 0320 C parfrac :: fraction of Qsw that is PAR 0321 C k0 :: light attentuation coefficient of water [1/m] 0322 C kchl :: light attentuation coefficient of chlorophyll [m2/mg] 0323 C alphamax,alphamin :: not used (legacy adjoint param) 0324 C calpha :: not used (legacy adjoint param) 0325 C crain_ratio :: not used (legacy adjoint param) 0326 C cInputFe :: not used (legacy adjoint param) 0327 C calpfe :: not used (legacy adjoint param) 0328 C feload :: not used (legacy adjoint param) 0329 C cfeload :: not used (legacy adjoint param) 0330 C QSW_underice :: is Qsw is masked by ice fraction? 0331 0332 COMMON /BIOTIC_NEEDS/ 2e3e8c330d Jona*0333 & par, alpha, rain_ratio, InputFe, CHL, 1507ef23b7 Davi*0334 & Kpo4, DOPfraction, zcrit, KRemin, 2e3e8c330d Jona*0335 & KDOPremin, R_op, R_cp, R_NP, R_FeP, 1507ef23b7 Davi*0336 & O2crit, alpfe, KScav, ligand_stab, ligand_tot, KFE, 13adbda1b1 Jean*0337 & freefemax, fesedflux_pcm, FeIntSec, 0338 & parfrac, k0, kchl, lit0, 7dccb58cca Jean*0339 & alphaUniform, rainRatioUniform, 1507ef23b7 Davi*0340 & nlev, QSW_underice 6acab690ae Jona*0341 c & alphamax, alphamin, 0342 c & calpha, crain_ratio, cInputFe, calpfe, feload, cfeload, 1507ef23b7 Davi*0343 0344 _RL par(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0345 _RL alpha(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0346 _RL rain_ratio(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0347 _RL InputFe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 13adbda1b1 Jean*0348 _RL CHL(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 946d3aa393 Jean*0349 _RL Kpo4 0350 _RL DOPfraction 1507ef23b7 Davi*0351 _RL zcrit 0352 _RL KRemin 0353 _RL KDOPremin 0354 _RL R_op 0355 _RL R_cp 0356 _RL R_NP 0357 _RL R_FeP 0358 _RL O2crit 0359 _RL alpfe 13adbda1b1 Jean*0360 _RL fesedflux_pcm 0361 _RL FeIntSec 1507ef23b7 Davi*0362 _RL KScav 0363 _RL ligand_stab 0364 _RL ligand_tot 2e3e8c330d Jona*0365 _RL KFe 1507ef23b7 Davi*0366 _RL freefemax 13adbda1b1 Jean*0367 _RL k0, kchl, parfrac, lit0 7dccb58cca Jean*0368 _RL alphaUniform 0369 _RL rainRatioUniform 6acab690ae Jona*0370 c _RL alphamax, alphamin 0371 c _RL calpha 0372 c _RL crain_ratio 0373 c _RL cInputFe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0374 c _RL calpfe 0375 c _RL cfeload 0376 c _RL feload(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 1507ef23b7 Davi*0377 6acab690ae Jona*0378 INTEGER nlev 1507ef23b7 Davi*0379 LOGICAL QSW_underice 0380 2e3e8c330d Jona*0381 #ifdef ALLOW_TIMEAVE 0382 C For time-averages output using TIMEAVE pkg: 0383 C Note: Although some of the Time-Ave arrays are for A-Biotic fields, all these 0384 C arrays are only used + filled within ifdef DIC_BIOTIC code, so we keep 0385 C them only if DIC_BIOTIC is defined. 0386 C BIOave :: biological productivity [mol P/m3/s] 0387 C CARave :: carbonate changes due to biological productivity 0388 C and remineralization [mol C/m3/s] 0389 C SURave :: tendency of DIC due to air-sea exchange 0390 C and virtual flux [mol C/m3/s] 0391 C SUROave :: tendency of O2 due to air-sea exchange [mol O2/m3/s] 0392 C pCO2ave :: surface ocean pCO2 [uatm] 0393 C pHave :: surface ocean pH 0394 C fluxCO2ave :: Air-sea flux of CO2 [mol C/m2/s] 0395 C pfluxave :: changes to PO4 due to flux and remineralization [mol P/m3/s] 0396 C epfluxave :: export flux of PO4 through each layer [mol P/m2/s] 0397 C cfluxave :: carbonate changes due to flux and remineralization [mol C/m3/s] 0398 C DIC_timeAve :: period over which DIC averages are calculated [s] 0399 0400 COMMON /BIOTIC_TAVE/ 0401 & BIOave, CARave, SURave, SUROave, pCO2ave, pHave, 0402 & fluxCO2ave, pfluxave, epfluxave, cfluxave, 0403 & DIC_timeAve 0404 0405 _RL BIOave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) 0406 _RL CARave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) 0407 _RL SURave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0408 _RL SUROave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0409 _RL pCO2ave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0410 _RL pHave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0411 _RL fluxCO2ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) 0412 _RL pfluxave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) 0413 _RL epfluxave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) 0414 _RL cfluxave (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy) 0415 _RL DIC_timeAve(nSx,nSy) 0416 #endif /* ALLOW_TIMEAVE */ 51e381e9c9 Jean*0417 2e3e8c330d Jona*0418 #endif /* DIC_BIOTIC */
[ Source navigation ] | [ Diff markup ] | [ Identifier search ] | [ general search ] |
This page was automatically generated from https://github.com/MITgcm/MITgcm by the 2.2.1-MITgcm-0.1 LXR engine. The LXR team |
![]() ![]() |