OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
airbaga1.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!|| airbaga1 ../engine/source/airbag/airbaga1.F
25!||--- called by ------------------------------------------------------
26!|| monvol0 ../engine/source/airbag/monvol0.F
27!||--- uses -----------------------------------------------------
28!|| finter_mixed_mod ../engine/source/tools/finter_mixed.F90
29!|| python_funct_mod ../common_source/modules/python_mod.F90
30!|| sensor_mod ../common_source/modules/sensor_mod.F90
31!||====================================================================
32 SUBROUTINE airbaga1(IVOLU ,NJET ,IBAGJET ,NVENT ,IBAGHOL ,
33 2 ICBAG ,RVOLU ,RBAGJET ,RBAGHOL ,RCBAG ,
34 3 RVOLUV ,NPC ,TF ,NSENSOR ,SENSOR_TAB ,
35 4 VOL ,PMAIN ,GEO ,IGEO ,PM ,
36 5 WFEXT, PYTHON )
37C-----------------------------------------------
38C M o d u l e s
39C-----------------------------------------------
40 USE python_funct_mod
41 USE sensor_mod
42 USE finter_mixed_mod
43C-----------------------------------------------
44C I m p l i c i t T y p e s
45C-----------------------------------------------
46#include "implicit_f.inc"
47C-----------------------------------------------
48C C o m m o n B l o c k s
49C-----------------------------------------------
50#include "com04_c.inc"
51#include "com06_c.inc"
52#include "com08_c.inc"
53#include "units_c.inc"
54#include "scr02_c.inc"
55#include "scr07_c.inc"
56#include "scr18_c.inc"
57#include "task_c.inc"
58#include "param_c.inc"
59C-----------------------------------------------
60C D u m m y A r g u m e n t s
61C-----------------------------------------------
62 INTEGER ,INTENT(IN) :: NSENSOR
63 INTEGER NPC(*),IVOLU(*),NJET,IBAGJET(NIBJET,*),
64 . NVENT,IBAGHOL(NIBHOL,*),ICBAG(NICBAG,*),PMAIN,
65 . IGEO(NPROPGI,*)
66 my_real
67 . tf(*),rvolu(*),rbagjet(nrbjet,*),
68 . rbaghol(nrbhol,*),rcbag(nrcbag,*),rvoluv(nrvolu,*),vol,
69 . geo(npropg,*),pm(npropm,*)
70 TYPE (SENSOR_STR_) ,DIMENSION(NSENSOR) ,INTENT(IN) :: SENSOR_TAB
71 DOUBLE PRECISION,INTENT(INOUT) :: WFEXT
72 TYPE(PYTHON_), intent(inout) :: PYTHON
73C-----------------------------------------------
74C L o c a l V a r i a b l e s
75C-----------------------------------------------
76 INTEGER I, IMASS, ITEMP, ISENS, IFLU, IEQUI,
77 . IDEF, NAV, II, IS, IDTPDEF,IDPDEF,
78 . IV, IJ,ITTF
79 INTEGER IK, KK, I_GAS, I_INJ, IMOL, NGASES, I_TYPINJ
80 my_real
81 . pdef, pext, pvois, dtpdefc,
82 . gama, amu, tstart, tsg,tsg2,
83 . amtot, p, pold, dv,
84 . rot, ftemp, fmass, dt, area,
85 . cv, cp, cvg, cpg, cpa, cpb, cpc, cvi, cpi, cpai, cpbi, cpci,
86 . rmwi, rmwg,
87 . rnm_old, rnm, rnmg,
88 . gmtot, gmtot_old, gmass, gmass_old, dgmass, dgmout,
89 . gmi, gmi_old,
90 . efac, dgeout, right, left,
91 . vold, veps, vmin, amtot_old, tbag_old, tbag,
92 . qold,q,cx,qx,temp
93 my_real
94 . mw, cpd, cpe, cpf, cpdi, cpei, cpfi, r_igc1, dgmin, astime,
95 . dgein, hconv, text, molfr, mw_mix
96C-----------------------------------------------
97 ftemp = zero
98 pext =rvolu(3)
99 pold =rvolu(12)
100 r_igc1 =pm(27,ivolu(66))
101 tbag_old = zero
102 fmass = zero
103 area = zero
104 dv = zero
105 i_gas = 0
106 qold = zero
107 vold = zero
108 rot = zero
109C------------------------
110 DO iv=1,nvent
111 idef = ibaghol(1,iv)
112 pdef = rbaghol(1,iv)
113 dtpdefc= rbaghol(5,iv)
114 idtpdef= ibaghol(11,iv)
115 idpdef = ibaghol(12,iv)
116C
117 IF (idtpdef==0) THEN
118 IF(idef==0.AND.pold>pdef+pext) THEN
119 dtpdefc = dtpdefc+dt1
120 END IF
121 ELSE IF (idtpdef==1) THEN
122 IF (pold>pdef+pext) THEN
123 idpdef = 1
124 END IF
125 IF (idpdef==1) THEN
126 dtpdefc = dtpdefc+dt1
127 END IF
128 END IF
129 ibaghol(12,iv) = idpdef
130 rbaghol(5,iv) = dtpdefc
131 ENDDO
132 IF(ispmd+1==pmain) THEN
133C Mass loss & Enthalpy Loss
134 DO iv=1,nvent
135 rbaghol(19,iv)= rbaghol(19,iv) + rbaghol(21,iv)*dt1
136 rbaghol(20,iv)= rbaghol(20,iv) + rbaghol(22,iv)*dt1
137 ENDDO
138C Heat Loss
139 rvolu(54)=rvolu(54)+rvolu(55)*dt1
140 END IF
141C-------------------------
142C AIRBAG COMMUNIQUANTS
143C-------------------------
144 nav = ivolu(3)
145 DO i=1,nav
146 ii = icbag(1,i)
147 idef = icbag(3,i)
148 pdef = rcbag(1,i)
149 dtpdefc= rcbag(5,i)
150c PVOIS =RVOLUV(12,II)
151c afin d'eviter probleme de meme adressage memoire
152c de RVOLU(12) et RVOLUV(12) on utilise RVOLU(31)
153c qui etait dispo
154 pvois =rvoluv(31,ii)
155 IF(idef==0.AND.pold>pdef+pvois)
156 . rcbag(5,i)=dtpdefc+dt1
157 ENDDO
158C---------------------
159 amu = rvolu(2)
160C
161 iequi=ivolu(15)
162C changement d'etat ?
163 is = 0
164 ittf=ivolu(17)
165 DO ij=1,njet
166 isens=ibagjet(4,ij)
167 IF(isens==0)THEN
168 tstart=zero
169 ELSE
170 tstart=sensor_tab(isens)%TSTART
171 ENDIF
172 IF(tt>=tstart)is=1
173 IF(iequi==0)THEN
174 dgmout=rbagjet(9,ij)
175 IF(dgmout/=0.)is=1
176 END IF
177 IF(is==1.AND.(ittf==1.OR.ittf==2.OR.ittf==3)) THEN
178 ittf=ittf+10
179C store in rvolu(60) airbag time to fire
180 rvolu(60)=tstart
181 ivolu(17)=ittf
182 END IF
183 ENDDO
184 IF(iequi==0)THEN
185 tbag_old =rvolu(13)
186 vold =rvolu(16)
187 veps =rvolu(17)
188 area =rvolu(18)
189 rot =rvolu(21)
190 qold =rvolu(23)
191 vol = vol + veps
192C
193 dgmout=rvolu(24)
194 IF(dgmout/=zero)is=1
195C
196 dv = vol-vold
197 IF(dv/=zero)is=1
198 ELSE
199 IF(is==1)THEN
200 tbag_old =rvolu(13)
201 area =rvolu(18)
202 IF(iequi==1)THEN
203 vmin =em4*area**three_half
204 veps =max(zero,vmin-vol)
205 rvolu(17)=veps
206 rmwi = rvolu(10)
207 gmi = pext*(vol+veps)/(rmwi*tbag_old)
208 rvolu(11)= gmi
209 rvolu(14)= rmwi*gmi
210 rvolu(20)= gmi
211 IF(ispmd+1==pmain) THEN
212 WRITE(iout,*)
213 . ' *** MONITORED VOLUME : INITIAL EQUILIBRIUM IS SET ***'
214 WRITE(iout,'(A,I10,A,G20.13,A)')
215 . ' *** MONITORED VOLUME ',ivolu(1),' VOLUME ',vol,' ***'
216 WRITE(istdo,*)
217 . ' *** MONITORED VOLUME : INITIAL EQUILIBRIUM IS SET ***'
218 ENDIF
219 ivolu(15)= -1
220C
221 vol =vol + veps
222 vold =vol
223 qold =zero
224 ELSE
225 vold =rvolu(16)
226 veps =rvolu(17)
227 qold =rvolu(23)
228 vol = vol + veps
229 ENDIF
230 rot =rvolu(21)
231 dv =vol-vold
232 ENDIF
233 ENDIF
234C--------------------------
235 dgmin =rvolu(64)*dt1
236 IF(dgmin > zero) is=1
237C
238 IF (is/=0)THEN
239C-----------------------------
240C MASSE MOLAIRE DU MELANGE
241C-----------------------------
242 DO ij=1,njet
243 i_inj = ibagjet(13,ij)
244 IF(i_inj <= 0) cycle
245 i_typinj = igeo(22,i_inj)
246 IF(i_typinj /= 2) cycle
247 isens=ibagjet(4,ij)
248 IF(isens==0)THEN
249 tstart=zero
250 ELSE
251 tstart=sensor_tab(isens)%TSTART
252 ENDIF
253 IF(tt < tstart) cycle
254C
255 astime = geo(201,i_inj)
256 tsg = (tt - tstart) / astime
257 ngases = igeo(23,i_inj)
258 mw_mix = zero
259C
260 DO ik=1,ngases
261 molfr = geo(211+(ik-1)+1,i_inj)
262 i_gas = igeo(100+(ik-1)*2+1,i_inj)
263 imol = igeo(100+(ik-1)*2+2,i_inj)
264 mw = pm(20,i_gas)
265 IF(imol == 0) THEN
266 mw_mix = mw_mix + mw*molfr
267 ELSE
268 mw_mix = mw_mix + mw*finter_mixed(python,nfunct,imol,tsg,npc,tf)
269 ENDIF
270 ENDDO
271 rbagjet(19,ij) = mw_mix
272 ENDDO
273C---------------------------------------------------------
274C FRACTION MOLAIRE
275C inclus caracteristiques du gaz a l'initial.
276C---------------------------------------------------------
277C RNM=R*NM=R*M/MW
278 rnm_old =rvolu(14)
279 amtot_old=rvolu(20)
280C
281 left =zero
282 right=zero
283 amtot=amtot_old
284C---------------------
285C INJECTEURS
286C---------------------
287 DO ij=1,njet
288 i_inj = ibagjet(13,ij)
289 IF(i_inj <= 0) cycle
290C
291 i_typinj = igeo(22,i_inj)
292 ngases = igeo(23,i_inj)
293 iflu = igeo(24,i_inj)
294 astime = geo(201,i_inj)
295C
296 isens=ibagjet(4,ij)
297 IF(isens==0)THEN
298 tstart=zero
299 ELSE
300 tstart=sensor_tab(isens)%TSTART
301 ENDIF
302C
303 DO ik=1,ngases
304 IF (i_typinj==1) THEN
305 fmass = geo(208+(ik-1)*2+1,i_inj)
306 ftemp = geo(208+(ik-1)*2+2,i_inj)
307 i_gas = igeo(100+(ik-1)*3+1,i_inj)
308 imass = igeo(100+(ik-1)*3+2,i_inj)
309 itemp = igeo(100+(ik-1)*3+3,i_inj)
310 ELSE IF (i_typinj==2) THEN
311 fmass = geo(209,i_inj)
312 ftemp = geo(210,i_inj)
313 molfr = geo(211+(ik-1)+1,i_inj)
314 imass = igeo(25,i_inj)
315 itemp = igeo(26,i_inj)
316 i_gas = igeo(100+(ik-1)*2+1,i_inj)
317 imol = igeo(100+(ik-1)*2+2,i_inj)
318 END IF
319C
320 mw =pm(20,i_gas)
321 cpa =pm(21,i_gas)
322 cpb =pm(22,i_gas)
323 cpc =pm(23,i_gas)
324 cpd =pm(24,i_gas)
325 cpe =pm(25,i_gas)
326 cpf =pm(26,i_gas)
327C
328 IF(tt>=tstart)THEN
329 tsg = (tt - tstart) / astime
330 tsg2 = (tt - dt1 - tstart) / astime
331 IF (tsg2 <= zero) tsg2 = tsg
332C
333 IF(itemp>0) THEN
334 temp=ftemp*finter_mixed(python,nfunct,itemp,tsg,npc,tf)
335 ELSE
336 temp=ftemp
337 ENDIF
338 efac= temp*(cpa+half*cpb*temp+third*cpc*temp*temp
339 . +fourth*cpd*temp*temp*temp
340 . -cpe/(temp*temp)
341 . +one_fifth*cpf*temp*temp*temp*temp)
342 ELSE
343 efac=zero
344 ENDIF
345C
346 kk=20+(ik-1)*4
347 gmtot_old=rbagjet(kk+1,ij)
348 gmass_old=rbagjet(kk+2,ij)
349 IF(tt>=tstart)THEN
350 IF(imass>0) THEN
351 gmass=half * fmass * (finter_mixed(python,nfunct,imass,tsg,npc,tf) +
352 . finter_mixed(python,nfunct,imass,tsg2,npc,tf))
353 IF(i_typinj == 2) THEN
354 mw_mix = rbagjet(19,ij)
355 IF(imol == 0) THEN
356 gmass=gmass*mw*molfr/mw_mix
357 ELSE
358 gmass=gmass*mw*finter_mixed(python,nfunct,imol,tsg,npc,tf)/mw_mix
359 ENDIF
360 ENDIF
361 IF(iflu==1)gmass = gmass*dt1 + gmass_old
362 dgmass=max(zero,gmass-gmass_old)
363 right =right+dgmass*efac
364 rbagjet(11,ij) = rbagjet(11,ij)+dgmass*efac
365 ELSE
366 gmass=zero
367 dgmass=zero
368 ENDIF
369 ELSE
370 dgmass=zero
371 gmass =zero
372 ENDIF
373C
374 rmwg = r_igc1/mw
375 gmtot = gmtot_old + dgmass
376 right = right -dgmass*tbag_old*
377 . (cpa+half*cpb*tbag_old
378 . +third*cpc*tbag_old*tbag_old
379 . +fourth*cpd*tbag_old*tbag_old*tbag_old
380 . -cpe/(tbag_old*tbag_old)
381 . +one_fifth*cpf*tbag_old*tbag_old*tbag_old*tbag_old
382 . -rmwg)
383 kk=20+(ik-1)*4
384 rbagjet(kk+1,ij) = gmtot
385 rbagjet(kk+2,ij) = gmass
386 amtot = amtot + dgmass
387 ENDDO !fin de boucle gaz
388 ENDDO !fin de boucle injecteurs
389C-------------------------
390C FUITES ET COMMUNICATIONS
391C-------------------------
392 DO ij=1,njet
393 i_inj = iabs(ibagjet(13,ij))
394 i_typinj = igeo(22,i_inj)
395 ngases = igeo(23,i_inj)
396 DO ik=1,ngases
397 IF (i_typinj==1) THEN
398 i_gas = igeo(100+(ik-1)*3+1,i_inj)
399 ELSE IF (i_typinj==2) THEN
400 i_gas = igeo(100+(ik-1)*2+1,i_inj)
401 END IF
402
403 mw =pm(20,i_gas)
404 cpa =pm(21,i_gas)
405 cpb =pm(22,i_gas)
406 cpc =pm(23,i_gas)
407 cpd =pm(24,i_gas)
408 cpe =pm(25,i_gas)
409 cpf =pm(26,i_gas)
410 rmwg = r_igc1/mw
411
412 kk=20+(ik-1)*4
413 dgmout = rbagjet(kk+3,ij)*dt1
414 dgeout = rbagjet(kk+4,ij)*dt1
415C
416 right = right -dgeout
417 right = right +dgmout*tbag_old*
418 . (cpa+half*cpb*tbag_old
419 . +third*cpc*tbag_old*tbag_old
420 . +fourth*cpd*tbag_old*tbag_old*tbag_old
421 . -cpe/(tbag_old*tbag_old)
422 . +one_fifth*cpf*tbag_old*tbag_old*tbag_old*tbag_old
423 . -rmwg)
424 cvg = cpa+cpb*tbag_old+cpc*tbag_old*tbag_old
425 . +cpd*tbag_old*tbag_old*tbag_old
426 . +cpe/(tbag_old*tbag_old)
427 . +cpf*tbag_old*tbag_old*tbag_old*tbag_old-rmwg
428
429 gmtot = rbagjet(kk+1,ij)-dgmout
430 left = left +gmtot*cvg
431 right = right+gmtot*cvg*tbag_old
432
433 rbagjet(kk+1,ij) = gmtot
434 amtot = amtot - dgmout
435 ENDDO !fin de boucle gaz
436 ENDDO !fin de boucle injecteurs
437C------------------------------------
438C CONTRIBUTION DU GAZ INITIAL...
439C------------------------------------
440 dgmout=rvolu(24)*dt1
441 dgmin =rvolu(64)*dt1
442C
443 cpai =rvolu(7)
444 cpbi =rvolu(8)
445 cpci =rvolu(9)
446 cpdi =rvolu(56)
447 cpei =rvolu(57)
448 cpfi =rvolu(58)
449 rmwi =rvolu(10)
450 gmi_old=rvolu(11)
451 cvi = cpai+cpbi*tbag_old+cpci*tbag_old*tbag_old
452 . +cpdi*tbag_old*tbag_old*tbag_old
453 . +cpei/(tbag_old*tbag_old)
454 . +cpfi*tbag_old*tbag_old*tbag_old*tbag_old-rmwi
455 gmi=gmi_old-dgmout+dgmin
456
457 dgeout=rvolu(22)*dt1
458 dgein =rvolu(65)*dt1
459C
460 right = right+(dgmout-dgmin)*tbag_old
461 . *(cpai+half*cpbi*tbag_old
462 . +third*cpci*tbag_old*tbag_old
463 . +fourth*cpdi*tbag_old*tbag_old*tbag_old
464 . -cpei/(tbag_old*tbag_old)
465 . +one_fifth*cpfi*tbag_old*tbag_old*tbag_old*tbag_old
466 . -rmwi)
467 right = right+dgein-dgeout
468 left = left+gmi*cvi
469 right = right+gmi*cvi*tbag_old
470C
471 amtot = amtot-dgmout+dgmin
472 rvolu(11)=gmi
473C------------------------------------------
474C VARIATION DE TEMPERATURE et PRESSION
475C------------------------------------------
476 gmi =rvolu(11)
477 rnm =gmi*rmwi
478 DO ij=1,njet
479 i_inj = iabs(ibagjet(13,ij))
480 i_typinj = igeo(22,i_inj)
481 ngases = igeo(23,i_inj)
482 DO ik=1,ngases
483 IF (i_typinj==1) THEN
484 i_gas = igeo(100+(ik-1)*3+1,i_inj)
485 ELSE IF (i_typinj==2) THEN
486 i_gas = igeo(100+(ik-1)*2+1,i_inj)
487 END IF
488 mw = pm(20,i_gas)
489 rmwg = r_igc1/mw
490 kk=20+(ik-1)*4
491 gmtot= rbagjet(kk+1,ij)
492 rnmg = gmtot*rmwg
493 rnm = rnm+rnmg
494 ENDDO
495 ENDDO
496C---------------------
497 hconv = rvolu(19)
498 text = rvolu(25)
499 left = left + half*(rnm*dv/vol+hconv*area*dt1)
500 right = right- half*tbag_old*(rnm_old*dv/vold+hconv*area*dt1)
501 . + hconv*area*text*dt1
502C---------------------------------------------------------------------
503 tbag = right/left
504 tbag = max(tbag,zero)
505 p=rnm*tbag/vol
506C--------------------
507C RATE OF HEAT LOSS
508C--------------------
509 rvolu(55) = hconv*area*max((tbag-text),zero)
510C-------------------------------------------------------------------
511C GAMA DU MELANGE (inclus caracteristiques du gaz a l'initial).
512C-------------------------------------------------------------------
513 cpi=cpai+cpbi*tbag+cpci*tbag*tbag
514 . +cpdi*tbag*tbag*tbag
515 . +cpei/(tbag*tbag)
516 . +cpfi*tbag*tbag*tbag*tbag
517 cvi=cpi-rmwi
518 cp =gmi*cpi
519 cv =gmi*cvi
520 DO ij=1,njet
521 i_inj = iabs(ibagjet(13,ij))
522 i_typinj = igeo(22,i_inj)
523 ngases = igeo(23,i_inj)
524 DO ik=1,ngases
525 IF (i_typinj==1) THEN
526 i_gas = igeo(100+(ik-1)*3+1,i_inj)
527 ELSE IF (i_typinj==2) THEN
528 i_gas = igeo(100+(ik-1)*2+1,i_inj)
529 END IF
530 mw = pm(20,i_gas)
531 rmwg = r_igc1/mw
532 cpa =pm(21,i_gas)
533 cpb =pm(22,i_gas)
534 cpc =pm(23,i_gas)
535 cpd =pm(24,i_gas)
536 cpe =pm(25,i_gas)
537 cpf =pm(26,i_gas)
538 cpg =cpa+cpb*tbag+cpc*tbag*tbag
539 . +cpd*tbag*tbag*tbag
540 . +cpe/(tbag*tbag)
541 . +cpf*tbag*tbag*tbag*tbag
542 cvg =cpg-rmwg
543 kk=20+(ik-1)*4
544 gmtot= rbagjet(kk+1,ij)
545 cp =cp+gmtot*cpg
546 cv =cv+gmtot*cvg
547 ENDDO
548 ENDDO
549C CP =CP/AMTOT
550C CV =CV/AMTOT
551 gama=cp/cv
552C------
553 rvolu(1) =gama
554 rvolu(14)=rnm
555C------
556 IF(dt1==zero.OR.dv>zero)THEN
557 q=zero
558 ELSE
559 q=-amu*sqrt(p*area*rot/vol)*dv/area/dt1
560 ENDIF
561 IF (ispmd+1==pmain) THEN
562 wfext=wfext+(half*(p+pold+q+qold)-pext)*dv
563 rvolu(32)=rvolu(32) + (half*(p+pold)-pext)*dv
564 ENDIF
565C------
566 rvolu(12)=p
567 rvolu(13)=tbag
568 rvolu(16)=vol
569 rvolu(20)=amtot
570 rvolu(23)=q
571
572 rvolu(22)=zero
573 rvolu(24)=zero
574 DO ij=1,njet
575 i_inj = iabs(ibagjet(13,ij))
576 ngases = igeo(23,i_inj)
577 DO ik=1,ngases
578 kk=20+(ik-1)*4
579 rbagjet(kk+3,ij)=zero
580 rbagjet(kk+4,ij)=zero
581 ENDDO
582 rbagjet( 9,ij)=zero
583 rbagjet(10,ij)=zero
584 ENDDO
585C------
586C DP/DRO GAZ PARFAIT+ADIABATIQUE=GAMA*R/MW*TBAG
587 cx = sqrt(2.*gama*p*vol/(gama-1.)/(amtot+area*rot))
588 qx = amu*cx
589 dt = dtfac1(9)*vol/area/max(em20,qx+sqrt(qx*qx+cx*cx))
590 IF(dt<dt2)THEN
591 dt2=dt
592 nelts =ivolu(1)
593 itypts=9
594 ENDIF
595 IF(idtmin(9)==1.AND.dt<dtmin1(9)) THEN
596 tstop = tt
597 IF (ispmd+1==pmain) THEN
598 WRITE(iout,*) '-- MINIMUM MONITORED VOLUME TIME STEP '
599 WRITE(istdo,*) '-- MINIMUM MONITORED VOLUME TIME STEP '
600 ENDIF
601 ELSEIF(idtmin(9)==5.AND.dt<dtmin1(9)) THEN
602 mstop = 2
603 IF (ispmd+1==pmain) THEN
604 WRITE(iout,*) '-- MINIMUM MONITORED VOLUME TIME STEP '
605 WRITE(istdo,*) '-- MINIMUM MONITORED VOLUME TIME STEP '
606 ENDIF
607 ENDIF
608 ENDIF
609C
610 RETURN
611 END
airbaga1
subroutine airbaga1(ivolu, njet, ibagjet, nvent, ibaghol, icbag, rvolu, rbagjet, rbaghol, rcbag, rvoluv, npc, tf, nsensor, sensor_tab, vol, pmain, geo, igeo, pm, wfext, python)
Definition airbaga1.F:37
my_real
#define my_real
Definition cppsort.cpp:32
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
max
#define max(a, b)
Definition macros.h:21