/[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.12 by stephd, Tue Nov 28 21:16:03 2006 UTC revision 1.13 by stephd, Mon May 7 15:58:20 2007 UTC
# Line 63  cccccccccccccccccccccccccccccccccccccccc Line 63  cccccccccccccccccccccccccccccccccccccccc
63    
64        kLev=1        kLev=1
65    
66    c if coupled to atmsopheric model, use the
67    c Co2 value passed from the coupler
68    #ifndef USE_ATMOSCO2
69  C PRE-INDUSTRIAL STEADY STATE pCO2 = 278.0 ppmv  C PRE-INDUSTRIAL STEADY STATE pCO2 = 278.0 ppmv
70         DO j=1-OLy,sNy+OLy         DO j=1-OLy,sNy+OLy
71          DO i=1-OLx,sNx+OLx          DO i=1-OLx,sNx+OLx
72             AtmospCO2(i,j,bi,bj)=278.0d-6             AtmospCO2(i,j,bi,bj)=278.0d-6
73          ENDDO          ENDDO
74         ENDDO         ENDDO
75    #endif
76    
77    
78  C =================================================================  C =================================================================
# Line 132  C calculate SCHMIDT NO. for CO2 Line 136  C calculate SCHMIDT NO. for CO2
136       &          + sca4 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj)       &          + sca4 * theta(i,j,kLev,bi,bj)*theta(i,j,kLev,bi,bj)
137       &                *theta(i,j,kLev,bi,bj)       &                *theta(i,j,kLev,bi,bj)
138    
139    c
140    #ifdef USE_PLOAD
141    c convert from Pa to atmos
142                   AtmosP(i,j,bi,bj)=pLoad(i,j,bi,bj)/1.01295e5
143    #endif
144    
145  C Determine surface flux (FDIC)  C Determine surface flux (FDIC)
146  C first correct pCO2at for surface atmos pressure  C first correct pCO2at for surface atmos pressure
147                pCO2sat(i,j) =                pCO2sat(i,j) =
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22