/[MITgcm]/MITgcm/pkg/dic/dic_surfforcing.F
ViewVC logotype

Diff of /MITgcm/pkg/dic/dic_surfforcing.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.27 by stephd, Thu May 5 22:23:27 2011 UTC revision 1.28 by jmc, Fri Jun 24 01:35:32 2011 UTC
# Line 186  C calculate SCHMIDT NO. for CO2 Line 186  C calculate SCHMIDT NO. for CO2
186       &          + 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)
187       &          + 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)
188       &                *theta(i,j,kLev,bi,bj)       &                *theta(i,j,kLev,bi,bj)
189  c make sure Schmidt number isn't negative (will happen if temp>39C)  c make sure Schmidt number is not negative (will happen if temp>39C)
190               SchmidtNoDIC(i,j)=max(1.0 _d -2, SchmidtNoDIC(i,j))               SchmidtNoDIC(i,j)=max(1.0 _d -2, SchmidtNoDIC(i,j))
191    
192  C Determine surface flux (FDIC)  C Determine surface flux (FDIC)
# Line 209  C Converting pCO2 to [CO2] using ff, as Line 209  C Converting pCO2 to [CO2] using ff, as
209       &         ff(i,j,bi,bj)*pCO2(i,j,bi,bj)       &         ff(i,j,bi,bj)*pCO2(i,j,bi,bj)
210       &         )       &         )
211  #else  #else
212  C Corrected by Val Bennington Nov 2010 per G.A. McKinley's finding  C Corrected by Val Bennington Nov 2010 per G.A. McKinley s finding
213  C of error in application of water vapor correction  C of error in application of water vapor correction
214  c Flux = kw*rho*(ff*pCO2atm-k0*FugFac*pCO2ocean)  c Flux = kw*rho*(ff*pCO2atm-k0*FugFac*pCO2ocean)
215                 FluxCO2(i,j,bi,bj) =                 FluxCO2(i,j,bi,bj) =
Legend:
Removed from v.1.27  
changed lines
  Added in v.1.28
  ViewVC Help
Powered by ViewVC 1.1.22