OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_monvol_type3.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_read_monvol_type3_mod ../starter/source/airbag/hm_read_monvol_type3.F
25!||--- called by ------------------------------------------------------
26!|| read_monvol ../starter/source/airbag/read_monvol.F
27!||====================================================================
28 MODULE hm_read_monvol_type3_mod
29 CONTAINS
30!||====================================================================
31!|| hm_read_monvol_type3 ../starter/source/airbag/hm_read_monvol_type3.F
32!||--- called by ------------------------------------------------------
33!|| read_monvol ../starter/source/airbag/read_monvol.F
34!||--- calls -----------------------------------------------------
35!|| ancmsg ../starter/source/output/message/message.F
36!|| freerr ../starter/source/starter/freform.F
37!|| hm_get_float_array_index ../starter/source/devtools/hm_reader/hm_get_float_array_index.F
38!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
39!|| hm_get_floatv_dim ../starter/source/devtools/hm_reader/hm_get_floatv_dim.F
40!|| hm_get_int_array_index ../starter/source/devtools/hm_reader/hm_get_int_array_index.F
41!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
42!|| monvol_check_surfclose ../starter/share/modules1/monvol_struct_mod.F
43!|| monvol_compute_volume ../starter/share/modules1/monvol_struct_mod.F
44!|| monvol_orient_surf ../starter/share/modules1/monvol_struct_mod.F
45!|| monvol_reverse_normals ../starter/share/modules1/monvol_struct_mod.F
46!||--- uses -----------------------------------------------------
47!|| message_mod ../starter/share/message_module/message_mod.F
48!|| monvol_struct_mod ../starter/share/modules1/monvol_struct_mod.F
49!|| submodel_mod ../starter/share/modules1/submodel_mod.F
50!||====================================================================
51 SUBROUTINE hm_read_monvol_type3(T_MONVOLN,
52 . UNITAB, LUID, NPC, IGRSURF,
53 . ITAB, X, PM, GEO, IXC, IXTG, LSUBMODEL)
54C-----------------------------------------------
55C M o d u l e s
56C-----------------------------------------------
57 USE unitab_mod
58 USE groupdef_mod
59 USE message_mod
60 USE monvol_struct_mod
61 USE submodel_mod
62 use element_mod , only : nixc,nixtg
63C-----------------------------------------------
64C I m p l i c i t T y p e s
65C-----------------------------------------------
66#include "implicit_f.inc"
67C-----------------------------------------------
68C C o m m o n B l o c k s
69C-----------------------------------------------
70C NSURF, NFUNCT
71#include "com04_c.inc"
72C KMONVO, IREC
73C NIMV, NRVOLU
74#include "param_c.inc"
75C IIN
76#include "units_c.inc"
77C-----------------------------------------------
78C D u m m y A r g u m e n t s
79C-----------------------------------------------
80 TYPE(unit_type_), INTENT(IN) :: UNITAB
81 INTEGER, INTENT(IN) :: LUID
82 INTEGER, INTENT(IN) :: NPC(*), ITAB(*), IXC(NIXC, *), IXTG(NIXTG, *)
83 my_real, INTENT(IN) :: x(3, *), geo(npropg, *), pm(npropm, *)
84 TYPE (SURF_), INTENT(INOUT), DIMENSION(NSURF) :: IGRSURF
85 TYPE(monvol_struct_), INTENT(INOUT) :: T_MONVOLN
86 TYPE(submodel_data), DIMENSION(NSUBMOD), INTENT(IN) :: LSUBMODEL
87C-----------------------------------------------
88C L o c a l V a r i a b l e s
89C-----------------------------------------------
90 INTEGER :: II, J
91 INTEGER :: SURFID
92 my_real :: fac_m, fac_l, fac_t, fac_c, fac_gen
93 my_real :: scal_t, scal_p, scal_s, scal_a, scal_d
94 LOGICAL :: FOUND
95 my_real :: sa, rot, vol, vmin, veps, amu, sv, rhoi
96 my_real :: gamma, t_relax, pext, pini, pmax, vinc, mi, ei, avent
97 my_real :: bvent, tvent, dpdef, dtpdef, fport, fporp, fpora, ti
98 INTEGER :: NVENT, IPVENT, IDEF, IPORT, IPORP, IPORA, IEQUI
99 LOGICAL :: IS_AVAILABLE
100C-----------------------------------------------
101C B e g i n n i n g o f s o u r c e
102C-----------------------------------------------
103C ==================
104C Conversion factors
105C ==================
106 fac_m = unitab%FAC_M(luid)
107 fac_l = unitab%FAC_L(luid)
108 fac_t = unitab%FAC_T(luid)
109 fac_c = fac_m / (fac_l * fac_t * fac_t)
110C =======
111C Reading
112C =======
113! Line 1
114 CALL hm_get_intv('surf_IDex', surfid, is_available, lsubmodel)
115 CALL hm_get_intv('I_equi', iequi, is_available, lsubmodel)
116
117! line 2
118 CALL hm_get_floatv('Ascalet', scal_t, is_available, lsubmodel, unitab)
119 CALL hm_get_floatv('AscaleP', scal_p, is_available, lsubmodel, unitab)
120 CALL hm_get_floatv('AscaleS', scal_s, is_available, lsubmodel, unitab)
121 CALL hm_get_floatv('AscaleA', scal_a, is_available, lsubmodel, unitab)
122 CALL hm_get_floatv('AscaleD', scal_d, is_available, lsubmodel, unitab)
123! Line 3
124 CALL hm_get_floatv('Gamma', gamma, is_available, lsubmodel, unitab)
125 CALL hm_get_floatv('Mu', amu, is_available, lsubmodel, unitab)
126 CALL hm_get_floatv('Trelax', t_relax, is_available, lsubmodel, unitab)
127 CALL hm_get_floatv('Tini', ti, is_available, lsubmodel, unitab)
128 CALL hm_get_floatv('Rhoi', rhoi, is_available, lsubmodel, unitab)
129! Line 4
130 CALL hm_get_floatv('Pext', pext, is_available, lsubmodel, unitab)
131 CALL hm_get_floatv('Pini', pini, is_available, lsubmodel, unitab)
132 CALL hm_get_floatv('Pmax', pmax, is_available, lsubmodel, unitab)
133 CALL hm_get_floatv('Vinc', vinc, is_available, lsubmodel, unitab)
134 CALL hm_get_floatv('Mini', mi, is_available, lsubmodel, unitab)
135! Ventholes
136 CALL hm_get_intv('Nvent', nvent, is_available, lsubmodel)
137 t_monvoln%IVOLU(11) = nvent
138 t_monvoln%NVENT = nvent
139 IF (nvent > 0) THEN
140 ALLOCATE(t_monvoln%IBAGHOL(nibhol, nvent))
141 t_monvoln%IBAGHOL(1:nibhol, 1:nvent) = 0
142 ALLOCATE(t_monvoln%RBAGHOL(nrbhol, nvent))
143 t_monvoln%RBAGHOL(1:nrbhol, 1:nvent) = zero
144 ENDIF
145C **********
146C Vent-holes
147C **********
148C NVENT dependant number of lines
149 DO ii = 1, nvent
150 CALL hm_get_int_array_index('surf_IDv', ipvent, ii,is_available, lsubmodel)
151 CALL hm_get_float_array_index('Avent', avent, ii, is_available, lsubmodel, unitab)
152 CALL hm_get_int_array_index('Ideleted', idef, ii,is_available, lsubmodel)
153 bvent = zero
154 CALL hm_get_float_array_index('Tvent', tvent, ii, is_available, lsubmodel, unitab)
155 CALL hm_get_float_array_index('DeltaPdef', dpdef, ii, is_available, lsubmodel, unitab)
156 CALL hm_get_float_array_index('DeltatPdef', dtpdef, ii, is_available, lsubmodel, unitab)
157
158 CALL hm_get_int_array_index('fct_IDt', iport, ii,is_available, lsubmodel)
159 CALL hm_get_int_array_index('fct_IDP', iporp, ii,is_available, lsubmodel)
160 CALL hm_get_int_array_index('fct_IDA', ipora, ii,is_available, lsubmodel)
161 CALL hm_get_float_array_index('Fscalet', fport, ii, is_available, lsubmodel, unitab)
162 CALL hm_get_float_array_index('FscaleP', fporp, ii, is_available, lsubmodel, unitab)
163 CALL hm_get_float_array_index('FscaleA', fpora, ii, is_available, lsubmodel, unitab)
164
165 IF (fport == zero) fport = one
166 IF (fporp == zero) fporp = one
167 IF (fpora == zero) fpora = one
168 IF(ipvent == 0 .OR. idef == 0)THEN
169 idef = -1
170 ELSE
171 idef = -2
172 ENDIF
173 t_monvoln%IBAGHOL(1, ii) = idef
174 IF(ipvent == 0) THEN
175 t_monvoln%IBAGHOL(2, ii) = 0
176 ELSE
177 t_monvoln%IBAGHOL(2, ii) = 0
178 DO j = 1, nsurf
179 IF (ipvent == igrsurf(j)%ID) t_monvoln%IBAGHOL(2, ii) = j
180 ENDDO
181 IF(avent == zero) avent = one
182 ENDIF
183 IF(avent == zero) dpdef = infinity
184 IF(avent == zero) tvent = infinity
185 IF(dpdef == zero .AND. dtpdef == zero .AND. tvent == zero) THEN
186 t_monvoln%IBAGHOL(1, ii) = -t_monvoln%IBAGHOL(1, ii)
187 ENDIF
188 t_monvoln%RBAGHOL(1, ii) = dpdef
189 t_monvoln%RBAGHOL(2, ii) = avent
190 t_monvoln%RBAGHOL(3, ii) = tvent
191 t_monvoln%RBAGHOL(4, ii) = dtpdef
192 t_monvoln%RBAGHOL(6, ii) = bvent
193 t_monvoln%RBAGHOL(7, ii) = fport
194 t_monvoln%RBAGHOL(8, ii) = fporp
195 t_monvoln%RBAGHOL(9, ii) = fpora
196C
197 t_monvoln%IBAGHOL(3, ii) = -1
198 t_monvoln%IBAGHOL(4, ii) = -1
199 t_monvoln%IBAGHOL(5, ii) = -1
200 DO j = 1, nfunct
201 IF(iport == npc(j)) t_monvoln%IBAGHOL(3, ii) = j
202 IF(iporp == npc(j)) t_monvoln%IBAGHOL(4, ii) = j
203 IF(ipora == npc(j)) t_monvoln%IBAGHOL(5, ii) = j
204 ENDDO
205 IF(iport == 0) t_monvoln%IBAGHOL(3, ii) = 0
206 IF(iporp == 0) t_monvoln%IBAGHOL(4, ii) = 0
207 IF(ipora == 0) t_monvoln%IBAGHOL(5, ii) = 0
208 t_monvoln%IBAGHOL(6, ii) = 0
209 t_monvoln%IBAGHOL(7, ii) = 0
210 t_monvoln%IBAGHOL(8, ii) = 0
211 IF(t_monvoln%IBAGHOL(3, ii) == -1)THEN
212 t_monvoln%IBAGHOL(3, ii) = 0
213 CALL ancmsg(msgid = 331, anmode = aninfo, msgtype = msgerror,
214 . i1 = t_monvoln%ID, c1 = t_monvoln%TITLE, i2 = iport)
215 ENDIF
216 IF(t_monvoln%IBAGHOL(4, ii) == -1)THEN
217 t_monvoln%IBAGHOL(4, ii) = 0
218 CALL ancmsg(msgid = 332, anmode = aninfo, msgtype = msgerror,
219 . i1 = t_monvoln%ID, c1 = t_monvoln%TITLE, i2 = iporp)
220 ENDIF
221 IF(t_monvoln%IBAGHOL(5, ii) == -1)THEN
222 t_monvoln%IBAGHOL(5, ii) = 0
223 CALL ancmsg(msgid = 333, anmode = aninfo, msgtype = msgerror,
224 . i1 = t_monvoln%ID, c1 = t_monvoln%TITLE, i2 = ipora)
225 ENDIF
226 ENDDO
227C ================
228C Check operations
229C ================
230C External surface check
231 t_monvoln%IVOLU(4) = 0
232 found = .false.
233 DO ii = 1, nsurf
234 IF (surfid == igrsurf(ii)%ID) THEN
235 t_monvoln%IVOLU(4) = ii
236 t_monvoln%EXT_SURFID = ii
237 found = .true.
238 EXIT
239 ENDIF
240 ENDDO
241 IF (.NOT. found) THEN
242 CALL freerr(3)
243 ELSEIF (igrsurf(t_monvoln%EXT_SURFID)%ISH4N3N == 0) THEN
244 CALL ancmsg(msgid = 18, anmode = aninfo, msgtype = msgerror,
245 . i1 = t_monvoln%ID, c1 = t_monvoln%TITLE, i2 = surfid)
246 CALL freerr(3)
247 ENDIF
248
249C Check surface closure
250 CALL monvol_check_surfclose(t_monvoln, itab, igrsurf(t_monvoln%EXT_SURFID), x)
251C Set all normal on same side
252 CALL monvol_orient_surf(t_monvoln, t_monvoln%TITLE, t_monvoln%IVOLU, itab,
253 . igrsurf(t_monvoln%EXT_SURFID),ixc, ixtg, x, 3)
254C Compute Monvon volume
255 CALL monvol_compute_volume(t_monvoln, t_monvoln%TITLE, t_monvoln%IVOLU, igrsurf(t_monvoln%EXT_SURFID),
256 . itab, x, pm, geo, ixc, ixtg,
257 . sa, rot, vol, vmin, veps, sv)
258C Reverse all normals to ensure positive volume
259 CALL monvol_reverse_normals(t_monvoln, t_monvoln%TITLE, t_monvoln%IVOLU, itab,
260 . igrsurf(t_monvoln%EXT_SURFID),ixc,ixtg,vol, x, 3)
261
262C Parameter checks
263 IF (gamma == one) THEN
264 CALL ancmsg(msgid = 641, msgtype = msgerror, anmode = aninfo_blind_1, i1 = t_monvoln%ID, c1 = t_monvoln%TITLE)
265 ENDIF
266 IF (amu == zero) THEN
267 amu = em02
268 ENDIF
269
270C =====
271C Units
272C =====
273C Default value for time scale factor
274 IF (scal_t == zero) THEN
275 CALL hm_get_floatv_dim('Ascalet', fac_gen, is_available, lsubmodel, unitab)
276 scal_t = one * fac_gen
277 ENDIF
278 IF (scal_p == zero) THEN
279 CALL hm_get_floatv_dim('AscaleP', fac_gen, is_available, lsubmodel, unitab)
280 scal_p = one * fac_gen
281 ENDIF
282 IF (scal_s == zero) THEN
283 CALL hm_get_floatv_dim('AscaleS', fac_gen, is_available, lsubmodel, unitab)
284 scal_s = one * fac_gen
285 ENDIF
286 IF (scal_a == zero) THEN
287 CALL hm_get_floatv_dim('AscaleA', fac_gen, is_available, lsubmodel, unitab)
288 scal_a = one * fac_gen
289 ENDIF
290 IF (scal_d == zero) THEN
291 CALL hm_get_floatv_dim('AscaleD', fac_gen, is_available, lsubmodel, unitab)
292 scal_d = one * fac_gen
293 ENDIF
294 IF (pmax == zero) THEN
295 CALL hm_get_floatv_dim('Pmax', fac_gen, is_available, lsubmodel, unitab)
296 pmax = infinity * fac_gen
297 ENDIF
298! Initial temperature
299 IF (ti == zero) ti = twohundred95
300 IF (iequi > 0) THEN
301 IF (rhoi /= zero .AND. mi /= zero) THEN
302 mi = rhoi / vol
303 CALL ancmsg(msgid = 1718, anmode = aninfo, msgtype = msgwarning,
304 . i1 = t_monvoln%ID, c1 = t_monvoln%TITLE)
305 ENDIF
306 IF (rhoi == zero .AND. mi == zero) THEN
307 CALL ancmsg(msgid = 1717, anmode = aninfo, msgtype = msgerror,
308 . i1 = t_monvoln%ID, c1 = t_monvoln%TITLE)
309 ENDIF
310 IF (rhoi /= zero .AND. mi == zero) THEN
311 mi = rhoi * vol
312 ENDIF
313 IF (mi /= zero .AND. rhoi == zero) THEN
314 rhoi = mi / vol
315 ENDIF
316 ENDIF
317C =====
318C Store
319C =====
320C
321! Store CV
322 t_monvoln%RVOLU(19) = zero
323 IF (rhoi /= zero) THEN
324 t_monvoln%RVOLU(19) = pext / ((gamma - one) * rhoi * ti)
325 ENDIF
326 t_monvoln%IVOLU(15) = iequi
327 t_monvoln%IVOLU(11) = nvent
328 t_monvoln%RVOLU(1) = gamma
329 t_monvoln%RVOLU(3) = pext
330 veps = max(zero, vmin - vol + vinc)
331 t_monvoln%RVOLU(4) = vol + veps
332 t_monvoln%RVOLU(5) = vinc
333 t_monvoln%RVOLU(6) = pmax
334 t_monvoln%RVOLU(17) = veps
335 IF (t_relax > zero) THEN
336 ei = pext * (vol + veps - vinc) / (gamma - one)
337 t_monvoln%RVOLU(12) = pext
338 t_monvoln%RVOLU(52) = (pini - pext) * (vol + veps - vinc) / t_relax / (gamma-one)
339 ELSE
340 ei = pini * (vol + veps - vinc) / (gamma - one)
341 t_monvoln%RVOLU(12) = pini
342 t_monvoln%RVOLU(52) = zero
343 ENDIF
344 t_monvoln%RVOLU(56) = pini
345 t_monvoln%RVOLU(13) = ei
346
347 t_monvoln%RVOLU(20) = mi
348 t_monvoln%RVOLU(57) = mi
349 t_monvoln%RVOLU(48) = t_relax
350
351C Store in data structure
352 t_monvoln%RVOLU(26) = one / scal_t
353 t_monvoln%RVOLU(27) = one / scal_p
354 t_monvoln%RVOLU(28) = one / scal_s
355 t_monvoln%RVOLU(29) = one / scal_a
356 t_monvoln%RVOLU(30) = one / scal_d
357C
358
359C
360 t_monvoln%RVOLU(2) = amu
361 t_monvoln%RVOLU(16) = vol + veps
362 t_monvoln%RVOLU(18) = sa
363 t_monvoln%RVOLU(21) = rot
364 t_monvoln%RVOLU(22:24) = zero
365
366C =========
367C Print out
368C =========
369 WRITE(iout, 1005) surfid
370 WRITE(iout, 1003) scal_t, scal_p, scal_s, scal_a, scal_d
371 WRITE(iout, 1002) sa, sv, vol
372 WRITE(iout, 1300) gamma, amu, pext, pini, pmax, vinc, mi, t_relax, rhoi, ti
373 IF (iequi == 1) THEN
374 WRITE(iout, 1301)
375 ELSE IF (iequi == 2) THEN
376 WRITE(iout, 1302)
377 ENDIF
378 WRITE(iout,1370)
379 WRITE(iout,1371) nvent
380 DO ii = 1, nvent
381 idef = t_monvoln%IBAGHOL(1, ii)
382 IF (ipvent /= 0) ipvent = igrsurf(t_monvoln%IBAGHOL(2, ii))%ID
383 iport = t_monvoln%IBAGHOL(3, ii)
384 iporp = t_monvoln%IBAGHOL(4, ii)
385 ipora = t_monvoln%IBAGHOL(5, ii)
386 IF(iport /= 0) iport = npc(iport)
387 IF(iporp /= 0) iporp = npc(iporp)
388 IF(ipora /= 0) ipora = npc(ipora)
389 avent = t_monvoln%RBAGHOL(2, ii)
390 tvent = t_monvoln%RBAGHOL(3, ii)
391 dpdef = t_monvoln%RBAGHOL(1, ii)
392 IF(abs(idef) == 1)THEN
393 WRITE(iout,1380) ii, ipvent, avent, iport, iporp, ipora,tvent, dpdef, dtpdef, fport, fporp, fpora
394 ELSE
395 WRITE(iout,1381) ii, ipvent, avent, iport, iporp, ipora,tvent, dpdef, dtpdef, fport, fporp, fpora
396 ENDIF
397 ENDDO
398C-----------------------------------------------
399C E n d o f s o u r c e
400C-----------------------------------------------
401
402 RETURN
403 1002 FORMAT(
404 . /5x,'INITIAL SURFACE OF MONITORED VOLUME . .=',1pg20.13,
405 . /5x,'SURFACE ERROR(NE.0 FOR NON CLOSED SURF)=',1pg20.13,
406 . /5x,'INITIAL VOLUME OF MONITORED VOLUME. . .=',1pg20.13)
407 1003 FORMAT(
408 . 5x,'UNIT SCALE FOR TIME FUNCTIONS =',1pg20.13,
409 . /5x,'UNIT SCALE FOR PRESSURE FUNCTIONS =',1pg20.13,
410 . /5x,'UNIT SCALE FOR AREA FUNCTIONS =',1pg20.13,
411 . /5x,'UNIT SCALE FOR ANGLE FUNCTIONS =',1pg20.13,
412 . /5x,'UNIT SCALE FOR DISTANCE FUNCTIONS =',1pg20.13)
413 1005 FORMAT( 5x,'EXTERNAL SURFACE ID . . . . . . . . . .=',i10)
414 1300 FORMAT(
415 . 5x,'GAMMA . . . . . . . . . . . . . . . . .=',1pg20.13,
416 . /5x,'VOLUMIC VISCOSITY . . . . . . . . . . .=',1pg20.13,
417 . /5x,'EXTERNAL PRESSURE . . . . . . . . . . .=',1pg20.13,
418 . /5x,'INITIAL PRESSURE. . . . . . . . . . . .=',1pg20.13,
419 . /5x,'BURSTING PRESSURE . . . . . . . . . . .=',1pg20.13,
420 . /5x,'INCOMPRESSIBLE VOLUME.. . . . . . . . .=',1pg20.13,
421 . /5x,'INITIAL MASS. . . . . . . . . . . . . .=',1pg20.13,
422 . /5x,'RELAXATION TIME . . . . . . . . . . . .=',1pg20.13,
423 . /5x,'INITIAL MASS DENSITY. . . . . . . . . .=',1pg20.13,
424 . /5x,'INITIAL TEMPERATURE . . . . . . . . . .=',1pg20.13)
425 1370 FORMAT(
426 . /5x,'VENT HOLES ',
427 . /5x,'---------- ')
428 1301 FORMAT(
429 . /5x,'ISOTHERMAL EQUILIBRIUM . . . . . . . . .')
430 1302 FORMAT(
431 . /5x,'ISENTROPIC EQUILIBRIUM . . . . . . . . .')
432 1371 FORMAT(
433 . 5x,'NUMBER OF VENT HOLES. . . . . . . . . .=',i10/)
434 1380 FORMAT(
435 . 5x,'VENT HOLE NUMBER . . . . . . . . . . . =',i10,
436 . /15x,'VENT HOLE SURFACE ID. . . . . . . . . .=',i10,
437 . /15x,' (TOTAL AREA IS CONSIDERED) ',
438 . /15x,'AVENT (VENT HOLE AREA OR SCALE FACTOR) =',1pg20.13,
439 . /15x,' IF IPDEF=0 : VENT HOLE AREA '
440 . /15x,' IF IPDEF/=0 : VENT HOLE AREA SCALE FACTOR'
441 . /15x,'POROSITY FUNCTION / TIME. . . . . . . .=',i10,
442 . /15x,'POROSITY FUNCTION / PRESSURE. . . . . .=',i10,
443 . /15x,'POROSITY FUNCTION / AREA. . . . . . . .=',i10,
444 . /15x,'START TIME FOR VENTING. . . . . . . . .=',1pg20.13,
445 . /15x,'RELATIVE PRES. FOR MEMBRANE DEFLATION .=',1pg20.13,
446 . /15x,' (DPDEF = PDEF - PEXT) ',
447 . /15x,'MIMINUM DURATION PRES. EXCEEDS PDEF ',
448 . /15x,' FOR MEMBRANE DEFLATION . . . .=',1pg20.13,
449 . /15x,'POROSITY TIME FUNCTION SCALE FACTOR =',1pg20.13,
450 . /15x,'POROSITY PRESSURE FUNCTION SCALE FACTOR=',1pg20.13,
451 . /15x,'POROSITY AREA FUNCTION SCALE FACTOR . .=',1pg20.13)
452 1381 FORMAT(
453 . 5x,'VENT HOLE NUMBER . . . . . . . . . . . =',i10,
454 . /15x,'VENT HOLE SURFACE ID. . . . . . . . . .=',i10,
455 . /15x,' (DELETED ELEMENTS AREA IS CONSIDERED) ',
456 . /15x,'AVENT (VENT HOLE AREA OR SCALE FACTOR) =',1pg20.13,
457 . /15x,' IF IPDEF=0 : VENT HOLE AREA '
458 . /15x,' IF IPDEF/=0 : VENT HOLE AREA SCALE FACTOR'
459 . /15x,'POROSITY FUNCTION / TIME. . . . . . . .=',i10,
460 . /15x,'POROSITY FUNCTION / PRESSURE. . . . . .=',i10,
461 . /15x,'POROSITY FUNCTION / AREA. . . . . . . .=',i10,
462 . /15x,'START TIME FOR VENTING. . . . . . . . .=',1pg20.13,
463 . /15x,'RELATIVE PRES. FOR MEMBRANE DEFLATION .=',1pg20.13,
464 . /15x,' (DPDEF = PDEF - PEXT) ',
465 . /15x,'MIMINUM DURATION PRES. EXCEEDS PDEF ',
466 . /15x,' FOR MEMBRANE DEFLATION . . . .=',1pg20.13,
467 . /15x,'POROSITY TIME FUNCTION SCALE FACTOR =',1pg20.13,
468 . /15x,'POROSITY PRESSURE FUNCTION SCALE FACTOR=',1pg20.13,
469 . /15x,'POROSITY AREA FUNCTION SCALE FACTOR . .=',1pg20.13)
470 END SUBROUTINE hm_read_monvol_type3
471 END MODULE hm_read_monvol_type3_mod
my_real
#define my_real
Definition cppsort.cpp:32
hm_get_float_array_index
subroutine hm_get_float_array_index(name, rval, index, is_available, lsubmodel, unitab)
Definition hm_get_float_array_index.F:105
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:418
hm_get_floatv_dim
subroutine hm_get_floatv_dim(name, dim_fac, is_available, lsubmodel, unitab)
Definition hm_get_floatv_dim.F:161
hm_get_int_array_index
subroutine hm_get_int_array_index(name, ival, index, is_available, lsubmodel)
Definition hm_get_int_array_index.F:164
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:528
max
#define max(a, b)
Definition macros.h:21
groupdef_mod
Definition groupdef_mod.F:668
hm_read_monvol_type3_mod
Definition hm_read_monvol_type3.F:28
hm_read_monvol_type3_mod::hm_read_monvol_type3
subroutine hm_read_monvol_type3(t_monvoln, unitab, luid, npc, igrsurf, itab, x, pm, geo, ixc, ixtg, lsubmodel)
Definition hm_read_monvol_type3.F:54
message_mod
Definition message_mod.F:1257
monvol_struct_mod
Definition monvol_struct_mod.F:66
monvol_struct_mod::monvol_compute_volume
subroutine monvol_compute_volume(t_monvoln, title, ivolu, surf, itab, node_coord, pm, geo, ixc, ixtg, sa, rot, vol, vmin, veps, sv)
Definition monvol_struct_mod.F:656
monvol_struct_mod::monvol_check_surfclose
subroutine monvol_check_surfclose(t_monvoln, itab, surf, x)
Definition monvol_struct_mod.F:421
submodel_mod
Definition submodel_mod.F:717
unitab_mod
Definition unitab_mod.F:517
monvol_orient_surf
subroutine monvol_orient_surf(t_monvoln, title, ivolu, itab, surf, ixc, ixtg, x, itype)
Definition monvol_struct_mod.F:1112
monvol_reverse_normals
subroutine monvol_reverse_normals(t_monvoln, title, ivolu, itab, surf, ixc, ixtg, vol, x, itype)
Definition monvol_struct_mod.F:1783
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
freerr
subroutine freerr(it)
Definition freform.F:501
monvol_struct_mod::monvol_struct_
Definition monvol_struct_mod.F:84
submodel_mod::submodel_data
Definition submodel_mod.F:728
unitab_mod::unit_type_
Definition unitab_mod.F:524