File indexing completed on 2024-07-04 05:10:17 UTC
view on githubraw file Latest commit e37161e0 on 2024-07-03 22:04:10 UTC
c0d1c06c15 Matt*0001
0002
0003
0004
0005 COMMON /CARBON_NEEDS/
a284455135 Matt*0006 & apCO2, AtmosP, pH, pCO2, FluxCO2,
c0d1c06c15 Matt*0007 & wind, FIce, Silica
e0f9a7ba0b Matt*0008 #ifdef ALLOW_EXF
0009 & , apco20, apco21
be72e7ae9e Jean*0010 #endif
a284455135 Matt*0011 _RL apCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
c0d1c06c15 Matt*0012 _RL AtmosP(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0013 _RL pH(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
0014 _RL pCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0015 _RL FluxCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0016 _RL wind(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0017 _RL FIce(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0018 _RL Silica(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
e0f9a7ba0b Matt*0019 #ifdef ALLOW_EXF
0020 _RL apco20(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0021 _RL apco21(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
be72e7ae9e Jean*0022 #endif
c0d1c06c15 Matt*0023
0024
0025
0026
0027
0028 COMMON /CARBON_CHEM/
0029 & ak0,ak1,ak2,akw,akb,aks,akf,
e0f9a7ba0b Matt*0030 & ak1p,ak2p,ak3p,aksi,fugf,
be72e7ae9e Jean*0031 & ff,ft,st,bt,
c0d1c06c15 Matt*0032 & Ksp_TP_Calc,Ksp_TP_Arag
0033 _RL ak0(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0034 _RL ak1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0035 _RL ak2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0036 _RL akw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0037 _RL akb(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0038 _RL aks(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0039 _RL akf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0040 _RL ak1p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0041 _RL ak2p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0042 _RL ak3p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0043 _RL aksi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0044 _RL ff(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0045 _RL fugf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0046 _RL ft(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0047 _RL st(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0048 _RL bt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0049 _RL Ksp_TP_Calc(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0050 _RL Ksp_TP_Arag(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0051
e0f9a7ba0b Matt*0052 #ifdef CARBONCHEM_SOLVESAPHE
0053
0054
0055
0056
0057
0058
0059
0060
0061
0062
0063
0064
0065
0066
0067
0068
0069
0070
0071
0072
0073
0074
0075
0076
0077
0078
0079
0080
0081
0082
0083
0084
0085
0086
0087
0088 COMMON /CARBONCHEM_SOLVESAPHE_ARIANE/
0089
0090 & cat, akn, akhs, aphscale,
0091 & at_maxniter,
0092 & selectBTconst,selectFTconst,
0093 & selectHFconst,selectK1K2const,
0094 & selectPHsolver
0095
0096 _RL cat(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0097 _RL akn(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0098 _RL akhs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0099 _RL aphscale(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0100
0101
0102 INTEGER at_maxniter
0103 INTEGER selectBTconst
0104 INTEGER selectFTconst
0105 INTEGER selectHFconst
0106 INTEGER selectK1K2const
0107 INTEGER selectPHsolver
0108 #endif /* CARBONCHEM_SOLVESAPHE */
0109
c0d1c06c15 Matt*0110 COMMON /OXYGEN_CHEM/
0111 & oA0,oA1,oA2,oA3,oA4,oA5,
0112 & oB0,oB1,oB2,oB3,
0113 & oC0
0114 _RL oA0,oA1,oA2,oA3,oA4,oA5
0115 _RL oB0,oB1,oB2,oB3
0116 _RL oC0
0117
0118 COMMON /GLOBAL_SURF_MEAN/
e0f9a7ba0b Matt*0119 & permil,Pa2Atm,epsln
c0d1c06c15 Matt*0120
0121
0122
0123
0124
0125
0126 _RL permil
0127 _RL Pa2Atm
e0f9a7ba0b Matt*0128 _RL epsln
c0d1c06c15 Matt*0129
0130 COMMON /SCHMIDT_NO/
0131 & sca1, sca2, sca3, sca4,
0132 & sox1, sox2, sox3, sox4
0133
0134
0135
0136 _RL sca1,sca2,sca3,sca4
0137 _RL sox1,sox2,sox3,sox4
0138
e0f9a7ba0b Matt*0139 #ifdef ALLOW_EXF
079948e6a6 Matt*0140 integer apco2startdate1
0141 integer apco2startdate2
be72e7ae9e Jean*0142 _RL apco2StartTime
079948e6a6 Matt*0143 _RL apco2period
041f4605e9 Jean*0144 _RL apco2RepCycle
079948e6a6 Matt*0145 _RL apco2const
0146 _RL apco2_exfremo_intercept
0147 _RL apco2_exfremo_slope
e0f9a7ba0b Matt*0148 CHARACTER*1 apco2mask
079948e6a6 Matt*0149 #endif
c0d1c06c15 Matt*0150
0151
0152
0153
0154
0155 COMMON /BLING_INPUTS/
0156 & bling_windFile, bling_atmospFile, bling_iceFile,
54b6d6aa91 Matt*0157 & bling_ironFile, bling_silicaFile,
be72e7ae9e Jean*0158 & bling_psmFile, bling_plgFile, bling_pdiazFile,
c0d1c06c15 Matt*0159 & bling_forcingPeriod, bling_forcingCycle,
be72e7ae9e Jean*0160 & bling_pCO2,
e0f9a7ba0b Matt*0161 & river_conc_po4, river_dom_to_nut,
47de7c1e0e Matt*0162 & bling_Pc_2dFile, bling_Pc_2d_diazFile,
af8f3b9ffc Matt*0163 & bling_alpha_photo2dFile,bling_phi_DOM2dFile,
0164 & bling_k_Fe2dFile, bling_k_Fe_diaz2dFile,
0165 & bling_gamma_POM2dFile, bling_wsink0_2dFile,
0166 & bling_phi_lg2dFile, bling_phi_sm2dFile
e0f9a7ba0b Matt*0167 #ifdef ALLOW_EXF
079948e6a6 Matt*0168 & ,apco2startdate1,apco2startdate2,
041f4605e9 Jean*0169 & apco2StartTime, apco2period, apco2RepCycle,
079948e6a6 Matt*0170 & apco2const,
0171 & apco2_exfremo_intercept,
0172 & apco2_exfremo_slope,
0173 & apco2file, apco2mask
0174 #endif
c0d1c06c15 Matt*0175
0176
0177
0178
0179
0180
54b6d6aa91 Matt*0181
0182
0183
c0d1c06c15 Matt*0184
be72e7ae9e Jean*0185
c0d1c06c15 Matt*0186
0187
079948e6a6 Matt*0188
c0d1c06c15 Matt*0189
0190 CHARACTER*(MAX_LEN_FNAM) bling_windFile
0191 CHARACTER*(MAX_LEN_FNAM) bling_atmospFile
0192 CHARACTER*(MAX_LEN_FNAM) bling_iceFile
0193 CHARACTER*(MAX_LEN_FNAM) bling_ironFile
0194 CHARACTER*(MAX_LEN_FNAM) bling_silicaFile
54b6d6aa91 Matt*0195 CHARACTER*(MAX_LEN_FNAM) bling_psmFile
0196 CHARACTER*(MAX_LEN_FNAM) bling_plgFile
0197 CHARACTER*(MAX_LEN_FNAM) bling_pdiazFile
be72e7ae9e Jean*0198 CHARACTER*(MAX_LEN_FNAM) bling_Pc_2dFile
47de7c1e0e Matt*0199 CHARACTER*(MAX_LEN_FNAM) bling_Pc_2d_diazFile
0200 CHARACTER*(MAX_LEN_FNAM) bling_alpha_photo2dFile
0201 CHARACTER*(MAX_LEN_FNAM) bling_k_Fe2dFile
0202 CHARACTER*(MAX_LEN_FNAM) bling_k_Fe_diaz2dFile
af8f3b9ffc Matt*0203 CHARACTER*(MAX_LEN_FNAM) bling_gamma_POM2dFile
0204 CHARACTER*(MAX_LEN_FNAM) bling_wsink0_2dFile
0205 CHARACTER*(MAX_LEN_FNAM) bling_phi_DOM2dFile
0206 CHARACTER*(MAX_LEN_FNAM) bling_phi_lg2dFile
0207 CHARACTER*(MAX_LEN_FNAM) bling_phi_sm2dFile
e0f9a7ba0b Matt*0208 #ifdef ALLOW_EXF
079948e6a6 Matt*0209 CHARACTER*(MAX_LEN_FNAM) apco2file
0210 #endif
c0d1c06c15 Matt*0211 _RL bling_forcingPeriod
0212 _RL bling_forcingCycle
0213 _RL bling_pCO2
e0f9a7ba0b Matt*0214 _RL river_conc_po4
0215 _RL river_dom_to_nut
c0d1c06c15 Matt*0216
079948e6a6 Matt*0217
0218
0219
e0f9a7ba0b Matt*0220 #ifdef ALLOW_EXF
079948e6a6 Matt*0221 _RL exf_inscal_apco2
0222 _RL exf_outscal_apco2
0223 COMMON /BLG_PARAM_SCAL/
0224 & exf_inscal_apco2,
0225 & exf_outscal_apco2
0226 #endif
0227
0228
be72e7ae9e Jean*0229
079948e6a6 Matt*0230
e0f9a7ba0b Matt*0231 #ifdef ALLOW_EXF
079948e6a6 Matt*0232 #ifdef USE_EXF_INTERPOLATION
0233 _RL apco2_lon0, apco2_lon_inc
0234 _RL apco2_lat0, apco2_lat_inc(MAX_LAT_INC)
0235 INTEGER apco2_nlon, apco2_nlat, apco2_interpMethod
0236
0237 COMMON /BLG_EXF_INTERPOLATION/
0238 & apco2_lon0, apco2_lon_inc,
0239 & apco2_lat0, apco2_lat_inc,
0240 & apco2_nlon, apco2_nlat,apco2_interpMethod
0241 #endif
0242 #endif
c0d1c06c15 Matt*0243
0244
0245
0246
e37161e05a Jean*0247
0248
0249
82e538d851 aver*0250
0251
c0d1c06c15 Matt*0252
0253 COMMON /BIOTIC_NEEDS/
0254 & InputFe,
be72e7ae9e Jean*0255 & omegaC,
0256 & omegaAr,
c0d1c06c15 Matt*0257 & irr_mem,
4ac06494d5 Matt*0258 & phyto_lg,
0259 & phyto_sm,
c0d1c06c15 Matt*0260 & chl,
82e538d851 aver*0261 & chl_sat,
3b75ee2a42 Matt*0262 & poc,
e0f9a7ba0b Matt*0263 & Pc_0_2d,
0264 & k_Fe_2d,
0265 & wsink0_2d,
0266 & gamma_POM_2d,
0267 & phi_DOM_2d,
0268 & phi_sm_2d,
0269 & phi_lg_2d,
0270 #ifndef USE_BLING_V1
0271 & phyto_diaz,
0272 & Pc_0_diaz_2d,
0273 & k_Fe_diaz_2d,
0274 & alpha_photo_2d,
c0d1c06c15 Matt*0275 & Pc_0_diaz,
e0f9a7ba0b Matt*0276 & alpha_photo,
0277 & gamma_DON,
0278 & k_Fe_diaz,
0279 & k_NO3,
0280 & k_PtoN,
0281 & k_FetoN,
0282 & PtoN_min,
0283 & PtoN_max,
0284 & FetoN_min,
0285 & FetoN_max,
0286 & kappa_eppley_diaz,
0287 & phi_dvm,
0288 & sigma_dvm,
0289 #ifdef USE_SIBLING
0290 & k_Si,
0291 & gamma_Si_0,
0292 & kappa_remin_Si,
0293 & wsink_Si,
00fa2d4ddd mmaz*0294 & SitoN_uptake_min,
0295 & SitoN_uptake_max,
0296 & SitoN_uptake_scale,
0297 & SitoN_uptake_exp,
e0f9a7ba0b Matt*0298 & q_SitoN_diss,
0299 #endif
0300 #else
0301 & alpha_max,
0302 & alpha_min,
0303 & gamma_biomass,
0304 & k_FetoP,
0305 & FetoP_max,
0306 & Fe_lim_min,
0307 #endif
c0d1c06c15 Matt*0308 & CtoN,
e0f9a7ba0b Matt*0309 & CtoP,
0310 & NtoP,
c0d1c06c15 Matt*0311 & HtoC,
e0f9a7ba0b Matt*0312 & NO3toN,
c0d1c06c15 Matt*0313 & O2toN,
e0f9a7ba0b Matt*0314 & O2toP,
c0d1c06c15 Matt*0315 & CatoN,
e0f9a7ba0b Matt*0316 & CatoP,
c0d1c06c15 Matt*0317 & masstoN,
e0f9a7ba0b Matt*0318 & pivotal,
0319 & Pc_0,
0320 & lambda_0,
0321 & resp_frac,
0322 & chl_min,
c0d1c06c15 Matt*0323 & theta_Fe_max_hi,
0324 & theta_Fe_max_lo,
0325 & gamma_irr_mem,
0326 & gamma_DOP,
0327 & gamma_POM,
0328 & k_O2,
e0f9a7ba0b Matt*0329 & k_Fe,
c0d1c06c15 Matt*0330 & k_PO4,
0331 & kFe_eq_lig_max,
0332 & kFe_eq_lig_min,
0333 & kFe_eq_lig_Femin,
be72e7ae9e Jean*0334 & kFe_eq_lig_irr,
c0d1c06c15 Matt*0335 & kFe_org,
0336 & kFe_inorg,
0337 & FetoC_sed,
0338 & remin_min,
0339 & oxic_min,
0340 & ligand,
0341 & kappa_eppley,
0342 & kappa_remin,
0343 & ca_remin_depth,
0344 & phi_DOM,
0345 & phi_sm,
0346 & phi_lg,
0347 & wsink0,
e0f9a7ba0b Matt*0348 & wsink0z,
c0d1c06c15 Matt*0349 & wsinkacc,
0350 & parfrac,
0351 & alpfe,
0352 & k0,
82e538d851 aver*0353 & MLmix_max,
0354 & chlsat_locTimWindow
c0d1c06c15 Matt*0355
0356 _RL InputFe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0357 _RL omegaC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
0358 _RL omegaAr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
0359 _RL irr_mem(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
4ac06494d5 Matt*0360 _RL phyto_lg(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
0361 _RL phyto_sm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
c0d1c06c15 Matt*0362 _RL chl(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
82e538d851 aver*0363 _RL chl_sat(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
3b75ee2a42 Matt*0364 _RL poc(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
e0f9a7ba0b Matt*0365 _RL Pc_0_2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0366 _RL k_Fe_2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
af8f3b9ffc Matt*0367 _RL wsink0_2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
e0f9a7ba0b Matt*0368 _RL gamma_POM_2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0369 _RL phi_DOM_2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0370 _RL phi_sm_2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0371 _RL phi_lg_2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0372 #ifndef USE_BLING_V1
0373 _RL phyto_diaz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
0374 _RL Pc_0_diaz_2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0375 _RL k_Fe_diaz_2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
0376 _RL alpha_photo_2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
c0d1c06c15 Matt*0377 _RL Pc_0_diaz
e0f9a7ba0b Matt*0378 _RL alpha_photo
0379 _RL gamma_DON
0380 _RL k_Fe_diaz
0381 _RL k_NO3
0382 _RL k_PtoN
0383 _RL k_FetoN
0384 _RL PtoN_min
0385 _RL PtoN_max
0386 _RL FetoN_min
0387 _RL FetoN_max
0388 _RL kappa_eppley_diaz
0389 _RL phi_dvm
0390 _RL sigma_dvm
0391 #ifdef USE_SIBLING
0392 _RL k_Si
0393 _RL gamma_Si_0
0394 _RL kappa_remin_Si
0395 _RL wsink_Si
00fa2d4ddd mmaz*0396 _RL SitoN_uptake_min
0397 _RL SitoN_uptake_max
0398 _RL SitoN_uptake_scale
0399 _RL SitoN_uptake_exp
e0f9a7ba0b Matt*0400 _RL q_SitoN_diss
0401 #endif
0402 #else
0403 _RL alpha_max
0404 _RL alpha_min
0405 _RL gamma_biomass
0406 _RL k_FetoP
0407 _RL FetoP_max
0408 _RL Fe_lim_min
0409 #endif
c0d1c06c15 Matt*0410 _RL CtoN
e0f9a7ba0b Matt*0411 _RL CtoP
0412 _RL NtoP
c0d1c06c15 Matt*0413 _RL HtoC
e0f9a7ba0b Matt*0414 _RL NO3toN
c0d1c06c15 Matt*0415 _RL O2toN
e0f9a7ba0b Matt*0416 _RL O2toP
be72e7ae9e Jean*0417 _RL CatoN
e0f9a7ba0b Matt*0418 _RL CatoP
be72e7ae9e Jean*0419 _RL masstoN
e0f9a7ba0b Matt*0420 _RL pivotal
0421 _RL Pc_0
0422 _RL lambda_0
0423 _RL resp_frac
0424 _RL chl_min
c0d1c06c15 Matt*0425 _RL theta_Fe_max_hi
0426 _RL theta_Fe_max_lo
0427 _RL gamma_irr_mem
0428 _RL gamma_DOP
0429 _RL gamma_POM
0430 _RL k_O2
e0f9a7ba0b Matt*0431 _RL k_Fe
c0d1c06c15 Matt*0432 _RL k_PO4
0433 _RL kFe_eq_lig_max
0434 _RL kFe_eq_lig_min
0435 _RL kFe_eq_lig_Femin
be72e7ae9e Jean*0436 _RL kFe_eq_lig_irr
c0d1c06c15 Matt*0437 _RL kFe_org
0438 _RL kFe_inorg
0439 _RL FetoC_sed
0440 _RL remin_min
0441 _RL oxic_min
0442 _RL ligand
0443 _RL kappa_eppley
0444 _RL kappa_remin
0445 _RL ca_remin_depth
0446 _RL phi_DOM
0447 _RL phi_sm
0448 _RL phi_lg
0449 _RL wsink0
e0f9a7ba0b Matt*0450 _RL wsink0z
c0d1c06c15 Matt*0451 _RL wsinkacc
0452 _RL parfrac
0453 _RL alpfe
0454 _RL k0
e0f9a7ba0b Matt*0455 _RL MLmix_max
82e538d851 aver*0456 _RL chlsat_locTimWindow(2)
c0d1c06c15 Matt*0457
0458
0459
0460