162 |
& + sca3 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj) |
& + sca3 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj) |
163 |
& + sca4 * 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) |
164 |
& *theta(i,j,kLev,bi,bj) |
& *theta(i,j,kLev,bi,bj) |
165 |
|
c make sure Schmidt number isn't negative (will happen if temp>39C) |
166 |
|
SchmidtNoDIC(i,j)=max(1.0 _d -2, SchmidtNoDIC(i,j)) |
167 |
|
|
168 |
C Determine surface flux (FDIC) |
C Determine surface flux (FDIC) |
169 |
C first correct pCO2at for surface atmos pressure |
C first correct pCO2at for surface atmos pressure |