/[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.8 by dfer, Mon Apr 7 20:31:16 2008 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
8    
9  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
10        SUBROUTINE CALCITE_SATURATION(        SUBROUTINE CALCITE_SATURATION(PTR_DIC, PTR_ALK, PTR_PO4,
11       I           bi,bj,imin,imax,jmin,jmax,       I           bi,bj,imin,imax,jmin,jmax,
12       I           myIter,myTime,myThid)       I           myIter,myTime,myThid)
13    
# 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"  
 #ifdef ALLOW_PTRACERS  
 #include "PTRACERS_SIZE.h"  
 #include "PTRACERS.h"  
 #endif  
25    
26  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
27  C  myThid               :: thread number  C  myThid               :: thread number
28  C  myIter               :: current timestep  C  myIter               :: current timestep
29  C  myTime               :: current time  C  myTime               :: current time
30  C  bioac                :: biological productivity  C  bioac                :: biological productivity
31           _RL  PTR_DIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
32           _RL  PTR_ALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
33           _RL  PTR_PO4(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
34    
35        INTEGER myIter        INTEGER myIter
36        _RL myTime        _RL myTime
37        INTEGER myThid        INTEGER myThid
# 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,
# Line 80  c determine calcite saturation state Line 81  c determine calcite saturation state
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 = Ptracer(i,j,k,bi,bj,3)               po4local = PTR_PO4(i,j,k)
91               diclocal = Ptracer(i,j,k,bi,bj,1)               diclocal = PTR_DIC(i,j,k)
92               alklocal = Ptracer(i,j,k,bi,bj,2)               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 152  c     &                bt(i,j,bi,bj),st( Line 152  c     &                bt(i,j,bi,bj),st(
152  c             end if  c             end if
153  cmick....................................................  cmick....................................................
154             else             else
155               omegaC(i,j,k,bi,bj) = 0.0d0               omegaC(i,j,k,bi,bj) = 0. _d 0
156             endif             endif
157    
            pHlast(i,j,k,bi,bj) = pHlocal  
   
158           ENDDO           ENDDO
159          ENDDO          ENDDO
160    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22