/[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.6 by dfer, Thu Apr 3 22:49:53 2008 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)               pHlocal = 7.9 _d 0
              pHlocal = 7.9d0  
96    
97  CMICK - TEMPORARY!!!!!  CMICK - TEMPORARY!!!!!
98  CMICK silica = fixed  CMICK silica = fixed
# Line 155  c     &                bt(i,j,bi,bj),st( Line 154  c     &                bt(i,j,bi,bj),st(
154  c             end if  c             end if
155  cmick....................................................  cmick....................................................
156             else             else
157               omegaC(i,j,k,bi,bj) = 0.0d0               omegaC(i,j,k,bi,bj) = 0. _d 0
158             endif             endif
159    
            pHlast(i,j,k,bi,bj) = pHlocal  
   
160           ENDDO           ENDDO
161          ENDDO          ENDDO
162    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22