/[MITgcm]/MITgcm_contrib/darwin2/pkg/darwin/DARWIN_FLUX.h
ViewVC logotype

Annotation of /MITgcm_contrib/darwin2/pkg/darwin/DARWIN_FLUX.h

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


Revision 1.3 - (hide annotations) (download)
Wed Oct 9 17:14:37 2013 UTC (11 years, 9 months ago) by stephd
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt64r_20131210, ctrb_darwin2_ckpt64o_20131024, ctrb_darwin2_ckpt64p_20131118, ctrb_darwin2_ckpt64q_20131118, ctrb_darwin2_ckpt64p_20131024
Changes since 1.2: +10 -0 lines
File MIME type: text/plain
o changes so that pH and pCO2 can be calculated for full water column
  by defining pH_3D in DARWIN_OPTIONS.h. Includes 3D diags and pickup

1 jahn 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 stephd 1.3 #ifdef pH_3D
13     _RL pH(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
14     _RL pCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
15     #else
16 jahn 1.1 _RL pH(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
17     _RL pCO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
18 stephd 1.3 #endif
19 jahn 1.1 _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 Kwexch_Pre(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
22    
23     C permil : is conversion factor for mol/m3 to mol/kg
24     C assumes uniform (surface) density
25     C Pa2Atm : for conversion of atmospheric pressure
26     C when coming from atmospheric model
27     COMMON /GLOBAL_SURF_MEAN/
28     & gsm_alk,gsm_s,gsm_t,gsm_dic,
29     & gsm_c14,permil,Pa2Atm
30     _RL gsm_alk
31     _RL gsm_s
32     _RL gsm_t
33     _RL gsm_DIC
34     _RL gsm_C14
35     _RL permil
36     _RL Pa2Atm
37    
38     COMMON /CARBON_CHEM/
39     & ak0,ak1,ak2,akw,akb,aks,akf,
40 stephd 1.2 & ak1p,ak2p,ak3p,aksi, fugf,
41 jahn 1.1 & ff,ft,st,bt, Ksp_TP_Calc
42     _RL ak0(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
43     _RL ak1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
44     _RL ak2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
45     _RL akw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
46     _RL akb(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
47     _RL aks(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
48     _RL akf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
49     _RL ak1p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
50     _RL ak2p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
51     _RL ak3p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
52     _RL aksi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
53 stephd 1.2 C Fugacity Factor added following Val Bennington
54     _RL fugf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
55 jahn 1.1 _RL ff(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
56     _RL ft(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
57     _RL st(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
58     _RL bt(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
59     _RL Ksp_TP_Calc(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
60    
61     COMMON /OXYGEN_CHEM/
62     & oA0,oA1,oA2,oA3,oA4,oA5,
63     & oB0,oB1,oB2,oB3,
64     & oC0
65     _RL oA0,oA1,oA2,oA3,oA4,oA5
66     _RL oB0,oB1,oB2,oB3
67     _RL oC0
68    
69     C schmidt number coefficients
70     COMMON /DIC_SCHMIDT_NO/
71     & sca1, sca2, sca3, sca4,
72     & sox1, sox2, sox3, sox4
73     _RL sca1
74     _RL sca2
75     _RL sca3
76     _RL sca4
77     _RL sox1
78     _RL sox2
79     _RL sox3
80     _RL sox4
81    
82     COMMON /DIC_DIAGS_NEEDS/
83     & SURave, SURCave, SUROave, pCO2ave, pHave,
84     & dic_timeave
85     _RL SURave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
86     _RL SURCave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
87     _RL SUROave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
88 stephd 1.3 #ifdef pH_3D
89     _RL pCO2ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
90     _RL pHave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
91     #else
92 jahn 1.1 _RL pCO2ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
93     _RL pHave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
94 stephd 1.3 #endif
95 jahn 1.1 _RL DIC_timeave(nSx,nSy,nR)
96    
97     #endif

  ViewVC Help
Powered by ViewVC 1.1.22