75 |
|
|
76 |
C ================================================================= |
C ================================================================= |
77 |
C determine inorganic carbon chem coefficients |
C determine inorganic carbon chem coefficients |
78 |
DO j=1-OLy,sNy+OLy |
DO j=jmin,jmax |
79 |
DO i=1-OLx,sNx+OLx |
DO i=imin,imax |
80 |
|
|
81 |
#ifdef DIC_BIOTIC |
#ifdef DIC_BIOTIC |
82 |
cQQQQ check ptracer numbers |
cQQQQ check ptracer numbers |
101 |
|
|
102 |
c pCO2 solver... |
c pCO2 solver... |
103 |
C$TAF LOOP = parallel |
C$TAF LOOP = parallel |
104 |
DO j=1-OLy,sNy+OLy |
DO j=jmin,jmax |
105 |
C$TAF LOOP = parallel |
C$TAF LOOP = parallel |
106 |
DO i=1-OLx,sNx+OLx |
DO i=imin,imax |
107 |
|
|
108 |
IF(maskC(i,j,kLev,bi,bj) .NE. 0.)THEN |
IF(maskC(i,j,kLev,bi,bj) .NE. 0.)THEN |
109 |
CALL CALC_PCO2_APPROX( |
CALL CALC_PCO2_APPROX( |
122 |
ENDDO |
ENDDO |
123 |
ENDDO |
ENDDO |
124 |
|
|
125 |
DO j=1-OLy,sNy+OLy |
DO j=jmin,jmax |
126 |
DO i=1-OLx,sNx+OLx |
DO i=imin,imax |
127 |
|
|
128 |
IF (maskC(i,j,kLev,bi,bj).NE.0.) THEN |
IF (maskC(i,j,kLev,bi,bj).NE.0.) THEN |
129 |
C calculate SCHMIDT NO. for CO2 |
C calculate SCHMIDT NO. for CO2 |
184 |
ENDDO |
ENDDO |
185 |
|
|
186 |
C update tendency |
C update tendency |
187 |
DO j=1-OLy,sNy+OLy |
DO j=jmin,jmax |
188 |
DO i=1-OLx,sNx+OLx |
DO i=imin,imax |
189 |
GDC(i,j)= maskC(i,j,kLev,bi,bj)*recip_drF(kLev)* |
GDC(i,j)= maskC(i,j,kLev,bi,bj)*recip_drF(kLev)* |
190 |
& recip_hFacC(i,j,kLev,bi,bj)*( |
& recip_hFacC(i,j,kLev,bi,bj)*( |
191 |
& FluxCO2(i,j,bi,bj) + VirtualFlux(i,j) |
& FluxCO2(i,j,bi,bj) + VirtualFlux(i,j) |