38
39
40
41 USE intbufdef_mod
43 USE output_mod
44
45
46
47#include "implicit_f.inc"
48
49
50
51 TYPE(OUTPUT_) :: OUTPUT
52 INTEGER IPARI(*), ICODT(*), ISKY(*)
53
55 . x(*), a(*), fsav(*),fskyi(*),
56 . fcont(3,*),fncont(3,*), ftcont(3,*)
57
58 TYPE(INTBUF_STRUCT_) INTBUF_TAB
59 TYPE(H3D_DATABASE) :: H3D_DATA
60
61
62
63#include "com08_c.inc"
64
65
66
67
69 . startt,ptmax
70
71 INTEGER IBC,NSN,NMN,NTY
72
73
74 nsn =ipari(5)
75 nmn =ipari(6)
76 nty =ipari(7)
77 ibc =ipari(11)
78
79
80 startt=intbuf_tab%VARIABLES(3)
81 IF(startt>tt) RETURN
82
83 ptmax=intbuf_tab%VARIABLES(16)
84 IF(nty==3) THEN
85 ptmax = ptmax*half
86 ENDIF
87 CALL invoi2(x ,intbuf_tab%IRECTM,intbuf_tab%LMSR,intbuf_tab%MSR,
88 1 intbuf_tab%NSV,intbuf_tab%ILOCS,intbuf_tab%IRTLM,intbuf_tab%NSEGM,nsn )
89 IF(nty==3)THEN
90 CALL invoi2(x ,intbuf_tab%IRECTS,intbuf_tab%LNSV,intbuf_tab%NSV,intbuf_tab%MSR,
91 1 intbuf_tab%ILOCM,intbuf_tab%IRTLS,intbuf_tab%NSEGS,nmn )
92 ENDIF
93 CALL i3for2(output, x ,a ,intbuf_tab%IRECTM,intbuf_tab%LMSR,intbuf_tab%MSR,
94 1 intbuf_tab%NSV,intbuf_tab%ILOCS,intbuf_tab%IRTLM,intbuf_tab%STFM,nsn ,
95 2 nmn ,nty ,intbuf_tab%CSTS,intbuf_tab%IRTLOM,intbuf_tab%FRICOS,
96 3 intbuf_tab%VARIABLES(1),intbuf_tab%STFNS,ibc ,icodt ,1 ,
97 4 fsav ,fskyi ,isky ,ptmax , intbuf_tab%AREAS ,
98 5 fcont ,fncont , ftcont, h3d_data )
99 IF(nty==3)THEN
100 CALL i3for2(output, x ,a ,intbuf_tab%IRECTS,intbuf_tab%LNSV,intbuf_tab%NSV,
101 1 intbuf_tab%MSR,intbuf_tab%ILOCM,intbuf_tab%IRTLS,intbuf_tab%STFS,nmn ,
102 2 nsn ,nty ,intbuf_tab%CSTM,intbuf_tab%IRTLOS,intbuf_tab%FRICOM,
103 3 intbuf_tab%VARIABLES(1),intbuf_tab%STFNM,0 ,icodt ,-1 ,
104 4 fsav ,fskyi ,isky ,ptmax ,intbuf_tab%AREAM ,
105 5 fcont ,fncont , ftcont, h3d_data )
106 ENDIF
107
108 RETURN
subroutine i3for2(output, x, e, irect, lmsr, msr, nsv, iloc, irtl, stf, nsn, nmn, nty, cst, irtlo, fric0, frigap, stfn, ibc, icodt, imast, fsav, fskyi, isky, ptmax, areas, fcont, fncont, ftcont, h3d_data)
subroutine invoi2(x, irect, lmsr, msr, nsv, iloc, irtl, nseg, nsn, nrt)