75 |
CALL CARBON_COEFFS_PRESSURE_DEP( |
CALL CARBON_COEFFS_PRESSURE_DEP( |
76 |
I theta,salt, |
I theta,salt, |
77 |
I bi,bj,iMin,iMax,jMin,jMax, |
I bi,bj,iMin,iMax,jMin,jMax, |
78 |
I k) |
I k,myThid) |
79 |
|
|
80 |
|
|
81 |
DO j=jmin,jmax |
DO j=jmin,jmax |
114 |
I aks(i,j,bi,bj),akb(i,j,bi,bj),akw(i,j,bi,bj), |
I aks(i,j,bi,bj),akb(i,j,bi,bj),akw(i,j,bi,bj), |
115 |
I aksi(i,j,bi,bj),akf(i,j,bi,bj),ff(i,j,bi,bj), |
I aksi(i,j,bi,bj),akf(i,j,bi,bj),ff(i,j,bi,bj), |
116 |
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), |
117 |
U pHlocal,pCO2local, |
U pHlocal,pCO2local,carbonate, |
118 |
U carbonate ) |
I myThid ) |
119 |
c........................................................ |
c........................................................ |
120 |
c if(i .eq. 76 .and. j .eq. 36 .and. k .eq. 15) then |
c if(i .eq. 76 .and. j .eq. 36 .and. k .eq. 15) then |
121 |
c write(6,*)'Iteration, pH = ',CO3iter,pHlocal |
c write(6,*)'Iteration, pH = ',CO3iter,pHlocal |