45
46
47
48 USE elbufdef_mod
50 USE sensor_mod
51 use glob_therm_mod
52 use element_mod , only : nixs,nixq,nixc,nixt,nixr,nixp,nixtg
53
54
55
56#include "implicit_f.inc"
57
58
59
60#include "param_c.inc"
61#include "com01_c.inc"
62#include "com04_c.inc"
63#include "com08_c.inc"
64#include "parit_c.inc"
65
66
67
68 INTEGER ,INTENT(IN) :: NSENSOR
69 INTEGER IACTIV(LACTIV,*),IPARG(NPARG,*),
70 . IXS(NIXS,*), IXQ(NIXQ,*), IXC(NIXC,*),
71 . IXT(NIXT,*),IXP(NIXP,*), IXR(NIXR,*), IXTG(NIXTG,*)
72 INTEGER IBCV(*), IBCR(*), IGROUPS(*)
73 my_real fsky(*), fconv(*), fradia(*)
74 my_real x(3,*), factiv(lractiv,*), temp(*), mcp(*), pm(npropm,*)
76 TYPE(ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
77 TYPE (SENSOR_STR_) ,DIMENSION(NSENSOR) ,INTENT(IN) :: SENSOR_TAB
78 type (glob_therm_) ,intent(inout) :: glob_therm
79
80 TYPE (GROUP_) , DIMENSION(NGRBRIC) :: IGRBRIC
81 TYPE (GROUP_) , DIMENSION(NGRQUAD) :: IGRQUAD
82 TYPE (GROUP_) , DIMENSION(NGRSHEL) :: IGRSH4N
83 TYPE (GROUP_) , DIMENSION(NGRSH3N) :: IGRSH3N
84 TYPE (GROUP_) , DIMENSION(NGRTRUS) :: IGRTRUSS
85 TYPE (GROUP_) , DIMENSION(NGRBEAM) :: IGRBEAM
86 TYPE (GROUP_) , DIMENSION(NGRSPRI) :: IGRSPRING
87
88
89
90 INTEGER I,N,IFLAG,ISENS,OFFCHANGE
91 INTEGER ITHERM_FE
93
94 offchange = 0
95 itherm_fe = glob_therm%ITHERM_FE
96
97 DO n = 1,nactiv
98 isens = iactiv(2,n)
99 iflag = iactiv(1,n)
100 IF (iactiv(10,n) == 1) THEN
101 IF (isens > 0) THEN
102
103
104
105 IF (iflag == 1 .AND. tt > sensor_tab(isens)%TSTART) THEN
106 CALL eloff(ixs ,ixq ,ixc ,ixp ,ixt ,
107 . ixr ,ixtg ,iparg ,
108 . iactiv ,tt ,iflag ,n ,elbuf_tab,
109 . x ,temp ,mcp ,pm ,igroups ,
110 . mcp_off ,igrbric ,igrquad ,igrsh4n ,igrsh3n ,
111 . igrtruss,igrbeam ,igrspring,itherm_fe)
112
113 iactiv(1,n) = 0
114 offchange=1
115 IF (iparit /= 0) THEN
116 DO i=1,8*lsky
117 fsky(i)=zero
118 ENDDO
119 ENDIF
120
121
122
123 ELSEIF (iflag == 0 .AND. tt <= sensor_tab(isens)%TSTART) THEN
124 CALL eloff(ixs ,ixq ,ixc ,ixp ,ixt ,
125 . ixr ,ixtg ,iparg ,
126 . iactiv ,tt ,iflag ,n ,elbuf_tab,
127 . x ,temp ,mcp ,pm ,igroups ,
128 . mcp_off ,igrbric ,igrquad ,igrsh4n ,igrsh3n ,
129 . igrtruss,igrbeam ,igrspring,itherm_fe)
130
131 iactiv(1,n) = 1
132 offchange=1
133 ENDIF
134 ENDIF
135
136 ELSEIF(iactiv(10,n) == 2) THEN
137
138
139
140 IF(tt == zero) THEN
141 CALL eloff(ixs ,ixq ,ixc ,ixp ,ixt ,
142 . ixr ,ixtg ,iparg ,
143 . iactiv ,tt ,1 ,n ,elbuf_tab,
144 . x ,temp ,mcp ,pm ,igroups ,
145 . mcp_off ,igrbric ,igrquad ,igrsh4n ,igrsh3n ,
146 . igrtruss,igrbeam ,igrspring,itherm_fe)
147
148 iactiv(1,n) = 0
149 offchange=1
150 ENDIF
151
152
153
154 startt = factiv(1,n)
155 stopt = factiv(2,n)
156 iflag = iactiv(1,n)
157 IF(itherm_fe > 0) THEN
158 startt = startt / glob_therm%THEACCFACT
159 stopt = stopt / glob_therm%THEACCFACT
160 ENDIF
161 IF(iflag == 0 .AND. tt >= startt .AND. tt < stopt) THEN
162 CALL eloff(ixs ,ixq ,ixc ,ixp ,ixt ,
163 . ixr ,ixtg ,iparg ,
164 . iactiv ,tt ,iflag ,n ,elbuf_tab,
165 . x ,temp ,mcp ,pm ,igroups ,
166 . mcp_off ,igrbric ,igrquad ,igrsh4n ,igrsh3n ,
167 . igrtruss,igrbeam ,igrspring,itherm_fe)
168 iactiv(1,n) = 1
169 offchange=1
170 ENDIF
171
172
173
174 IF(iflag == 1 .AND. tt > stopt) THEN
175 CALL eloff(ixs ,ixq ,ixc ,ixp ,ixt ,
176 . ixr ,ixtg ,iparg ,
177 . iactiv ,tt ,1 ,n ,elbuf_tab,
178 . x ,temp ,mcp ,pm ,igroups ,
179 . mcp_off ,igrbric ,igrquad ,igrsh4n ,igrsh3n ,
180 . igrtruss,igrbeam ,igrspring,itherm_fe)
181 iactiv(1,n) = 0
182 offchange=1
183 ENDIF
184 ENDIF
185 ENDDO
186
187 IF(itherm_fe > 0 .AND. offchange == 1) THEN
188 IF(glob_therm%NUMCONV > 0)
CALL convecoff(ibcv, fconv, iparg, igroups, elbuf_tab,glob_therm)
189 IF(glob_therm%NUMRADIA > 0)
CALL radiatoff(ibcr, fradia,iparg, igroups, ixs, elbuf_tab,glob_therm)
190 ENDIF
191
192 IF(itherm_fe > 0 .AND. tt == zero) THEN
193
194 iflag=-1
195 CALL eloff(ixs ,ixq ,ixc ,ixp ,ixt ,
196 . ixr ,ixtg ,iparg ,
197 . iactiv ,tt ,iflag ,n ,elbuf_tab,
198 . x ,temp ,mcp ,pm ,igroups ,
199 . mcp_off ,igrbric ,igrquad ,igrsh4n ,igrsh3n ,
200 . igrtruss,igrbeam ,igrspring,itherm_fe)
201 ENDIF
202
203
204
205 RETURN
subroutine convecoff(ibcv, fconv, iparg, igroups, elbuf_tab, glob_therm)
subroutine eloff(ixs, ixq, ixc, ixp, ixt, ixr, ixtg, iparg, iactiv, time, iflag, nn, elbuf_tab, x, temp, mcp, pm, igroups, mcp_off, igrbric, igrquad, igrsh4n, igrsh3n, igrtruss, igrbeam, igrspring, itherm_fe)
subroutine radiatoff(ibcr, fradia, iparg, igroups, ixs, elbuf_tab, glob_therm)