--- MITgcm/pkg/dic/dic_surfforcing.F 2004/08/06 19:50:52 1.8 +++ 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.8 2004/08/06 19:50:52 stephd 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" @@ -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,7 +186,8 @@ C update tendency DO j=1-OLy,sNy+OLy DO i=1-OLx,sNx+OLx - GDC(i,j)= maskC(i,j,kLev,bi,bj)*recip_drF(kLev)*( + 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