9 |
|
|
10 |
C !INTERFACE: ========================================================== |
C !INTERFACE: ========================================================== |
11 |
SUBROUTINE DIC_SURFFORCING( PTR_CO2 , PTR_ALK, PTR_PO4, GDC, |
SUBROUTINE DIC_SURFFORCING( PTR_CO2 , PTR_ALK, PTR_PO4, GDC, |
12 |
I bi,bj,imin,imax,jmin,jmax, |
I bi,bj,iMin,iMax,jMin,jMax, |
13 |
I myIter,myTime,myThid) |
I myIter,myTime,myThid) |
14 |
|
|
15 |
C !DESCRIPTION: |
C !DESCRIPTION: |
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) |
80 |
c ENDDO |
c ENDDO |
81 |
c#endif |
c#endif |
82 |
|
|
|
|
|
83 |
C ================================================================= |
C ================================================================= |
84 |
C determine inorganic carbon chem coefficients |
C determine inorganic carbon chem coefficients |
85 |
DO j=jmin,jmax |
DO j=jMin,jMax |
86 |
DO i=imin,imax |
DO i=iMin,iMax |
87 |
|
|
88 |
#ifdef DIC_BIOTIC |
#ifdef DIC_BIOTIC |
89 |
cQQQQ check ptracer numbers |
cQQQQ check ptracer numbers |
90 |
#ifdef DIC_BOUNDS |
#ifdef DIC_BOUNDS |
91 |
surfalk(i,j) = max(0.4 _d 0, |
surfalk(i,j) = max(0.4 _d 0, |
92 |
& min(10. _d 0,PTR_ALK(i,j,klev))) |
& min(10. _d 0,PTR_ALK(i,j,klev))) |
93 |
& * maskC(i,j,kLev,bi,bj) |
& * maskC(i,j,kLev,bi,bj) |
94 |
surfphos(i,j) = max(1.0 _d -11, |
surfphos(i,j) = max(1.0 _d -11, |
95 |
& min(1._d -1, PTR_PO4(i,j,klev))) |
& min(1. _d -1,PTR_PO4(i,j,klev))) |
96 |
& * maskC(i,j,kLev,bi,bj) |
& * maskC(i,j,kLev,bi,bj) |
97 |
#else |
#else |
98 |
surfalk(i,j) = PTR_ALK(i,j,klev) |
surfalk(i,j) = PTR_ALK(i,j,klev) |
108 |
C FOR NON-INTERACTIVE Si |
C FOR NON-INTERACTIVE Si |
109 |
surfsi(i,j) = SILICA(i,j,bi,bj) * maskC(i,j,kLev,bi,bj) |
surfsi(i,j) = SILICA(i,j,bi,bj) * maskC(i,j,kLev,bi,bj) |
110 |
#ifdef DIC_BOUNDS |
#ifdef DIC_BOUNDS |
111 |
surftemp(i,j) = max(-4. _d 0, |
surftemp(i,j) = max(-4. _d 0, |
112 |
& min(50. _d 0, theta(i,j,kLev,bi,bj))) |
& min(50. _d 0, theta(i,j,kLev,bi,bj))) |
113 |
surfsalt(i,j) = max(4. _d 0, |
surfsalt(i,j) = max(4. _d 0, |
114 |
& min(50. _d 0, salt(i,j,kLev,bi,bj))) |
& min(50. _d 0, salt(i,j,kLev,bi,bj))) |
115 |
surfdic(i,j) = max(0.4 _d 0, |
surfdic(i,j) = max(0.4 _d 0, |
116 |
& min(10. _d 0, PTR_CO2(i,j,kLev))) |
& min(10. _d 0, PTR_CO2(i,j,kLev))) |
117 |
#else |
#else |
118 |
surftemp(i,j) = theta(i,j,kLev,bi,bj) |
surftemp(i,j) = theta(i,j,kLev,bi,bj) |
127 |
I bi,bj,iMin,iMax,jMin,jMax,myThid) |
I bi,bj,iMin,iMax,jMin,jMax,myThid) |
128 |
C==================================================================== |
C==================================================================== |
129 |
|
|
130 |
DO j=jmin,jmax |
DO j=jMin,jMax |
131 |
DO i=imin,imax |
DO i=iMin,iMax |
132 |
C Compute AtmosP and Kwexch_Pre which are re-used for flux of O2 |
C Compute AtmosP and Kwexch_Pre which are re-used for flux of O2 |
133 |
|
|
134 |
#ifdef USE_PLOAD |
#ifdef USE_PLOAD |
151 |
|
|
152 |
c pCO2 solver... |
c pCO2 solver... |
153 |
C$TAF LOOP = parallel |
C$TAF LOOP = parallel |
154 |
DO j=jmin,jmax |
DO j=jMin,jMax |
155 |
C$TAF LOOP = parallel |
C$TAF LOOP = parallel |
156 |
DO i=imin,imax |
DO i=iMin,iMax |
157 |
|
|
158 |
IF ( maskC(i,j,kLev,bi,bj).NE.0. _d 0 ) THEN |
IF ( maskC(i,j,kLev,bi,bj).NE.0. _d 0 ) THEN |
159 |
CALL CALC_PCO2_APPROX( |
CALL CALC_PCO2_APPROX( |
167 |
I ak0(i,j,bi,bj), fugf(i,j,bi,bj), |
I ak0(i,j,bi,bj), fugf(i,j,bi,bj), |
168 |
I ff(i,j,bi,bj), |
I ff(i,j,bi,bj), |
169 |
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), |
170 |
U pH(i,j,bi,bj),pCO2(i,j,bi,bj), |
U pH(i,j,bi,bj),pCO2(i,j,bi,bj),co3dummy, |
171 |
I i,j,kLev,bi,bj,myIter,myThid ) |
I i,j,kLev,bi,bj,myIter,myThid ) |
172 |
ELSE |
ELSE |
173 |
pCO2(i,j,bi,bj)=0. _d 0 |
pCO2(i,j,bi,bj)=0. _d 0 |
175 |
ENDDO |
ENDDO |
176 |
ENDDO |
ENDDO |
177 |
|
|
178 |
DO j=jmin,jmax |
DO j=jMin,jMax |
179 |
DO i=imin,imax |
DO i=iMin,iMax |
180 |
|
|
181 |
IF ( maskC(i,j,kLev,bi,bj).NE.0. _d 0 ) THEN |
IF ( maskC(i,j,kLev,bi,bj).NE.0. _d 0 ) THEN |
182 |
C calculate SCHMIDT NO. for CO2 |
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) |
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) = |
247 |
ENDDO |
ENDDO |
248 |
|
|
249 |
C update tendency |
C update tendency |
250 |
DO j=jmin,jmax |
DO j=jMin,jMax |
251 |
DO i=imin,imax |
DO i=iMin,iMax |
252 |
GDC(i,j)= recip_drF(kLev)*recip_hFacC(i,j,kLev,bi,bj) |
GDC(i,j)= recip_drF(kLev)*recip_hFacC(i,j,kLev,bi,bj) |
253 |
& *(FluxCO2(i,j,bi,bj) |
& *(FluxCO2(i,j,bi,bj) |
254 |
#ifdef ALLOW_OLD_VIRTUALFLUX |
#ifdef ALLOW_OLD_VIRTUALFLUX |