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

Annotation of /MITgcm/pkg/dic/o2_surfforcing.F

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

Revision 1.22 - (hide annotations) (download)
Mon Mar 22 02:25:20 2010 UTC (14 years, 2 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint64, checkpoint65, checkpoint63, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, HEAD
Changes since 1.21: +12 -12 lines 
change a comment
1 jmc 1.22 C $Header: /u/gcmpack/MITgcm/pkg/dic/o2_surfforcing.F,v 1.21 2008/04/07 20:31:17 dfer Exp $
2 jmc 1.6 C $Name: $
3    
4 edhill 1.3 #include "DIC_OPTIONS.h"
5 stephd 1.1
6 stephd 1.4 CBOP
7     C !ROUTINE: O2_SURFFORCING
8    
9     C !INTERFACE: ==========================================================
10 stephd 1.5 SUBROUTINE O2_SURFFORCING( PTR_O2, SGO2,
11 stephd 1.1 I bi,bj,iMin,iMax,jMin,jMax,
12     I myIter, myTime, myThid )
13 stephd 1.4
14     C !DESCRIPTION:
15 jmc 1.22 C Calculate the oxygen air-sea flux terms
16 stephd 1.4
17     C !USES: ===============================================================
18 stephd 1.1 IMPLICIT NONE
19     #include "SIZE.h"
20     #include "DYNVARS.h"
21     #include "EEPARAMS.h"
22     #include "PARAMS.h"
23     #include "GRID.h"
24     #include "FFIELDS.h"
25 dfer 1.20 #include "DIC_VARS.h"
26 stephd 1.1
27 stephd 1.4 c !INPUT PARAMETERS: ===================================================
28     C myThid :: thread number
29     C myIter :: current timestep
30     C myTime :: current time
31     C PTR_O2 :: oxygen tracer field
32 stephd 1.1 _RL myTime
33     _RL PTR_O2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
34     INTEGER iMin,iMax,jMin,jMax, bi, bj
35 stephd 1.5 INTEGER myIter, myThid
36 stephd 1.1
37 stephd 1.4 c !OUTPUT PARAMETERS: ===================================================
38 stephd 1.5 C SGO2 :: air-sea exchange of oxygen
39     _RL SGO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
40 stephd 1.4
41 stephd 1.1 #ifdef ALLOW_PTRACERS
42    
43 stephd 1.12 #ifdef ALLOW_O2
44 stephd 1.1
45 stephd 1.4 C !LOCAL VARIABLES: ===================================================
46 stephd 1.1 C I, J, K - Loop counters
47     INTEGER I,J,K
48     C Solubility relation coefficients
49     _RL SchmidtNoO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
50     _RL O2sat(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
51     _RL Kwexch(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
52     _RL FluxO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
53 jmc 1.22 _RL aTT
54 stephd 1.1 _RL aTK
55     _RL aTS
56     _RL aTS2
57     _RL aTS3
58     _RL aTS4
59     _RL aTS5
60     _RL o2s
61     _RL ttemp
62     _RL stemp
63     _RL oCnew
64 stephd 1.4 CEOP
65 stephd 1.1
66    
67     K=1
68    
69     C calculate SCHMIDT NO. for O2
70 stephd 1.9 DO j=jmin,jmax
71     DO i=imin,imax
72 dfer 1.16 IF (maskC(i,j,k,bi,bj).NE.0.) THEN
73     ttemp = theta(i,j,k,bi,bj)
74     stemp = salt(i,j,k,bi,bj)
75    
76 jmc 1.22 SchmidtNoO2(i,j) =
77     & sox1
78 dfer 1.16 & + sox2 * ttemp
79 jmc 1.22 & + sox3 * ttemp*ttemp
80 dfer 1.16 & + sox4 * ttemp*ttemp*ttemp
81 stephd 1.1
82 jmc 1.22 C Determine surface flux of O2
83 dfer 1.17 C exchange coeff accounting for ice cover and Schmidt no.
84     C Kwexch_Pre= pisvel*(1-fice): previously computed in dic_surfforcing.F
85    
86     Kwexch(i,j) = Kwexch_Pre(i,j,bi,bj)
87     & / sqrt(SchmidtNoO2(i,j)/660.0 _d 0)
88 stephd 1.1
89     C determine saturation O2
90 jmc 1.22 C using Garcia and Gordon (1992), L&O (mistake in original ?)
91 dfer 1.15 aTT = 298.15 _d 0 -ttemp
92     aTK = 273.15 _d 0 +ttemp
93 stephd 1.1 aTS = log(aTT/aTK)
94     aTS2 = aTS*aTS
95     aTS3 = aTS2*aTS
96     aTS4 = aTS3*aTS
97     aTS5 = aTS4*aTS
98    
99 jmc 1.22 oCnew = oA0 + oA1*aTS + oA2*aTS2 + oA3*aTS3 +
100 stephd 1.1 & oA4*aTS4 + oA5*aTS5
101     & + stemp*(oB0 + oB1*aTS + oB2*aTS2 + oB3*aTS3)
102     & + oC0*(stemp*stemp)
103    
104     o2s = EXP(oCnew)
105    
106     c Convert from ml/l to mol/m^3
107 dfer 1.15 O2sat(i,j) = o2s/22391.6 _d 0 * 1. _d 3
108 stephd 1.1
109 dfer 1.17 C Determine flux, inc. correction for local atmos surface pressure
110 dfer 1.16 FluxO2(i,j) = Kwexch(i,j)*
111 jmc 1.22 & (AtmosP(i,j,bi,bj)*O2sat(i,j)
112     & - PTR_O2(i,j,K))
113 stephd 1.1 ELSE
114 dfer 1.15 FluxO2(i,j) = 0. _d 0
115 stephd 1.1 ENDIF
116    
117    
118     END DO
119     END DO
120    
121 jmc 1.22 C update surface tendencies
122 stephd 1.9 DO j=jmin,jmax
123     DO i=imin,imax
124 dfer 1.16 SGO2(i,j)= FluxO2(i,j)
125     & *recip_drF(K) * recip_hFacC(i,j,K,bi,bj)
126 stephd 1.1 ENDDO
127     ENDDO
128     #endif
129 stephd 1.12 #endif
130 stephd 1.1
131    
132     RETURN
133     END
134    
  ViewVC Help
Powered by ViewVC 1.1.22