pkg/darwin/DARWIN_FIELDS.h

8fbfd1f382 Oliv*0001 #ifdef ALLOW_DARWIN
                0002
                0003 CBOP
                0004 C     !ROUTINE: DARWIN_FIELDS.h
                0005 C     !INTERFACE:
                0006 C #include DARWIN_FIELDS.h
                0007
                0008 C     !DESCRIPTION:
                0009 C Contains fields for darwin package
                0010 C
                0011 C Requires: SIZE.h
                0012 C
                0013
                0014 #ifndef DARWIN_ALLOW_CHLQUOTA
                0015       COMMON /DARWIN_CHL_STORE/ chlPrev
                0016 #ifdef ALLOW_RADTRANS
                0017       _RL chlPrev(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy,nPhoto)
                0018 #else
                0019       _RL chlPrev(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0020 #endif
                0021 #endif
                0022
                0023
                0024 #ifdef DARWIN_ALLOW_CONS
                0025       COMMON /DARWIN_CHECK_CONS/
                0026      &      DAR_cons_C_unit,
                0027      &      DAR_cons_P_unit,
                0028      &      DAR_cons_N_unit,
                0029      &      DAR_cons_Fe_unit,
                0030      &      DAR_cons_Si_unit,
                0031      &      DAR_cons_A_unit,
                0032      &      DAR_cons_O_unit
                0033       INTEGER DAR_cons_C_unit
                0034       INTEGER DAR_cons_P_unit
                0035       INTEGER DAR_cons_N_unit
                0036       INTEGER DAR_cons_Fe_unit
                0037       INTEGER DAR_cons_Si_unit
                0038       INTEGER DAR_cons_A_unit
                0039       INTEGER DAR_cons_O_unit
                0040
a59ebcf0a3 Oliv*0041       COMMON /DARWIN_CONS_3D/
                0042      &      DARWIN_Nfix,
                0043      &      DARWIN_Ndenit,
                0044      &      DARWIN_partScav,
5e411acc9e Oliv*0045      &      DARWIN_minFeLoss,
                0046      &      DARWIN_O2prod,
                0047      &      DARWIN_O2cons,
                0048      &      DARWIN_AlkSrc
8fbfd1f382 Oliv*0049       _RL DARWIN_Nfix(sNx,sNy,Nr,nSx,nSy)
                0050       _RL DARWIN_Ndenit(sNx,sNy,Nr,nSx,nSy)
a59ebcf0a3 Oliv*0051       _RL DARWIN_partScav(sNx,sNy,Nr,nSx,nSy)
                0052       _RL DARWIN_minFeLoss(sNx,sNy,Nr,nSx,nSy)
5e411acc9e Oliv*0053       _RL DARWIN_O2prod(sNx,sNy,Nr,nSx,nSy)
                0054       _RL DARWIN_O2cons(sNx,sNy,Nr,nSx,nSy)
                0055       _RL DARWIN_AlkSrc(sNx,sNy,Nr,nSx,nSy)
a59ebcf0a3 Oliv*0056
5e411acc9e Oliv*0057       COMMON /DARWIN_CONS_2D/
                0058      &      ironSedFlx,
16faaabac2 Oliv*0059 # ifdef DARWIN_ALLOW_HYDROTHERMAL_VENTS
2c93eb88ef Oliv*0060      &      ironVentFlx,
16faaabac2 Oliv*0061 # endif
5e411acc9e Oliv*0062      &      carbSfcFlx,
                0063      &      carbVirFlx,
                0064      &      oxySfcFlx,
                0065      &      alkVirFlx
                0066       _RL ironSedFlx(sNx,sNy,nSx,nSy)
16faaabac2 Oliv*0067 # ifdef DARWIN_ALLOW_HYDROTHERMAL_VENTS
2c93eb88ef Oliv*0068       _RL ironVentFlx(sNx,sNy,nSx,nSy)
16faaabac2 Oliv*0069 # endif
5e411acc9e Oliv*0070       _RL carbSfcFlx(sNx,sNy,nSx,nSy)
                0071       _RL carbVirFlx(sNx,sNy,nSx,nSy)
                0072       _RL oxySfcFlx(sNx,sNy,nSx,nSy)
                0073       _RL alkVirFlx(sNx,sNy,nSx,nSy)
248275f1c4 Oliv*0074
                0075 # ifdef DARWIN_BOTTOM_SINK
                0076       COMMON /DARWIN_CONS_BOTTOM/
                0077      &      botSnkC,
                0078      &      botSnkN,
                0079      &      botSnkP,
                0080      &      botSnkFe,
                0081      &      botSnkSi
                0082       _RL botSnkC(sNx,sNy,nSx,nSy)
                0083       _RL botSnkN(sNx,sNy,nSx,nSy)
                0084       _RL botSnkP(sNx,sNy,nSx,nSy)
                0085       _RL botSnkFe(sNx,sNy,nSx,nSy)
                0086       _RL botSnkSi(sNx,sNy,nSx,nSy)
                0087 # endif
407760706f Oliv*0088
2a97638e3d Oliv*0089 #ifdef DARWIN_ALLOW_CARBON
8a10163480 Oliv*0090 # ifdef DARWIN_ALLOW_RADIv1
2a97638e3d Oliv*0091 C     Sediment fluxes
                0092       COMMON /DARWIN_CONS_RADI/
                0093      &      radiFLuxC, radiFLuxA, radiFLuxO
                0094       _RL radiFluxC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0095       _RL radiFluxA(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0096       _RL radiFluxO(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
8a10163480 Oliv*0097 # endif
                0098 # ifdef DARWIN_ALLOW_RADIv2
                0099 C     Sediment fluxes
                0100       COMMON /DARWIN_CONS_RADI/
                0101      &      radiFLuxC, radiFLuxA, radiFLuxO, radiFluxP, radiFluxN
                0102       _RL radiFluxC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0103       _RL radiFluxA(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0104       _RL radiFluxO(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0105       _RL radiFluxP(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0106       _RL radiFluxN(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0107 # endif
2a97638e3d Oliv*0108 #endif
                0109
16faaabac2 Oliv*0110       COMMON /DARWIN_CONS_IMPFS/
                0111      &      totFS,
                0112      &      totFSC,
                0113      &      totFSN,
                0114      &      totFSP,
                0115      &      totFSFe,
                0116      &      totFSSi,
                0117      &      totFSA,
                0118      &      totFSO
                0119       _RL totFS
                0120       _RL totFSC
                0121       _RL totFSN
                0122       _RL totFSP
                0123       _RL totFSFe
                0124       _RL totFSSi
                0125       _RL totFSA
                0126       _RL totFSO
                0127
407760706f Oliv*0128       COMMON /DARWIN_CONS_EPR/
16faaabac2 Oliv*0129      &      totPER,
407760706f Oliv*0130      &      totEPRC,
                0131      &      totEPRN,
                0132      &      totEPRP,
                0133      &      totEPRFe,
                0134      &      totEPRSi,
                0135      &      totEPRA,
                0136      &      totEPRO
16faaabac2 Oliv*0137       _RL totPER
407760706f Oliv*0138       _RL totEPRC
                0139       _RL totEPRN
                0140       _RL totEPRP
                0141       _RL totEPRFe
                0142       _RL totEPRSi
                0143       _RL totEPRA
                0144       _RL totEPRO
895c6145db Oliv*0145 #endif /* DARWIN_ALLOW_CONS */
8fbfd1f382 Oliv*0146
                0147 C Carbon Variables
                0148
895c6145db Oliv*0149 #ifdef DARWIN_ALLOW_CARBON
1a75f0d7fd Oliv*0150 C-     COMMON /CARBON_NEEDS/ main carbon chem arrays
                0151 C  omegaC           :: Local saturation state with respect to calcite
8fbfd1f382 Oliv*0152        COMMON /CARBON_NEEDS/
1a75f0d7fd Oliv*0153      &   pH, omegaC, Atmosp
8fbfd1f382 Oliv*0154       _RL  pH(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
1a75f0d7fd Oliv*0155       _RL  omegaC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
8fbfd1f382 Oliv*0156       _RL  AtmosP(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0157
                0158        COMMON /CARBON_CHEM/
                0159      &                     ak0,ak1,ak2,akw,akb,aks,akf,
                0160      &                     ak1p,ak2p,ak3p,aksi, fugf,
                0161      &                     ff,ft,st,bt, Ksp_TP_Calc
                0162       _RL  ak0(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0163       _RL  ak1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0164       _RL  ak2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0165       _RL  akw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0166       _RL  akb(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0167       _RL  aks(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0168       _RL  akf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0169       _RL  ak1p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0170       _RL  ak2p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0171       _RL  ak3p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0172       _RL  aksi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0173       _RL  ff(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0174 C Fugacity Factor added by Val Bennington Nov. 2010
                0175       _RL  fugf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0176       _RL  ft(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0177       _RL  st(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0178       _RL  bt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0179       _RL  Ksp_TP_Calc(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0180
f0a72e2151 Oliv*0181 #ifdef DARWIN_SOLVESAPHE
                0182 C If using Solvesaphe routines (Munhoven, 2013) then in adittion,
                0183 C  cat  = total calcium concentration, Ref.: Culkin (1965)
                0184 C   akn = the dissociation constant of ammonium [H][NH3]/[NH4]
                0185 C           Ref.: Yao and Millero (1995)
                0186 C   akhs = the dissociation constant of hydrogen sulfide [H][HS]/[H2S]
                0187 C           Ref.: Millero et al. (1988)
                0188 C          (cited by Millero (1995) and Yao and Millero (1995))
                0189 C  aphscale = pH scale conversion factor ; convert from the total to the free
                0190 C          scale for solvesaphe calculations;  Ref.: Munhoven, 2013
                0191 C   Ksp_TP_Arag = solubility product for aragonite, Ref.: Mucci (1983)
                0192
                0193        COMMON /DARWIN_SOLVESAPHE/
                0194      &                     cat, akn, akhs, aphscale, Ksp_TP_Arag
                0195
                0196       _RL  cat(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0197       _RL  akn(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0198       _RL  akhs(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0199       _RL  aphscale(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0200       _RL  Ksp_TP_Arag(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0201 #endif /* DARWIN_SOLVESAPHE */
895c6145db Oliv*0202 #endif /* DARWIN_ALLOW_CARBON */
f0a72e2151 Oliv*0203
8fbfd1f382 Oliv*0204 CEOP
                0205 #endif /* ALLOW_DARWIN */