/[MITgcm]/MITgcm/pkg/bling/BLING_VARS.h
ViewVC logotype

Annotation of /MITgcm/pkg/bling/BLING_VARS.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.6 - (hide annotations) (download)
Fri Jan 27 23:11:31 2017 UTC (7 years, 4 months ago) by mmazloff
Branch: MAIN
Changes since 1.5: +20 -3 lines
File MIME type: text/plain
Make a few more BGC params 2D.

1 mmazloff 1.6 C $Header: /u/gcmpack/MITgcm/pkg/bling/BLING_VARS.h,v 1.5 2017/01/20 00:31:49 mmazloff Exp $
2 mmazloff 1.1 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 mmazloff 1.3 #ifdef USE_EXFCO2
13     & ,apco2, apco20, apco21
14     #endif
15 mmazloff 1.1 _RL AtmospCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
16     _RL AtmosP(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
17     _RL pH(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
18     _RL pCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
19     _RL FluxCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
20     _RL wind(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
21     _RL FIce(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
22     _RL Silica(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
23 mmazloff 1.3 #ifdef USE_EXFCO2
24     _RL apco2 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
25     _RL apco20 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
26     _RL apco21 (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
27     #endif
28 mmazloff 1.1
29     C ==========================================================
30     C Carbon and oxygen chemistry parameters
31     C ==========================================================
32    
33     COMMON /CARBON_CHEM/
34     & ak0,ak1,ak2,akw,akb,aks,akf,
35     & ak1p,ak2p,ak3p,aksi, fugf,
36     & ff,ft,st,bt,
37     & Ksp_TP_Calc,Ksp_TP_Arag
38     _RL ak0(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
39     _RL ak1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
40     _RL ak2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
41     _RL akw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
42     _RL akb(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
43     _RL aks(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
44     _RL akf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
45     _RL ak1p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
46     _RL ak2p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
47     _RL ak3p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
48     _RL aksi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
49     _RL ff(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
50     _RL fugf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
51     _RL ft(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
52     _RL st(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
53     _RL bt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
54     _RL Ksp_TP_Calc(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
55     _RL Ksp_TP_Arag(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
56    
57    
58     COMMON /OXYGEN_CHEM/
59     & oA0,oA1,oA2,oA3,oA4,oA5,
60     & oB0,oB1,oB2,oB3,
61     & oC0
62     _RL oA0,oA1,oA2,oA3,oA4,oA5
63     _RL oB0,oB1,oB2,oB3
64     _RL oC0
65    
66    
67     COMMON /GLOBAL_SURF_MEAN/
68     & permil,Pa2Atm
69    
70     C permil : is conversion factor for mol/m3 to mol/kg
71     C assumes uniform (surface) density
72     C Pa2Atm : for conversion of atmospheric pressure
73     C when coming from atmospheric model
74    
75     _RL permil
76     _RL Pa2Atm
77    
78    
79     COMMON /SCHMIDT_NO/
80     & sca1, sca2, sca3, sca4,
81     & sox1, sox2, sox3, sox4
82    
83     C Schmidt number coefficients
84    
85     _RL sca1,sca2,sca3,sca4
86     _RL sox1,sox2,sox3,sox4
87    
88 mmazloff 1.3 #ifdef USE_EXFCO2
89     integer apco2startdate1
90     integer apco2startdate2
91     _RL apco2startdate
92     _RL apco2period
93     _RL apco2const
94     _RL apco2_exfremo_intercept
95     _RL apco2_exfremo_slope
96     character*1 apco2mask
97     #endif
98 mmazloff 1.1
99     C ==========================================================
100     C Bling inputs (specified in data.bling)
101     C ==========================================================
102    
103     COMMON /BLING_INPUTS/
104     & bling_windFile, bling_atmospFile, bling_iceFile,
105 mmazloff 1.4 & bling_ironFile, bling_silicaFile,
106     & bling_psmFile, bling_plgFile, bling_pdiazFile,
107 mmazloff 1.1 & bling_forcingPeriod, bling_forcingCycle,
108     & bling_pCO2,
109 mmazloff 1.5 & river_conc_trac,
110     & bling_Pc_2dFile, bling_Pc_2d_diazFile,
111 mmazloff 1.6 & bling_alpha_photo2dFile,bling_phi_DOM2dFile,
112     & bling_k_Fe2dFile, bling_k_Fe_diaz2dFile,
113     & bling_gamma_POM2dFile, bling_wsink0_2dFile,
114     & bling_phi_lg2dFile, bling_phi_sm2dFile
115 mmazloff 1.3 #ifdef USE_EXFCO2
116     & ,apco2startdate1,apco2startdate2,
117     & apco2period, apco2startdate,
118     & apco2const,
119     & apco2_exfremo_intercept,
120     & apco2_exfremo_slope,
121     & apco2file, apco2mask
122     #endif
123 mmazloff 1.1
124     C bling_windFile :: file name of wind speeds
125     C bling_atmospFile :: file name of atmospheric pressure
126     C bling_iceFile :: file name of sea ice fraction
127     C bling_ironFile :: file name of aeolian iron flux
128     C bling_silicaFile :: file name of surface silica
129 mmazloff 1.4 C bling_psmFile :: file name of init small phyto biomass
130     C bling_plgFile :: file name of init lg phyto biomass
131     C bling_pdiazFile :: file name of init diaz biomass
132 mmazloff 1.1 C bling_forcingPeriod :: period of forcing for biogeochemistry (seconds)
133     C bling_forcingCycle :: periodic forcing parameter for biogeochemistry
134     C bling_pCO2 :: Atmospheric pCO2 to be read in data.bling
135     C river_conc_trac :: River concentration, bgc tracers
136 mmazloff 1.3 C apco2 :: Atmospheric pCO2 to be read in with exf pkg
137 mmazloff 1.1
138     CHARACTER*(MAX_LEN_FNAM) bling_windFile
139     CHARACTER*(MAX_LEN_FNAM) bling_atmospFile
140     CHARACTER*(MAX_LEN_FNAM) bling_iceFile
141     CHARACTER*(MAX_LEN_FNAM) bling_ironFile
142     CHARACTER*(MAX_LEN_FNAM) bling_silicaFile
143 mmazloff 1.4 CHARACTER*(MAX_LEN_FNAM) bling_psmFile
144     CHARACTER*(MAX_LEN_FNAM) bling_plgFile
145     CHARACTER*(MAX_LEN_FNAM) bling_pdiazFile
146 mmazloff 1.5 CHARACTER*(MAX_LEN_FNAM) bling_Pc_2dFile
147     CHARACTER*(MAX_LEN_FNAM) bling_Pc_2d_diazFile
148     CHARACTER*(MAX_LEN_FNAM) bling_alpha_photo2dFile
149     CHARACTER*(MAX_LEN_FNAM) bling_k_Fe2dFile
150     CHARACTER*(MAX_LEN_FNAM) bling_k_Fe_diaz2dFile
151 mmazloff 1.6 CHARACTER*(MAX_LEN_FNAM) bling_gamma_POM2dFile
152     CHARACTER*(MAX_LEN_FNAM) bling_wsink0_2dFile
153     CHARACTER*(MAX_LEN_FNAM) bling_phi_DOM2dFile
154     CHARACTER*(MAX_LEN_FNAM) bling_phi_lg2dFile
155     CHARACTER*(MAX_LEN_FNAM) bling_phi_sm2dFile
156 mmazloff 1.3 #ifdef USE_EXFCO2
157     CHARACTER*(MAX_LEN_FNAM) apco2file
158     #endif
159 mmazloff 1.1 _RL bling_forcingPeriod
160     _RL bling_forcingCycle
161     _RL bling_pCO2
162     c _RL river_conc_trac(PTRACERS_num)
163     _RL river_conc_trac(8)
164    
165 mmazloff 1.3 C ==========================================================
166     C EXF input/output scaling factors for unit conversion if needed
167     C ==========================================================
168     #ifdef USE_EXFCO2
169     _RL exf_inscal_apco2
170     _RL exf_outscal_apco2
171     COMMON /BLG_PARAM_SCAL/
172     & exf_inscal_apco2,
173     & exf_outscal_apco2
174     #endif
175    
176     C ==========================================================
177     C EXF interpolation needs
178     C ==========================================================
179     #ifdef USE_EXFCO2
180     #ifdef USE_EXF_INTERPOLATION
181     _RL apco2_lon0, apco2_lon_inc
182     _RL apco2_lat0, apco2_lat_inc(MAX_LAT_INC)
183     INTEGER apco2_nlon, apco2_nlat, apco2_interpMethod
184    
185     COMMON /BLG_EXF_INTERPOLATION/
186     & apco2_lon0, apco2_lon_inc,
187     & apco2_lat0, apco2_lat_inc,
188     & apco2_nlon, apco2_nlat,apco2_interpMethod
189     #endif
190     #endif
191 mmazloff 1.1
192     C ==========================================================
193     C Ecosystem variables and parameters
194     C ==========================================================
195    
196     COMMON /BIOTIC_NEEDS/
197     & InputFe,
198     & omegaC,
199     & omegaAr,
200     & irr_mem,
201 mmazloff 1.2 & phyto_lg,
202     & phyto_sm,
203     & phyto_diaz,
204 mmazloff 1.1 & chl,
205     & pivotal,
206     & Pc_0,
207     & Pc_0_diaz,
208 mmazloff 1.5 & Pc_2d,
209     & Pc_2d_diaz,
210 mmazloff 1.1 & lambda_0,
211     & chl_min,
212     & CtoN,
213     & NO3toN,
214     & HtoC,
215     & O2toN,
216     & CatoN,
217     & masstoN,
218     & alpha_photo,
219 mmazloff 1.5 & alpha_photo2d,
220 mmazloff 1.1 & theta_Fe_max_hi,
221     & theta_Fe_max_lo,
222     & gamma_irr_mem,
223     & gamma_DON,
224     & gamma_DOP,
225     & gamma_POM,
226 mmazloff 1.6 & gamma_POM2d,
227 mmazloff 1.1 & k_Fe,
228     & k_Fe_diaz,
229 mmazloff 1.5 & k_Fe2d,
230     & k_Fe_diaz2d,
231 mmazloff 1.1 & k_O2,
232     & k_NO3,
233     & k_PO4,
234     & k_PtoN,
235     & k_FetoN,
236     & kFe_eq_lig_max,
237     & kFe_eq_lig_min,
238     & kFe_eq_lig_Femin,
239     & kFe_eq_lig_irr,
240     & kFe_org,
241     & kFe_inorg,
242     & PtoN_min,
243     & PtoN_max,
244     & FetoN_min,
245     & FetoN_max,
246     & FetoC_sed,
247     & remin_min,
248     & oxic_min,
249     & ligand,
250     & kappa_eppley,
251     & kappa_eppley_diaz,
252     & kappa_remin,
253     & ca_remin_depth,
254     & phi_DOM,
255     & phi_sm,
256     & phi_lg,
257 mmazloff 1.6 & phi_DOM2d,
258     & phi_sm2d,
259     & phi_lg2d,
260 mmazloff 1.1 & phi_dvm,
261     & sigma_dvm,
262     & wsink0z,
263     & wsink0,
264 mmazloff 1.6 & wsink0_2d,
265 mmazloff 1.1 & wsinkacc,
266     & parfrac,
267     & alpfe,
268     & k0,
269     & epsln,
270     & QSW_underice
271    
272     _RL InputFe(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
273     _RL omegaC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
274     _RL omegaAr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
275     _RL irr_mem(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
276 mmazloff 1.2 _RL phyto_lg(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
277     _RL phyto_sm(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
278     _RL phyto_diaz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
279 mmazloff 1.1 _RL chl(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
280 mmazloff 1.5 _RL Pc_2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
281     _RL Pc_2d_diaz(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
282     _RL alpha_photo2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
283     _RL k_Fe2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
284     _RL k_Fe_diaz2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
285 mmazloff 1.6 _RL wsink0_2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
286     _RL gamma_POM2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
287     _RL phi_DOM2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
288     _RL phi_sm2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
289     _RL phi_lg2d(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
290 mmazloff 1.1 _RL pivotal
291     _RL Pc_0
292     _RL Pc_0_diaz
293     _RL lambda_0
294     _RL chl_min
295     _RL CtoN
296     _RL NO3toN
297     _RL HtoC
298     _RL O2toN
299     _RL CatoN
300     _RL masstoN
301     _RL alpha_photo
302     _RL theta_Fe_max_hi
303     _RL theta_Fe_max_lo
304     _RL gamma_irr_mem
305     _RL gamma_DON
306     _RL gamma_DOP
307     _RL gamma_POM
308     _RL k_Fe
309     _RL k_Fe_diaz
310     _RL k_O2
311     _RL k_NO3
312     _RL k_PO4
313     _RL k_PtoN
314     _RL k_FetoN
315     _RL kFe_eq_lig_max
316     _RL kFe_eq_lig_min
317     _RL kFe_eq_lig_Femin
318     _RL kFe_eq_lig_irr
319     _RL kFe_org
320     _RL kFe_inorg
321     _RL PtoN_min
322     _RL PtoN_max
323     _RL FetoN_min
324     _RL FetoN_max
325     _RL FetoC_sed
326     _RL remin_min
327     _RL oxic_min
328     _RL ligand
329     _RL kappa_eppley
330     _RL kappa_eppley_diaz
331     _RL kappa_remin
332     _RL ca_remin_depth
333     _RL phi_DOM
334     _RL phi_sm
335     _RL phi_lg
336     _RL phi_dvm
337     _RL sigma_dvm
338     _RL wsink0z
339     _RL wsink0
340     _RL wsinkacc
341     _RL parfrac
342     _RL alpfe
343     _RL k0
344     _RL epsln
345     LOGICAL QSW_underice
346    
347    
348     CEH3 ;;; Local Variables: ***
349     CEH3 ;;; mode:fortran ***
350     CEH3 ;;; End: ***

  ViewVC Help
Powered by ViewVC 1.1.22