/[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.30 by jmc, Sun Jan 11 20:02:41 2015 UTC
# Line 9  C !ROUTINE: DIC_SURFFORCING Line 9  C !ROUTINE: DIC_SURFFORCING
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:
# 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 79  c        ENDDO Line 80  c        ENDDO
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)
# Line 108  cQQQQ check ptracer numbers Line 108  cQQQQ check ptracer numbers
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)
# Line 127  C FOR NON-INTERACTIVE Si Line 127  C FOR NON-INTERACTIVE Si
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
# Line 151  C Schmidt number is accounted for later Line 151  C Schmidt number is accounted for later
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(
# Line 167  C$TAF LOOP = parallel Line 167  C$TAF LOOP = parallel
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
# Line 175  C$TAF LOOP = parallel Line 175  C$TAF LOOP = parallel
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
# 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) =
# Line 247  c Line 247  c
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
Legend:
Removed from v.1.27  
changed lines
  Added in v.1.30
  ViewVC Help
Powered by ViewVC 1.1.22