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