56 |
_RL surfalk(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL surfalk(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
57 |
_RL surfphos(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL surfphos(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
58 |
_RL surfsi(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL surfsi(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
59 |
|
_RL surftemp(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
60 |
|
_RL surfsalt(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
61 |
|
_RL surfdic(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
62 |
#ifdef ALLOW_OLD_VIRTUALFLUX |
#ifdef ALLOW_OLD_VIRTUALFLUX |
63 |
_RL VirtualFlux(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
_RL VirtualFlux(1-OLx:sNx+OLx,1-OLy:sNy+OLy) |
64 |
#endif |
#endif |
87 |
|
|
88 |
#ifdef DIC_BIOTIC |
#ifdef DIC_BIOTIC |
89 |
cQQQQ check ptracer numbers |
cQQQQ check ptracer numbers |
90 |
|
#ifdef DIC_BOUNDS |
91 |
|
surfalk(i,j) = max(0.4 _d 0, |
92 |
|
& min(10. _d 0,PTR_ALK(i,j,klev))) |
93 |
|
& * maskC(i,j,kLev,bi,bj) |
94 |
|
surfphos(i,j) = max(1.0 _d -11, |
95 |
|
& min(1._d -1, PTR_PO4(i,j,klev))) |
96 |
|
& * maskC(i,j,kLev,bi,bj) |
97 |
|
#else |
98 |
surfalk(i,j) = PTR_ALK(i,j,klev) |
surfalk(i,j) = PTR_ALK(i,j,klev) |
99 |
& * maskC(i,j,kLev,bi,bj) |
& * maskC(i,j,kLev,bi,bj) |
100 |
surfphos(i,j) = PTR_PO4(i,j,klev) |
surfphos(i,j) = PTR_PO4(i,j,klev) |
101 |
& * maskC(i,j,kLev,bi,bj) |
& * maskC(i,j,kLev,bi,bj) |
102 |
|
#endif |
103 |
#else |
#else |
104 |
surfalk(i,j) = 2.366595 _d 0 * salt(i,j,kLev,bi,bj)/gsm_s |
surfalk(i,j) = 2.366595 _d 0 * salt(i,j,kLev,bi,bj)/gsm_s |
105 |
& * maskC(i,j,kLev,bi,bj) |
& * maskC(i,j,kLev,bi,bj) |
107 |
#endif |
#endif |
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 |
111 |
|
surftemp(i,j) = max(-4. _d 0, |
112 |
|
& min(50. _d 0, theta(i,j,kLev,bi,bj))) |
113 |
|
surfsalt(i,j) = max(4. _d 0, |
114 |
|
& min(50. _d 0, salt(i,j,kLev,bi,bj))) |
115 |
|
surfdic(i,j) = max(0.4 _d 0, |
116 |
|
& min(10. _d 0, PTR_CO2(i,j,kLev))) |
117 |
|
#else |
118 |
|
surftemp(i,j) = theta(i,j,kLev,bi,bj) |
119 |
|
surfsalt(i,j) = salt(i,j,kLev,bi,bj) |
120 |
|
surfdic(i,j) = PTR_CO2(i,j,kLev) |
121 |
|
#endif |
122 |
ENDDO |
ENDDO |
123 |
ENDDO |
ENDDO |
124 |
|
|
125 |
CALL CARBON_COEFFS( |
CALL CARBON_COEFFS( |
126 |
I theta,salt, |
I surftemp,surfsalt, |
127 |
I bi,bj,iMin,iMax,jMin,jMax,myThid) |
I bi,bj,iMin,iMax,jMin,jMax,myThid) |
128 |
C==================================================================== |
C==================================================================== |
129 |
|
|
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( |
160 |
I theta(i,j,kLev,bi,bj),salt(i,j,kLev,bi,bj), |
I surftemp(i,j),surfsalt(i,j), |
161 |
I PTR_CO2(i,j,kLev), surfphos(i,j), |
I surfdic(i,j), surfphos(i,j), |
162 |
I surfsi(i,j),surfalk(i,j), |
I surfsi(i,j),surfalk(i,j), |
163 |
I ak1(i,j,bi,bj),ak2(i,j,bi,bj), |
I ak1(i,j,bi,bj),ak2(i,j,bi,bj), |
164 |
I ak1p(i,j,bi,bj),ak2p(i,j,bi,bj),ak3p(i,j,bi,bj), |
I ak1p(i,j,bi,bj),ak2p(i,j,bi,bj),ak3p(i,j,bi,bj), |