/[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.2 by stephd, Fri Dec 16 21:07:53 2005 UTC revision 1.7 by dfer, Fri Apr 4 21:32:36 2008 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "DIC_OPTIONS.h"  #include "DIC_OPTIONS.h"
5  #include "GCHEM_OPTIONS.h"  #include "GCHEM_OPTIONS.h"
6    
# Line 5  CBOP Line 8  CBOP
8  C !ROUTINE: CAR_FLUX  C !ROUTINE: CAR_FLUX
9    
10  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
11        SUBROUTINE CALCITE_SATURATION(        SUBROUTINE CALCITE_SATURATION(PTR_DIC, PTR_ALK, PTR_PO4,
12       I           bi,bj,imin,imax,jmin,jmax,       I           bi,bj,imin,imax,jmin,jmax,
13       I           myIter,myTime,myThid)       I           myIter,myTime,myThid)
14    
# Line 19  C !USES: =============================== Line 22  C !USES: ===============================
22  #include "EEPARAMS.h"  #include "EEPARAMS.h"
23  #include "PARAMS.h"  #include "PARAMS.h"
24  #include "GRID.h"  #include "GRID.h"
25  #include "DIC_BIOTIC.h"  #include "DIC_VARS.h"
 #include "DIC_ABIOTIC.h"  
 #ifdef ALLOW_PTRACERS  
 #include "PTRACERS_SIZE.h"  
 #include "PTRACERS.h"  
 #endif  
26    
27  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
28  C  myThid               :: thread number  C  myThid               :: thread number
29  C  myIter               :: current timestep  C  myIter               :: current timestep
30  C  myTime               :: current time  C  myTime               :: current time
31  C  bioac                :: biological productivity  C  bioac                :: biological productivity
32           _RL  PTR_DIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
33           _RL  PTR_ALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
34           _RL  PTR_PO4(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
35    
36        INTEGER myIter        INTEGER myIter
37        _RL myTime        _RL myTime
38        INTEGER myThid        INTEGER myThid
# Line 69  cmick................................... Line 71  cmick...................................
71    
72  c determine carbonate ion concentration through full domain  c determine carbonate ion concentration through full domain
73  c determine calcite saturation state  c determine calcite saturation state
74         DO k=1,nR         DO k=1,Nr
75    
76          CALL CARBON_COEFFS_PRESSURE_DEP(          CALL CARBON_COEFFS_PRESSURE_DEP(
77       I                       theta,salt,       I                       theta,salt,
# Line 80  c determine calcite saturation state Line 82  c determine calcite saturation state
82          DO j=jmin,jmax          DO j=jmin,jmax
83           DO i=imin,imax           DO i=imin,imax
84    
85             if(hFacC(i,j,k,bi,bj) .gt. 0.0d0)then             if ( hFacC(i,j,k,bi,bj) .gt. 0. _d 0 ) then
86    
87               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
88    
89  c 30 micromol = 0.03 mol m-3  c 30 micromol = 0.03 mol m-3
90               silicaTEST = 0.03d0               silicaTEST = 0.03 _d 0
91               po4local = Ptracer(i,j,k,bi,bj,3)               po4local = PTR_PO4(i,j,k)
92               diclocal = Ptracer(i,j,k,bi,bj,1)               diclocal = PTR_DIC(i,j,k)
93               alklocal = Ptracer(i,j,k,bi,bj,2)               alklocal = PTR_ALK(i,j,k)
94  c            pHlocal = pHlast(i,j,k,bi,bj)               pHlocal = 7.9 _d 0
              pHlocal = 7.9d0  
95    
96  CMICK - TEMPORARY!!!!!  CMICK - TEMPORARY!!!!!
97  CMICK silica = fixed  CMICK silica = fixed
# 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.2  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22