/[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.3 by stephd, Tue Nov 28 21:16:03 2006 UTC revision 1.11 by dfer, Fri Oct 7 21:36:39 2011 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "DIC_OPTIONS.h"  #include "DIC_OPTIONS.h"
 #include "GCHEM_OPTIONS.h"  
5    
6  CBOP  CBOP
7  C !ROUTINE: CAR_FLUX  C !ROUTINE: CAR_FLUX
# Line 19  C !USES: =============================== Line 21  C !USES: ===============================
21  #include "EEPARAMS.h"  #include "EEPARAMS.h"
22  #include "PARAMS.h"  #include "PARAMS.h"
23  #include "GRID.h"  #include "GRID.h"
24  #include "DIC_BIOTIC.h"  #include "DIC_VARS.h"
 #include "DIC_ABIOTIC.h"  
25    
26  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
27  C  myThid               :: thread number  C  myThid               :: thread number
# Line 69  cmick................................... Line 70  cmick...................................
70    
71  c determine carbonate ion concentration through full domain  c determine carbonate ion concentration through full domain
72  c determine calcite saturation state  c determine calcite saturation state
73         DO k=1,nR         DO k=1,Nr
74    
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
82           DO i=imin,imax           DO i=imin,imax
83    
84             if(hFacC(i,j,k,bi,bj) .gt. 0.0d0)then             if ( hFacC(i,j,k,bi,bj) .gt. 0. _d 0 ) then
85    
86               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
87    
88  c 30 micromol = 0.03 mol m-3  c 30 micromol = 0.03 mol m-3
89               silicaTEST = 0.03d0               silicaTEST = 0.03 _d 0
90               po4local = PTR_PO4(i,j,k)               po4local = PTR_PO4(i,j,k)
91               diclocal = PTR_DIC(i,j,k)               diclocal = PTR_DIC(i,j,k)
92               alklocal = PTR_ALK(i,j,k)               alklocal = PTR_ALK(i,j,k)
93  c            pHlocal = pHlast(i,j,k,bi,bj)               pHlocal = 7.9 _d 0
              pHlocal = 7.9d0  
94    
95  CMICK - TEMPORARY!!!!!  CMICK - TEMPORARY!!!!!
96  CMICK silica = fixed  CMICK silica = fixed
# 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,       U          pHlocal,pCO2local,carbonate,
119       U          carbonate )       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
# Line 152  c     &                bt(i,j,bi,bj),st( Line 153  c     &                bt(i,j,bi,bj),st(
153  c             end if  c             end if
154  cmick....................................................  cmick....................................................
155             else             else
156               omegaC(i,j,k,bi,bj) = 0.0d0               omegaC(i,j,k,bi,bj) = 0. _d 0
157             endif             endif
158    
            pHlast(i,j,k,bi,bj) = pHlocal  
   
159           ENDDO           ENDDO
160          ENDDO          ENDDO
161    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22