/[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.4 by jmc, Tue Oct 9 00:01:42 2007 UTC revision 1.5 by dfer, Fri Oct 26 21:08:12 2007 UTC
# Line 72  cmick................................... Line 72  cmick...................................
72    
73  c determine carbonate ion concentration through full domain  c determine carbonate ion concentration through full domain
74  c determine calcite saturation state  c determine calcite saturation state
75         DO k=1,nR         DO k=1,Nr
76    
77          CALL CARBON_COEFFS_PRESSURE_DEP(          CALL CARBON_COEFFS_PRESSURE_DEP(
78       I                       theta,salt,       I                       theta,salt,
# Line 83  c determine calcite saturation state Line 83  c determine calcite saturation state
83          DO j=jmin,jmax          DO j=jmin,jmax
84           DO i=imin,imax           DO i=imin,imax
85    
86             if(hFacC(i,j,k,bi,bj) .gt. 0.0d0)then             if ( hFacC(i,j,k,bi,bj) .gt. 0. _d 0 ) then
87    
88               calcium = 1.028d-2*salt(i,j,k,bi,bj)/35.0               calcium = 1.028 _d -2*salt(i,j,k,bi,bj)/35. _d 0
89    
90  c 30 micromol = 0.03 mol m-3  c 30 micromol = 0.03 mol m-3
91               silicaTEST = 0.03d0               silicaTEST = 0.03 _d 0
92               po4local = PTR_PO4(i,j,k)               po4local = PTR_PO4(i,j,k)
93               diclocal = PTR_DIC(i,j,k)               diclocal = PTR_DIC(i,j,k)
94               alklocal = PTR_ALK(i,j,k)               alklocal = PTR_ALK(i,j,k)
95  c            pHlocal = pHlast(i,j,k,bi,bj)  c            pHlocal = pHlast(i,j,k,bi,bj)
96               pHlocal = 7.9d0               pHlocal = 7.9 _d 0
97    
98  CMICK - TEMPORARY!!!!!  CMICK - TEMPORARY!!!!!
99  CMICK silica = fixed  CMICK silica = fixed
# Line 155  c     &                bt(i,j,bi,bj),st( Line 155  c     &                bt(i,j,bi,bj),st(
155  c             end if  c             end if
156  cmick....................................................  cmick....................................................
157             else             else
158               omegaC(i,j,k,bi,bj) = 0.0d0               omegaC(i,j,k,bi,bj) = 0. _d 0
159             endif             endif
160    
161             pHlast(i,j,k,bi,bj) = pHlocal             pHlast(i,j,k,bi,bj) = pHlocal
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22