/[MITgcm]/MITgcm_contrib/dcarroll/highres_darwin/code/DARWIN_FLUX.h
ViewVC logotype

Annotation of /MITgcm_contrib/dcarroll/highres_darwin/code/DARWIN_FLUX.h

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


Revision 1.1 - (hide annotations) (download)
Sun Sep 22 21:23:45 2019 UTC (5 years, 10 months ago) by dcarroll
Branch: MAIN
CVS Tags: HEAD
File MIME type: text/plain
Initial check in of high resolution Darwin simulation code

1 dcarroll 1.1 #ifdef ALLOW_CARBON
2     c air-sea flux info
3     C permil : is conversion factor for mol/m3 to mol/kg
4     C assumes uniform (surface) density
5     C Pa2Atm : for conversion of atmospheric pressure
6     C when coming from atmospheric model
7     COMMON /CARBON_NEEDS/
8     & AtmospCO2, AtmosP, pH, pCO2, FluxCO2,
9     & wind, Kwexch_Pre
10     _RL AtmospCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
11     _RL AtmosP(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
12     _RL pH(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
13     _RL pCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
14     _RL FluxCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
15     _RL wind(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
16     _RL Kwexch_Pre(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
17    
18     C permil : is conversion factor for mol/m3 to mol/kg
19     C assumes uniform (surface) density
20     C Pa2Atm : for conversion of atmospheric pressure
21     C when coming from atmospheric model
22     COMMON /GLOBAL_SURF_MEAN/
23     & gsm_alk,gsm_s,gsm_t,gsm_dic,
24     & gsm_c14,permil,Pa2Atm
25     _RL gsm_alk
26     _RL gsm_s
27     _RL gsm_t
28     _RL gsm_DIC
29     _RL gsm_C14
30     _RL permil
31     _RL Pa2Atm
32    
33     COMMON /CARBON_CHEM/
34     & ak0,ak1,ak2,akw,akb,aks,akf,
35     & ak1p,ak2p,ak3p,aksi, fugf,
36     & ff,ft,st,bt, Ksp_TP_Calc,CO3,
37     & omegaC,disscPIC
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     C Fugacity Factor added following Val Bennington
50     _RL fugf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
51     _RL ff(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
52     _RL ft(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
53     _RL st(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
54     _RL bt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
55     _RL Ksp_TP_Calc(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
56     _RL CO3(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
57     _RL omegaC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
58     _RL disscPIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nR,nSx,nSy)
59    
60     COMMON /CO2_FLUX_BUDGET/
61     & fluxCO2_1,Kwexch1,budgetTemp1,budgetSalt1,budgetAlk1,
62     & budgetDic1,budgetPV1,budgetApCO21,budgetBio1,
63     & budgetPhos1,budgetSi1,
64     & dCO2Flux,
65     & dFluxCO2Temp,dFluxCO2Salt,dFluxCO2Alk,dFluxCO2Dic,
66     & dFluxCO2Pv,dFluxCO2ApCO2,dFluxCO2Sf,
67     & budgetTStep1,
68     & budgetPert,
69     & budgetConsumpDIC,budgetConsumpDIC_PIC,
70     & budgetDOCRemin,budgetPReminC
71    
72     _RL fluxCO2_1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
73     _RL Kwexch1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
74     _RL budgetTemp1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
75     _RL budgetSalt1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
76     _RL budgetAlk1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
77     _RL budgetDic1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
78     _RL budgetPV1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
79     _RL budgetApCO21(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
80     _RL budgetBio1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
81     _RL budgetPhos1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
82     _RL budgetSi1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
83    
84     _RL dCO2Flux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
85    
86     _RL dFluxCO2Temp(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
87     _RL dFluxCO2Salt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
88     _RL dFluxCO2Alk(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
89     _RL dFluxCO2Dic(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
90     _RL dFluxCO2Pv(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
91     _RL dFluxCO2ApCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
92     _RL dFluxCO2Sf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
93    
94     _RL budgetTStep1
95     _RL budgetPert
96    
97     _RL budgetConsumpDIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
98     _RL budgetConsumpDIC_PIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,
99     & Nr,nSx,nSy)
100    
101     _RL budgetDOCRemin(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
102     _RL budgetPReminC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
103    
104     COMMON /OXYGEN_CHEM/
105     & oA0,oA1,oA2,oA3,oA4,oA5,
106     & oB0,oB1,oB2,oB3,
107     & oC0
108     _RL oA0,oA1,oA2,oA3,oA4,oA5
109     _RL oB0,oB1,oB2,oB3
110     _RL oC0
111    
112     C schmidt number coefficients
113     COMMON /DIC_SCHMIDT_NO/
114     & sca1, sca2, sca3, sca4,
115     & sox1, sox2, sox3, sox4
116     _RL sca1
117     _RL sca2
118     _RL sca3
119     _RL sca4
120     _RL sox1
121     _RL sox2
122     _RL sox3
123     _RL sox4
124    
125     COMMON /DIC_DIAGS_NEEDS/
126     & SURave, SURCave, SUROave, pCO2ave, pHave,
127     & dic_timeave
128     _RL SURave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
129     _RL SURCave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
130     _RL SUROave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
131     _RL pCO2ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
132     _RL pHave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
133     _RL DIC_timeave(nSx,nSy)
134    
135     #endif

  ViewVC Help
Powered by ViewVC 1.1.22