/[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.3 - (hide annotations) (download)
Mon Apr 18 23:47:25 2011 UTC (14 years, 8 months ago) by mmazloff
Branch: MAIN
Changes since 1.2: +7 -1 lines
File MIME type: text/plain
Updating

1 mmazloff 1.3 C $Header: /u/gcmpack/MITgcm/pkg/obcs/OBCS.h,v 1.29 2011/04/17 21:12:24 jmc 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 mmazloff 1.3 #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 mmazloff 1.1 #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 mmazloff 1.2 #ifdef ALLOW_OBCS_CONTROL_MODES
377 mmazloff 1.1 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