-
Notifications
You must be signed in to change notification settings - Fork 1
/
DARWIN_FLUX.h
87 lines (81 loc) · 3.19 KB
/
DARWIN_FLUX.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
#ifdef ALLOW_CARBON
c air-sea flux info
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
COMMON /CARBON_NEEDS/
& AtmospCO2, AtmosP, pH, pCO2, FluxCO2,
& wind, Kwexch_Pre
_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,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 Kwexch_Pre(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
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
COMMON /GLOBAL_SURF_MEAN/
& gsm_alk,gsm_s,gsm_t,gsm_dic,
& gsm_c14,permil,Pa2Atm
_RL gsm_alk
_RL gsm_s
_RL gsm_t
_RL gsm_DIC
_RL gsm_C14
_RL permil
_RL Pa2Atm
COMMON /CARBON_CHEM/
& ak0,ak1,ak2,akw,akb,aks,akf,
& ak1p,ak2p,ak3p,aksi, fugf,
& ff,ft,st,bt, Ksp_TP_Calc
_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)
C Fugacity Factor added following Val Bennington
_RL fugf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
_RL ff(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)
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
C schmidt number coefficients
COMMON /DIC_SCHMIDT_NO/
& sca1, sca2, sca3, sca4,
& sox1, sox2, sox3, sox4
_RL sca1
_RL sca2
_RL sca3
_RL sca4
_RL sox1
_RL sox2
_RL sox3
_RL sox4
COMMON /DIC_DIAGS_NEEDS/
& SURave, SURCave, SUROave, pCO2ave, pHave,
& dic_timeave
_RL SURave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
_RL SURCave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
_RL SUROave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
_RL pCO2ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
_RL pHave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
_RL DIC_timeave(nSx,nSy,nR)
#endif