/[MITgcm]/MITgcm/pkg/dic/dic_surfforcing.F
ViewVC logotype

Diff of /MITgcm/pkg/dic/dic_surfforcing.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.18 by dfer, Fri Oct 26 21:08:13 2007 UTC revision 1.23 by jmc, Sun Apr 11 22:03:53 2010 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "DIC_OPTIONS.h"  #include "DIC_OPTIONS.h"
5  #include "PTRACERS_OPTIONS.h"  #include "PTRACERS_OPTIONS.h"
 #include "GCHEM_OPTIONS.h"  
6    
7  CBOP  CBOP
8  C !ROUTINE: DIC_SURFFORCING  C !ROUTINE: DIC_SURFFORCING
# Line 25  C !USES: =============================== Line 24  C !USES: ===============================
24  #include "PARAMS.h"  #include "PARAMS.h"
25  #include "GRID.h"  #include "GRID.h"
26  #include "FFIELDS.h"  #include "FFIELDS.h"
27  #include "DIC_ABIOTIC.h"  #include "DIC_VARS.h"
28    
29  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
30  C  myThid               :: thread number  C  myThid               :: thread number
# Line 46  c GDC                   :: tendency due Line 45  c GDC                   :: tendency due
45  #ifdef ALLOW_PTRACERS  #ifdef ALLOW_PTRACERS
46    
47  C !LOCAL VARIABLES: ====================================================  C !LOCAL VARIABLES: ====================================================
48         INTEGER I,J, kLev, it         INTEGER i,j, kLev
49  C Number of iterations for pCO2 solvers...  C Number of iterations for pCO2 solvers...
50  C Solubility relation coefficients  C Solubility relation coefficients
51        _RL SchmidtNoDIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL SchmidtNoDIC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
# Line 66  cccccccccccccccccccccccccccccccccccccccc Line 65  cccccccccccccccccccccccccccccccccccccccc
65    
66        kLev=1        kLev=1
67    
68  c if coupled to atmsopheric model, use the  cc if coupled to atmsopheric model, use the
69  c Co2 value passed from the coupler  cc Co2 value passed from the coupler
70  #ifndef USE_ATMOSCO2  c#ifndef USE_ATMOSCO2
71  C PRE-INDUSTRIAL STEADY STATE pCO2 = 278.0 ppmv  cC PRE-INDUSTRIAL STEADY STATE pCO2 = 278.0 ppmv
72         DO j=1-OLy,sNy+OLy  c       DO j=1-OLy,sNy+OLy
73          DO i=1-OLx,sNx+OLx  c        DO i=1-OLx,sNx+OLx
74             AtmospCO2(i,j,bi,bj)=278.0 _d -6  c           AtmospCO2(i,j,bi,bj)=278.0 _d -6
75          ENDDO  c        ENDDO
76         ENDDO  c       ENDDO
77  #endif  c#endif
78    
79    
80  C =================================================================  C =================================================================
# Line 101  C FOR NON-INTERACTIVE Si Line 100  C FOR NON-INTERACTIVE Si
100    
101           CALL CARBON_COEFFS(           CALL CARBON_COEFFS(
102       I                       theta,salt,       I                       theta,salt,
103       I                       bi,bj,iMin,iMax,jMin,jMax)       I                       bi,bj,iMin,iMax,jMin,jMax,myThid)
104  C====================================================================  C====================================================================
105    
106         DO j=jmin,jmax         DO j=jmin,jmax
# Line 142  C$TAF LOOP = parallel Line 141  C$TAF LOOP = parallel
141       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),
142       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),ff(i,j,bi,bj),
143       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),
144       U        pH(i,j,bi,bj),pCO2(i,j,bi,bj) )       U        pH(i,j,bi,bj),pCO2(i,j,bi,bj),
145         I        myThid )
146            ELSE            ELSE
147              pCO2(i,j,bi,bj)=0. _d 0              pCO2(i,j,bi,bj)=0. _d 0
148            ENDIF            ENDIF
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.23
  ViewVC Help
Powered by ViewVC 1.1.22