/[MITgcm]/MITgcm_contrib/bling/pkg/bling_carbon_chem.F
ViewVC logotype

Diff of /MITgcm_contrib/bling/pkg/bling_carbon_chem.F

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

revision 1.1 by mmazloff, Fri May 23 17:33:42 2014 UTC revision 1.2 by mmazloff, Sun Feb 28 21:49:24 2016 UTC
# Line 786  C LOCAL VARIABLES Line 786  C LOCAL VARIABLES
786          _RL  invtk          _RL  invtk
787          _RL  is          _RL  is
788          _RL  is2          _RL  is2
789    c add Bennington
790            _RL  P1atm
791            _RL  Rgas
792            _RL  RT
793            _RL  delta
794            _RL  B1
795            _RL  B
796          INTEGER i          INTEGER i
797          INTEGER j          INTEGER j
798          INTEGER k          INTEGER k
# Line 853  C terms used more than once Line 860  C terms used more than once
860             s15=s**1.5             s15=s**1.5
861             scl=s/1.80655             scl=s/1.80655
862    
863    C-----------------------------------------------------------------------
864    C Ariane adding this - copied from non-pressure-dep case
865    C added by Val Bennington Nov 2010
866    C Fugacity Factor needed for non-ideality in ocean
867    C ff used for atmospheric correction for water vapor and pressure
868    C Weiss (1974) Marine Chemistry
869               P1atm = 1.01325 _d 0 ! bars
870               Rgas = 83.1451 _d 0 ! bar*cm3/(mol*K)
871               RT = Rgas*tk
872               delta = (57.7 _d 0 - 0.118 _d 0*tk)
873               B1 = -1636.75 _d 0 + 12.0408 _d 0*tk - 0.0327957 _d 0*tk*tk
874               B = B1 + 3.16528 _d 0*tk*tk*tk*(0.00001 _d 0)
875               fugf(i,j,bi,bj) = exp( (B+2. _d 0*delta) * P1atm / RT)
876    
877  C------------------------------------------------------------------------  C------------------------------------------------------------------------
878  C f = k0(1-pH2O)*correction term for non-ideality  C f = k0(1-pH2O)*correction term for non-ideality
879  C Weiss & Price (1980, Mar. Chem., 8, 347-359; Eq 13 with table 6 values)  C Weiss & Price (1980, Mar. Chem., 8, 347-359; Eq 13 with table 6 values)
# Line 1021  c alternative pressure dependence from I Line 1042  c alternative pressure dependence from I
1042    
1043  C------------------------------------------------------------------------  C------------------------------------------------------------------------
1044           else           else
1045                fugf(i,j,bi,bj)=0.d0
1046              ff(i,j,bi,bj)=0.d0              ff(i,j,bi,bj)=0.d0
1047              ak0(i,j,bi,bj)= 0.d0              ak0(i,j,bi,bj)= 0.d0
1048              ak1(i,j,bi,bj)= 0.d0              ak1(i,j,bi,bj)= 0.d0
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22