/[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.3 by stephd, Tue Nov 28 21:16:03 2006 UTC
# Line 5  CBOP Line 5  CBOP
5  C !ROUTINE: CAR_FLUX  C !ROUTINE: CAR_FLUX
6    
7  C !INTERFACE: ==========================================================  C !INTERFACE: ==========================================================
8        SUBROUTINE CALCITE_SATURATION(        SUBROUTINE CALCITE_SATURATION(PTR_DIC, PTR_ALK, PTR_PO4,
9       I           bi,bj,imin,imax,jmin,jmax,       I           bi,bj,imin,imax,jmin,jmax,
10       I           myIter,myTime,myThid)       I           myIter,myTime,myThid)
11    
# Line 21  C !USES: =============================== Line 21  C !USES: ===============================
21  #include "GRID.h"  #include "GRID.h"
22  #include "DIC_BIOTIC.h"  #include "DIC_BIOTIC.h"
23  #include "DIC_ABIOTIC.h"  #include "DIC_ABIOTIC.h"
 #ifdef ALLOW_PTRACERS  
 #include "PTRACERS_SIZE.h"  
 #include "PTRACERS.h"  
 #endif  
24    
25  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
26  C  myThid               :: thread number  C  myThid               :: thread number
27  C  myIter               :: current timestep  C  myIter               :: current timestep
28  C  myTime               :: current time  C  myTime               :: current time
29  C  bioac                :: biological productivity  C  bioac                :: biological productivity
30           _RL  PTR_DIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
31           _RL  PTR_ALK(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
32           _RL  PTR_PO4(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
33    
34        INTEGER myIter        INTEGER myIter
35        _RL myTime        _RL myTime
36        INTEGER myThid        INTEGER myThid
# Line 86  c determine calcite saturation state Line 86  c determine calcite saturation state
86    
87  c 30 micromol = 0.03 mol m-3  c 30 micromol = 0.03 mol m-3
88               silicaTEST = 0.03d0               silicaTEST = 0.03d0
89               po4local = Ptracer(i,j,k,bi,bj,3)               po4local = PTR_PO4(i,j,k)
90               diclocal = Ptracer(i,j,k,bi,bj,1)               diclocal = PTR_DIC(i,j,k)
91               alklocal = Ptracer(i,j,k,bi,bj,2)               alklocal = PTR_ALK(i,j,k)
92  c            pHlocal = pHlast(i,j,k,bi,bj)  c            pHlocal = pHlast(i,j,k,bi,bj)
93               pHlocal = 7.9d0               pHlocal = 7.9d0
94    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22