146 |
c OR use a constant coeff |
c OR use a constant coeff |
147 |
c Kwexch(i,j) = 5e-5 |
c Kwexch(i,j) = 5e-5 |
148 |
c ice influence |
c ice influence |
149 |
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) |
150 |
|
|
151 |
|
|
152 |
C Calculate flux in terms of DIC units using K0, solubility |
C Calculate flux in terms of DIC units using K0, solubility |
186 |
C update tendency |
C update tendency |
187 |
DO j=1-OLy,sNy+OLy |
DO j=1-OLy,sNy+OLy |
188 |
DO i=1-OLx,sNx+OLx |
DO i=1-OLx,sNx+OLx |
189 |
GDC(i,j)= maskC(i,j,kLev,bi,bj)*recip_drF(kLev)*( |
GDC(i,j)= maskC(i,j,kLev,bi,bj)*recip_drF(kLev)* |
190 |
|
& recip_hFacC(i,j,kLev,bi,bj)*( |
191 |
& FluxCO2(i,j,bi,bj) + VirtualFlux(i,j) |
& FluxCO2(i,j,bi,bj) + VirtualFlux(i,j) |
192 |
& ) |
& ) |
193 |
ENDDO |
ENDDO |