--- MITgcm/pkg/dic/dic_surfforcing.F	2004/07/18 01:13:50	1.7
+++ MITgcm/pkg/dic/dic_surfforcing.F	2005/08/25 22:32:44	1.9
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/dic_surfforcing.F,v 1.7 2004/07/18 01:13:50 jmc Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/dic_surfforcing.F,v 1.9 2005/08/25 22:32:44 stephd Exp $
 C $Name:  $
 
 #include "DIC_OPTIONS.h"
@@ -42,7 +42,7 @@
       INTEGER iMin,iMax,jMin,jMax, bi, bj
 
 C !OUTPUT PARAMETERS: ===================================================
-c GDC                   :: surface flux due to air-sea exchange
+c GDC                   :: tendency due to air-sea exchange
       _RL  GDC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
 
 #ifdef ALLOW_PTRACERS
@@ -146,7 +146,7 @@
 c OR use a constant  coeff
 c             Kwexch(i,j) = 5e-5
 c ice influence
-cQQ           Kwexch(i,j)  =(1.d0-Fice(i,j,bi,bj))*Kwexch(i,j)
+              Kwexch(i,j)  =(1.d0-Fice(i,j,bi,bj))*Kwexch(i,j)
 
 
 C Calculate flux in terms of DIC units using K0, solubility
@@ -186,9 +186,9 @@
 C update tendency      
          DO j=1-OLy,sNy+OLy
           DO i=1-OLx,sNx+OLx
-           GDC(i,j)= maskC(i,j,kLev,bi,bj)*(
-     &                    FluxCO2(i,j,bi,bj)
-     &                    + VirtualFlux(i,j)
+           GDC(i,j)= maskC(i,j,kLev,bi,bj)*recip_drF(kLev)*
+     &                     recip_hFacC(i,j,kLev,bi,bj)*(
+     &                    FluxCO2(i,j,bi,bj) + VirtualFlux(i,j)
      &                                              )
           ENDDO
          ENDDO