/[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.2 - (hide annotations) (download)
Wed Apr 20 19:19:27 2011 UTC (14 years, 3 months ago) by stephd
Branch: MAIN
CVS Tags: ctrb_darwin2_ckpt64k_20130723, ctrb_darwin2_ckpt63l_20120405, ctrb_darwin2_ckpt64h_20130528, ctrb_darwin2_ckpt64m_20130820, ctrb_darwin2_ckpt64f_20130405, ctrb_darwin2_ckpt63f_20111201, ctrb_darwin2_ckpt64a_20121116, ctrb_darwin2_ckpt64n_20130826, ctrb_darwin2_ckpt62y_20110526, ctrb_darwin2_ckpt64i_20130622, ctrb_darwin2_ckpt62x_20110513, ctrb_darwin2_ckpt62w_20110426, ctrb_darwin2_ckpt63o_20120629, ctrb_darwin2_ckpt64e_20130305, ctrb_darwin2_ckpt63c_20111011, ctrb_darwin2_ckpt63i_20120124, ctrb_darwin2_ckpt63m_20120506, ctrb_darwin2_ckpt63s_20120908, ctrb_darwin2_ckpt63e_20111107, ctrb_darwin2_ckpt63b_20110830, ctrb_darwin2_ckpt63j_20120217, ctrb_darwin2_ckpt63r_20120817, ctrb_darwin2_ckpt64g_20130503, ctrb_darwin2_ckpt64l_20130806, ctrb_darwin2_ckpt63g_20111220, ctrb_darwin2_ckpt64c_20130120, ctrb_darwin2_ckpt63a_20110804, ctrb_darwin2_ckpt64j_20130704, ctrb_darwin2_ckpt63h_20111230, ctrb_darwin2_ckpt63p_20120707, ctrb_darwin2_ckpt63d_20111107, ctrb_darwin2_ckpt63q_20120731, ctrb_darwin2_ckpt63_20110728, ctrb_darwin2_ckpt64b_20121224, ctrb_darwin2_ckpt64d_20130219, ctrb_darwin2_ckpt64_20121012, ctrb_darwin2_ckpt63n_20120604, ctrb_darwin2_ckpt63k_20120317, ctrb_darwin2_ckpt62z_20110622
Changes since 1.1: +3 -1 lines
File MIME type: text/plain
o add fix to two bugs in handling water vapour pressure in carbonate chemisty
  and air-sea fluxes. Found by Val Bennington and Galen McKinley

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     _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 stephd 1.2 & ak1p,ak2p,ak3p,aksi, fugf,
36 jahn 1.1 & ff,ft,st,bt, Ksp_TP_Calc
37     _RL ak0(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
38     _RL ak1(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
39     _RL ak2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
40     _RL akw(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
41     _RL akb(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
42     _RL aks(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
43     _RL akf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
44     _RL ak1p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
45     _RL ak2p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
46     _RL ak3p(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
47     _RL aksi(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
48 stephd 1.2 C Fugacity Factor added following Val Bennington
49     _RL fugf(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
50 jahn 1.1 _RL ff(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    
56     COMMON /OXYGEN_CHEM/
57     & oA0,oA1,oA2,oA3,oA4,oA5,
58     & oB0,oB1,oB2,oB3,
59     & oC0
60     _RL oA0,oA1,oA2,oA3,oA4,oA5
61     _RL oB0,oB1,oB2,oB3
62     _RL oC0
63    
64     C schmidt number coefficients
65     COMMON /DIC_SCHMIDT_NO/
66     & sca1, sca2, sca3, sca4,
67     & sox1, sox2, sox3, sox4
68     _RL sca1
69     _RL sca2
70     _RL sca3
71     _RL sca4
72     _RL sox1
73     _RL sox2
74     _RL sox3
75     _RL sox4
76    
77     COMMON /DIC_DIAGS_NEEDS/
78     & SURave, SURCave, SUROave, pCO2ave, pHave,
79     & dic_timeave
80     _RL SURave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
81     _RL SURCave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
82     _RL SUROave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
83     _RL pCO2ave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
84     _RL pHave(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
85     _RL DIC_timeave(nSx,nSy,nR)
86    
87     #endif

  ViewVC Help
Powered by ViewVC 1.1.22