/[MITgcm]/MITgcm/pkg/dic/calcite_saturation.F
ViewVC logotype

Diff of /MITgcm/pkg/dic/calcite_saturation.F

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

revision 1.9 by jmc, Wed Apr 9 22:13:15 2008 UTC revision 1.11 by dfer, Fri Oct 7 21:36:39 2011 UTC
# Line 105  CMICK - SO NOW WE ITERATE, UPDATING THE Line 105  CMICK - SO NOW WE ITERATE, UPDATING THE
105  CMICK - SINCE WE CALL THIS FOR DEEP OCEAN INFREQUENTLY (MONTHLY?)  CMICK - SINCE WE CALL THIS FOR DEEP OCEAN INFREQUENTLY (MONTHLY?)
106  CMIKC - CAN AFFORD TO MAKE SEVERAL ITERATIONS...  CMIKC - CAN AFFORD TO MAKE SEVERAL ITERATIONS...
107               DO CO3iter = 1, CO3itermax               DO CO3iter = 1, CO3itermax
108                 CALL CALC_PCO2_APPROX_CO3(                 CALL CALC_PCO2_APPROX(
109       I          theta(i,j,k,bi,bj),salt(i,j,k,bi,bj),       I          theta(i,j,k,bi,bj),salt(i,j,k,bi,bj),
110       I          diclocal, po4local,       I          diclocal, po4local,
111       I          silicaTEST,alklocal,       I          silicaTEST,alklocal,
112       I          ak1(i,j,bi,bj),ak2(i,j,bi,bj),       I          ak1(i,j,bi,bj),ak2(i,j,bi,bj),
113       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),
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),
116         I          ak0(i,j,bi,bj), fugf(i,j,bi,bj), ff(i,j,bi,bj),
117       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),
118       U          pHlocal,pCO2local,carbonate,       U          pHlocal,pCO2local,carbonate,
119       I          myThid )       I          i,j,k,bi,bj,myIter,myThid )
120  c........................................................  c........................................................
121  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
122  c                 write(6,*)'Iteration, pH = ',CO3iter,pHlocal  c                 write(6,*)'Iteration, pH = ',CO3iter,pHlocal
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22