/[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.4 by jmc, Tue Oct 9 00:01:42 2007 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 21  C !USES: =============================== Line 24  C !USES: ===============================
24  #include "GRID.h"  #include "GRID.h"
25  #include "DIC_BIOTIC.h"  #include "DIC_BIOTIC.h"
26  #include "DIC_ABIOTIC.h"  #include "DIC_ABIOTIC.h"
 #ifdef ALLOW_PTRACERS  
 #include "PTRACERS_SIZE.h"  
 #include "PTRACERS.h"  
 #endif  
27    
28  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
29  C  myThid               :: thread number  C  myThid               :: thread number
30  C  myIter               :: current timestep  C  myIter               :: current timestep
31  C  myTime               :: current time  C  myTime               :: current time
32  C  bioac                :: biological productivity  C  bioac                :: biological productivity
33           _RL  PTR_DIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
34           _RL  PTR_ALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
35           _RL  PTR_PO4(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
36    
37        INTEGER myIter        INTEGER myIter
38        _RL myTime        _RL myTime
39        INTEGER myThid        INTEGER myThid
# Line 86  c determine calcite saturation state Line 89  c determine calcite saturation state
89    
90  c 30 micromol = 0.03 mol m-3  c 30 micromol = 0.03 mol m-3
91               silicaTEST = 0.03d0               silicaTEST = 0.03d0
92               po4local = Ptracer(i,j,k,bi,bj,3)               po4local = PTR_PO4(i,j,k)
93               diclocal = Ptracer(i,j,k,bi,bj,1)               diclocal = PTR_DIC(i,j,k)
94               alklocal = Ptracer(i,j,k,bi,bj,2)               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.9d0
97    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22