/[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.14 by dfer, Mon Aug 13 02:29:40 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 anomalous pressure pLoad (in Pa) from atmospheric model
142    C to total pressure (in Atm)
143    C Note: it is assumed the reference atmospheric pressure is 1Atm=1013mb
144    C       rather than the actual ref. pressure from Atm. model so that on
145    C       average AtmosP is about 1 Atm.
146                   AtmosP(i,j,bi,bj)= 1. _d 0 + pLoad(i,j,bi,bj)/Pa2Atm
147    #endif
148    
149  C Determine surface flux (FDIC)  C Determine surface flux (FDIC)
150  C first correct pCO2at for surface atmos pressure  C first correct pCO2at for surface atmos pressure
151                pCO2sat(i,j) =                pCO2sat(i,j) =

Legend:
Removed from v.1.12  
changed lines
  Added in v.1.14

  ViewVC Help
Powered by ViewVC 1.1.22