OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_bem.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| hm_preread_bem ../starter/source/loads/bem/hm_read_bem.F
25!||--- called by ------------------------------------------------------
26!|| lectur ../starter/source/starter/lectur.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| arret ../starter/source/system/arret.F
30!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.f
31!|| hm_option_count ../starter/source/devtools/hm_reader/hm_option_count.F
32!|| hm_option_read_key ../starter/source/devtools/hm_reader/hm_option_read_key.F
33!|| hm_option_start ../starter/source/devtools/hm_reader/hm_option_start.F
34!||--- uses -----------------------------------------------------
35!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
36!|| message_mod ../starter/share/message_module/message_mod.F
37!|| submodel_mod ../starter/share/modules1/submodel_mod.F
38!||====================================================================
39 SUBROUTINE hm_preread_bem(IGRSURF,IGRNOD, NNFT ,
40 . UNITAB, NOM_OPT, LSUBMODEL)
41C-----------------------------------------------
42C M o d u l e s
43C-----------------------------------------------
44 USE unitab_mod
45 USE message_mod
46 USE groupdef_mod
47 USE submodel_mod
48 USE hm_option_read_mod
49 USE names_and_titles_mod , ONLY : nchartitle
50C-----------------------------------------------
51C I m p l i c i t T y p e s
52C-----------------------------------------------
53#include "implicit_f.inc"
54C-----------------------------------------------
55C C o m m o n B l o c k s
56C-----------------------------------------------
57#include "com01_c.inc"
58#include "com04_c.inc"
59#include "scr17_c.inc"
60#include "flowcom.inc"
61C-----------------------------------------------
62C D u m m y A r g u m e n t s
63C-----------------------------------------------
64 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
65 INTEGER NNFT
66 INTEGER NOM_OPT(LNOPT1,*)
67C-----------------------------------------------
68 TYPE (GROUP_) , DIMENSION(NGRNOD) :: IGRNOD
69 TYPE (SURF_) , DIMENSION(NSURF) :: IGRSURF
70 TYPE(submodel_data), DIMENSION(NSUBMOD), INTENT(IN) :: LSUBMODEL
71C-----------------------------------------------
72C L o c a l V a r i a b l e s
73C-----------------------------------------------
74 INTEGER I, ID, ITYP, II, NINOUT, J, ISU, NN, IAD, ITAG(NUMNOD),
75 . iad1, nel, j1, j2, j3, j4, ish34, nno, l, ibid,
76 . iii, igr, nnn, uid, iflagunit, iunit
77 INTEGER JFORM, FREESURF, NELMAX, HG
78 my_real
79 . rbid
80 CHARACTER(LEN=NCHARTITLE)::TITR
81 LOGICAL :: IS_AVAILABLE
82 INTEGER :: HM_NDAA, HM_NFLOW
83
84 hm_ndaa = 0
85 hm_nflow = 0
86 CALL hm_option_count('/BEM/DAA', hm_ndaa)
87 CALL hm_option_count('/BEM/FLOW', hm_nflow)
88C
89 liflow=0
90 lrflow=0
91 nnft=0
92 CALL hm_option_start('/BEM/FLOW')
93 DO i = 1, hm_nflow
94 CALL hm_option_read_key(lsubmodel, option_titr = titr, option_id = id)
95
96 ityp = 1
97C
98 CALL hm_get_intv('surf_IDex', ii, is_available, lsubmodel)
99 CALL hm_get_intv('Nio', ninout, is_available, lsubmodel)
100 isu=0
101 DO j=1,nsurf
102 IF (ii==igrsurf(j)%ID) isu=j
103 ENDDO
104 IF (isu==0) THEN
105 CALL ancmsg(msgid=621,msgtype=msgerror,anmode=aninfo,
106 . i1=id,c1=titr,c2='SURFACE',i2=ii)
107 ENDIF
108C
109 nn=igrsurf(isu)%NSEG
110 DO j=1,numnod
111 itag(j)=0
112 ENDDO
113 nel=0
114 DO j=1,nn
115 j1=igrsurf(isu)%NODES(j,1)
116 j2=igrsurf(isu)%NODES(j,2)
117 j3=igrsurf(isu)%NODES(j,3)
118 j4=igrsurf(isu)%NODES(j,4)
119 ish34=igrsurf(isu)%ELTYP(j)
120 itag(j1)=1
121 itag(j2)=1
122 itag(j3)=1
123 IF (ish34==3) THEN
124 nel=nel+2
125 itag(j4)=1
126 ELSEIF (ish34==7) THEN
127 nel=nel+1
128 ENDIF
129 ENDDO
130 nno=0
131 DO j=1,numnod
132 IF (itag(j)==1) nno=nno+1
133 ENDDO
134 nnft=nnft+nno
135C
136 CALL hm_get_intv('grn_IDaux', iii, is_available, lsubmodel)
137
138 nnn=0
139 IF (iii/=0) THEN
140 igr=0
141 DO j=1,ngrnod
142 IF (igrnod(j)%ID==iii) igr=j
143 ENDDO
144 IF (igr==0) THEN
145 CALL ancmsg(msgid=621,
146 . msgtype=msgerror,
147 . anmode=aninfo,
148 . i1=id,
149 . c1=titr,
150 . c2='NODE GROUP',
151 . i2=iii)
152 ENDIF
153 nnn=igrnod(igr)%NENTITY
154 ENDIF
155C
156 IF (nspmd == 1) THEN
157 liflow=liflow+niflow+nno+3*nel+ninout*niioflow+nnn+nel+nno+nnn
158 ELSE
159 liflow=liflow+niflow+nno+3*nel+ninout*niioflow+nnn+nel+4*nno+2*nnn+2*nel
160 ENDIF
161 lrflow=lrflow+nrflow+5*(nno+nnn)+ninout*nrioflow
162
163 ENDDO
164C
165
166
167 CALL hm_option_start('/BEM/DAA')
168 DO i = 1, hm_ndaa
169 CALL hm_option_read_key(lsubmodel, option_titr = titr, option_id = id)
170
171 ityp = 3
172C DAA SURFACES
173 nelmax=0
174 CALL hm_get_intv('surf_ID', ii, is_available, lsubmodel)
175 CALL hm_get_intv('Freesurf', freesurf, is_available, lsubmodel)
176 jform=2
177 IF(freesurf == 0) freesurf=1
178 isu=0
179 DO j=1,nsurf
180 IF (ii==igrsurf(j)%ID) isu=j
181 ENDDO
182 IF (isu==0) THEN
183 CALL ancmsg(msgid=1603,msgtype=msgerror,anmode=aninfo,
184 . i1=id,c1=titr,c2='SURFACE NUMBER NOT FOUND')
185 ENDIF
186C
187 nn= igrsurf(isu)%NSEG
188 DO j=1,numnod
189 itag(j)=0
190 ENDDO
191 nel=0
192 DO j=1,nn
193 j1=igrsurf(isu)%NODES(j,1)
194 j2=igrsurf(isu)%NODES(j,2)
195 j3=igrsurf(isu)%NODES(j,3)
196 j4=igrsurf(isu)%NODES(j,4)
197 ish34=igrsurf(isu)%ELTYP(j)
198 itag(j1)=1
199 itag(j2)=1
200 itag(j3)=1
201 IF (ish34==3) THEN
202 itag(j4)=1
203 IF(jform == 1) THEN
204 nel=nel+2
205 ELSEIF(jform == 2) THEN
206 nel=nel+1
207 ENDIF
208 ELSEIF (ish34==7) THEN
209 nel=nel+1
210 ENDIF
211 ENDDO
212 nno=0
213 DO j=1,numnod
214 IF (itag(j)==1) nno=nno+1
215 ENDDO
216 IF(nspmd == 1) THEN
217 IF(jform == 1) THEN
218 liflow=liflow+niflow+nno+3*nel+nno
219 ELSEIF(jform == 2) THEN
220 liflow=liflow+niflow+nno+5*nel+nno+nbgauge
221 ENDIF
222 ELSE
223 IF(jform == 1) THEN
224 liflow=liflow+niflow+nno+3*nel+nno+nno+2*nel
225 ELSEIF(jform == 2) THEN
226 liflow=liflow+niflow+nno+5*nel+nno+nbgauge+nno+2*nel
227 ENDIF
228 ENDIF
229 hg = huge(nel)
230 IF (nel > int(sqrt(real(hg)))) THEN
231 CALL ancmsg(msgid = 1711, anmode=aninfo, msgtype = msgerror,
232 . i1 = int(sqrt(real(hg))))
233 CALL arret(2)
234 ENDIF
235 IF(nel < nelmax) THEN
236 lrflow=lrflow+nrflow+7*nel+nel*nel+3*nel
237 ELSE
238 liflow=liflow+nel
239 lrflow=lrflow+nrflow+7*nel+2*nel*nel+3*nel
240 ENDIF
241 IF(freesurf == 2) lrflow=lrflow+3*nel
242 ENDDO
243C
244 RETURN
245 END
246!||====================================================================
247!|| hm_read_bem ../starter/source/loads/bem/hm_read_bem.F
248!||--- called by ------------------------------------------------------
249!|| lectur ../starter/source/starter/lectur.F
250!||--- calls -----------------------------------------------------
251!|| ancmsg ../starter/source/output/message/message.F
252!|| arret ../starter/source/system/arret.F
253!|| hm_get_float_array_index ../starter/source/devtools/hm_reader/hm_get_float_array_index.F
254!|| hm_get_float_array_index_dim ../starter/source/devtools/hm_reader/hm_get_float_array_index_dim.F
255!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
256!|| hm_get_floatv_dim ../starter/source/devtools/hm_reader/hm_get_floatv_dim.F
257!|| hm_get_int_array_index ../starter/source/devtools/hm_reader/hm_get_int_array_index.F
258!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
259!|| hm_option_count ../starter/source/devtools/hm_reader/hm_option_count.F
260!|| hm_option_read_key ../starter/source/devtools/hm_reader/hm_option_read_key.F
261!|| hm_option_start ../starter/source/devtools/hm_reader/hm_option_start.F
262!|| init_qd ../starter/source/fluid/init_qd.F
263!|| init_tg ../starter/source/fluid/init_tg.F
264!|| mass_fluid_qd ../starter/source/fluid/mass-fluid_qd.F
265!|| mass_fluid_tg ../starter/source/fluid/mass-fluid_tg.F
266!||--- uses -----------------------------------------------------
267!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
268!|| message_mod ../starter/share/message_module/message_mod.F
269!|| submodel_mod ../starter/share/modules1/submodel_mod.F
270!||====================================================================
271 SUBROUTINE hm_read_bem(IGRSURF, IFLOW , RFLOW ,
272 . NPC , IGRNOD , MEMFLOW,UNITAB,
273 . X , NOM_OPT , LGAUGE ,IGRV, LSUBMODEL,IRESP)
274C-----------------------------------------------
275C M o d u l e s
276C-----------------------------------------------
277 USE unitab_mod
278 USE message_mod
279 USE groupdef_mod
280 USE submodel_mod
281 USE hm_option_read_mod
282 USE names_and_titles_mod , ONLY : nchartitle, ncharkey
283C-----------------------------------------------
284C I m p l i c i t T y p e s
285C-----------------------------------------------
286#include "implicit_f.inc"
287C-----------------------------------------------
288C C o m m o n B l o c k s
289C-----------------------------------------------
290#include "com01_c.inc"
291#include "com04_c.inc"
292#include "scr17_c.inc"
293#include "param_c.inc"
294#include "units_c.inc"
295#include "flowcom.inc"
296C-----------------------------------------------
297C D u m m y A r g u m e n t s
298C-----------------------------------------------
299 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
300 INTEGER IFLOW(*),NPC(*)
301 INTEGER NOM_OPT(LNOPT1,*), LGAUGE(3,*), IGRV(NIGRV,*)
302 INTEGER(KIND=8) MEMFLOW(*)
303 my_real rflow(*), x(3,*)
304 TYPE(submodel_data), DIMENSION(NSUBMOD), INTENT(IN) :: LSUBMODEL
305 INTEGER,INTENT(IN) :: IRESP
306C-----------------------------------------------
307 TYPE (GROUP_) , DIMENSION(NGRNOD) :: IGRNOD
308 TYPE (SURF_) , DIMENSION(NSURF) :: IGRSURF
309C-----------------------------------------------
310C L o c a l V a r i a b l e s
311C-----------------------------------------------
312 INTEGER IADI, IADR, I, ID, ITYP, II, NINOUT, ISU, J, IPIMP,
313 . NN, IAD, ITAG(NUMNOD), IAD1, IAD2, J1, J2, J3, J4, ISH34,
314 . NNO, ITABINV(NUMNOD), NEL, IVFREE, IFVEL, IFPRES,
315 . k, ifunc, l, ifpa, iadmati, iadmatr, iii, igr, nnn,
316 . nbloc, iform, ilvout, itagio(numnod), ng1, ng2, ng3, ng4,
317 . n1, n2, n3, n4, prod, nprow, npcol, itest, ivini, ifvini,
318 . iinside, nrmax, nr, isuio,nflow0, nblocmax, uid,
319 . iflagunit,iunit
320 INTEGER JFORM, KFORM, II1, II2, II3, II4, II5, II6
321 INTEGER IR1, IR2, IR3,IR4, IR5, IR6, IR7, IR8, IR9, IR10, IR11
322 INTEGER IPRES, IWAVE, INTEGR, FREESURF, AFTERFLOW, GRAV_ID, NELMAX, IBID, HG
323 INTEGER, DIMENSION(:), ALLOCATABLE :: N_SHELL
324 my_real sfpa, sfvel, sfpres, scalt, dtsub, rho, tole, scalt_pa,
325 . sfvini, dirx, diry, dirz, norm, rnspmd, scalt_vi
326 my_real xc, yc, zc, xs, ys, zs, ssp, pmax, theta
327 my_real apmax, atheta, pmin
328 my_real xa, ya, za, xd, yd, zd, tt
329 my_real, DIMENSION(:,:), ALLOCATABLE :: cbem
330 my_real :: fac_gen
331 CHARACTER(LEN=NCHARKEY) :: KEY
332 CHARACTER(LEN=NCHARTITLE) :: TITR
333 LOGICAL :: IS_AVAILABLE
334 INTEGER :: HM_NDAA, HM_NFLOW
335C-----------------------------------------------
336 hm_ndaa = 0
337 hm_nflow = 0
338 CALL hm_option_count('/BEM/DAA', hm_ndaa)
339 CALL hm_option_count('/BEM/FLOW', hm_nflow)
340C
341 iadi=0
342 iadr=0
343 iadmati=1
344 iadmatr=1
345 nflow0=nflow
346 ifunc = 0
347
348 CALL hm_option_start('/BEM/FLOW')
349
350 DO i=1,hm_nflow
351 CALL hm_option_read_key(lsubmodel, option_titr = titr, option_id = id)
352
353 ityp = 1
354 CALL hm_get_intv('surf_idex', II, IS_AVAILABLE, LSUBMODEL)
355 CALL HM_GET_INTV('nio', NINOUT, IS_AVAILABLE, LSUBMODEL)
356 CALL HM_GET_INTV('iinside', IINSIDE, IS_AVAILABLE, LSUBMODEL)
357 CALL HM_GET_INTV('ifsp', IFPA, IS_AVAILABLE, LSUBMODEL)
358 CALL HM_GET_FLOATV('fscalesp', SFPA, IS_AVAILABLE, LSUBMODEL, UNITAB)
359 CALL HM_GET_FLOATV('ascalesp', SCALT, IS_AVAILABLE, LSUBMODEL, UNITAB)
360
361 ISU=0
362 DO J=1,NSURF
363 IF (II==IGRSURF(J)%ID) ISU=J
364 ENDDO
365 SCALT_PA=SCALT
366C
367 NN=IGRSURF(ISU)%NSEG
368 DO J=1,NUMNOD
369 ITAG(J)=0
370 ENDDO
371 DO J=1,NN
372 J1=IGRSURF(ISU)%NODES(J,1)
373 J2=IGRSURF(ISU)%NODES(J,2)
374 J3=IGRSURF(ISU)%NODES(J,3)
375 J4=IGRSURF(ISU)%NODES(J,4)
376 ISH34=IGRSURF(ISU)%ELTYP(J)
377 ITAG(J1)=1
378 ITAG(J2)=1
379 ITAG(J3)=1
380 IF (ISH34==3) ITAG(J4)=1
381 ENDDO
382 NNO=0
383 DO J=1,NUMNOD
384 IF (ITAG(J)==1) THEN
385 NNO=NNO+1
386 IFLOW(IADI+NIFLOW+NNO)=J
387 ITABINV(J)=NNO
388 ENDIF
389 ENDDO
390 NEL=0
391 DO J=1,NN
392 J1=IGRSURF(ISU)%NODES(J,1)
393 J2=IGRSURF(ISU)%NODES(J,2)
394 J3=IGRSURF(ISU)%NODES(J,3)
395 J4=IGRSURF(ISU)%NODES(J,4)
396 ISH34=IGRSURF(ISU)%ELTYP(J)
397 IF (ISH34==7) THEN
398 NEL=NEL+1
399 IFLOW(IADI+NIFLOW+NNO+3*(NEL-1)+1)=ITABINV(J1)
400 IFLOW(IADI+NIFLOW+NNO+3*(NEL-1)+2)=ITABINV(J2)
401 IFLOW(IADI+NIFLOW+NNO+3*(NEL-1)+3)=ITABINV(J3)
402 ELSEIF (ISH34==3) THEN
403 NEL=NEL+1
404 IFLOW(IADI+NIFLOW+NNO+3*(NEL-1)+1)=ITABINV(J1)
405 IFLOW(IADI+NIFLOW+NNO+3*(NEL-1)+2)=ITABINV(J2)
406 IFLOW(IADI+NIFLOW+NNO+3*(NEL-1)+3)=ITABINV(J4)
407 NEL=NEL+1
408 IFLOW(IADI+NIFLOW+NNO+3*(NEL-1)+1)=ITABINV(J2)
409 IFLOW(IADI+NIFLOW+NNO+3*(NEL-1)+2)=ITABINV(J3)
410 IFLOW(IADI+NIFLOW+NNO+3*(NEL-1)+3)=ITABINV(J4)
411 ENDIF
412 ENDDO
413C
414 IF (IINSIDE/=2) IINSIDE=1
415C
416 IPIMP=0
417 IF (IFPA/=0) THEN
418 IFUNC=0
419 DO J=1,NFUNCT
420 IF (IFPA==NPC(J)) IFUNC=J
421 ENDDO
422 IF (IFUNC==0) THEN
423 CALL ANCMSG(MSGID=621,
424 . MSGTYPE=MSGERROR,
425 . ANMODE=ANINFO,
426 . I1=ID,
427 . C1=TITR,
428 . C2='function',
429 . I2=IFPA)
430 ENDIF
431 IFPA=IFUNC
432 IPIMP=IPIMP+1
433 ENDIF
434C
435 CALL HM_GET_INTV('grn_idaux', III, IS_AVAILABLE, LSUBMODEL)
436 CALL HM_GET_INTV('itest', ITEST, IS_AVAILABLE, LSUBMODEL)
437 CALL HM_GET_FLOATV('tole', TOLE, IS_AVAILABLE, LSUBMODEL, UNITAB)
438
439.AND. IF (IINSIDE==2ITEST==1) ITEST=2
440 IF (TOLE==ZERO) TOLE=EM5
441C
442 CALL HM_GET_FLOATV('rho', RHO, IS_AVAILABLE, LSUBMODEL, UNITAB)
443 CALL HM_GET_INTV('ivinf', IVINI, IS_AVAILABLE, LSUBMODEL)
444C
445 NNN=0
446 IF (III/=0) THEN
447 DO J=1,NGRNOD
448 IF (IGRNOD(J)%ID==III) IGR=J
449 ENDDO
450 NNN=IGRNOD(IGR)%NENTITY
451 DO J=1,NNN
452 IFLOW(IADI+NIFLOW+NNO+3*NEL+NINOUT*NIIOFLOW+J)=
453 . IGRNOD(IGR)%ENTITY(J)
454 ENDDO
455 ENDIF
456C
457 IVFREE=0
458 DO J = 1, NINOUT
459 CALL HM_GET_INT_ARRAY_INDEX('surf_idio', II, J, IS_AVAILABLE, LSUBMODEL)
460 CALL HM_GET_INT_ARRAY_INDEX('fct_idvel', IFVEL, J, IS_AVAILABLE, LSUBMODEL)
461 CALL HM_GET_INT_ARRAY_INDEX('fct_idpr', IFPRES, J, IS_AVAILABLE, LSUBMODEL)
462 CALL HM_GET_FLOAT_ARRAY_INDEX('fscalenv', SFVEL, J, IS_AVAILABLE, LSUBMODEL, UNITAB)
463 CALL HM_GET_FLOAT_ARRAY_INDEX('fscalepres', SFPRES, J, IS_AVAILABLE, LSUBMODEL, UNITAB)
464 CALL HM_GET_FLOAT_ARRAY_INDEX('ascalet', SCALT, J, IS_AVAILABLE, LSUBMODEL, UNITAB)
465
466
467 IF(SFVEL == ZERO) THEN
468 CALL HM_GET_FLOAT_ARRAY_INDEX_DIM('sfvel', FAC_GEN, J, IS_AVAILABLE, LSUBMODEL, UNITAB)
469 SFVEL = ONE * FAC_GEN
470 ENDIF
471 IF(SFPRES == ZERO) THEN
472 CALL HM_GET_FLOAT_ARRAY_INDEX_DIM('sfpres', FAC_GEN, J, IS_AVAILABLE, LSUBMODEL, UNITAB)
473 SFPRES = ONE * FAC_GEN
474 ENDIF
475 IF(SCALT == ZERO) THEN
476 CALL HM_GET_FLOAT_ARRAY_INDEX_DIM('scal_t', FAC_GEN, J, IS_AVAILABLE, LSUBMODEL, UNITAB)
477 SCALT = ONE * FAC_GEN
478 ENDIF
479C
480 ISUIO=0
481 DO K=1,NSURF
482 IF (II==IGRSURF(K)%ID) ISUIO=K
483 ENDDO
484 IF (ISUIO==0) THEN
485 CALL ANCMSG(MSGID=621,
486 . MSGTYPE=MSGERROR,
487 . ANMODE=ANINFO,
488 . I1=ID,
489 . C1=TITR,
490 . C2='surface',
491 . I2=II)
492 ENDIF
493 IFLOW(IADI+NIFLOW+NNO+3*NEL+NIIOFLOW*(J-1)+1)=ISUIO
494C
495 IFLOW(IADI+NIFLOW+NNO+3*NEL+NIIOFLOW*(J-1)+2)=0
496C
497 IF (IFVEL/=0) THEN
498 IFUNC=0
499 DO K=1,NFUNCT
500 IF (IFVEL==NPC(K)) IFUNC=K
501 ENDDO
502 IF (IFUNC==0) THEN
503 CALL ANCMSG(MSGID=621,
504 . MSGTYPE=MSGERROR,
505 . ANMODE=ANINFO,
506 . I1=ID,
507 . C1=TITR,
508 . C2='function',
509 . I2=IFVEL)
510 ENDIF
511 IFVEL=IFUNC
512 ELSE
513 IVFREE=IVFREE+1
514 ENDIF
515 IFLOW(IADI+NIFLOW+NNO+3*NEL+NIIOFLOW*(J-1)+3)=IFVEL
516 RFLOW(IADR+NRFLOW+5*(NNO+NNN)+NRIOFLOW*(J-1)+1)=SFVEL
517C
518 IF (IFPRES/=0) THEN
519 IPIMP=IPIMP+1
520 IFUNC=0
521 DO K=1,NFUNCT
522 IF (IFPRES==NPC(K)) IFUNC=K
523 ENDDO
524 IF (IFUNC==0) THEN
525 CALL ANCMSG(MSGID=621,
526 . MSGTYPE=MSGERROR,
527 . ANMODE=ANINFO,
528 . I1=ID,
529 . C1=TITR,
530 . C2='function',
531 . I2=IFPRES)
532 ENDIF
533 IFPRES=IFUNC
534 ENDIF
535 IFLOW(IADI+NIFLOW+NNO+3*NEL+NIIOFLOW*(J-1)+4)=IFPRES
536 RFLOW(IADR+NRFLOW+5*(NNO+NNN)+NRIOFLOW*(J-1)+2)=SFPRES
537C
538 RFLOW(IADR+NRFLOW+5*(NNO+NNN)+NRIOFLOW*(J-1)+3)=SCALT
539 ENDDO
540 IF (NINOUT>=1) THEN
541 IF (IPIMP/=1) THEN
542 CALL ANCMSG(MSGID=622,
543 . MSGTYPE=MSGERROR,
544 . ANMODE=ANINFO,
545 . I1=ID,
546 . C1=TITR,
547 . C2='one and only one pressure must be imposed')
548.AND..OR. ELSEIF (IVFREE/=1(IINSIDE==0
549.AND. . (IINSIDE==1IVINI==0))) THEN
550 CALL ANCMSG(MSGID=622,
551 . MSGTYPE=MSGERROR,
552 . ANMODE=ANINFO,
553 . I1=ID,
554 . C1=TITR,
555 . C2='one and only one velocity must be left free')
556.AND..AND. ELSEIF (IVFREE/=0IINSIDE==1IVINI==1) THEN
557 CALL ANCMSG(MSGID=622,
558 . MSGTYPE=MSGERROR,
559 . ANMODE=ANINFO,
560 . I1=ID,
561 . C1=TITR,
562 . C2='no free velocity allowed')
563 ENDIF
564 ELSE
565 IF (IFPA==0) THEN
566 CALL ANCMSG(MSGID=622,
567 . MSGTYPE=MSGERROR,
568 . ANMODE=ANINFO,
569 . I1=ID,
570 . C1=TITR,
571 . C2='no imposed pressure')
572 ENDIF
573 ENDIF
574C tag of input and output elements
575 DO J=1,NNO
576 ITAGIO(J)=0
577 ENDDO
578 DO J=1,NINOUT
579 ISUIO=IFLOW(IADI+NIFLOW+NNO+3*NEL+NIIOFLOW*(J-1)+1)
580 NN=IGRSURF(ISUIO)%NSEG
581 DO K=1,NN
582 NG1=IGRSURF(ISUIO)%NODES(K,1)
583 NG2=IGRSURF(ISUIO)%NODES(K,2)
584 NG3=IGRSURF(ISUIO)%NODES(K,3)
585 NG4=IGRSURF(ISUIO)%NODES(K,4)
586 ISH34=IGRSURF(ISUIO)%ELTYP(K)
587 N1=ITABINV(NG1)
588 N2=ITABINV(NG2)
589 N3=ITABINV(NG3)
590 ITAGIO(N1)=J
591 ITAGIO(N2)=J
592 ITAGIO(N3)=J
593 IF (ISH34==3) THEN
594 N4=ITABINV(NG4)
595 ITAGIO(N4)=J
596 ENDIF
597 ENDDO
598 ENDDO
599C
600 DO J=1,NEL
601 N1=IFLOW(IADI+NIFLOW+NNO+3*(J-1)+1)
602 N2=IFLOW(IADI+NIFLOW+NNO+3*(J-1)+2)
603 N3=IFLOW(IADI+NIFLOW+NNO+3*(J-1)+3)
604 PROD=ITAGIO(N1)*ITAGIO(N2)*ITAGIO(N3)
605 IFLOW(IADI+NIFLOW+NNO+3*NEL+NINOUT*NIIOFLOW+NNN+J)=0
606 IF (PROD/=0)
607 . IFLOW(IADI+NIFLOW+NNO+3*NEL+NINOUT*NIIOFLOW+NNN+J)=
608 . MAX(ITAGIO(N1),ITAGIO(N2),ITAGIO(N3))
609 ENDDO
610C
611 CALL HM_GET_INTV('iform', IFORM, IS_AVAILABLE, LSUBMODEL)
612 CALL HM_GET_INTV('ipri', ILVOUT, IS_AVAILABLE, LSUBMODEL)
613 CALL HM_GET_FLOATV('dtflow', DTSUB, IS_AVAILABLE, LSUBMODEL, UNITAB)
614
615
616 IF (IFORM==0) IFORM=1
617C
618 NBLOC=0
619 NPROW=1
620 NPCOL=1
621 IF (NSPMD > 1) THEN
622 RNSPMD=NSPMD
623 NRMAX=INT(SQRT(RNSPMD))
624 IF (IFORM==1) THEN
625 DO NR=1,NRMAX
626 IF (MOD(NSPMD,NR)==0) THEN
627 NPROW=NR
628 NPCOL=NSPMD/NR
629 ENDIF
630 ENDDO
631 ELSEIF (IFORM==2) THEN
632 NPROW=NSPMD
633 NPCOL=1
634 ENDIF
635 IF (NNO<1000) THEN
636 NBLOCMAX=32
637 ELSE
638 NBLOCMAX=64
639 ENDIF
640 NBLOC=MIN(NNO/NPROW, NNO/NPCOL)
641 NBLOC=MIN(NBLOCMAX,NBLOC)
642 NBLOC=MAX(1,NBLOC)
643 ENDIF
644C
645 IFVINI=0
646 SFVINI=ZERO
647 SCALT_VI=ZERO
648 DIRX=ZERO
649 DIRY=ZERO
650 DIRZ=ZERO
651 IF (IVINI==1) THEN
652 CALL HM_GET_INTV('ifvinf', IFVINI, IS_AVAILABLE, LSUBMODEL)
653 CALL HM_GET_FLOATV('fscalevel', SFVINI, IS_AVAILABLE, LSUBMODEL, UNITAB)
654 CALL HM_GET_FLOATV('ascalevel', SCALT_VI, IS_AVAILABLE, LSUBMODEL, UNITAB)
655
656 CALL HM_GET_FLOATV('dirx', DIRX, IS_AVAILABLE, LSUBMODEL, UNITAB)
657 CALL HM_GET_FLOATV('diry', DIRY, IS_AVAILABLE, LSUBMODEL, UNITAB)
658 CALL HM_GET_FLOATV('dirz', DIRZ, IS_AVAILABLE, LSUBMODEL, UNITAB)
659
660 IFUNC=0
661 DO J=1,NFUNCT
662 IF (IFVINI==NPC(J)) IFUNC=J
663 ENDDO
664 IF (IFUNC==0) THEN
665 CALL ANCMSG(MSGID=621,
666 . MSGTYPE=MSGERROR,
667 . ANMODE=ANINFO,
668 . I1=ID,
669 . C1=TITR,
670 . C2='function',
671 . I2=IFVINI)
672 ENDIF
673 IFVINI=IFUNC
674 NORM=SQRT(DIRX**2+DIRY**2+DIRZ**2)
675 IF (NORM==ZERO) THEN
676 CALL ANCMSG(MSGID=622,
677 . MSGTYPE=MSGERROR,
678 . ANMODE=ANINFO,
679 . I1=ID,
680 . C1=TITR,
681 . C2='null velocity direction vector')
682 ENDIF
683 DIRX=DIRX/NORM
684 DIRY=DIRY/NORM
685 DIRZ=DIRZ/NORM
686 ENDIF
687C
688 IFLOW(IADI+1)=ID
689 IFLOW(IADI+2)=ITYP
690 IFLOW(IADI+3)=ISU
691 IFLOW(IADI+4)=NINOUT
692 IFLOW(IADI+5)=NNO
693 IFLOW(IADI+6)=NEL
694 IFLOW(IADI+7)=NNN
695 IFLOW(IADI+8)=NNO
696 IFLOW(IADI+9)=NNO*NNO+NNO*(NEL+1)
697 IFLOW(IADI+10)=IADMATI
698 IFLOW(IADI+11)=IADMATR
699 IFLOW(IADI+12)=NBLOC
700 IFLOW(IADI+13)=IFORM
701 IF (NSPMD == 1) THEN
702 IFLOW(IADI+14)=NIFLOW+NNO+3*NEL+NINOUT*NIIOFLOW+NNN+NEL+NNO+NNN
703 ELSE
704 IFLOW(IADI+14)=NIFLOW+NNO+3*NEL+NINOUT*NIIOFLOW+NNN+NEL+4*NNO+2*NNN+2*NEL
705 ENDIF
706 IFLOW(IADI+15)=NRFLOW+5*(NNO+NNN)+NINOUT*NRIOFLOW
707 IFLOW(IADI+17)=ILVOUT
708 IFLOW(IADI+18)=NPROW
709 IFLOW(IADI+19)=NPCOL
710 IFLOW(IADI+20)=IADMATR+NNO*NNO
711 IFLOW(IADI+21)=ITEST
712 IFLOW(IADI+23)=IFPA
713 IFLOW(IADI+24)=IFVINI
714 RFLOW(IADR+1)=SFPA
715 RFLOW(IADR+2)=SCALT_PA
716 RFLOW(IADR+3)=DTSUB
717 RFLOW(IADR+4)=ZERO
718 RFLOW(IADR+5)=RHO
719 RFLOW(IADR+6)=TOLE
720 RFLOW(IADR+7)=SFVINI
721 RFLOW(IADR+8)=SCALT_VI
722 RFLOW(IADR+9)=DIRX
723 RFLOW(IADR+10)=DIRY
724 RFLOW(IADR+11)=DIRZ
725 RFLOW(IADR+12)=ZERO
726C
727 MEMFLOW(1)=MEMFLOW(1)+IFLOW(IADI+8)
728 MEMFLOW(2)=MEMFLOW(2)+IFLOW(IADI+9)
729C PRINTOUTS
730 WRITE(IOUT,1000)
731 WRITE(IOUT,1100) I, ID, TRIM(TITR),IGRSURF(ISU)%ID, NNO, NEL
732 IF (IINSIDE==1) THEN
733 WRITE(IOUT,1110)
734 ELSEIF (IINSIDE==2) THEN
735 WRITE(IOUT,1120)
736 ENDIF
737 WRITE(IOUT,1200) IFPA, SFPA, SCALT_PA, III, NNN
738 IF (ITEST==1) THEN
739 WRITE(IOUT,1210)
740 WRITE(IOUT,1225) TOLE
741 ELSEIF (ITEST==2) THEN
742 WRITE(IOUT,1220)
743 WRITE(IOUT,1225) TOLE
744 ELSEIF (ITEST==0) THEN
745 WRITE(IOUT,1230)
746 ENDIF
747 WRITE(IOUT,1300) RHO
748 WRITE(IOUT,1400) NINOUT
749 DO J=1,NINOUT
750 ISUIO=IFLOW(IADI+NIFLOW+NNO+3*NEL+NIIOFLOW*(J-1)+1)
751 IFVEL=IFLOW(IADI+NIFLOW+NNO+3*NEL+NIIOFLOW*(J-1)+3)
752 SFVEL=RFLOW(IADR+NRFLOW+5*(NNO+NNN)+NRIOFLOW*(J-1)+1)
753 IFPRES=IFLOW(IADI+NIFLOW+NNO+3*NEL+NIIOFLOW*(J-1)+4)
754 SFPRES=RFLOW(IADR+NRFLOW+5*(NNO+NNN)+NRIOFLOW*(J-1)+2)
755 SCALT=RFLOW(IADR+NRFLOW+5*(NNO+NNN)+NRIOFLOW*(J-1)+3)
756 WRITE(IOUT,1410) J, IGRSURF(ISUIO)%ID
757 IF (IFVEL>0) WRITE(IOUT,1420) IFVEL, SFVEL
758 IF (IFPRES>0) WRITE(IOUT,1430) IFPRES, SFPRES
759 WRITE(IOUT,1440) SCALT
760 ENDDO
761 WRITE(IOUT,1500) IFORM, ILVOUT, DTSUB
762 IF (IVINI==1) WRITE(IOUT,1600) IFVINI, SFVINI, SCALT_VI,
763 . DIRX, DIRY, DIRZ
764 IF (NSPMD > 1) WRITE(IOUT,1700) NPROW, NPCOL, NBLOC
765C
766 IADR=IADR+IFLOW(IADI+15)
767 IADI=IADI+IFLOW(IADI+14)
768 IADMATI=IADMATI+NNO
769 IADMATR=IADMATR+NNO*NNO+NNO*(NEL+1)
770
771 ENDDO
772
773 CALL HM_OPTION_START('/bem/daa')
774
775 DO I=1,HM_NDAA
776 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_TITR = TITR, OPTION_ID = ID)
777C
778
779 CALL HM_GET_INTV('surf_id', II, IS_AVAILABLE, LSUBMODEL)
780 CALL HM_GET_INTV('grav_id', GRAV_ID, IS_AVAILABLE, LSUBMODEL)
781 NELMAX=0
782 CALL HM_GET_FLOATV('rho', RHO, IS_AVAILABLE, LSUBMODEL, UNITAB)
783 CALL HM_GET_FLOATV('c', SSP, IS_AVAILABLE, LSUBMODEL, UNITAB)
784 CALL HM_GET_FLOATV('pmin', PMIN, IS_AVAILABLE, LSUBMODEL, UNITAB)
785
786 IF(PMIN == ZERO) PMIN=-EP30
787
788 CALL HM_GET_FLOATV('xs', XS, IS_AVAILABLE, LSUBMODEL, UNITAB)
789 CALL HM_GET_FLOATV('ys', YS, IS_AVAILABLE, LSUBMODEL, UNITAB)
790 CALL HM_GET_FLOATV('zs', ZS, IS_AVAILABLE, LSUBMODEL, UNITAB)
791
792 CALL HM_GET_INTV('iform', IFORM, IS_AVAILABLE, LSUBMODEL)
793 CALL HM_GET_INTV('ipri', ILVOUT, IS_AVAILABLE, LSUBMODEL)
794 CALL HM_GET_INTV('ipres', IPRES, IS_AVAILABLE, LSUBMODEL)
795 CALL HM_GET_INTV('kform', KFORM, IS_AVAILABLE, LSUBMODEL)
796 CALL HM_GET_INTV('freesurf', FREESURF, IS_AVAILABLE, LSUBMODEL)
797 CALL HM_GET_INTV('afterflow', AFTERFLOW, IS_AVAILABLE, LSUBMODEL)
798 CALL HM_GET_INTV('integr', INTEGR, IS_AVAILABLE, LSUBMODEL)
799 ITYP = 3
800 IF(IFORM == 0) IFORM=1
801 IF(ILVOUT == 0) ILVOUT=1
802 IWAVE=1
803 JFORM=2
804 IF(KFORM == 0) KFORM=1
805 IF(INTEGR == 0) INTEGR=2
806 IF(FREESURF == 0) FREESURF=1
807 IF(AFTERFLOW == 0) AFTERFLOW=2
808.AND. IF (IWAVE ==2 FREESURF == 2) THEN
809 CALL ANCMSG(MSGID=1603,
810 . MSGTYPE=MSGERROR,
811 . ANMODE=ANINFO,
812 . I1=ID,
813 . C1=TITR,
814 . C2='free surface is not compatible with plane wave')
815 ENDIF
816 IF(KFORM ==2) INTEGR=1
817C
818.AND. IF(NBGAUGE > 0 JFORM == 2) ALLOCATE(N_SHELL(NUMELC))
819C
820 ISU=0
821 DO J=1,NSURF
822 IF (II==IGRSURF(J)%ID) ISU=J
823 ENDDO
824C
825 NN =IGRSURF(ISU)%NSEG
826 DO J=1,NUMNOD
827 ITAG(J)=0
828 ENDDO
829 DO J=1,NN
830 J1=IGRSURF(ISU)%NODES(J,1)
831 J2=IGRSURF(ISU)%NODES(J,2)
832 J3=IGRSURF(ISU)%NODES(J,3)
833 J4=IGRSURF(ISU)%NODES(J,4)
834 ISH34=IGRSURF(ISU)%ELTYP(J)
835 ITAG(J1)=1
836 ITAG(J2)=1
837 ITAG(J3)=1
838 IF (ISH34==3) ITAG(J4)=1
839 ENDDO
840 NNO=0
841 DO J=1,NUMNOD
842 IF (ITAG(J)==1) THEN
843 NNO=NNO+1
844 IFLOW(IADI+NIFLOW+NNO)=J
845 ITABINV(J)=NNO
846 ENDIF
847 ENDDO
848 NEL=0
849 DO J=1,NN
850 J1=IGRSURF(ISU)%NODES(J,1)
851 J2=IGRSURF(ISU)%NODES(J,2)
852 J3=IGRSURF(ISU)%NODES(J,3)
853 J4=IGRSURF(ISU)%NODES(J,4)
854 ISH34=IGRSURF(ISU)%ELTYP(J)
855 IF(JFORM == 1) THEN
856 IF (ISH34==7) THEN
857 NEL=NEL+1
858 IAD2=IADI+NIFLOW+NNO+3*(NEL-1)
859 IFLOW(IAD2+1)=ITABINV(J1)
860 IFLOW(IAD2+2)=ITABINV(J2)
861 IFLOW(IAD2+3)=ITABINV(J3)
862 ELSEIF (ISH34==3) THEN
863 NEL=NEL+1
864 IAD2=IADI+NIFLOW+NNO+3*(NEL-1)
865 IFLOW(IAD2+1)=ITABINV(J1)
866 IFLOW(IAD2+2)=ITABINV(J2)
867 IFLOW(IAD2+3)=ITABINV(J4)
868 NEL=NEL+1
869 IAD2=IADI+NIFLOW+NNO+3*(NEL-1)
870 IFLOW(IAD2+1)=ITABINV(J2)
871 IFLOW(IAD2+2)=ITABINV(J3)
872 IFLOW(IAD2+3)=ITABINV(J4)
873 ENDIF
874 ELSEIF(JFORM == 2) THEN
875 IF (ISH34==7) THEN
876 NEL =NEL+1
877 IAD2=IADI+NIFLOW+NNO+5*(NEL-1)
878 IFLOW(IAD2+1)=ITABINV(J1)
879 IFLOW(IAD2+2)=ITABINV(J2)
880 IFLOW(IAD2+3)=ITABINV(J3)
881 IFLOW(IAD2+4)=ITABINV(J3)
882 IFLOW(IAD2+5)=2
883 IF(NBGAUGE > 0) N_SHELL(NEL)=0
884 ELSEIF (ISH34==3) THEN
885 NEL =NEL+1
886 IAD2=IADI+NIFLOW+NNO+5*(NEL-1)
887 IFLOW(IAD2+1)=ITABINV(J1)
888 IFLOW(IAD2+2)=ITABINV(J2)
889 IFLOW(IAD2+3)=ITABINV(J3)
890 IFLOW(IAD2+4)=ITABINV(J4)
891 IFLOW(IAD2+5)=1
892 IF(NBGAUGE > 0) N_SHELL(NEL)=IGRSURF(ISU)%ELEM(J)
893 ENDIF
894 ENDIF
895 ENDDO
896C
897 PMAX = ZERO
898 THETA = ZERO
899 APMAX = ZERO
900 ATHETA = ZERO
901 IFPRES = 0
902 SFPRES = ZERO
903 IF(IPRES == 1)THEN
904 CALL HM_GET_FLOATV('pm', PMAX, IS_AVAILABLE, LSUBMODEL, UNITAB)
905 CALL HM_GET_FLOATV('theta', THETA, IS_AVAILABLE, LSUBMODEL, UNITAB)
906 CALL HM_GET_FLOATV('a', APMAX, IS_AVAILABLE, LSUBMODEL, UNITAB)
907 CALL HM_GET_FLOATV('atheta', ATHETA, IS_AVAILABLE, LSUBMODEL, UNITAB)
908
909 IF(APMAX == ZERO) APMAX = ONE
910 ELSEIF(IPRES == 2) THEN
911C Initialization
912 PMAX = EP30
913 THETA = EP30
914 CALL HM_GET_INTV('fct_idp', IFPRES, IS_AVAILABLE, LSUBMODEL)
915 CALL HM_GET_FLOATV('fscalep', SFPRES, IS_AVAILABLE, LSUBMODEL, UNITAB)
916 IF(SFPRES == ZERO) THEN
917 CALL HM_GET_FLOATV_DIM('fscalep', FAC_GEN, IS_AVAILABLE, LSUBMODEL, UNITAB)
918 SFPRES = ONE * FAC_GEN
919 ENDIF
920 IF (IFPRES/=0) THEN
921 IFUNC=0
922 DO K=1,NFUNCT
923 IF (IFPRES==NPC(K)) IFUNC=K
924 ENDDO
925 IF (IFUNC==0) THEN
926 CALL ANCMSG(MSGID=621,MSGTYPE=MSGERROR,ANMODE=ANINFO,
927 . I1=ID,C1=TITR,C2='function',I2=IFPRES)
928 ENDIF
929 ENDIF
930 ELSEIF(IPRES == 3) THEN
931
932 ELSEIF(IPRES == 4) THEN
933 ENDIF
934
935 CALL HM_GET_FLOATV('xc', XC, IS_AVAILABLE, LSUBMODEL, UNITAB)
936 CALL HM_GET_FLOATV('yc', YC, IS_AVAILABLE, LSUBMODEL, UNITAB)
937 CALL HM_GET_FLOATV('zc', ZC, IS_AVAILABLE, LSUBMODEL, UNITAB)
938
939
940
941
942 XD=ZERO
943 YD=ZERO
944 ZD=ZERO
945 XA=ZERO
946 YA=ZERO
947 ZA=ZERO
948 DIRX=ZERO
949 DIRY=ZERO
950 DIRZ=ZERO
951.OR. IF(FREESURF == 2 GRAV_ID > 0)THEN
952 CALL HM_GET_FLOATV('xa', XA, IS_AVAILABLE, LSUBMODEL, UNITAB)
953 CALL HM_GET_FLOATV('ya', YA, IS_AVAILABLE, LSUBMODEL, UNITAB)
954 CALL HM_GET_FLOATV('za', ZA, IS_AVAILABLE, LSUBMODEL, UNITAB)
955 CALL HM_GET_FLOATV('dir-x', DIRX, IS_AVAILABLE, LSUBMODEL, UNITAB)
956 CALL HM_GET_FLOATV('dir-y', DIRY, IS_AVAILABLE, LSUBMODEL, UNITAB)
957 CALL HM_GET_FLOATV('dir-z', DIRZ, IS_AVAILABLE, LSUBMODEL, UNITAB)
958 NORM=SQRT(DIRX**2+DIRY**2+DIRZ**2)
959 IF (NORM==ZERO) THEN
960 CALL ANCMSG(MSGID=622,
961 . MSGTYPE=MSGERROR,ANMODE=ANINFO,
962 . I1=ID,C1=TITR,C2='null free surface normal')
963 ENDIF
964 DIRX=DIRX/NORM
965 DIRY=DIRY/NORM
966 DIRZ=DIRZ/NORM
967C Compute Image D of Charge C
968 IF(FREESURF == 2) THEN
969 TT = DIRX*(XC-XA)+DIRY*(YC-YA)+DIRZ*(ZC-ZA)
970 XD = XC - TWO*TT*DIRX
971 YD = YC - TWO*TT*DIRY
972 ZD = ZC - TWO*TT*DIRZ
973 ENDIF
974 ENDIF
975
976 NBLOC=0
977 NPROW=1
978 NPCOL=1
979 IF (NSPMD > 1) THEN
980 RNSPMD=NSPMD
981 NRMAX=INT(SQRT(RNSPMD))
982 DO NR=1,NRMAX
983 IF (MOD(NSPMD,NR)==0) THEN
984 NPROW=NR
985 NPCOL=NSPMD/NR
986 ENDIF
987 ENDDO
988 IF (NEL<1000) THEN
989 NBLOCMAX=32
990 ELSE
991 NBLOCMAX=64
992 ENDIF
993 NBLOC=MIN(NEL/NPROW, NEL/NPCOL)
994 NBLOC=MIN(NBLOCMAX,NBLOC)
995 NBLOC=MAX(1,NBLOC)
996 ENDIF
997C
998 IFLOW(IADI+1)=ID
999 IFLOW(IADI+2)=ITYP
1000 IFLOW(IADI+3)=ISU
1001 IFLOW(IADI+4)=JFORM
1002 IFLOW(IADI+5)=NNO
1003 IFLOW(IADI+6)=NEL
1004 IFLOW(IADI+7)=IFUNC
1005 IFLOW(IADI+8)=0
1006 IFLOW(IADI+9)=0
1007 IFLOW(IADI+10)=IADMATI
1008 IFLOW(IADI+11)=IADMATR
1009 IFLOW(IADI+12)=NBLOC
1010 IFLOW(IADI+13)=IFORM
1011 IF(NSPMD == 1) THEN
1012 IF(JFORM == 1) THEN
1013 IFLOW(IADI+14)=NIFLOW+NNO+3*NEL+NNO
1014 ELSEIF(JFORM == 2) THEN
1015 IFLOW(IADI+14)=NIFLOW+NNO+5*NEL+NNO+NBGAUGE
1016 ENDIF
1017 ELSE
1018 IF(JFORM == 1) THEN
1019 IFLOW(IADI+14)=NIFLOW+NNO+3*NEL+NNO+NNO+NEL
1020 ELSEIF(JFORM == 2) THEN
1021 IFLOW(IADI+14)=NIFLOW+NNO+5*NEL+NNO+NBGAUGE+NNO+NEL
1022 ENDIF
1023 ENDIF
1024 HG = HUGE(NEL)
1025 IF (NEL > INT(SQRT(REAL(HG)))) THEN
1026 CALL ANCMSG(MSGID = 1711, ANMODE=ANINFO, MSGTYPE = MSGERROR,
1027 . I1 = INT(SQRT(REAL(HG))))
1028 CALL ARRET(2)
1029 ENDIF
1030 IF(FREESURF == 1)IFLOW(IADI+15)=NRFLOW+ 7*NEL+NEL*NEL+3*NEL
1031 IF(FREESURF == 2)IFLOW(IADI+15)=NRFLOW+10*NEL+NEL*NEL+3*NEL
1032 IFLOW(IADI+17)=ILVOUT
1033 IFLOW(IADI+18)=NPROW
1034 IFLOW(IADI+19)=NPCOL
1035 IFLOW(IADI+20)=0
1036 IFLOW(IADI+21)=IPRES
1037 IFLOW(IADI+22)=IWAVE
1038 IFLOW(IADI+23)=KFORM
1039 IFLOW(IADI+24)=INTEGR
1040 IFLOW(IADI+25)=FREESURF
1041 IFLOW(IADI+26)=AFTERFLOW
1042 K=0
1043 IF(GRAV_ID > 0) THEN
1044 DO J=1,NGRAV
1045 IF(IGRV(5,J) == GRAV_ID) K=J
1046 ENDDO
1047 ENDIF
1048 IFLOW(IADI+27)=K
1049 IFLOW(IADI+28)=NELMAX
1050
1051 RFLOW(IADR+1)=RHO*SSP
1052 RFLOW(IADR+2)=SSP
1053 RFLOW(IADR+3)=SFPRES
1054 RFLOW(IADR+5)=RHO
1055 RFLOW(IADR+6)=PMAX
1056 RFLOW(IADR+7)=THETA
1057 RFLOW(IADR+8)=-ONE
1058 RFLOW(IADR+9) =XC
1059 RFLOW(IADR+10)=YC
1060 RFLOW(IADR+11)=ZC
1061 RFLOW(IADR+12)=SQRT((XS-XC)**2+(YS-YC)**2+(ZS-ZC)**2)
1062 RFLOW(IADR+13)=XD
1063 RFLOW(IADR+14)=YD
1064 RFLOW(IADR+15)=ZD
1065
1066 RFLOW(IADR+16)=XA
1067 RFLOW(IADR+17)=YA
1068 RFLOW(IADR+18)=ZA
1069 RFLOW(IADR+19)=DIRX
1070 RFLOW(IADR+20)=DIRY
1071 RFLOW(IADR+21)=DIRZ
1072 RFLOW(IADR+22)=PMIN
1073 RFLOW(IADR+23)=APMAX
1074 RFLOW(IADR+24)=ATHETA
1075
1076 MEMFLOW(1)=MEMFLOW(1)+IFLOW(IADI+8)
1077 MEMFLOW(2)=MEMFLOW(2)+IFLOW(IADI+9)
1078C
1079 WRITE(IOUT,2000)
1080 IF(JFORM == 1) WRITE(IOUT,2100) I, ID, TRIM(TITR), IGRSURF(ISU)%ID, NNO, NEL, GRAV_ID
1081 IF(JFORM == 2) WRITE(IOUT,2200) I, ID, TRIM(TITR), IGRSURF(ISU)%ID, NNO, NEL, GRAV_ID
1082 WRITE(IOUT,2300) RHO, SSP, PMIN
1083 WRITE(IOUT,2400) XS, YS, ZS
1084 WRITE(IOUT,2500) IFORM, ILVOUT, IPRES, KFORM, FREESURF, AFTERFLOW, INTEGR
1085 IF(IPRES == 1)THEN
1086 WRITE(IOUT,2600) PMAX, THETA, APMAX, ATHETA
1087 ELSEIF(IPRES == 2) THEN
1088 WRITE(IOUT,2700) IFPRES,SFPRES
1089 ENDIF
1090 IF(IWAVE == 1) WRITE(IOUT,3000) XC, YC, ZC
1091.OR. IF(GRAV_ID > 0 FREESURF == 2) WRITE(IOUT,3500) XA,YA,ZA,DIRX,DIRY,DIRZ
1092 IF(FREESURF == 2) WRITE(IOUT,3600) XD,YD,ZD
1093
1094C-----------------------------------------------------------------------------------
1095C Compute arrival time of incident wave, area, direction cosine, distance to charge
1096C Compute fluid mass matrix
1097C-----------------------------------------------------------------------------------
1098C description of arrays
1099C II1 -> IFlow Niflow: entiers parameters
1100C II2 -> IBUF NNO : local-global correspondences surface nodes
1101C II3 -> Elem 3 or 5*nel: local connectivity of triangles/quadrangles+flag
1102C II4 -> IBUF_L NNO : local-global correspondences surface nodes
1103C II5 -> SHELL_GA NBGAUGE : local-gauge correspondences 4-node shell
1104C II6 -> CNP NNO: Number of processor for each NSPMD node> 1
1105C engine
1106C II7 -> IPIV NEL : lapack resolution nel > nelmax
1107C ii8 -> ibufelr nel : spmd row of elements in the process grid (0 if absent from current processor)
1108C ii9 -> ibufelc nel : spmd column of elements in the process grid (0 if absent from current processor)
1109C
1110C IR1 -> RFLOW : parametres reels
1111C IR2 -> NORMAL : normale
1112C IR3 -> TA : arrival time
1113C IR4 -> AREA : element area
1114C IR5 -> COSG : direction cosine
1115C IR6 -> DIST : distance charge
1116C IR7 -> MFLE : fluide mass matrix (inverse) if nel < nelmax
1117C IR7 -> MFLE : C**t B + B**t C if nel >= nelmax
1118C engine
1119C IR8 -> ACCF : acceleration point fluide
1120C IR9 -> PS : scattered pressure
1121C IR10 -> PTI : incident pressure time integral
1122C IR11 -> CBEM : C matrix if nel >= nelmax
1123C-----------------------------------------------------------------------------------
1124 II1=IADI+1
1125 II2=II1+NIFLOW
1126 II3=II2+NNO
1127 IF(JFORM==1) THEN
1128 II4=II3+3*NEL
1129 II5=II4+NNO
1130 II6=II5
1131 ELSEIF(JFORM==2) THEN
1132 II4=II3+5*NEL
1133 II5=II4+NNO
1134 II6=II5+NBGAUGE
1135 ENDIF
1136 IR1=IADR+1
1137 IR2=IR1+NRFLOW
1138 IR3=IR2+NEL*3
1139 IR4=IR3+NEL*FREESURF
1140 IR5=IR4+NEL
1141 IR6=IR5+NEL*FREESURF
1142 IR7=IR6+NEL*FREESURF
1143 HG = HUGE(NEL)
1144 IF (NEL > INT(SQRT(REAL(HG)))) THEN
1145 CALL ANCMSG(MSGID = 1711, ANMODE=ANINFO, MSGTYPE = MSGERROR,
1146 . I1 = INT(SQRT(REAL(HG))))
1147 CALL ARRET(2)
1148 ENDIF
1149 IR8=IR7+NEL*NEL
1150 IR9=IR8+NEL
1151 IR10=IR9+NEL
1152 IR11=IR10+NEL
1153 IF(JFORM == 1) THEN
1154 CALL INIT_TG(IFLOW(II1), IFLOW(II2), IFLOW(II3), X, XS, YS, ZS, XD, YD, ZD,
1155 . RFLOW(IR1), RFLOW(IR2), RFLOW(IR3), RFLOW(IR4), RFLOW(IR5), RFLOW(IR6))
1156 CALL MASS_FLUID_TG(IFORM, ILVOUT, NNO, NEL, IFLOW(II2), IFLOW(II3), X,
1157 . RFLOW(IR2), RFLOW(IR4), RFLOW(IR7), RHO)
1158 ELSEIF(JFORM == 2) THEN
1159 CALL INIT_QD(IFLOW(II1), IFLOW(II2), IFLOW(II3), X, XS, YS, ZS, XD, YD, ZD,
1160 . RFLOW(IR1), RFLOW(IR2), RFLOW(IR3), RFLOW(IR4), RFLOW(IR5), RFLOW(IR6))
1161 IF(NEL < NELMAX) THEN
1162 ALLOCATE(CBEM(NEL,NEL))
1163 CALL MASS_FLUID_QD(NNO, NEL, IFLOW(II1), IFLOW(II2), IFLOW(II3), X,
1164 . RFLOW(IR2), RFLOW(IR4), RFLOW(IR7), CBEM, RHO,IRESP)
1165 DEALLOCATE(CBEM)
1166 ELSE
1167 CALL MASS_FLUID_QD(NNO, NEL, IFLOW(II1), IFLOW(II2), IFLOW(II3), X,
1168 . RFLOW(IR2), RFLOW(IR4), RFLOW(IR7), RFLOW(IR11), RHO,IRESP)
1169 ENDIF
1170 IF(NBGAUGE > 0) THEN
1171 WRITE (IOUT,'(/5x,a)') 'gauge element element'
1172 DO J=1,NBGAUGE
1173 IFLOW(II5+J-1)=0
1174 IF(LGAUGE(1,J) /= 0) CYCLE
1175C Node or Shell input
1176 N1=LGAUGE(3,J)
1177 IF(N1 < 0) THEN
1178 J1=-N1
1179 DO K=1,NEL
1180 IF(J1 /= N_SHELL(K)) CYCLE
1181 IFLOW(II5+J-1)=K
1182 GO TO 100
1183 ENDDO
1184 ENDIF
1185 100 CONTINUE
1186 WRITE(IOUT,'(3i10)') J,-LGAUGE(3,J),IFLOW(II5+J-1)
1187 ENDDO
1188 DEALLOCATE(N_SHELL)
1189 ENDIF
1190 ENDIF
1191C
1192 IADR=IADR+IFLOW(IADI+15)
1193 IADI=IADI+IFLOW(IADI+14)
1194C
1195 ENDDO
1196C
1197 RETURN
1198C
1199 1000 FORMAT(/
1200 . /
1201 . ' incompressible flow(boundary elements method)'/
1202 . ' ----------------------------------------------'/)
1203 1100 FORMAT( 5X,'bem problem number ',I10
1204 . /10X,'flow id ',I10,1X,A,
1205 . /10X,'EXTERNAL surface id ',I10
1206 . /10X,'number of surface nodes ',I10
1207 . /10X,'number of triangular boundary elements ',I10)
1208 1110 FORMAT( 10X,'flow inside the surface')
1209 1120 FORMAT( 10X,'flow outside the surface')
1210 1200 FORMAT( 10X,'stagnation pressure curve ',I10
1211 . /10X,'stagnation pressure scale factor ',1PE10.3
1212 . /10X,'time scale factor for stag. pres. curve ',1PE10.3
1213 . /10X,'auxiliary node group id ',I10
1214 . /10X,'number of auxiliary nodes ',I10)
1215 1210 FORMAT( 10X,'point-inside-surface test for aux. nodes')
1216 1220 FORMAT( 10X,'point-outside-surface test for aux. nodes')
1217 1225 FORMAT( 10X,'adimensional tolerance for testing ',1PE10.3)
1218 1230 FORMAT( 10X,'no test for aux. nodes')
1219 1300 FORMAT( 10X,'fluid density ',1PE10.3)
1220 1400 FORMAT(/10X,'inflow-outflow'
1221 . /10X,'--------------'
1222 . /10X,'number of inflow-outflow surfaces ',I10)
1223 1410 FORMAT(/10X,'surface number ',I10
1224 . /10X,'surface id ',I10)
1225 1420 FORMAT( 10X,'imposed velocity curve ',I10
1226 . /10X,'imposed velocity scale factor ',1PE10.3)
1227 1430 FORMAT( 10X,'imposed pressure curve ',I10
1228 . /10X,'imposed pressure scale factor ',1PE10.3)
1229 1440 FORMAT( 10X,'time scale factor for curves ',1PE10.3)
1230 1500 FORMAT(/10X,'bem parameters'
1231 . /10X,'--------------'
1232 . /10X,'bem formulation flag ',I10
1233 . /10X,'bem solver output level ',I10
1234 . /10X,'time step for matrices assembly ',1PE10.3)
1235 1600 FORMAT(/10X,'velocity field at infinity'
1236 . /10X,'--------------------------'
1237 . /10X,'velocity curve ',I10
1238 . /10X,'velocity scale factor ',1PE10.3
1239 . /10X,'time scale for velocity curve ',1PE10.3
1240 . /10X,'x component of velocity vector ',1PE10.3
1241 . /10X,'y component of velocity vector ',1PE10.3
1242 . /10X,'z component of velocity vector ',1PE10.3)
1243 1700 FORMAT(/10X,'parallel solver parameters(scalapack)'
1244 . /10X,'--------------------------------------'
1245 . /10X,'number of row of process grid ',I10
1246 . /10X,'number of columns of process grid ',I10
1247 . /10X,'2d-cyclic decomposition block-SIZE ',I10)
1248
1249 2000 FORMAT(//
1250 . ' daa surface(boundary element method) '/
1251 . ' -------------------------------------- '/)
1252 2100 FORMAT( 5X,'daa surface number ',I10
1253 . /10X,'daa id ',I10,1X,A,
1254 . /10X,'wet surface id ',I10
1255 . /10X,'number of surface nodes ',I10
1256 . /10X,'number of triangular elements ',I10
1257 . /10X,'gravity id(/grav) ',I10)
1258 2200 FORMAT( 5X,'daa surface number ',I10
1259 . /10X,'daa id ',I10,1X,A,
1260 . /10X,'wet surface id ',I10
1261 . /10X,'number of surface nodes ',I10
1262 . /10X,'number of shell elements ',I10
1263 . /10X,'gravity id(/grav) ',I10)
1264 2300 FORMAT( 10X,'fluid density ',1PE13.6
1265 . /10X,'fluid sound speed ',1PE13.6
1266 . /10X,'minimum pressure ',1PE13.6)
1267 2400 FORMAT( 10X,'x-coordinate of standoff point ',1PE13.6
1268 . /10X,'y-coordinate of standoff point ',1PE13.6
1269 . /10X,'z-coordinate of standoff point ',1PE13.6)
1270 2500 FORMAT(/10X,'bem formulation flag iform ',I10
1271 . /10X,'daa solver output level ',I10
1272 . /10X,'incident pressure input flag ',I10
1273 . /10X,'daa formulation flag kform ',I10
1274 . /10X,'free surface flag ',I10
1275 . /10X,'afterflow velocity flag ',I10
1276 . /10X,'integration flag ',I10)
1277 2600 FORMAT(/10X,'maximum pressure at standoff point ',1PE13.6
1278 . /10X,'decay time at standoff point ',1PE13.6
1279 . /10X,'exponent for pmax(apmax) ',1PE13.6
1280 . /10X,'exponent for decay time(atheta) ',1PE13.6)
1281 2700 FORMAT(/10X,'incident pressure FUNCTION ',I10
1282 . /10X,'pressure scale factor ',1PE13.6)
1283 3000 FORMAT( 10X,'x-coordinate of explosive charge ',1PE13.6
1284 . /10X,'y-coordinate of explosive charge ',1PE13.6
1285 . /10X,'z-coordinate of explosive charge ',1PE13.6)
1286 3100 FORMAT(/10X,'plane wave direction '
1287 . /10X,'x-direction ',1PE13.6
1288 . /10X,'y-direction ',1PE13.6
1289 . /10X,'z-direction ',1PE13.6)
1290 3500 FORMAT(/10X,'free surface '
1291 . /10X,'x-coordinate of surface point a ',1PE13.6
1292 . /10X,'y-coordinate of surface point a ',1PE13.6
1293 . /10X,'z-coordinate of surface point a ',1PE13.6
1294 . /10X,'surface normal x-component ',1PE13.6
1295 . /10X,'surface normal y-component ',1PE13.6
1296 . /10X,'surface normal z-component ',1PE13.6)
1297 3600 FORMAT(/10X,'x-coordinate of charge image ',1PE13.6
1298 . /10X,'y-coordinate of charge image ',1PE13.6
1299 . /10X,'z-coordinate of charge image ',1PE13.6)
1300 END
1301
my_real
#define my_real
Definition cppsort.cpp:32
the
end diagonal values have been computed in the(sparse) matrix id.SOL
norm
norm(diag(diag(diag(inv(mat))) -id.SOL), 2) % destroy mumps instance id.JOB
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:528
hm_option_count
subroutine hm_option_count(entity_type, hm_option_number)
Definition hm_option_count.F:105
hm_option_start
subroutine hm_option_start(entity_type)
Definition hm_option_start.F:253
hm_read_bem
subroutine hm_read_bem(igrsurf, iflow, rflow, npc, igrnod, memflow, unitab, x, nom_opt, lgauge, igrv, lsubmodel, iresp)
Definition hm_read_bem.F:274
hm_preread_bem
subroutine hm_preread_bem(igrsurf, igrnod, nnft, unitab, nom_opt, lsubmodel)
Definition hm_read_bem.F:41
hm_option_read_mod::hm_option_read_key
Definition hm_option_read_mod.F:414
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
groupdef_mod
Definition groupdef_mod.F:668
hm_option_read_mod
Definition hm_option_read_mod.F:412
message_mod
Definition message_mod.F:1257
names_and_titles_mod
Definition names_and_titles_mod.F:1002
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1008
names_and_titles_mod::ncharkey
integer, parameter ncharkey
Definition names_and_titles_mod.F:1005
submodel_mod
Definition submodel_mod.F:717
unitab_mod
Definition unitab_mod.F:517
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:895
arret
subroutine arret(nn)
Definition arret.F:86
starter
program starter
Definition starter.F:39
submodel_mod::submodel_data
Definition submodel_mod.F:728
velocity
subroutine velocity(a, ar, v, vr, fzero, itab, nale)
Definition velocity.F:29