/[MITgcm]/MITgcm_contrib/ecco_darwin/v4_llc270/code_darwin/DARWIN_FLUX.h
ViewVC logotype

Annotation of /MITgcm_contrib/ecco_darwin/v4_llc270/code_darwin/DARWIN_FLUX.h

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


Revision 1.1 - (hide annotations) (download)
Wed Jun 5 02:36:32 2019 UTC (7 years, 3 months ago) by dimitri
Branch: MAIN
File MIME type: text/plain
Adding Naviaux (2019, Marine Chemistry) dissolution rate.

1 dimitri 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 CO3(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
15     _RL FluxCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
16     _RL wind(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
17     _RL Kwexch_Pre(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
18    
19     C permil : is conversion factor for mol/m3 to mol/kg
20     C assumes uniform (surface) density
21     C Pa2Atm : for conversion of atmospheric pressure
22     C when coming from atmospheric model
23     COMMON /GLOBAL_SURF_MEAN/
24     & gsm_alk,gsm_s,gsm_t,gsm_dic,
25     & gsm_c14,permil,Pa2Atm
26     _RL gsm_alk
27     _RL gsm_s
28     _RL gsm_t
29     _RL gsm_DIC
30     _RL gsm_C14
31     _RL permil
32     _RL Pa2Atm
33    
34     COMMON /CARBON_CHEM/
35     & ak0,ak1,ak2,akw,akb,aks,akf,
36     & ak1p,ak2p,ak3p,aksi, fugf,
37     & ff,ft,st,bt, Ksp_TP_Calc
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    
57     COMMON /OXYGEN_CHEM/
58     & oA0,oA1,oA2,oA3,oA4,oA5,
59     & oB0,oB1,oB2,oB3,
60     & oC0
61     _RL oA0,oA1,oA2,oA3,oA4,oA5
62     _RL oB0,oB1,oB2,oB3
63     _RL oC0
64    
65     C schmidt number coefficients
66     COMMON /DIC_SCHMIDT_NO/
67     & sca1, sca2, sca3, sca4,
68     & sox1, sox2, sox3, sox4
69     _RL sca1
70     _RL sca2
71     _RL sca3
72     _RL sca4
73     _RL sox1
74     _RL sox2
75     _RL sox3
76     _RL sox4
77    
78     COMMON /DIC_DIAGS_NEEDS/
79     & SURave, SURCave, SUROave, pCO2ave, pHave,
80     & dic_timeave
81     _RL SURave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
82     _RL SURCave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
83     _RL SUROave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
84     _RL pCO2ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
85     _RL pHave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
86     _RL DIC_timeave(nSx,nSy)
87    
88     #endif

  ViewVC Help
Powered by ViewVC 1.1.22