1 | include(dom.inc)
2 | ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
3 | !
4 | ! module: avbp_exBound
5 | ! read *.exBound files
6 | !
7 | ! author: Andre Kaufmann, 03.08.2000
8 | ! status: untested
9 | ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
10 |
11 | module avbp_exBound
12 |
13 | implicit none
14 |
15 | character*80 exBoundfile
16 | DOM_INT :: exBound_npart = 1
17 | DOM_INT :: exBound_mnbpatch = 0
18 | DOM_INT :: exBound_ntbound = 0
19 | DOM_INT :: exBound_ntbface = 0
20 | DOM_INT :: exBound_npbound = 0
21 | DOM_INT :: exBound_nbbegin = 0
22 | DOM_INT :: exBound_nblength = 0
23 | DOM_INT :: exBound_nbfbegin = 0
24 | DOM_INT :: exBound_nbflength = 0
25 | DOM_INT :: exBound_nbaxipoints = 0
26 | DOM_INT :: exBound_nbbndyconn = 0
27 | DOM_INT :: exBound_ndim = 0
28 | logical :: exBound_readaxiinfo = .false.
29 | logical :: exBound_readconninfo = .false.
30 |
31 | logical :: exBound_debug = .false.
32 | logical :: exBound_itexists = .false.
33 | logical :: exBound_show = .true.
34 |
35 | DOM_INT,dimension(:),allocatable :: exBound_ibound
36 | DOM_INT,dimension(:),allocatable :: exBound_ibound1
37 | DOM_INT,dimension(:),allocatable :: exBound_ibound2
38 | DOM_INT,dimension(:),allocatable :: exBound_ibfael
39 | DOM_INT,dimension(:),allocatable :: exBound_axipoints
40 | DOM_INT,dimension(:),allocatable :: exBound_bndyconn
41 | DOM_INT,dimension(:,:),allocatable :: exBound_ifbound
42 |
43 | DOM_REAL,dimension(:),allocatable :: exBound_snbound
44 |
45 | contains
46 | ! ================================================================
47 | !
48 | ! deallocate memory
49 | !
50 | ! ================================================================
51 | subroutine destroyexBoundmemory
52 | implicit none
53 | if(exBound_debug) then
54 | write(*,*) ' deallocating exBound memory ... '
55 | endif
56 | if(allocated(exBound_ibound)) deallocate(exBound_ibound)
57 | if(allocated(exBound_ibound1)) deallocate(exBound_ibound1)
58 | if(allocated(exBound_ibound2)) deallocate(exBound_ibound2)
59 | if(allocated(exBound_snbound)) deallocate(exBound_snbound)
60 | if(allocated(exBound_ibfael)) deallocate(exBound_ibfael)
61 | if(allocated(exBound_axipoints)) deallocate(exBound_axipoints)
62 | if(allocated(exBound_bndyconn)) deallocate(exBound_bndyconn)
63 | if(allocated(exBound_ifbound)) deallocate(exBound_ifbound)
64 |
65 | return
66 | end subroutine destroyexBoundmemory
67 | ! ================================================================
68 | !
69 | ! get the boundary info
70 | !
71 | ! ================================================================
72 | subroutine getexBoundinfo
73 | implicit none
74 | DOM_INT n,m,n1,n2,n3,n4,n5,n6,k
75 | DOM_INT read_iostatus
76 |
77 | if(exBound_debug) then
78 | write(*,*) 'getexBoundinfo ... '
79 | endif
80 |
81 | inquire(file=exBoundfile,exist=exBound_itexists)
82 |
83 | if(exBound_itexists) then
84 |
85 | n = len_trim(exBoundfile)
86 |
87 | open(unit=2,file=exBoundfile,form='unformatted')
88 |
89 | read(2,IOSTAT=read_iostatus) &
90 | & exBound_npart,exBound_mnbpatch, &
91 | & exBound_ntbound,exBound_ntbface
92 |
93 | if(exBound_debug.or.(read_iostatus.ne.0)) then
94 | write(*,*) '==================================='
95 | write(*,*) ' npart = ',exBound_npart
96 | write(*,*) ' mnbpatch = ',exBound_mnbpatch
97 | write(*,*) ' ntbound = ',exBound_ntbound
98 | write(*,*) ' ntbface = ',exBound_ntbface
99 | write(*,*) ' IOSTATUS = ',read_iostatus
100 | endif
101 |
102 | read(2,IOSTAT=read_iostatus) &
103 | & exBound_npbound,exBound_nbbegin, &
104 | & exBound_nblength,exBound_nbfbegin, &
105 | & exBound_nbflength
106 |
107 | if(exBound_debug.or.(read_iostatus.ne.0)) then
108 | write(*,*) '==================================='
109 | write(*,*) ' npbound = ',exBound_npbound
110 | write(*,*) ' nbbegin = ',exBound_nbbegin
111 | write(*,*) ' nblength = ',exBound_nblength
112 | write(*,*) ' nbfbegin = ',exBound_nbfbegin
113 | write(*,*) ' nbflength = ',exBound_nbflength
114 | write(*,*) ' IOSTATUS = ',read_iostatus
115 | endif
116 |
117 | else
118 | write(*,*) ' could not read exBoundfile ',exBoundfile
119 | endif
120 |
121 | return
122 | end subroutine getexBoundinfo
123 | ! ================================================================
124 | !
125 | ! read exBound file
126 | !
127 | ! ================================================================
128 | subroutine readexBound
129 | implicit none
130 | !
131 | ! local varialbes
132 | !
133 | DOM_INT :: n,m,n1,n2,n3,n4,n5,n6,n7,k
134 | DOM_INT :: read_iostatus,allocstat
135 |
136 | if(exBound_debug) then
137 | write(*,*) 'readexBound ... '
138 | endif
139 |
140 |
141 | inquire(file=exBoundfile,exist=exBound_itexists)
142 |
143 | if(exBound_itexists) then
144 |
145 | n = len_trim(exBoundfile)
146 | if (exBound_show) then
147 | write(*,*) ' >>>> Reading the exBound file ',exBoundfile(1:n)
148 | end if
149 |
150 | open(unit=2,file=exBoundfile,form='unformatted')
151 | read(2,IOSTAT=read_iostatus) &
152 | & exBound_npart,exBound_mnbpatch, &
153 | & exBound_ntbound,exBound_ntbface
154 |
155 | if(exBound_debug.or.(read_iostatus.ne.0)) then
156 | write(*,*) '==================================='
157 | write(*,*) ' npart = ',exBound_npart
158 | write(*,*) ' mnbpatch = ',exBound_mnbpatch
159 | write(*,*) ' ntbound = ',exBound_ntbound
160 | write(*,*) ' ntbface = ',exBound_ntbface
161 | write(*,*) ' IOSTATUS = ',read_iostatus
162 | endif
163 | read(2,IOSTAT=read_iostatus) &
164 | & exBound_npbound,exBound_nbbegin, &
165 | & exBound_nblength,exBound_nbfbegin, &
166 | & exBound_nbflength
167 |
168 |
169 | if(exBound_debug.or.(read_iostatus.ne.0)) then
170 | write(*,*) '==================================='
171 | write(*,*) ' npbound = ',exBound_npbound
172 | write(*,*) ' nbbegin = ',exBound_nbbegin
173 | write(*,*) ' nblength = ',exBound_nblength
174 | write(*,*) ' nbfbegin = ',exBound_nbfbegin
175 | write(*,*) ' nbflength = ',exBound_nbflength
176 | write(*,*) ' ndim = ',exBound_ndim
177 | write(*,*) ' IOSTATUS = ',read_iostatus
178 | endif
179 |
180 | !
181 | ! allocate memory for this data
182 | !
183 | if(allocated(exBound_ibound)) deallocate(exBound_ibound)
184 | if(allocated(exBound_ibound1)) deallocate(exBound_ibound1)
185 | if(allocated(exBound_ibound2)) deallocate(exBound_ibound2)
186 | if(allocated(exBound_snbound)) deallocate(exBound_snbound)
187 | if(allocated(exBound_ibfael)) deallocate(exBound_ibfael)
188 | if(allocated(exBound_axipoints)) deallocate(exBound_axipoints)
189 | if(allocated(exBound_bndyconn)) deallocate(exBound_bndyconn)
190 | if(allocated(exBound_ifbound)) deallocate(exBound_ifbound)
191 |
192 | allocate(exBound_ibound(1:exBound_nblength))
193 | allocate(exBound_ibound1(1:exBound_npbound))
194 | allocate(exBound_ibound2(1:exBound_npbound))
195 | allocate(exBound_snbound(1:exBound_nblength*exBound_ndim))
196 | allocate(exBound_ifbound(1:4,1:exBound_npbound))
197 | allocate(exBound_ibfael(1:exBound_nbflength*2))
198 |
199 | !
200 | ! calculate the running index
201 | !
202 | n1 = exBound_nbbegin
203 | n2 = exBound_nbbegin+exBound_nblength - 1
204 | n3 = 1+(exBound_nbbegin-1)*exBound_ndim
205 | n4 = (exBound_nbbegin+exBound_nblength-1)*exBound_ndim
206 | n5 = 1+(exBound_nbfbegin -1)*2
207 | n6 = (exBound_nbfbegin+exBound_nbflength-1)*2
208 | !
209 | ! read start adress and length for each patch (nodes)
210 | !
211 | read(2,IOSTAT=read_iostatus) &
212 | & (exBound_ibound1(n), &
213 | & exBound_ibound2(n),n=1,exBound_npbound)
214 |
215 | if(exBound_debug.or.(read_iostatus.ne.0)) then
216 | write(*,*) '==================================='
217 | write(*,*) 'exBound_ibound1, ibound2 read '
218 | write(*,*) ' IOSTATUS = ',read_iostatus
219 | endif
220 |
221 | !
222 | ! pointer array from boundary nodes to global nodes
223 | !
224 | read(2,IOSTAT=read_iostatus) &
225 | & (exBound_ibound(n),n=n1,n2)
226 |
227 | if(exBound_debug.or.(read_iostatus.ne.0)) then
228 | write(*,*) '==================================='
229 | write(*,*) 'exBound_ibound read '
230 | write(*,*) ' IOSTATUS = ',read_iostatus
231 | endif
232 |
233 | !
234 | ! boundary normal vectors
235 | !
236 | read(2,IOSTAT=read_iostatus) (exBound_snbound(n),n=n3,n4)
237 |
238 | if(exBound_debug.or.(read_iostatus.ne.0)) then
239 | write(*,*) '==================================='
240 | write(*,*) 'exBound_snbound read '
241 | write(*,*) ' IOSTATUS = ',read_iostatus
242 | endif
243 | !
244 | ! for each patch start adress and length of boundary faces
245 | ! (quad and tri see hip doc)
246 | !
247 | read(2,IOSTAT=read_iostatus) &
248 | & ((exBound_ifbound(k,n),k=1,4),n=1,exBound_npbound)
249 |
250 | if(exBound_debug.or.(read_iostatus.ne.0)) then
251 | write(*,*) '==================================='
252 | write(*,*) 'exBound_snbound read '
253 | write(*,*) ' IOSTATUS = ',read_iostatus
254 | endif
255 | !
256 | ! boundary face to element pointer
257 | !
258 | read(2,IOSTAT=read_iostatus) (exBound_ibfael(n),n=n5,n6)
259 |
260 | if(exBound_debug.or.(read_iostatus.ne.0)) then
261 | write(*,*) '==================================='
262 | write(*,*) 'exBound_ibfael read '
263 | write(*,*) ' IOSTATUS = ',read_iostatus
264 | endif
265 |
266 | !
267 | ! read the axi information if present
268 | ! a) read the record length, and check if existance
269 | ! b) allocate mem
270 | ! c) read the record
271 | !
272 |
273 | if(exBound_readaxiinfo) then
274 |
275 | read(2,IOSTAT=read_iostatus) exBound_nbaxipoints
276 |
277 | if(exBound_debug.or.(read_iostatus.ne.0)) then
278 | write(*,*) '==================================='
279 | write(*,*) 'nbaxipoints = ',exbound_nbaxipoints
280 | write(*,*) ' IOSTATUS = ',read_iostatus
281 | endif
282 |
283 | if(read_iostatus.lt.0) then
284 | write(*,*) '==================================='
285 | write(*,*) 'nbaxipoints = ',exbound_nbaxipoints
286 | write(*,*) ' IOSTATUS = ',read_iostatus
287 | write(*,*) ' no axisymmetric information present', &
288 | & ' in this exBound; check exBound and ', &
289 | & ' run new hip ! '
290 | endif
291 |
292 | !
293 | ! does it exist ?
294 | !
295 | if(read_iostatus.ge.0) then
296 | !
297 | ! axi symmetric grid ?
298 | !
299 | if(exBound_nbaxipoints.gt.0) then
300 | !
301 | ! debug info
302 | !
303 | if(exBound_debug) then
304 | write(*,*) '==================================='
305 | write(*,*) ' reading axi information ... '
306 | endif
307 | !
308 | ! b) allocate mem
309 | !
310 | allocate(exBound_axipoints(1:exBound_nbaxipoints))
311 | !
312 | ! c) read the record
313 | !
314 | read(2,IOSTAT=read_iostatus) (exBound_axipoints(n), &
315 | & n=1,exBound_nbaxipoints)
316 |
317 | if(read_iostatus.ne.0) then
318 | write(*,*) 'avbp_exBound : could not read axisymmetric', &
319 | & 'boundary nodes; check exBound file '
320 | write(*,*) '(IOSTATUS = ',read_iostatus,' )'
321 | endif
322 |
323 | else
324 | !
325 | ! debug info
326 | !
327 | if(exBound_debug) then
328 | write(*,*) '==================================='
329 | write(*,*) ' no axi information empty read statement ... '
330 | endif
331 |
332 | !
333 | ! no axi points ==> empty read statement
334 | !
335 | read(2,IOSTAT=read_iostatus)
336 | !
337 | !
338 | endif
339 | !
340 | !
341 | endif
342 |
343 | if(exBound_readaxiinfo.and.exBound_readconninfo) then
344 | !
345 | ! read the boundary connectivity (same index as ibfael ...)
346 | ! a) calculate the number of entries
347 | ! b) allocate mem
348 | ! c) read the record and check for endof file
349 | !
350 | !
351 | !
352 | ! a) calculate the number of entries and b) allocate mem
353 | ! number of tri + number of quads per patch
354 | !
355 | n7 = 0
356 | if(exBound_ndim.eq.2) then
357 | do n=1,exBound_npbound
358 | n7 = n7 + exBound_ifbound(4,n) * 2
359 | enddo
360 | else if(exBound_ndim.eq.3) then
361 | do n=1,exBound_npbound
362 | n7 = n7 + exBound_ifbound(2,n) * 3 + &
363 | & exBound_ifbound(4,n) * 4
364 | enddo
365 | endif
366 | !
367 | ! give a warning if boundary connectivity is zero
368 | !
369 | if(exBound_debug) then
370 | write(*,*) '==================================='
371 | write(*,*) 'boundary connectivity : ',n7
372 | write(*,*) 'allocating boundary info ... '
373 | endif
374 |
375 | exBound_nbbndyconn = n7
376 |
377 | allocate(exBound_bndyconn(1:exBound_nbbndyconn),stat=allocstat)
378 | !
379 | ! debug info
380 | !
381 | if(exBound_debug) then
382 | write(*,*) '==================================='
383 | write(*,*) 'allocation status: ',allocstat
384 | endif
385 | !
386 | !
387 | !
388 |
389 | !
390 | ! c) read the record
391 | !
392 | read(2,IOSTAT=read_iostatus) (exBound_bndyconn(n), &
393 | & n=1,exBound_nbbndyconn)
394 |
395 | if(read_iostatus.ne.0) then
396 | write(*,*) 'avbp_exBound : could not read boundary ', &
397 | & 'boundary connectivity; check exBound file '
398 | endif
399 | if((read_iostatus.ne.0).or.exBound_debug) then
400 | write(*,*) '(IOSTATUS = ',read_iostatus,' )'
401 | endif
402 | !
403 | !
404 | endif
405 | !
406 | !
407 | endif
408 | !
409 | !
410 | close(2)
411 | !
412 | !
413 | else
414 | n = len_trim(exBoundfile)
415 | write(*,*) ' could not read the exBound file ', &
416 | & exBoundfile(1:n)
417 | endif
418 |
419 | return
420 | end subroutine readexBound
421 | ! ================================================================
422 | !
423 | ! write the ex Bound variables
424 | !
425 | ! ================================================================
426 | subroutine writeexBoundvariables
427 | DOM_INT n,m,n1,n2,n3,n4,n5,n6,k
428 |
429 | if(exBound_debug) then
430 | write(*,*) 'write exBoundvariables ... '
431 | endif
432 |
433 |
434 | write(*,*) exBound_npart,exBound_mnbpatch, &
435 | & exBound_ntbound,exBound_ntbface
436 |
437 | write(*,*) exBound_npbound,exBound_nbbegin, &
438 | & exBound_nblength,exBound_nbfbegin, &
439 | & exBound_nbflength
440 |
441 | n1 = exBound_nbbegin
442 | n2 = exBound_nbbegin+exBound_nblength-1
443 | n3 = 1+(exBound_nbbegin-1)*exBound_ndim
444 | n4 = (exBound_nbbegin+exBound_nblength-1)*exBound_ndim
445 | n5 = 1+(exBound_nbfbegin -1)*2
446 | n6 = (exBound_nbfbegin+exBound_nbflength-1)*2
447 |
448 | write(*,*) (exBound_ibound1(n), &
449 | & exBound_ibound2(n),n=1,exBound_npbound)
450 |
451 | write(*,*) (exBound_ibound(n),n=n1,n2)
452 | write(*,*) (exBound_snbound(n),n=n3,n4)
453 | write(*,*) ((exBound_ifbound(k,n),k=1,4),n=1,exBound_npbound)
454 | write(*,*) (exBound_ibfael(n),n=n5,n6)
455 | write(*,*) exBound_nbaxipoints
456 | if(exBound_nbaxipoints.gt.0) then
457 | write(*,*) (exBound_axipoints(n),n=1,exBound_nbaxipoints)
458 | endif
459 | write(*,*) exBound_nbbndyconn
460 | if(exBound_nbbndyconn.gt.0) then
461 | write(*,*) (exBound_bndyconn(n),n=1,exBound_nbbndyconn)
462 | endif
463 |
464 |
465 | end subroutine writeexBoundvariables
466 | ! ================================================================
467 | end module avbp_exBound