139 |
I ak1(i,j,bi,bj),ak2(i,j,bi,bj), |
I ak1(i,j,bi,bj),ak2(i,j,bi,bj), |
140 |
I ak1p(i,j,bi,bj),ak2p(i,j,bi,bj),ak3p(i,j,bi,bj), |
I ak1p(i,j,bi,bj),ak2p(i,j,bi,bj),ak3p(i,j,bi,bj), |
141 |
I aks(i,j,bi,bj),akb(i,j,bi,bj),akw(i,j,bi,bj), |
I aks(i,j,bi,bj),akb(i,j,bi,bj),akw(i,j,bi,bj), |
142 |
I aksi(i,j,bi,bj),akf(i,j,bi,bj),ff(i,j,bi,bj), |
I aksi(i,j,bi,bj),akf(i,j,bi,bj), |
143 |
|
I ak0(i,j,bi,bj), fugf(i,j,bi,bj), |
144 |
|
I ff(i,j,bi,bj), |
145 |
I bt(i,j,bi,bj),st(i,j,bi,bj),ft(i,j,bi,bj), |
I bt(i,j,bi,bj),st(i,j,bi,bj),ft(i,j,bi,bj), |
146 |
U pH(i,j,bi,bj),pCO2(i,j,bi,bj), |
U pH(i,j,bi,bj),pCO2(i,j,bi,bj), |
147 |
I i,j,kLev,bi,bj,myIter,myThid ) |
I i,j,kLev,bi,bj,myIter,myThid ) |
172 |
Kwexch(i,j) = Kwexch_Pre(i,j,bi,bj) |
Kwexch(i,j) = Kwexch_Pre(i,j,bi,bj) |
173 |
& / sqrt(SchmidtNoDIC(i,j)/660.0 _d 0) |
& / sqrt(SchmidtNoDIC(i,j)/660.0 _d 0) |
174 |
|
|
175 |
|
#ifdef WATERVAP_BUG |
176 |
C Calculate flux in terms of DIC units using K0, solubility |
C Calculate flux in terms of DIC units using K0, solubility |
177 |
C Flux = Vp * ([CO2sat] - [CO2]) |
C Flux = Vp * ([CO2sat] - [CO2]) |
178 |
C CO2sat = K0*pCO2atmos*P/P0 |
C CO2sat = K0*pCO2atmos*P/P0 |
182 |
& ak0(i,j,bi,bj)*pCO2sat(i,j) - |
& ak0(i,j,bi,bj)*pCO2sat(i,j) - |
183 |
& ff(i,j,bi,bj)*pCO2(i,j,bi,bj) |
& ff(i,j,bi,bj)*pCO2(i,j,bi,bj) |
184 |
& ) |
& ) |
185 |
|
#else |
186 |
|
C Corrected by Val Bennington Nov 2010 per G.A. McKinley's finding |
187 |
|
C of error in application of water vapor correction |
188 |
|
c Flux = kw*rho*(ff*pCO2atm-k0*FugFac*pCO2ocean) |
189 |
|
FluxCO2(i,j,bi,bj) = |
190 |
|
& Kwexch(i,j)*( |
191 |
|
& ff(i,j,bi,bj)*pCO2sat(i,j) - |
192 |
|
& pCO2(i,j,bi,bj)*fugf(i,j,bi,bj) |
193 |
|
& *ak0(i,j,bi,bj) ) |
194 |
|
& |
195 |
|
#endif |
196 |
ELSE |
ELSE |
197 |
FluxCO2(i,j,bi,bj) = 0. _d 0 |
FluxCO2(i,j,bi,bj) = 0. _d 0 |
198 |
ENDIF |
ENDIF |