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

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

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


Revision 1.2 - (hide annotations) (download)
Fri Mar 18 18:47:49 2011 UTC (14 years, 9 months ago) by mmazloff
Branch: MAIN
Changes since 1.1: +21 -16 lines
File MIME type: text/plain
bringing ALLOW_OBCS_CONTROL_MODES up to date

1 mmazloff 1.2 C $Header: /u/gcmpack/MITgcm/pkg/obcs/OBCS.h,v 1.28 2011/03/14 17:31:07 mlosch Exp $
2 mmazloff 1.1 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 mmazloff 1.2 C useStevensPhaseVel
14 mmazloff 1.1 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 mmazloff 1.2 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 mmazloff 1.1 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 mmazloff 1.2 & OBCS_monSelect,
46 mmazloff 1.1 & spongeThickness,
47     & useOrlanskiNorth,useOrlanskiSouth,
48     & useOrlanskiEast,useOrlanskiWest,
49     & useStevensNorth,useStevensSouth,
50 mmazloff 1.2 & useStevensEast,useStevensWest,
51 mmazloff 1.1 & useStevensPhaseVel, useStevensAdvection,
52     & useOBCSsponge, useOBCSbalance, useOBCSprescribe,
53     & OBCSprintDiags, useOBCSYearlyFields,
54     & OBCSfixTopo
55 mmazloff 1.2 INTEGER OBCS_monSelect
56 mmazloff 1.1 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 mmazloff 1.2 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 mmazloff 1.1 COMMON /PARM_R_OB/
78 mmazloff 1.2 & OBCS_balanceFacN, OBCS_balanceFacS,
79     & OBCS_balanceFacE, OBCS_balanceFacW,
80     & OBCS_monitorFreq,
81 mmazloff 1.1 & Urelaxobcsinner,Urelaxobcsbound,
82     & Vrelaxobcsinner,Vrelaxobcsbound,
83     & TrelaxStevens, SrelaxStevens
84 mmazloff 1.2 _RL OBCS_balanceFacN, OBCS_balanceFacS
85     _RL OBCS_balanceFacE, OBCS_balanceFacW
86     _RL OBCS_monitorFreq
87 mmazloff 1.1 _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 mmazloff 1.2 & insideOBmaskFile,
110 mmazloff 1.1 & 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 mmazloff 1.2 & insideOBmaskFile,
122 mmazloff 1.1 & 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 mmazloff 1.2 C-- COMMON /GRID_OB/ Open boundary related stuff
134 mmazloff 1.1 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_NORTH
147     COMMON /GRID_N_OB/
148     & OBNu,OBNv,OBNt,OBNs
149     _RL OBNu (1-Olx:sNx+Olx,Nr,nSx,nSy)
150     _RL OBNv (1-Olx:sNx+Olx,Nr,nSx,nSy)
151     _RL OBNt (1-Olx:sNx+Olx,Nr,nSx,nSy)
152     _RL OBNs (1-Olx:sNx+Olx,Nr,nSx,nSy)
153     #ifdef ALLOW_OBCS_PRESCRIBE
154     COMMON /GRID_N_OB_AUX/
155     & OBNu0,OBNv0,OBNt0,OBNs0,
156     & OBNu1,OBNv1,OBNt1,OBNs1
157     _RL OBNu0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
158     _RL OBNv0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
159     _RL OBNt0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
160     _RL OBNs0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
161     _RL OBNu1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
162     _RL OBNv1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
163     _RL OBNt1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
164     _RL OBNs1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
165     #endif /* ALLOW_OBCS_PRESCRIBE */
166     #ifdef ALLOW_SEAICE
167     COMMON /SEAICE_N_OB/ OBNa,OBNh,OBNa0,OBNh0,OBNa1,OBNh1,
168     & OBNsl,OBNsn,OBNsl0,OBNsn0,OBNsl1,OBNsn1,
169     & OBNuice,OBNvice,OBNuice0,OBNvice0,OBNuice1,OBNvice1
170     _RL OBNa (1-Olx:sNx+Olx,nSx,nSy)
171     _RL OBNh (1-Olx:sNx+Olx,nSx,nSy)
172     _RL OBNa0 (1-Olx:sNx+Olx,nSx,nSy)
173     _RL OBNh0 (1-Olx:sNx+Olx,nSx,nSy)
174     _RL OBNa1 (1-Olx:sNx+Olx,nSx,nSy)
175     _RL OBNh1 (1-Olx:sNx+Olx,nSx,nSy)
176     _RL OBNsl (1-Olx:sNx+Olx,nSx,nSy)
177     _RL OBNsn (1-Olx:sNx+Olx,nSx,nSy)
178     _RL OBNsl0 (1-Olx:sNx+Olx,nSx,nSy)
179     _RL OBNsn0 (1-Olx:sNx+Olx,nSx,nSy)
180     _RL OBNsl1 (1-Olx:sNx+Olx,nSx,nSy)
181     _RL OBNsn1 (1-Olx:sNx+Olx,nSx,nSy)
182     _RL OBNuice (1-Olx:sNx+Olx,nSx,nSy)
183     _RL OBNvice (1-Olx:sNx+Olx,nSx,nSy)
184     _RL OBNuice0 (1-Olx:sNx+Olx,nSx,nSy)
185     _RL OBNvice0 (1-Olx:sNx+Olx,nSx,nSy)
186     _RL OBNuice1 (1-Olx:sNx+Olx,nSx,nSy)
187     _RL OBNvice1 (1-Olx:sNx+Olx,nSx,nSy)
188     #endif /* ALLOW_SEAICE */
189     #endif /* ALLOW_OBCS_NORTH */
190    
191     #ifdef ALLOW_OBCS_SOUTH
192     COMMON /GRID_S_OB/
193     & OBSu,OBSv,OBSt,OBSs
194     _RL OBSu (1-Olx:sNx+Olx,Nr,nSx,nSy)
195     _RL OBSv (1-Olx:sNx+Olx,Nr,nSx,nSy)
196     _RL OBSt (1-Olx:sNx+Olx,Nr,nSx,nSy)
197     _RL OBSs (1-Olx:sNx+Olx,Nr,nSx,nSy)
198     #ifdef ALLOW_OBCS_PRESCRIBE
199     COMMON /GRID_S_OB_AUX/
200     & OBSu0,OBSv0,OBSt0,OBSs0,
201     & OBSu1,OBSv1,OBSt1,OBSs1
202     _RL OBSu0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
203     _RL OBSv0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
204     _RL OBSt0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
205     _RL OBSs0 (1-Olx:sNx+Olx,Nr,nSx,nSy)
206     _RL OBSu1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
207     _RL OBSv1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
208     _RL OBSt1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
209     _RL OBSs1 (1-Olx:sNx+Olx,Nr,nSx,nSy)
210     #endif /* ALLOW_OBCS_PRESCRIBE */
211     #ifdef ALLOW_SEAICE
212     COMMON /SEAICE_S_OB/ OBSa,OBSh,OBSa0,OBSh0,OBSa1,OBSh1,
213     & OBSsl,OBSsn,OBSsl0,OBSsn0,OBSsl1,OBSsn1,
214     & OBSuice,OBSvice,OBSuice0,OBSvice0,OBSuice1,OBSvice1
215     _RL OBSa (1-Olx:sNx+Olx,nSx,nSy)
216     _RL OBSh (1-Olx:sNx+Olx,nSx,nSy)
217     _RL OBSa0 (1-Olx:sNx+Olx,nSx,nSy)
218     _RL OBSh0 (1-Olx:sNx+Olx,nSx,nSy)
219     _RL OBSa1 (1-Olx:sNx+Olx,nSx,nSy)
220     _RL OBSh1 (1-Olx:sNx+Olx,nSx,nSy)
221     _RL OBSsl (1-Olx:sNx+Olx,nSx,nSy)
222     _RL OBSsn (1-Olx:sNx+Olx,nSx,nSy)
223     _RL OBSsl0 (1-Olx:sNx+Olx,nSx,nSy)
224     _RL OBSsn0 (1-Olx:sNx+Olx,nSx,nSy)
225     _RL OBSsl1 (1-Olx:sNx+Olx,nSx,nSy)
226     _RL OBSsn1 (1-Olx:sNx+Olx,nSx,nSy)
227     _RL OBSuice (1-Olx:sNx+Olx,nSx,nSy)
228     _RL OBSvice (1-Olx:sNx+Olx,nSx,nSy)
229     _RL OBSuice0 (1-Olx:sNx+Olx,nSx,nSy)
230     _RL OBSvice0 (1-Olx:sNx+Olx,nSx,nSy)
231     _RL OBSuice1 (1-Olx:sNx+Olx,nSx,nSy)
232     _RL OBSvice1 (1-Olx:sNx+Olx,nSx,nSy)
233     #endif /* ALLOW_SEAICE */
234     #endif /* ALLOW_OBCS_SOUTH */
235    
236     #ifdef ALLOW_OBCS_EAST
237     COMMON /GRID_E_OB/
238     & OBEu,OBEv,OBEt,OBEs
239     _RL OBEu (1-Oly:sNy+Oly,Nr,nSx,nSy)
240     _RL OBEv (1-Oly:sNy+Oly,Nr,nSx,nSy)
241     _RL OBEt (1-Oly:sNy+Oly,Nr,nSx,nSy)
242     _RL OBEs (1-Oly:sNy+Oly,Nr,nSx,nSy)
243     #ifdef ALLOW_OBCS_PRESCRIBE
244     COMMON /GRID_E_OB_AUX/
245     & OBEu0,OBEv0,OBEt0,OBEs0,
246     & OBEu1,OBEv1,OBEt1,OBEs1
247     _RL OBEu0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
248     _RL OBEv0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
249     _RL OBEt0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
250     _RL OBEs0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
251     _RL OBEu1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
252     _RL OBEv1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
253     _RL OBEt1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
254     _RL OBEs1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
255     #endif /* ALLOW_OBCS_PRESCRIBE */
256     #ifdef ALLOW_SEAICE
257     COMMON /SEAICE_E_OB/ OBEa,OBEh,OBEa0,OBEh0,OBEa1,OBEh1,
258     & OBEsl,OBEsn,OBEsl0,OBEsn0,OBEsl1,OBEsn1,
259     & OBEuice,OBEvice,OBEuice0,OBEvice0,OBEuice1,OBEvice1
260     _RL OBEa (1-Oly:sNy+Oly,nSx,nSy)
261     _RL OBEh (1-Oly:sNy+Oly,nSx,nSy)
262     _RL OBEa0 (1-Oly:sNy+Oly,nSx,nSy)
263     _RL OBEh0 (1-Oly:sNy+Oly,nSx,nSy)
264     _RL OBEa1 (1-Oly:sNy+Oly,nSx,nSy)
265     _RL OBEh1 (1-Oly:sNy+Oly,nSx,nSy)
266     _RL OBEsl (1-Oly:sNy+Oly,nSx,nSy)
267     _RL OBEsn (1-Oly:sNy+Oly,nSx,nSy)
268     _RL OBEsl0 (1-Oly:sNy+Oly,nSx,nSy)
269     _RL OBEsn0 (1-Oly:sNy+Oly,nSx,nSy)
270     _RL OBEsl1 (1-Oly:sNy+Oly,nSx,nSy)
271     _RL OBEsn1 (1-Oly:sNy+Oly,nSx,nSy)
272     _RL OBEuice (1-Oly:sNy+Oly,nSx,nSy)
273     _RL OBEvice (1-Oly:sNy+Oly,nSx,nSy)
274     _RL OBEuice0 (1-Oly:sNy+Oly,nSx,nSy)
275     _RL OBEvice0 (1-Oly:sNy+Oly,nSx,nSy)
276     _RL OBEuice1 (1-Oly:sNy+Oly,nSx,nSy)
277     _RL OBEvice1 (1-Oly:sNy+Oly,nSx,nSy)
278     #endif /* ALLOW_SEAICE */
279     #endif /* ALLOW_OBCS_EAST */
280    
281     #ifdef ALLOW_OBCS_WEST
282     COMMON /GRID_W_OB/
283     & OBWu,OBWv,OBWt,OBWs
284     _RL OBWu (1-Oly:sNy+Oly,Nr,nSx,nSy)
285     _RL OBWv (1-Oly:sNy+Oly,Nr,nSx,nSy)
286     _RL OBWt (1-Oly:sNy+Oly,Nr,nSx,nSy)
287     _RL OBWs (1-Oly:sNy+Oly,Nr,nSx,nSy)
288     #ifdef ALLOW_OBCS_PRESCRIBE
289     COMMON /GRID_W_OB_AUX/
290     & OBWu0,OBWv0,OBWt0,OBWs0,
291     & OBWu1,OBWv1,OBWt1,OBWs1
292     _RL OBWu0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
293     _RL OBWv0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
294     _RL OBWt0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
295     _RL OBWs0 (1-Oly:sNy+Oly,Nr,nSx,nSy)
296     _RL OBWu1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
297     _RL OBWv1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
298     _RL OBWt1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
299     _RL OBWs1 (1-Oly:sNy+Oly,Nr,nSx,nSy)
300     #endif /* ALLOW_OBCS_PRESCRIBE */
301     #ifdef ALLOW_SEAICE
302     COMMON /SEAICE_W_OB/ OBWa,OBWh,OBWa0,OBWh0,OBWa1,OBWh1,
303     & OBWsl,OBWsn,OBWsl0,OBWsn0,OBWsl1,OBWsn1,
304     & OBWuice,OBWvice,OBWuice0,OBWvice0,OBWuice1,OBWvice1
305     _RL OBWa (1-Oly:sNy+Oly,nSx,nSy)
306     _RL OBWh (1-Oly:sNy+Oly,nSx,nSy)
307     _RL OBWa0 (1-Oly:sNy+Oly,nSx,nSy)
308     _RL OBWh0 (1-Oly:sNy+Oly,nSx,nSy)
309     _RL OBWa1 (1-Oly:sNy+Oly,nSx,nSy)
310     _RL OBWh1 (1-Oly:sNy+Oly,nSx,nSy)
311     _RL OBWsl (1-Oly:sNy+Oly,nSx,nSy)
312     _RL OBWsn (1-Oly:sNy+Oly,nSx,nSy)
313     _RL OBWsl0 (1-Oly:sNy+Oly,nSx,nSy)
314     _RL OBWsn0 (1-Oly:sNy+Oly,nSx,nSy)
315     _RL OBWsl1 (1-Oly:sNy+Oly,nSx,nSy)
316     _RL OBWsn1 (1-Oly:sNy+Oly,nSx,nSy)
317     _RL OBWuice (1-Oly:sNy+Oly,nSx,nSy)
318     _RL OBWvice (1-Oly:sNy+Oly,nSx,nSy)
319     _RL OBWuice0 (1-Oly:sNy+Oly,nSx,nSy)
320     _RL OBWvice0 (1-Oly:sNy+Oly,nSx,nSy)
321     _RL OBWuice1 (1-Oly:sNy+Oly,nSx,nSy)
322     _RL OBWvice1 (1-Oly:sNy+Oly,nSx,nSy)
323     #endif /* ALLOW_SEAICE */
324     #endif /* ALLOW_OBCS_WEST */
325    
326     #ifdef ALLOW_NONHYDROSTATIC
327     COMMON /GRID_OBNH/
328     & OBNw, OBSw, OBEw, OBWw
329     _RL OBNw (1-Olx:sNx+Olx,Nr,nSx,nSy)
330     _RL OBSw (1-Olx:sNx+Olx,Nr,nSx,nSy)
331     _RL OBEw (1-Oly:sNy+Oly,Nr,nSx,nSy)
332     _RL OBWw (1-Oly:sNy+Oly,Nr,nSx,nSy)
333     #ifdef ALLOW_OBCS_PRESCRIBE
334     COMMON /GRID_OBNH_AUX/
335     & OBNw0, OBSw0, OBEw0, OBWw0,
336     & OBNw1, OBSw1, OBEw1, OBWw1
337     _RL OBNw0(1-Olx:sNx+Olx,Nr,nSx,nSy)
338     _RL OBSw0(1-Olx:sNx+Olx,Nr,nSx,nSy)
339     _RL OBEw0(1-Oly:sNy+Oly,Nr,nSx,nSy)
340     _RL OBWw0(1-Oly:sNy+Oly,Nr,nSx,nSy)
341     _RL OBNw1(1-Olx:sNx+Olx,Nr,nSx,nSy)
342     _RL OBSw1(1-Olx:sNx+Olx,Nr,nSx,nSy)
343     _RL OBEw1(1-Oly:sNy+Oly,Nr,nSx,nSy)
344     _RL OBWw1(1-Oly:sNy+Oly,Nr,nSx,nSy)
345     #endif /* ALLOW_OBCS_PRESCRIBE */
346     #endif /* ALLOW_NONHYDROSTATIC */
347    
348     #ifdef NONLIN_FRSURF
349     COMMON /OB_NLFS/
350     & OBNeta, OBSeta, OBEeta, OBWeta
351     _RL OBNeta (1-Olx:sNx+Olx,nSx,nSy)
352     _RL OBSeta (1-Olx:sNx+Olx,nSx,nSy)
353     _RL OBEeta (1-Oly:sNy+Oly,nSx,nSy)
354     _RL OBWeta (1-Oly:sNy+Oly,nSx,nSy)
355     #ifdef ALLOW_OBCS_PRESCRIBE
356     COMMON /OB_NLFS_AUX/
357     & OBNeta0,OBSeta0,OBEeta0,OBWeta0,
358     & OBNeta1,OBSeta1,OBEeta1,OBWeta1
359     _RL OBNeta0(1-Olx:sNx+Olx,nSx,nSy)
360     _RL OBSeta0(1-Olx:sNx+Olx,nSx,nSy)
361     _RL OBEeta0(1-Oly:sNy+Oly,nSx,nSy)
362     _RL OBWeta0(1-Oly:sNy+Oly,nSx,nSy)
363     _RL OBNeta1(1-Olx:sNx+Olx,nSx,nSy)
364     _RL OBSeta1(1-Olx:sNx+Olx,nSx,nSy)
365     _RL OBEeta1(1-Oly:sNy+Oly,nSx,nSy)
366     _RL OBWeta1(1-Oly:sNy+Oly,nSx,nSy)
367     #endif /* ALLOW_OBCS_PRESCRIBE */
368     #endif /* NONLIN_FRSURF */
369     CMM(
370 mmazloff 1.2 #ifdef ALLOW_OBCS_CONTROL_MODES
371 mmazloff 1.1 common /ih_modes/ modesv
372     _RL modesv (nr,nr,nr)
373     #endif
374     CMM)
375    
376     #endif /* ALLOW_OBCS */

  ViewVC Help
Powered by ViewVC 1.1.22