bling/pkg/BLING_VARS.h

C $Header: $
C $Name: $


C ==========================================================
C   Carbon chemistry variables
C ==========================================================

       COMMON /CARBON_NEEDS/
     &                      AtmospCO2, AtmosP, pH, pCO2, FluxCO2,
     &                      wind, FIce, Silica
      _RL  AtmospCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  AtmosP(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  pH(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL  pCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  FluxCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  wind(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  FIce(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  Silica(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)


C ==========================================================
C   Carbon and oxygen chemistry parameters
C ==========================================================

       COMMON /CARBON_CHEM/
     &                     ak0,ak1,ak2,akw,akb,aks,akf,
     &                     ak1p,ak2p,ak3p,aksi, fugf, 
     &                     ff,ft,st,bt, 
     &                     Ksp_TP_Calc,Ksp_TP_Arag
      _RL  ak0(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  ak1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  ak2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  akw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  akb(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  aks(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  akf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  ak1p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  ak2p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  ak3p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  aksi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  ff(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  fugf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  ft(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  st(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  bt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  Ksp_TP_Calc(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL  Ksp_TP_Arag(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)


       COMMON /OXYGEN_CHEM/
     &                     oA0,oA1,oA2,oA3,oA4,oA5,
     &                     oB0,oB1,oB2,oB3,
     &                     oC0
      _RL oA0,oA1,oA2,oA3,oA4,oA5
      _RL oB0,oB1,oB2,oB3
      _RL oC0


       COMMON /GLOBAL_SURF_MEAN/
     &                          permil,Pa2Atm

C      permil : is conversion factor for mol/m3 to mol/kg
C               assumes uniform (surface) density
C      Pa2Atm : for conversion of atmospheric pressure
C               when coming from atmospheric model

      _RL  permil
      _RL  Pa2Atm


       COMMON /SCHMIDT_NO/
     &                    sca1, sca2, sca3, sca4,
     &                    sox1, sox2, sox3, sox4

C      Schmidt number coefficients

      _RL  sca1,sca2,sca3,sca4
      _RL  sox1,sox2,sox3,sox4


C ==========================================================
C   Bling inputs (specified in data.bling)
C ==========================================================

       COMMON /BLING_INPUTS/
     &        bling_windFile, bling_atmospFile, bling_iceFile,
     &        bling_ironFile, bling_silicaFile, 
     &        bling_forcingPeriod, bling_forcingCycle,
     &        bling_pCO2, 
     &        river_conc_trac

C      bling_windFile      :: file name of wind speeds
C      bling_atmospFile    :: file name of atmospheric pressure
C      bling_iceFile       :: file name of sea ice fraction
C      bling_ironFile      :: file name of aeolian iron flux
C      bling_silicaFile    :: file name of surface silica
C      bling_forcingPeriod :: period of forcing for biogeochemistry (seconds)
C      bling_forcingCycle  :: periodic forcing parameter for biogeochemistry 
C      bling_pCO2          :: Atmospheric pCO2 to be read in data.bling
C      river_conc_trac     :: River concentration, bgc tracers

      CHARACTER*(MAX_LEN_FNAM) bling_windFile
      CHARACTER*(MAX_LEN_FNAM) bling_atmospFile
      CHARACTER*(MAX_LEN_FNAM) bling_iceFile
      CHARACTER*(MAX_LEN_FNAM) bling_ironFile
      CHARACTER*(MAX_LEN_FNAM) bling_silicaFile
      _RL     bling_forcingPeriod
      _RL     bling_forcingCycle
      _RL     bling_pCO2
c      _RL     river_conc_trac(PTRACERS_num)
      _RL     river_conc_trac(8)


C ==========================================================
C   Ecosystem variables and parameters
C ==========================================================

      COMMON /BIOTIC_NEEDS/
     &                     InputFe,
     &                     omegaC, 
     &                     omegaAr, 
     &                     irr_mem,
     &                     P_sm,
     &                     P_lg,
     &                     P_diaz,
     &                     chl,
     &                     pivotal,
     &                     Pc_0,
     &                     Pc_0_diaz,
     &                     lambda_0,
     &                     chl_min,
     &                     CtoN,
     &                     NO3toN,
     &                     HtoC,
     &                     O2toN,
     &                     CatoN,
     &                     masstoN,
     &                     alpha_photo,
     &                     theta_Fe_max_hi,
     &                     theta_Fe_max_lo,
     &                     gamma_irr_mem,
     &                     gamma_DON,
     &                     gamma_DOP,
     &                     gamma_POM,
     &                     k_Fe,
     &                     k_Fe_diaz,
     &                     k_O2,
     &                     k_NO3,
     &                     k_PO4,
     &                     k_PtoN,
     &                     k_FetoN,
     &                     kFe_eq_lig_max,
     &                     kFe_eq_lig_min,
     &                     kFe_eq_lig_Femin,
     &                     kFe_eq_lig_irr,    
     &                     kFe_org,
     &                     kFe_inorg,
     &                     PtoN_min,
     &                     PtoN_max,
     &                     FetoN_min,
     &                     FetoN_max,
     &                     FetoC_sed,
     &                     remin_min,
     &                     oxic_min,
     &                     ligand,
     &                     kappa_eppley,
     &                     kappa_eppley_diaz,
     &                     kappa_remin,
     &                     ca_remin_depth,
     &                     phi_DOM,
     &                     phi_sm,
     &                     phi_lg,
     &                     phi_dvm,
     &                     sigma_dvm,
     &                     wsink0z,
     &                     wsink0,
     &                     wsinkacc,
     &                     parfrac,
     &                     alpfe,
     &                     k0,
     &                     epsln,
     &                     QSW_underice

      _RL InputFe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
      _RL omegaC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL omegaAr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL irr_mem(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL P_sm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL P_lg(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL P_diaz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL chl(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
      _RL pivotal
      _RL Pc_0
      _RL Pc_0_diaz
      _RL lambda_0
      _RL chl_min
      _RL CtoN
      _RL NO3toN
      _RL HtoC
      _RL O2toN
      _RL CatoN 
      _RL masstoN 
      _RL alpha_photo
      _RL theta_Fe_max_hi
      _RL theta_Fe_max_lo
      _RL gamma_irr_mem
      _RL gamma_DON
      _RL gamma_DOP
      _RL gamma_POM
      _RL k_Fe
      _RL k_Fe_diaz
      _RL k_O2
      _RL k_NO3
      _RL k_PO4
      _RL k_PtoN
      _RL k_FetoN
      _RL kFe_eq_lig_max
      _RL kFe_eq_lig_min
      _RL kFe_eq_lig_Femin
      _RL kFe_eq_lig_irr    
      _RL kFe_org
      _RL kFe_inorg
      _RL PtoN_min
      _RL PtoN_max
      _RL FetoN_min
      _RL FetoN_max
      _RL FetoC_sed
      _RL remin_min
      _RL oxic_min
      _RL ligand
      _RL kappa_eppley
      _RL kappa_eppley_diaz
      _RL kappa_remin
      _RL ca_remin_depth
      _RL phi_DOM
      _RL phi_sm
      _RL phi_lg
      _RL phi_dvm
      _RL sigma_dvm
      _RL wsink0z
      _RL wsink0
      _RL wsinkacc
      _RL parfrac
      _RL alpfe
      _RL k0
      _RL epsln
      LOGICAL QSW_underice


CEH3 ;;; Local Variables: ***
CEH3 ;;; mode:fortran ***
CEH3 ;;; End: ***
1	mmazloff	1.1	C $Header: $
2			C $Name: $
3
4
5			C ==========================================================
6			C Carbon chemistry variables
7			C ==========================================================
8
9			COMMON /CARBON_NEEDS/
10			& AtmospCO2, AtmosP, pH, pCO2, FluxCO2,
11			& wind, FIce, Silica
12			_RL AtmospCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
13			_RL AtmosP(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
14			_RL pH(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
15			_RL pCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
16			_RL FluxCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
17			_RL wind(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
18			_RL FIce(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
19			_RL Silica(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
20
21
22			C ==========================================================
23			C Carbon and oxygen chemistry parameters
24			C ==========================================================
25
26			COMMON /CARBON_CHEM/
27			& ak0,ak1,ak2,akw,akb,aks,akf,
28			& ak1p,ak2p,ak3p,aksi, fugf,
29			& ff,ft,st,bt,
30			& Ksp_TP_Calc,Ksp_TP_Arag
31			_RL ak0(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
32			_RL ak1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
33			_RL ak2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
34			_RL akw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
35			_RL akb(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
36			_RL aks(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
37			_RL akf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
38			_RL ak1p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
39			_RL ak2p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
40			_RL ak3p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
41			_RL aksi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
42			_RL ff(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
43			_RL fugf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
44			_RL ft(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
45			_RL st(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
46			_RL bt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
47			_RL Ksp_TP_Calc(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
48			_RL Ksp_TP_Arag(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
49
50
51			COMMON /OXYGEN_CHEM/
52			& oA0,oA1,oA2,oA3,oA4,oA5,
53			& oB0,oB1,oB2,oB3,
54			& oC0
55			_RL oA0,oA1,oA2,oA3,oA4,oA5
56			_RL oB0,oB1,oB2,oB3
57			_RL oC0
58
59
60			COMMON /GLOBAL_SURF_MEAN/
61			& permil,Pa2Atm
62
63			C permil : is conversion factor for mol/m3 to mol/kg
64			C assumes uniform (surface) density
65			C Pa2Atm : for conversion of atmospheric pressure
66			C when coming from atmospheric model
67
68			_RL permil
69			_RL Pa2Atm
70
71
72			COMMON /SCHMIDT_NO/
73			& sca1, sca2, sca3, sca4,
74			& sox1, sox2, sox3, sox4
75
76			C Schmidt number coefficients
77
78			_RL sca1,sca2,sca3,sca4
79			_RL sox1,sox2,sox3,sox4
80
81
82			C ==========================================================
83			C Bling inputs (specified in data.bling)
84			C ==========================================================
85
86			COMMON /BLING_INPUTS/
87			& bling_windFile, bling_atmospFile, bling_iceFile,
88			& bling_ironFile, bling_silicaFile,
89			& bling_forcingPeriod, bling_forcingCycle,
90			& bling_pCO2,
91	mmazloff	1.2	& river_conc_trac
92	mmazloff	1.1
93			C bling_windFile :: file name of wind speeds
94			C bling_atmospFile :: file name of atmospheric pressure
95			C bling_iceFile :: file name of sea ice fraction
96			C bling_ironFile :: file name of aeolian iron flux
97			C bling_silicaFile :: file name of surface silica
98			C bling_forcingPeriod :: period of forcing for biogeochemistry (seconds)
99			C bling_forcingCycle :: periodic forcing parameter for biogeochemistry
100			C bling_pCO2 :: Atmospheric pCO2 to be read in data.bling
101	mmazloff	1.2	C river_conc_trac :: River concentration, bgc tracers
102	mmazloff	1.1
103			CHARACTER*(MAX_LEN_FNAM) bling_windFile
104			CHARACTER*(MAX_LEN_FNAM) bling_atmospFile
105			CHARACTER*(MAX_LEN_FNAM) bling_iceFile
106			CHARACTER*(MAX_LEN_FNAM) bling_ironFile
107			CHARACTER*(MAX_LEN_FNAM) bling_silicaFile
108			_RL bling_forcingPeriod
109			_RL bling_forcingCycle
110			_RL bling_pCO2
111	mmazloff	1.2	c _RL river_conc_trac(PTRACERS_num)
112			_RL river_conc_trac(8)
113	mmazloff	1.1
114
115			C ==========================================================
116			C Ecosystem variables and parameters
117			C ==========================================================
118
119			COMMON /BIOTIC_NEEDS/
120			& InputFe,
121			& omegaC,
122			& omegaAr,
123			& irr_mem,
124			& P_sm,
125			& P_lg,
126	mmazloff	1.2	& P_diaz,
127			& chl,
128			& pivotal,
129	mmazloff	1.1	& Pc_0,
130	mmazloff	1.2	& Pc_0_diaz,
131	mmazloff	1.1	& lambda_0,
132			& chl_min,
133	mmazloff	1.2	& CtoN,
134			& NO3toN,
135			& HtoC,
136			& O2toN,
137			& CatoN,
138			& masstoN,
139			& alpha_photo,
140	mmazloff	1.1	& theta_Fe_max_hi,
141			& theta_Fe_max_lo,
142			& gamma_irr_mem,
143	mmazloff	1.2	& gamma_DON,
144			& gamma_DOP,
145	mmazloff	1.1	& gamma_POM,
146			& k_Fe,
147	mmazloff	1.2	& k_Fe_diaz,
148	mmazloff	1.1	& k_O2,
149	mmazloff	1.2	& k_NO3,
150			& k_PO4,
151			& k_PtoN,
152			& k_FetoN,
153			& kFe_eq_lig_max,
154			& kFe_eq_lig_min,
155			& kFe_eq_lig_Femin,
156			& kFe_eq_lig_irr,
157	mmazloff	1.1	& kFe_org,
158			& kFe_inorg,
159	mmazloff	1.2	& PtoN_min,
160			& PtoN_max,
161			& FetoN_min,
162			& FetoN_max,
163			& FetoC_sed,
164	mmazloff	1.1	& remin_min,
165	mmazloff	1.2	& oxic_min,
166	mmazloff	1.1	& ligand,
167			& kappa_eppley,
168	mmazloff	1.2	& kappa_eppley_diaz,
169	mmazloff	1.1	& kappa_remin,
170			& ca_remin_depth,
171			& phi_DOM,
172			& phi_sm,
173			& phi_lg,
174	mmazloff	1.2	& phi_dvm,
175			& sigma_dvm,
176	mmazloff	1.1	& wsink0z,
177			& wsink0,
178			& wsinkacc,
179			& parfrac,
180			& alpfe,
181			& k0,
182			& epsln,
183			& QSW_underice
184
185			_RL InputFe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
186			_RL omegaC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
187			_RL omegaAr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
188			_RL irr_mem(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
189			_RL P_sm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
190			_RL P_lg(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
191	mmazloff	1.2	_RL P_diaz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
192			_RL chl(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
193			_RL pivotal
194	mmazloff	1.1	_RL Pc_0
195	mmazloff	1.2	_RL Pc_0_diaz
196	mmazloff	1.1	_RL lambda_0
197			_RL chl_min
198	mmazloff	1.2	_RL CtoN
199			_RL NO3toN
200			_RL HtoC
201			_RL O2toN
202			_RL CatoN
203			_RL masstoN
204			_RL alpha_photo
205	mmazloff	1.1	_RL theta_Fe_max_hi
206			_RL theta_Fe_max_lo
207			_RL gamma_irr_mem
208	mmazloff	1.2	_RL gamma_DON
209			_RL gamma_DOP
210	mmazloff	1.1	_RL gamma_POM
211			_RL k_Fe
212	mmazloff	1.2	_RL k_Fe_diaz
213	mmazloff	1.1	_RL k_O2
214	mmazloff	1.2	_RL k_NO3
215			_RL k_PO4
216			_RL k_PtoN
217			_RL k_FetoN
218			_RL kFe_eq_lig_max
219			_RL kFe_eq_lig_min
220			_RL kFe_eq_lig_Femin
221			_RL kFe_eq_lig_irr
222	mmazloff	1.1	_RL kFe_org
223			_RL kFe_inorg
224	mmazloff	1.2	_RL PtoN_min
225			_RL PtoN_max
226			_RL FetoN_min
227			_RL FetoN_max
228			_RL FetoC_sed
229	mmazloff	1.1	_RL remin_min
230	mmazloff	1.2	_RL oxic_min
231	mmazloff	1.1	_RL ligand
232			_RL kappa_eppley
233	mmazloff	1.2	_RL kappa_eppley_diaz
234	mmazloff	1.1	_RL kappa_remin
235			_RL ca_remin_depth
236			_RL phi_DOM
237			_RL phi_sm
238			_RL phi_lg
239	mmazloff	1.2	_RL phi_dvm
240			_RL sigma_dvm
241	mmazloff	1.1	_RL wsink0z
242			_RL wsink0
243			_RL wsinkacc
244			_RL parfrac
245			_RL alpfe
246			_RL k0
247			_RL epsln
248			LOGICAL QSW_underice
249
250
251			CEH3 ;;; Local Variables: ***
252			CEH3 ;;; mode:fortran ***
253			CEH3 ;;; End: ***