/[MITgcm]/MITgcm_contrib/SOSE/BoxAdj/code_ad/OBCS.h
ViewVC logotype

Contents of /MITgcm_contrib/SOSE/BoxAdj/code_ad/OBCS.h

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


Revision 1.4 - (show annotations) (download)
Tue Apr 19 23:33:00 2011 UTC (14 years, 8 months ago) by mmazloff
Branch: MAIN
CVS Tags: HEAD
Changes since 1.3: +1 -1 lines
File MIME type: text/plain
FILE REMOVED
cleaning

1 C $Header: /u/gcmpack/MITgcm_contrib/SOSE/BoxAdj/code_ad/OBCS.h,v 1.3 2011/04/18 23:47:25 mmazloff Exp $
2 C $Name: $
3
4 #ifdef ALLOW_OBCS
5
6 C useOrlanskiNorth/South/East/West
7 C :: specify Orlanski boundary conditions for northern/
8 C southern/eastern/Western
9 C useStevensNorth/South/East/West
10 C :: use open boundary computations following Stevens (1990)
11 C T/SrelaxStevens :: relaxation time scale (in seconds) for T/S-points
12 C for Stevens boundary conditions
13 C useStevensPhaseVel
14 C :: use phase velocity contribution for open boundary
15 C computations following Stevens (1990), default = true
16 C useStevensAdvection
17 C :: use advective contribution for open boundary
18 C computations following Stevens (1990), default = true
19 C useOBCSsponge :: turns on sponge layer along boundary (def=false)
20 C spongeThickness :: number grid points that make up the sponge layer (def=0)
21 C U/Vrelaxobcsinner/bound :: relaxation time scale (in seconds) for U/V-points
22 C on the boundary (bound) and at the innermost grid point
23 C of the sponge layer (inner); relaxation time scales
24 C in-between are linearly interpolated from these values
25 C useOBCSbalance :: balance the volume flux through boundary
26 C at every time step
27 C useOBCSprescribe :: read boundary conditions from a file
28 C (overrides Orlanski and other boundary values)
29 C OBCSprintDiags :: print boundary values to STDOUT (def=true)
30 C useOBCSYearlyFields :: when reading boundary values by exf, assume yearly
31 C climatology (def=false)
32 C OBCSfixTopo :: check and adjust topography for problematic gradients
33 C across boundaries (def=true)
34 C tileHasOB[N,S,E,W] :: this tile has OB at Northern/Southern/Eastern/Western edge
35 C insideOBmaskFile :: File to specify Inside OB region mask (zero beyond OB).
36 C OB[N,S,E,W][u,v,t,s,a,h,sn,sl,uice,vice]File :: Files with boundary conditions,
37 C the letter combinations mean:
38 C N/S/E/W :: northern/southern/eastern/western boundary
39 C u/v/t/s :: ocean u/v velocities, temperature/salinity
40 C a/h :: sea ice concentration/effective thickness
41 C sn/sl :: effective snow thickness/sea ice salinity
42 C uice/vice :: sea ice u/v drift velocities
43
44 COMMON /PARM_IL_OB/
45 & OBCS_monSelect,
46 & spongeThickness,
47 & useOrlanskiNorth,useOrlanskiSouth,
48 & useOrlanskiEast,useOrlanskiWest,
49 & useStevensNorth,useStevensSouth,
50 & useStevensEast,useStevensWest,
51 & useStevensPhaseVel, useStevensAdvection,
52 & useOBCSsponge, useOBCSbalance, useOBCSprescribe,
53 & OBCSprintDiags, useOBCSYearlyFields,
54 & OBCSfixTopo
55 INTEGER OBCS_monSelect
56 INTEGER spongeThickness
57 LOGICAL useOrlanskiNorth
58 LOGICAL useOrlanskiSouth
59 LOGICAL useOrlanskiEast
60 LOGICAL useOrlanskiWest
61 LOGICAL useStevensNorth
62 LOGICAL useStevensSouth
63 LOGICAL useStevensEast
64 LOGICAL useStevensWest
65 LOGICAL useStevensPhaseVel
66 LOGICAL useStevensAdvection
67 LOGICAL useOBCSsponge
68 LOGICAL useOBCSbalance
69 LOGICAL useOBCSprescribe
70 LOGICAL OBCSprintDiags
71 LOGICAL useOBCSYearlyFields
72 LOGICAL OBCSfixTopo
73
74 C OBCS_balanceFacN/S/E/W :: weighting factor for balancing OB normal flow
75 C OBCS_monitorFreq :: monitor output frequency (s) for OB statistics
76 C OBCS_monSelect :: select group of variables to monitor
77 COMMON /PARM_R_OB/
78 & OBCS_balanceFacN, OBCS_balanceFacS,
79 & OBCS_balanceFacE, OBCS_balanceFacW,
80 & OBCS_monitorFreq,
81 & Urelaxobcsinner,Urelaxobcsbound,
82 & Vrelaxobcsinner,Vrelaxobcsbound,
83 & TrelaxStevens, SrelaxStevens
84 _RL OBCS_balanceFacN, OBCS_balanceFacS
85 _RL OBCS_balanceFacE, OBCS_balanceFacW
86 _RL OBCS_monitorFreq
87 _RS Urelaxobcsinner
88 _RS Urelaxobcsbound
89 _RS Vrelaxobcsinner
90 _RS Vrelaxobcsbound
91 _RS TrelaxStevens
92 _RS SrelaxStevens
93
94 COMMON /OBCS_ACTIVE_TILES/
95 & tileHasOBN, tileHasOBS, tileHasOBE, tileHasOBW
96 LOGICAL tileHasOBN(nSx,nSy)
97 LOGICAL tileHasOBS(nSx,nSy)
98 LOGICAL tileHasOBE(nSx,nSy)
99 LOGICAL tileHasOBW(nSx,nSy)
100
101 COMMON /GRID_IND_OB/
102 & OB_Jn,OB_Js,OB_Ie,OB_Iw
103 INTEGER OB_Jn(1-Olx:sNx+Olx,nSx,nSy)
104 INTEGER OB_Js(1-Olx:sNx+Olx,nSx,nSy)
105 INTEGER OB_Ie(1-Oly:sNy+Oly,nSx,nSy)
106 INTEGER OB_Iw(1-Oly:sNy+Oly,nSx,nSy)
107
108 COMMON /OB_FILES/
109 & insideOBmaskFile,
110 & OBNetaFile,OBSetaFile,OBEetaFile,OBWetaFile,
111 & OBNwFile, OBSwFile, OBEwFile, OBWwFile,
112 & OBNuFile,OBNvFile,OBNtFile,OBNsFile,OBNaFile,OBNhFile,
113 & OBSuFile,OBSvFile,OBStFile,OBSsFile,OBSaFile,OBShFile,
114 & OBEuFile,OBEvFile,OBEtFile,OBEsFile,OBEaFile,OBEhFile,
115 & OBWuFile,OBWvFile,OBWtFile,OBWsFile,OBWaFile,OBWhFile,
116 & OBNslFile,OBSslFile,OBEslFile,OBWslFile,
117 & OBNsnFile,OBSsnFile,OBEsnFile,OBWsnFile,
118 & OBNuiceFile,OBSuiceFile,OBEuiceFile,OBWuiceFile,
119 & OBNviceFile,OBSviceFile,OBEviceFile,OBWviceFile
120 CHARACTER*(MAX_LEN_FNAM)
121 & insideOBmaskFile,
122 & OBNetaFile,OBSetaFile,OBEetaFile,OBWetaFile,
123 & OBNwFile, OBSwFile, OBEwFile, OBWwFile,
124 & OBNuFile,OBNvFile,OBNtFile,OBNsFile,OBNaFile,OBNhFile,
125 & OBSuFile,OBSvFile,OBStFile,OBSsFile,OBSaFile,OBShFile,
126 & OBEuFile,OBEvFile,OBEtFile,OBEsFile,OBEaFile,OBEhFile,
127 & OBWuFile,OBWvFile,OBWtFile,OBWsFile,OBWaFile,OBWhFile,
128 & OBNslFile,OBSslFile,OBEslFile,OBWslFile,
129 & OBNsnFile,OBSsnFile,OBEsnFile,OBWsnFile,
130 & OBNuiceFile,OBSuiceFile,OBEuiceFile,OBWuiceFile,
131 & OBNviceFile,OBSviceFile,OBEviceFile,OBWviceFile
132
133 C-- COMMON /GRID_OB/ Open boundary related stuff
134 C OBNu is the U value imposed at the Northern OB
135 C OBNv is the V value imposed at the Northern OB
136 C OBNt is the T value imposed at the Northern OB
137 C OBNs is the S value imposed at the Northern OB
138 C OBNa is the ice AREA value imposed at the Northern OB
139 C OBNh is the ice HEFF value imposed at the Northern OB
140 C OBNsl is the ice HSALT value imposed at the Northern OB
141 C OBNsn is the ice HSNOW value imposed at the Northern OB
142 C OBNuice is the uice value imposed at the Northern OB
143 C OBNvice is the vice value imposed at the Northern OB
144 C etc
145
146 #ifdef ALLOW_OBCS_PRESCRIBE
147 C OBCS_ldRec :: time-record currently loaded (in temp arrays *[1])
148 COMMON /OBCS_LOAD_I/ OBCS_ldRec
149 INTEGER OBCS_ldRec(nSx,nSy)
150 #endif /* ALLOW_OBCS_PRESCRIBE */
151
152 #ifdef ALLOW_OBCS_NORTH
153 COMMON /GRID_N_OB/
154 & OBNu,OBNv,OBNt,OBNs
155 _RL OBNu (1-Olx:sNx+Olx,Nr,nSx,nSy)
156 _RL OBNv (1-Olx:sNx+Olx,Nr,nSx,nSy)
157 _RL OBNt (1-Olx:sNx+Olx,Nr,nSx,nSy)
158 _RL OBNs (1-Olx:sNx+Olx,Nr,nSx,nSy)
159 #ifdef ALLOW_OBCS_PRESCRIBE
160 COMMON /GRID_N_OB_AUX/
161 & OBNu0,OBNv0,OBNt0,OBNs0,
162 & OBNu1,OBNv1,OBNt1,OBNs1
163 _RL OBNu0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
164 _RL OBNv0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
165 _RL OBNt0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
166 _RL OBNs0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
167 _RL OBNu1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
168 _RL OBNv1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
169 _RL OBNt1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
170 _RL OBNs1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
171 #endif /* ALLOW_OBCS_PRESCRIBE */
172 #ifdef ALLOW_SEAICE
173 COMMON /SEAICE_N_OB/ OBNa,OBNh,OBNa0,OBNh0,OBNa1,OBNh1,
174 & OBNsl,OBNsn,OBNsl0,OBNsn0,OBNsl1,OBNsn1,
175 & OBNuice,OBNvice,OBNuice0,OBNvice0,OBNuice1,OBNvice1
176 _RL OBNa (1-Olx:sNx+Olx,nSx,nSy)
177 _RL OBNh (1-Olx:sNx+Olx,nSx,nSy)
178 _RL OBNa0 (1-Olx:sNx+Olx,nSx,nSy)
179 _RL OBNh0 (1-Olx:sNx+Olx,nSx,nSy)
180 _RL OBNa1 (1-Olx:sNx+Olx,nSx,nSy)
181 _RL OBNh1 (1-Olx:sNx+Olx,nSx,nSy)
182 _RL OBNsl (1-Olx:sNx+Olx,nSx,nSy)
183 _RL OBNsn (1-Olx:sNx+Olx,nSx,nSy)
184 _RL OBNsl0 (1-Olx:sNx+Olx,nSx,nSy)
185 _RL OBNsn0 (1-Olx:sNx+Olx,nSx,nSy)
186 _RL OBNsl1 (1-Olx:sNx+Olx,nSx,nSy)
187 _RL OBNsn1 (1-Olx:sNx+Olx,nSx,nSy)
188 _RL OBNuice (1-Olx:sNx+Olx,nSx,nSy)
189 _RL OBNvice (1-Olx:sNx+Olx,nSx,nSy)
190 _RL OBNuice0 (1-Olx:sNx+Olx,nSx,nSy)
191 _RL OBNvice0 (1-Olx:sNx+Olx,nSx,nSy)
192 _RL OBNuice1 (1-Olx:sNx+Olx,nSx,nSy)
193 _RL OBNvice1 (1-Olx:sNx+Olx,nSx,nSy)
194 #endif /* ALLOW_SEAICE */
195 #endif /* ALLOW_OBCS_NORTH */
196
197 #ifdef ALLOW_OBCS_SOUTH
198 COMMON /GRID_S_OB/
199 & OBSu,OBSv,OBSt,OBSs
200 _RL OBSu (1-Olx:sNx+Olx,Nr,nSx,nSy)
201 _RL OBSv (1-Olx:sNx+Olx,Nr,nSx,nSy)
202 _RL OBSt (1-Olx:sNx+Olx,Nr,nSx,nSy)
203 _RL OBSs (1-Olx:sNx+Olx,Nr,nSx,nSy)
204 #ifdef ALLOW_OBCS_PRESCRIBE
205 COMMON /GRID_S_OB_AUX/
206 & OBSu0,OBSv0,OBSt0,OBSs0,
207 & OBSu1,OBSv1,OBSt1,OBSs1
208 _RL OBSu0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
209 _RL OBSv0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
210 _RL OBSt0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
211 _RL OBSs0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
212 _RL OBSu1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
213 _RL OBSv1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
214 _RL OBSt1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
215 _RL OBSs1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
216 #endif /* ALLOW_OBCS_PRESCRIBE */
217 #ifdef ALLOW_SEAICE
218 COMMON /SEAICE_S_OB/ OBSa,OBSh,OBSa0,OBSh0,OBSa1,OBSh1,
219 & OBSsl,OBSsn,OBSsl0,OBSsn0,OBSsl1,OBSsn1,
220 & OBSuice,OBSvice,OBSuice0,OBSvice0,OBSuice1,OBSvice1
221 _RL OBSa (1-Olx:sNx+Olx,nSx,nSy)
222 _RL OBSh (1-Olx:sNx+Olx,nSx,nSy)
223 _RL OBSa0 (1-Olx:sNx+Olx,nSx,nSy)
224 _RL OBSh0 (1-Olx:sNx+Olx,nSx,nSy)
225 _RL OBSa1 (1-Olx:sNx+Olx,nSx,nSy)
226 _RL OBSh1 (1-Olx:sNx+Olx,nSx,nSy)
227 _RL OBSsl (1-Olx:sNx+Olx,nSx,nSy)
228 _RL OBSsn (1-Olx:sNx+Olx,nSx,nSy)
229 _RL OBSsl0 (1-Olx:sNx+Olx,nSx,nSy)
230 _RL OBSsn0 (1-Olx:sNx+Olx,nSx,nSy)
231 _RL OBSsl1 (1-Olx:sNx+Olx,nSx,nSy)
232 _RL OBSsn1 (1-Olx:sNx+Olx,nSx,nSy)
233 _RL OBSuice (1-Olx:sNx+Olx,nSx,nSy)
234 _RL OBSvice (1-Olx:sNx+Olx,nSx,nSy)
235 _RL OBSuice0 (1-Olx:sNx+Olx,nSx,nSy)
236 _RL OBSvice0 (1-Olx:sNx+Olx,nSx,nSy)
237 _RL OBSuice1 (1-Olx:sNx+Olx,nSx,nSy)
238 _RL OBSvice1 (1-Olx:sNx+Olx,nSx,nSy)
239 #endif /* ALLOW_SEAICE */
240 #endif /* ALLOW_OBCS_SOUTH */
241
242 #ifdef ALLOW_OBCS_EAST
243 COMMON /GRID_E_OB/
244 & OBEu,OBEv,OBEt,OBEs
245 _RL OBEu (1-Oly:sNy+Oly,Nr,nSx,nSy)
246 _RL OBEv (1-Oly:sNy+Oly,Nr,nSx,nSy)
247 _RL OBEt (1-Oly:sNy+Oly,Nr,nSx,nSy)
248 _RL OBEs (1-Oly:sNy+Oly,Nr,nSx,nSy)
249 #ifdef ALLOW_OBCS_PRESCRIBE
250 COMMON /GRID_E_OB_AUX/
251 & OBEu0,OBEv0,OBEt0,OBEs0,
252 & OBEu1,OBEv1,OBEt1,OBEs1
253 _RL OBEu0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
254 _RL OBEv0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
255 _RL OBEt0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
256 _RL OBEs0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
257 _RL OBEu1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
258 _RL OBEv1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
259 _RL OBEt1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
260 _RL OBEs1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
261 #endif /* ALLOW_OBCS_PRESCRIBE */
262 #ifdef ALLOW_SEAICE
263 COMMON /SEAICE_E_OB/ OBEa,OBEh,OBEa0,OBEh0,OBEa1,OBEh1,
264 & OBEsl,OBEsn,OBEsl0,OBEsn0,OBEsl1,OBEsn1,
265 & OBEuice,OBEvice,OBEuice0,OBEvice0,OBEuice1,OBEvice1
266 _RL OBEa (1-Oly:sNy+Oly,nSx,nSy)
267 _RL OBEh (1-Oly:sNy+Oly,nSx,nSy)
268 _RL OBEa0 (1-Oly:sNy+Oly,nSx,nSy)
269 _RL OBEh0 (1-Oly:sNy+Oly,nSx,nSy)
270 _RL OBEa1 (1-Oly:sNy+Oly,nSx,nSy)
271 _RL OBEh1 (1-Oly:sNy+Oly,nSx,nSy)
272 _RL OBEsl (1-Oly:sNy+Oly,nSx,nSy)
273 _RL OBEsn (1-Oly:sNy+Oly,nSx,nSy)
274 _RL OBEsl0 (1-Oly:sNy+Oly,nSx,nSy)
275 _RL OBEsn0 (1-Oly:sNy+Oly,nSx,nSy)
276 _RL OBEsl1 (1-Oly:sNy+Oly,nSx,nSy)
277 _RL OBEsn1 (1-Oly:sNy+Oly,nSx,nSy)
278 _RL OBEuice (1-Oly:sNy+Oly,nSx,nSy)
279 _RL OBEvice (1-Oly:sNy+Oly,nSx,nSy)
280 _RL OBEuice0 (1-Oly:sNy+Oly,nSx,nSy)
281 _RL OBEvice0 (1-Oly:sNy+Oly,nSx,nSy)
282 _RL OBEuice1 (1-Oly:sNy+Oly,nSx,nSy)
283 _RL OBEvice1 (1-Oly:sNy+Oly,nSx,nSy)
284 #endif /* ALLOW_SEAICE */
285 #endif /* ALLOW_OBCS_EAST */
286
287 #ifdef ALLOW_OBCS_WEST
288 COMMON /GRID_W_OB/
289 & OBWu,OBWv,OBWt,OBWs
290 _RL OBWu (1-Oly:sNy+Oly,Nr,nSx,nSy)
291 _RL OBWv (1-Oly:sNy+Oly,Nr,nSx,nSy)
292 _RL OBWt (1-Oly:sNy+Oly,Nr,nSx,nSy)
293 _RL OBWs (1-Oly:sNy+Oly,Nr,nSx,nSy)
294 #ifdef ALLOW_OBCS_PRESCRIBE
295 COMMON /GRID_W_OB_AUX/
296 & OBWu0,OBWv0,OBWt0,OBWs0,
297 & OBWu1,OBWv1,OBWt1,OBWs1
298 _RL OBWu0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
299 _RL OBWv0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
300 _RL OBWt0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
301 _RL OBWs0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
302 _RL OBWu1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
303 _RL OBWv1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
304 _RL OBWt1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
305 _RL OBWs1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
306 #endif /* ALLOW_OBCS_PRESCRIBE */
307 #ifdef ALLOW_SEAICE
308 COMMON /SEAICE_W_OB/ OBWa,OBWh,OBWa0,OBWh0,OBWa1,OBWh1,
309 & OBWsl,OBWsn,OBWsl0,OBWsn0,OBWsl1,OBWsn1,
310 & OBWuice,OBWvice,OBWuice0,OBWvice0,OBWuice1,OBWvice1
311 _RL OBWa (1-Oly:sNy+Oly,nSx,nSy)
312 _RL OBWh (1-Oly:sNy+Oly,nSx,nSy)
313 _RL OBWa0 (1-Oly:sNy+Oly,nSx,nSy)
314 _RL OBWh0 (1-Oly:sNy+Oly,nSx,nSy)
315 _RL OBWa1 (1-Oly:sNy+Oly,nSx,nSy)
316 _RL OBWh1 (1-Oly:sNy+Oly,nSx,nSy)
317 _RL OBWsl (1-Oly:sNy+Oly,nSx,nSy)
318 _RL OBWsn (1-Oly:sNy+Oly,nSx,nSy)
319 _RL OBWsl0 (1-Oly:sNy+Oly,nSx,nSy)
320 _RL OBWsn0 (1-Oly:sNy+Oly,nSx,nSy)
321 _RL OBWsl1 (1-Oly:sNy+Oly,nSx,nSy)
322 _RL OBWsn1 (1-Oly:sNy+Oly,nSx,nSy)
323 _RL OBWuice (1-Oly:sNy+Oly,nSx,nSy)
324 _RL OBWvice (1-Oly:sNy+Oly,nSx,nSy)
325 _RL OBWuice0 (1-Oly:sNy+Oly,nSx,nSy)
326 _RL OBWvice0 (1-Oly:sNy+Oly,nSx,nSy)
327 _RL OBWuice1 (1-Oly:sNy+Oly,nSx,nSy)
328 _RL OBWvice1 (1-Oly:sNy+Oly,nSx,nSy)
329 #endif /* ALLOW_SEAICE */
330 #endif /* ALLOW_OBCS_WEST */
331
332 #ifdef ALLOW_NONHYDROSTATIC
333 COMMON /GRID_OBNH/
334 & OBNw, OBSw, OBEw, OBWw
335 _RL OBNw (1-Olx:sNx+Olx,Nr,nSx,nSy)
336 _RL OBSw (1-Olx:sNx+Olx,Nr,nSx,nSy)
337 _RL OBEw (1-Oly:sNy+Oly,Nr,nSx,nSy)
338 _RL OBWw (1-Oly:sNy+Oly,Nr,nSx,nSy)
339 #ifdef ALLOW_OBCS_PRESCRIBE
340 COMMON /GRID_OBNH_AUX/
341 & OBNw0, OBSw0, OBEw0, OBWw0,
342 & OBNw1, OBSw1, OBEw1, OBWw1
343 _RL OBNw0(1-Olx:sNx+Olx,Nr,nSx,nSy)
344 _RL OBSw0(1-Olx:sNx+Olx,Nr,nSx,nSy)
345 _RL OBEw0(1-Oly:sNy+Oly,Nr,nSx,nSy)
346 _RL OBWw0(1-Oly:sNy+Oly,Nr,nSx,nSy)
347 _RL OBNw1(1-Olx:sNx+Olx,Nr,nSx,nSy)
348 _RL OBSw1(1-Olx:sNx+Olx,Nr,nSx,nSy)
349 _RL OBEw1(1-Oly:sNy+Oly,Nr,nSx,nSy)
350 _RL OBWw1(1-Oly:sNy+Oly,Nr,nSx,nSy)
351 #endif /* ALLOW_OBCS_PRESCRIBE */
352 #endif /* ALLOW_NONHYDROSTATIC */
353
354 #ifdef NONLIN_FRSURF
355 COMMON /OB_NLFS/
356 & OBNeta, OBSeta, OBEeta, OBWeta
357 _RL OBNeta (1-Olx:sNx+Olx,nSx,nSy)
358 _RL OBSeta (1-Olx:sNx+Olx,nSx,nSy)
359 _RL OBEeta (1-Oly:sNy+Oly,nSx,nSy)
360 _RL OBWeta (1-Oly:sNy+Oly,nSx,nSy)
361 #ifdef ALLOW_OBCS_PRESCRIBE
362 COMMON /OB_NLFS_AUX/
363 & OBNeta0,OBSeta0,OBEeta0,OBWeta0,
364 & OBNeta1,OBSeta1,OBEeta1,OBWeta1
365 _RL OBNeta0(1-Olx:sNx+Olx,nSx,nSy)
366 _RL OBSeta0(1-Olx:sNx+Olx,nSx,nSy)
367 _RL OBEeta0(1-Oly:sNy+Oly,nSx,nSy)
368 _RL OBWeta0(1-Oly:sNy+Oly,nSx,nSy)
369 _RL OBNeta1(1-Olx:sNx+Olx,nSx,nSy)
370 _RL OBSeta1(1-Olx:sNx+Olx,nSx,nSy)
371 _RL OBEeta1(1-Oly:sNy+Oly,nSx,nSy)
372 _RL OBWeta1(1-Oly:sNy+Oly,nSx,nSy)
373 #endif /* ALLOW_OBCS_PRESCRIBE */
374 #endif /* NONLIN_FRSURF */
375 CMM(
376 #ifdef ALLOW_OBCS_CONTROL_MODES
377 common /ih_modes/ modesv
378 _RL modesv (nr,nr,nr)
379 #endif
380 CMM)
381
382 #endif /* ALLOW_OBCS */

  ViewVC Help
Powered by ViewVC 1.1.22