/[MITgcm]/MITgcm_contrib/bling/pkg/bling_carbonate_init.F
ViewVC logotype

Diff of /MITgcm_contrib/bling/pkg/bling_carbonate_init.F

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

revision 1.1 by mmazloff, Sun Feb 28 21:53:37 2016 UTC revision 1.2 by mmazloff, Sun May 15 00:30:35 2016 UTC
# Line 163  C determine inorganic carbon chem coeffi Line 163  C determine inorganic carbon chem coeffi
163          IF ( .NOT.pH_isLoaded ) THEN          IF ( .NOT.pH_isLoaded ) THEN
164  C set guess of pH for first step here  C set guess of pH for first step here
165    
 C$TAF LOOP = parallel  
166          DO k=1,Nr          DO k=1,Nr
167          DO j=jMin,jMax          DO j=jMin,jMax
168           DO i=iMin,iMax           DO i=iMin,iMax
# Line 192  C======================================= Line 191  C=======================================
191    
192  c first approximation  c first approximation
193    
 C$TAF LOOP = parallel  
194            DO j=jMin,jMax            DO j=jMin,jMax
 C$TAF LOOP = parallel  
195             DO i=iMin,iMax             DO i=iMin,iMax
196              IF ( maskC(i,j,k,bi,bj) .NE. 0. _d 0) THEN              IF ( maskC(i,j,k,bi,bj) .NE. 0. _d 0) THEN
 C$TAF init dic_surf = static, 10  
197               DO it=1,10               DO it=1,10
 C$TAF STORE pH(i,j,k,bi,bj), dictmp(i,j)      = dic_surf  
 C$TAF STORE alktmp(i,j), phostmp(i,j), sitmp(i,j)   = dic_surf  
198                CALL CALC_PCO2_APPROX(                CALL CALC_PCO2_APPROX(
199       I        thetatmp(i,j),salttmp(i,j),       I        thetatmp(i,j),salttmp(i,j),
200       I        dictmp(i,j), phostmp(i,j),       I        dictmp(i,j), phostmp(i,j),
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22