Back to home page

MITgcm

 
 

    


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 C ==========================================================
                0002 C   Carbon chemistry variables
                0003 C ==========================================================
                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 C ==========================================================
                0025 C   Carbon and oxygen chemistry parameters
                0026 C ==========================================================
                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 C If using Solvesaphe routines (Munhoven, 2013) then in adittion,
                0054 C  cat  = total calcium concentration, Ref.: Culkin (1965)
                0055 C   akn = the dissociation constant of ammonium [H][NH3]/[NH4]
                0056 C           Ref.: Yao and Millero (1995)
                0057 C   akhs = the dissociation constant of hydrogen sulfide [H][HS]/[H2S]
                0058 C           Ref.: Millero et al. (1988)
                0059 C          (cited by Millero (1995) and Yao and Millero (1995))
                0060 C  aphscale = pH scale conversion factor ; convert from the total to the free
                0061 C          scale for solvesaphe calculations;  Ref.: Munhoven, 2013
                0062 C   Ksp_TP_Arag = solubility product for aragonite, Ref.: Mucci (1983)
                0063 C----
                0064 C  selectBTconst :: estimates borate concentration from salinity:
                0065 C     =1 :: use default formulation of Uppström (1974)(same as S/R CARBON_COEFFS)
                0066 C     =2 :: use new formulation from Lee et al (2010)
                0067 C
                0068 C  selectFTconst :: estimates fluoride concentration from salinity:
                0069 C     =1 :: use default formulation of Riley (1965) (same as S/R CARBON_COEFFS)
                0070 C     =2 :: use new formulation from Culkin (1965)
                0071 C
                0072 C  selectHFconst :: sets the first dissociation constant for hydrogen fluoride:
                0073 C     =1 :: use default  Dickson and Riley (1979) (same as S/R CARBON_COEFFS)
                0074 C     =2 :: use new formulation of Perez and Fraga (1987)
                0075 C
                0076 C  selectK1K2const :: sets the 1rst & 2nd dissociation constants of carbonic acid:
                0077 C     =1 :: use default formulation of Millero (1995) with data
                0078 C            from Mehrbach et al. (1973) (same as S/R CARBON_COEFFS)
                0079 C     =2 :: use formulation of Roy et al. (1993)
                0080 C     =3 :: use "combination" formulation of Millero (1995)
                0081 C     =4 :: use formulation of Luecker et al. (2000)
                0082 C     =5 :: use formulation of Millero (2010, Mar. Fresh Wat. Res.)
                0083 C     =6 :: use formulation of Waters, Millero, Woosley (2014, Mar. Chem.)
                0084 C  selectPHsolver :: sets the pH solver to use:
                0085 C     =1 :: use the GENERAL solver ;  =2 :: use SEC solver ;
                0086 C     =3 :: use FAST solver routine.
                0087 
                0088        COMMON /CARBONCHEM_SOLVESAPHE_ARIANE/
                0089 cav     &                     cat, akn, akhs, aphscale, Ksp_TP_Arag,
                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 cav      _RL  Ksp_TP_Arag(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                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 C      permil : is conversion factor for mol/m3 to mol/kg
                0122 C               assumes uniform (surface) density
                0123 C      Pa2Atm : for conversion of atmospheric pressure
                0124 C               when coming from atmospheric model
                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 C      Schmidt number coefficients
                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 C ==========================================================
                0152 C   Bling inputs (specified in data.bling)
                0153 C ==========================================================
                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 C      bling_windFile      :: file name of wind speeds
                0177 C      bling_atmospFile    :: file name of atmospheric pressure
                0178 C      bling_iceFile       :: file name of sea ice fraction
                0179 C      bling_ironFile      :: file name of aeolian iron flux
                0180 C      bling_silicaFile    :: file name of surface silica
54b6d6aa91 Matt*0181 C      bling_psmFile       :: file name of init small phyto biomass
                0182 C      bling_plgFile       :: file name of init lg phyto biomass
                0183 C      bling_pdiazFile     :: file name of init diaz biomass
c0d1c06c15 Matt*0184 C      bling_forcingPeriod :: period of forcing for biogeochemistry (seconds)
be72e7ae9e Jean*0185 C      bling_forcingCycle  :: periodic forcing parameter for biogeochemistry
c0d1c06c15 Matt*0186 C      bling_pCO2          :: Atmospheric pCO2 to be read in data.bling
                0187 C      river_conc_trac     :: River concentration, bgc tracers
079948e6a6 Matt*0188 C      apco2               :: Atmospheric pCO2 to be read in with exf pkg
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 C ==========================================================
                0218 C   EXF input/output scaling factors for unit conversion if needed
                0219 C ==========================================================
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 C ==========================================================
be72e7ae9e Jean*0229 C   EXF interpolation needs
079948e6a6 Matt*0230 C ==========================================================
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 C ==========================================================
                0245 C   Ecosystem variables and parameters
                0246 C ==========================================================
e37161e05a Jean*0247 C     irr_mem       :: Phyto irradiance memory
                0248 C          this is a temporally smoothed field carried between timesteps,
                0249 C          to represent photoadaptation.
82e538d851 aver*0250 C   chlsat_locTimWindow(1:2) :: local-time window (in h) for
                0251 C          satellite-equivalent chlorophyll diagnostic (and cost)
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 CEH3 ;;; Local Variables: ***
                0459 CEH3 ;;; mode:fortran ***
                0460 CEH3 ;;; End: ***