/[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.29 by dfer, Fri Oct 7 21:36:39 2011 UTC
# Line 46  c GDC                   :: tendency due Line 46  c GDC                   :: tendency due
46    
47  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
48         INTEGER i,j, kLev         INTEGER i,j, kLev
49           _RL co3dummy
50  C Number of iterations for pCO2 solvers...  C Number of iterations for pCO2 solvers...
51  C Solubility relation coefficients  C Solubility relation coefficients
52        _RL SchmidtNoDIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL SchmidtNoDIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 167  C$TAF LOOP = parallel Line 168  C$TAF LOOP = parallel
168       I        ak0(i,j,bi,bj), fugf(i,j,bi,bj),       I        ak0(i,j,bi,bj), fugf(i,j,bi,bj),
169       I        ff(i,j,bi,bj),       I        ff(i,j,bi,bj),
170       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),
171       U        pH(i,j,bi,bj),pCO2(i,j,bi,bj),       U        pH(i,j,bi,bj),pCO2(i,j,bi,bj),co3dummy,
172       I        i,j,kLev,bi,bj,myIter,myThid )       I        i,j,kLev,bi,bj,myIter,myThid )
173            ELSE            ELSE
174              pCO2(i,j,bi,bj)=0. _d 0              pCO2(i,j,bi,bj)=0. _d 0
# Line 186  C calculate SCHMIDT NO. for CO2 Line 187  C calculate SCHMIDT NO. for CO2
187       &          + 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)
188       &          + 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)
189       &                *theta(i,j,kLev,bi,bj)       &                *theta(i,j,kLev,bi,bj)
190  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)
191               SchmidtNoDIC(i,j)=max(1.0 _d -2, SchmidtNoDIC(i,j))               SchmidtNoDIC(i,j)=max(1.0 _d -2, SchmidtNoDIC(i,j))
192    
193  C Determine surface flux (FDIC)  C Determine surface flux (FDIC)
# Line 209  C Converting pCO2 to [CO2] using ff, as Line 210  C Converting pCO2 to [CO2] using ff, as
210       &         ff(i,j,bi,bj)*pCO2(i,j,bi,bj)       &         ff(i,j,bi,bj)*pCO2(i,j,bi,bj)
211       &         )       &         )
212  #else  #else
213  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
214  C of error in application of water vapor correction  C of error in application of water vapor correction
215  c Flux = kw*rho*(ff*pCO2atm-k0*FugFac*pCO2ocean)  c Flux = kw*rho*(ff*pCO2atm-k0*FugFac*pCO2ocean)
216                 FluxCO2(i,j,bi,bj) =                 FluxCO2(i,j,bi,bj) =
Legend:
Removed from v.1.27  
changed lines
  Added in v.1.29
  ViewVC Help
Powered by ViewVC 1.1.22