--- MITgcm/pkg/dic/dic_surfforcing.F	2011/05/05 22:23:27	1.27
+++ MITgcm/pkg/dic/dic_surfforcing.F	2011/10/07 21:36:39	1.29
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/dic_surfforcing.F,v 1.27 2011/05/05 22:23:27 stephd Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/dic_surfforcing.F,v 1.29 2011/10/07 21:36:39 dfer Exp $
 C $Name:  $
 
 #include "DIC_OPTIONS.h"
@@ -46,6 +46,7 @@
 
 C !LOCAL VARIABLES: ====================================================
        INTEGER i,j, kLev
+       _RL co3dummy
 C Number of iterations for pCO2 solvers...
 C Solubility relation coefficients
       _RL SchmidtNoDIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
@@ -167,7 +168,7 @@
      I        ak0(i,j,bi,bj), fugf(i,j,bi,bj),
      I        ff(i,j,bi,bj),
      I        bt(i,j,bi,bj),st(i,j,bi,bj),ft(i,j,bi,bj),
-     U        pH(i,j,bi,bj),pCO2(i,j,bi,bj),
+     U        pH(i,j,bi,bj),pCO2(i,j,bi,bj),co3dummy,
      I        i,j,kLev,bi,bj,myIter,myThid )
           ELSE
             pCO2(i,j,bi,bj)=0. _d 0
@@ -186,7 +187,7 @@
      &          + sca3 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj)
      &          + sca4 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj)
      &                *theta(i,j,kLev,bi,bj)
-c make sure Schmidt number isn't negative (will happen if temp>39C)
+c make sure Schmidt number is not negative (will happen if temp>39C)
              SchmidtNoDIC(i,j)=max(1.0 _d -2, SchmidtNoDIC(i,j))
 
 C Determine surface flux (FDIC)
@@ -209,7 +210,7 @@
      &         ff(i,j,bi,bj)*pCO2(i,j,bi,bj)
      &         )
 #else
-C Corrected by Val Bennington Nov 2010 per G.A. McKinley's finding
+C Corrected by Val Bennington Nov 2010 per G.A. McKinley s finding
 C of error in application of water vapor correction
 c Flux = kw*rho*(ff*pCO2atm-k0*FugFac*pCO2ocean)
                FluxCO2(i,j,bi,bj) =