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

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

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


Revision 1.6 - (show annotations) (download)
Tue Jul 13 18:03:31 2004 UTC (19 years, 11 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint57o_post, checkpoint57m_post, checkpoint55c_post, checkpoint54e_post, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint57b_post, checkpoint57c_pre, checkpoint55j_post, checkpoint56b_post, checkpoint57i_post, checkpoint57e_post, checkpoint55h_post, checkpoint57g_pre, checkpoint55b_post, checkpoint54d_post, checkpoint56c_post, checkpoint55, checkpoint57f_pre, checkpoint57a_post, checkpoint54f_post, checkpoint55g_post, checkpoint55f_post, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, eckpoint57e_pre, checkpoint57h_done, checkpoint57n_post, checkpoint57p_post, checkpoint57f_post, checkpoint57q_post, checkpoint57c_post, checkpoint55e_post, checkpoint55a_post, checkpoint54c_post, checkpoint57j_post, checkpoint57h_pre, checkpoint57l_post, checkpoint57h_post, checkpoint56a_post, checkpoint55d_post
Changes since 1.5: +4 -0 lines
max number of passive tracers is now defined in PTRACERS_SIZE.h

1 C $Header: $
2 C $Name: $
3
4 #include "DIC_OPTIONS.h"
5 #include "GCHEM_OPTIONS.h"
6
7 CBOP
8 C !ROUTINE: O2_SURFFORCING
9
10 C !INTERFACE: ==========================================================
11 SUBROUTINE O2_SURFFORCING( PTR_O2, SGO2,
12 I bi,bj,iMin,iMax,jMin,jMax,
13 I myIter, myTime, myThid )
14
15 C !DESCRIPTION:
16 C Calculate the oxygen air-sea flux terms
17
18 C !USES: ===============================================================
19 IMPLICIT NONE
20 #include "SIZE.h"
21 #include "DYNVARS.h"
22 #include "EEPARAMS.h"
23 #include "PARAMS.h"
24 #include "GRID.h"
25 #include "FFIELDS.h"
26 #ifdef DIC_BIOTIC
27 #include "DIC_ABIOTIC.h"
28 #include "PTRACERS_SIZE.h"
29 #include "PTRACERS.h"
30 #endif
31
32
33 c !INPUT PARAMETERS: ===================================================
34 C myThid :: thread number
35 C myIter :: current timestep
36 C myTime :: current time
37 C PTR_O2 :: oxygen tracer field
38 _RL myTime
39 _RL PTR_O2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
40 INTEGER iMin,iMax,jMin,jMax, bi, bj
41 INTEGER myIter, myThid
42
43 c !OUTPUT PARAMETERS: ===================================================
44 C SGO2 :: air-sea exchange of oxygen
45 _RL SGO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
46
47 #ifdef ALLOW_PTRACERS
48
49
50 C !LOCAL VARIABLES: ===================================================
51 C I, J, K - Loop counters
52 INTEGER I,J,K
53 C Solubility relation coefficients
54 _RL SchmidtNoO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
55 _RL O2sat(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
56 _RL Kwexch(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
57 _RL FluxO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
58 _RL AtmosO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
59 _RL aTT
60 _RL aTK
61 _RL aTS
62 _RL aTS2
63 _RL aTS3
64 _RL aTS4
65 _RL aTS5
66 _RL o2s
67 _RL ttemp
68 _RL stemp
69 _RL oCnew
70 CEOP
71
72
73 K=1
74
75 C calculate SCHMIDT NO. for O2
76 DO j=jMin,jMax
77 DO i=iMin,iMax
78 IF (hFacC(i,j,k,bi,bj).NE.0.) THEN
79 SchmidtNoO2(i,j) =
80 & sox1
81 & + sox2 * theta(i,j,k,bi,bj)
82 & + sox3 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj)
83 & + sox4 * theta(i,j,k,bi,bj)*theta(i,j,k,bi,bj)
84 & *theta(i,j,k,bi,bj)
85
86 C Determine surface flux of O2
87 C exchange coeff, accounting for ice cover and Schmidt no.
88 Kwexch(i,j) =
89 & pisvel(i,j,bi,bj)
90 & / sqrt(SchmidtNoO2(i,j)/660.0)
91 c ice influence
92 cQQ Kwexch(i,j) =(1.d0-Fice(i,j,bi,bj))*Kwexch(i,j)
93
94 ttemp = theta(i,j,k,bi,bj)
95 stemp = salt(i,j,k,bi,bj)
96 C determine saturation O2
97 C using Garcia and Gordon (1992), L&O (mistake in original???)
98 aTT = 298.15-ttemp
99 aTK = 273.15+ttemp
100 aTS = log(aTT/aTK)
101 aTS2 = aTS*aTS
102 aTS3 = aTS2*aTS
103 aTS4 = aTS3*aTS
104 aTS5 = aTS4*aTS
105
106 oCnew = oA0 + oA1*aTS + oA2*aTS2 + oA3*aTS3 +
107 & oA4*aTS4 + oA5*aTS5
108 & + stemp*(oB0 + oB1*aTS + oB2*aTS2 + oB3*aTS3)
109 & + oC0*(stemp*stemp)
110
111 o2s = EXP(oCnew)
112
113 c Convert from ml/l to mol/m^3
114 O2sat(i,j) = o2s/22391.6*1000.0
115
116 c Determine flux, inc. correction for local atmos surface pressure
117 cQQ PTR_O2?
118 FluxO2(i,j) = maskC(i,j,k,bi,bj)*Kwexch(i,j)*
119 & (atmosP(i,j,bi,bj)*O2sat(i,j)
120 & - PTR_O2(i,j,1))
121 ELSE
122 FluxO2(i,j) = 0.d0
123 ENDIF
124
125
126 END DO
127 END DO
128
129 C update surface tendencies
130 DO j=jMin,jMax
131 DO i=iMin,iMax
132 SGO2(i,j)= maskC(i,j,1,bi,bj)*FluxO2(i,j)*recip_drF(1)
133 ENDDO
134 ENDDO
135 #endif
136
137
138 RETURN
139 END
140

  ViewVC Help
Powered by ViewVC 1.1.22