OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_thgrou.F File Reference
#include "implicit_f.inc"
#include "scr23_c.inc"
#include "scr17_c.inc"
#include "com10_c.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "param_c.inc"
#include "sphcom.inc"
#include "fxbcom.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine hm_read_thgrou (ithgrp, ithbuf, itab, itabm1, ixtg, ixs, ixq, ixc, ixt, ixp, ixr, kxx, ixx, ipart, ifi, nthwa, kxsp, ixri, iskwn, iframe, nthgrp2, pathid, suthid, fxbipm, iparth, nparth, nvparth, nvsubth, imerge, ithvar, flagabf, nvarabf, nom_opt, ptr_nopt_fxby, ptr_nopt_inter, ptr_nopt_rwall, ptr_nopt_sect, ptr_nopt_joint, ptr_nopt_monv, ptr_nopt_acc, ptr_nopt_skw, ptr_nopt_gau, ptr_nopt_clus, ptr_nopt_sphio, isphio, rfi, t_monvol, igrsurf, subset, ithflag, npby, lsubmodel, iparg, iparts, ipartq, ipartc, ipartt, ipartp, ipartr, ipartg, ipartx, ipartsp, ipartig3d, lithbufmx, map_tables, iflag, ptr_nopt_slipring, ptr_nopt_retractor, sensors, interfaces, ipari, dump_thnms1_file, itherm_fe, checksum, nsubdom, ipri)

Function/Subroutine Documentation

◆ hm_read_thgrou()

subroutine hm_read_thgrou ( integer, dimension(nithgr,*) ithgrp,
integer, dimension(lithbufmx) ithbuf,
integer, dimension(numnod) itab,
integer, dimension(*) itabm1,
integer, dimension(nixtg,numeltg) ixtg,
integer, dimension(nixs,numels) ixs,
integer, dimension(nixq,numelq) ixq,
integer, dimension(nixc,numelc) ixc,
integer, dimension(nixt,numelt) ixt,
integer, dimension(nixp,numelp) ixp,
integer, dimension(nixr,numelr) ixr,
integer, dimension(nixx,*) kxx,
integer, dimension(*) ixx,
integer, dimension(lipart1,*) ipart,
integer ifi,
integer nthwa,
integer, dimension(nisp,*) kxsp,
integer, dimension(4,*) ixri,
integer, dimension(liskn,*) iskwn,
integer, dimension(liskn,numfram+1) iframe,
integer nthgrp2,
integer, dimension(*) pathid,
integer, dimension(*) suthid,
integer, dimension(nbipm,*) fxbipm,
integer, dimension(nparth,*) iparth,
integer nparth,
integer nvparth,
integer nvsubth,
integer, dimension(*) imerge,
integer, dimension(*) ithvar,
integer flagabf,
integer nvarabf,
integer, dimension(lnopt1,*) nom_opt,
integer ptr_nopt_fxby,
integer ptr_nopt_inter,
integer ptr_nopt_rwall,
integer ptr_nopt_sect,
integer ptr_nopt_joint,
integer ptr_nopt_monv,
integer ptr_nopt_acc,
integer ptr_nopt_skw,
integer ptr_nopt_gau,
integer ptr_nopt_clus,
integer ptr_nopt_sphio,
integer, dimension(nisphio,*) isphio,
integer rfi,
type(monvol_struct_), dimension(nvolu), intent(in) t_monvol,
type (surf_), dimension(nsurf) igrsurf,
type (subset_), dimension(nsubs) subset,
integer ithflag,
integer, dimension(*) npby,
type (submodel_data), dimension(nsubmod), intent(in) lsubmodel,
integer, dimension(nparg,*) iparg,
integer, dimension(numels), target iparts,
integer, dimension(numelq), target ipartq,
integer, dimension(numelc), target ipartc,
integer, dimension(numelt), target ipartt,
integer, dimension(numelp), target ipartp,
integer, dimension(numelr), target ipartr,
integer, dimension(numeltg), target ipartg,
integer, dimension(numelx), target ipartx,
integer, dimension(numsph), target ipartsp,
integer, dimension(numelig3d), target ipartig3d,
integer, intent(in) lithbufmx,
type(mapping_struct_), intent(in) map_tables,
integer iflag,
integer ptr_nopt_slipring,
integer ptr_nopt_retractor,
type (sensors_), intent(in) sensors,
type (interfaces_), intent(inout) interfaces,
integer, dimension(npari,ninter), intent(inout) ipari,
integer, intent(in) dump_thnms1_file,
integer, intent(in) itherm_fe,
type(checksum_option_), intent(in) checksum,
integer, intent(in) nsubdom,
integer, intent(in) ipri )

Definition at line 61 of file hm_read_thgrou.F.

78C-----------------------------------------------
79 USE my_alloc_mod
80 USE message_mod
81 USE groupdef_mod
82 USE monvol_struct_mod
83 USE submodel_mod
84 USE hm_option_read_mod
85 USE mapping_option_mod
86 USE sensor_mod
87 USE interfaces_mod
88 USE th_titles_mod
89 USE write_thnms1_mod
90 USE names_and_titles_mod , ONLY : nchartitle, ncharline
91 USE user_id_mod , ONLY : id_limit
92 USE checksum_output_option_mod
93 USE checksum_starter_option_mod
94 use hm_read_th_checksum_mod
95 use element_mod , only : nixs,nixq,nixc,nixt,nixp,nixr,nixtg
96C-----------------------------------------------
97C I m p l i c i t T y p e s
98C-----------------------------------------------
99#include "implicit_f.inc"
100C-----------------------------------------------
101C C o m m o n B l o c k s
102C-----------------------------------------------
103#include "scr23_c.inc"
104#include "scr17_c.inc"
105#include "com10_c.inc"
106#include "com01_c.inc"
107#include "com04_c.inc"
108#include "param_c.inc"
109#include "sphcom.inc"
110#include "fxbcom.inc"
111C-----------------------------------------------
112C D u m m y A r g u m e n t s
113C-----------------------------------------------
114 INTEGER,INTENT(IN) :: NSUBDOM
115 INTEGER,INTENT(IN) :: IPRI
116 INTEGER,INTENT(IN) :: LITHBUFMX
117 INTEGER,INTENT(IN) :: ITHERM_FE
118 INTEGER ITAB(NUMNOD),ITABM1(*), IFLAG,
119 . IXS(NIXS,NUMELS),IXQ(NIXQ,NUMELQ),IXC(NIXC,NUMELC),IXT(NIXT,NUMELT),
120 . IXP(NIXP,NUMELP),IXR(NIXR,NUMELR),IXTG(NIXTG,NUMELTG),KXX(NIXX,*),
121 . IXX(*),IPART(LIPART1,*),ITHGRP(NITHGR,*),ITHBUF(LITHBUFMX),
122 . IMERGE(*),
123 . IFI,NTHWA,KXSP(NISP,*),IXRI(4,*),
124 . IFRAME(LISKN,NUMFRAM+1), NTHGRP2,PATHID(*),SUTHID(*),
125 . FXBIPM(NBIPM,*),ISKWN(LISKN,*),
126 . NPARTH,IPARTH(NPARTH,*),NVPARTH,ITYP,
127 . NVSUBTH,ITHVAR(*),FLAGABF,NVARABF,
128 . RFI,ISPHIO(NISPHIO,*),ITHFLAG,IPARG(NPARG,*)
129 INTEGER NOM_OPT(LNOPT1,*),PTR_NOPT_FXBY,PTR_NOPT_INTER,
130 . PTR_NOPT_RWALL,PTR_NOPT_SECT,PTR_NOPT_JOINT,
131 . PTR_NOPT_MONV,PTR_NOPT_ACC,PTR_NOPT_SKW,PTR_NOPT_GAU,
132 . PTR_NOPT_CLUS,PTR_NOPT_SPHIO,PTR_NOPT_SURF,PTR_NOPT_RETRACTOR,
133 . PTR_NOPT_SLIPRING
134 INTEGER NPBY(*)
135 INTEGER, TARGET :: IPARTS(NUMELS),IPARTQ(NUMELQ),IPARTC(NUMELC),
136 . IPARTT(NUMELT),IPARTP(NUMELP),IPARTR(NUMELR),IPARTG(NUMELTG),
137 . IPARTX(NUMELX),IPARTSP(NUMSPH),IPARTIG3D(NUMELIG3D)
138 INTEGER,INTENT(INOUT) :: IPARI(NPARI,NINTER)
139 TYPE(MONVOL_STRUCT_), DIMENSION(NVOLU), INTENT(IN) :: T_MONVOL
140 TYPE (SUBMODEL_DATA) ,DIMENSION(NSUBMOD) ,INTENT(IN) :: LSUBMODEL
141 TYPE(MAPPING_STRUCT_),INTENT(IN) :: MAP_TABLES
142 TYPE (INTERFACES_) ,INTENT(INOUT):: INTERFACES
143 TYPE (SENSORS_) ,INTENT(IN) :: SENSORS
144C-----------------------------------------------
145 TYPE (SUBSET_) , DIMENSION(NSUBS) :: SUBSET
146 TYPE (SURF_) , DIMENSION(NSURF) :: IGRSURF
147 INTEGER, INTENT(IN):: DUMP_THNMS1_FILE
148 TYPE(checksum_option_) , INTENT(IN) :: checksum
149C-----------------------------------------------
150C L o c a l V a r i a b l e s
151C-----------------------------------------------
152 INTEGER J,JJ,II, I,ISU,ID,NOSYS,J10(10),NTOT,KK,IER,
153 . OK,IGS,IGRS,NSU,K,L,CONT,IAD0,IADV,NTRI,NL,
154 . IFITMP,IADFIN,NUMEL,NV,NVAR,M,N,IAD2,IAD,NSNE,NNE,ISK,
155 . NVALL,TAGP,TAGS,IFIX_TMP,IBID,NIN
156 CHARACTER(LEN=NCHARTITLE) :: TITR,TITR1
157 CHARACTER MESS*40
158 CHARACTER(LEN=NCHARLINE) :: KEY,COPT
159 INTEGER NVARN,NVARN1,NVARN1A,NVARN2,NVARNPINCH,NVARS,NVARC,NVART,NVARP,NVARR,NVARUR
160 INTEGER NVARS1,NVARS2,NVARS3,NVARS4,NVARS5,NVARS6,NVARS7,NVARS8,NVARS9,NVARS10,NVARSNLOC
161 INTEGER NVARNS,NVARSPH
162 INTEGER NVARIN,NVARRW,NVARRB,NVARAC,NVARSE,NVARJO,NVARFX,NVARFXM
163 INTEGER NVARAB,NVARMV4,NVARMV,NVARPA,NVARMVENT,NVARSENS,NVARSENSG
164 INTEGER NVARCHECKSUM,NVARCHECKSUMG
165 INTEGER NVARF1,NVARFR,NVARGAU,NVARCLUS,NVARCLUSG,NVARFLOW,NVARSURF
166 INTEGER NVARRIV,NVARRIVG,NVARSLIP,NVARRET,NSELRT,INOUT
167 INTEGER HM_NTHGRP,NTHACCEL,NTHINTER,NTHRWALL,NTHSECTIO,NTHCLUS
168 INTEGER NTHBEAM,NTHTRUS,NTHBRIC,NTHNODE,NTHSHEL,NTHSH3N,NTHSPRING,NTHRBODY
169 INTEGER NTHMONVOL,HM_NTHPART ,NUMTHPART, HM_NTHSUBS, HM_NTHSURF
170 INTEGER HM_NTHSPHCEL, HM_NTHQUAD, HM_NTHSPHFLOW
171 INTEGER HM_NTHGAUGE, HM_NTHFXBODY, HM_NTHFRAME, HM_NTHCYLJO, HM_NTHNSTRAND
172 INTEGER HM_NTHTRIA,HM_NTHSLIPRING,HM_NTHRETRACTOR,HM_NTHSENS,HM_NTHCHECKSUM
173 parameter(nvarriv = 10,nvarrivg = 1)
174 INTEGER IVARRIV(18,NVARRIVG)
175 CHARACTER*10 VARRIV(NVARRIV),VARRIVG(NVARRIVG)
176C-----------------------------------------------
177 parameter(nvarn1 = 19,nvarn1a = 600,nvarn2 = 10 ,nvarnpinch = 9)
178 parameter(nvars1 = 196,nvars2 = 567,nvars3 = 369,nvars4 = 492,nvars5 = 22)
179 parameter(nvars6 = 97200,nvars7 = 97200,nvars8 = 516,nvars9 = 6,nvars10 = 1,nvarsnloc = 2)
180 parameter(nvars =239555 ,nvarc = 37856,nvart = 6)
181 parameter(nvarp = 337,nvarr = 66,nvarur = 12)
182 parameter(nvarns = 4,nvarsph = 41)
183 parameter(nvarin = 29,nvarrw = 6,nvarrb =15,nvarfx =4)
184 parameter(nvarfxm = 3)
185 parameter(nvarac = 6,nvarse =39,nvarjo = 6)
186 parameter(nvarab = 7,nvarmv4= 9,nvarmv = 19,nvarmvent = 150)
187 parameter(nvarpa = 32)
188 PARAMETER (nvarf1 = 18,nvarfr = 24,nvargau = 8,nvarclus=11,nvarflow=1,nvarsurf = 6,nvarslip=6,nvarret=3)
189 CHARACTER*10, DIMENSION(:), ALLOCATABLE :: VARN
190 CHARACTER*10 :: VART(NVART)
191 CHARACTER*10 VARP(NVARP),VARR(NVARR),VARUR(NVARUR)
192 CHARACTER(LEN=10),DIMENSION(:),ALLOCATABLE::VARS
193 CHARACTER(LEN=10),DIMENSION(:),ALLOCATABLE::VARC
194 CHARACTER*10 VARS1(NVARS1),VARS2(NVARS2), VARS3(NVARS3), VARS4(NVARS4),
195 . VARS5(NVARS5),VARS8(NVARS8),VARS9(NVARS9),VARS10(NVARS10),
196 . VARN1(NVARN1),VARN1A(NVARN1A),
197 . VARN2(NVARN2),
198 . VARNPINCH(NVARNPINCH),
199 . VARSG1(202),VARSG4(1),
200 . VARSNLOC(NVARSNLOC)
201 CHARACTER(LEN=10),DIMENSION(:),ALLOCATABLE::VARS6
202 CHARACTER(LEN=10),DIMENSION(:),ALLOCATABLE::VARS7
203 CHARACTER(LEN=10),DIMENSION(:),ALLOCATABLE::VARSG2
204 CHARACTER(LEN=10),DIMENSION(:),ALLOCATABLE::VARSG3
205 CHARACTER*10 VARNS(NVARNS),VARSPH(NVARSPH)
206 CHARACTER*10 VARIN(NVARIN),VARRW(NVARRW),VARRB(NVARRB)
207 CHARACTER*10 VARAC(NVARAC),VARSE(NVARSE),VARJO(NVARJO)
208 CHARACTER*10 VARAB(NVARAB),VARMV4(NVARMV4),VARMV(NVARMV),
209 . VARMVENT(NVARMVENT)
210 CHARACTER*10 VARPA(NVARPA),VARFX(NVARFX),VARFXM(NVARFXM)
211 CHARACTER*10 VARF1(NVARF1),VARFR(NVARFR),VARGAU(NVARGAU),
212 . VARSLIP(NVARSLIP),VARRET(NVARRET)
213 CHARACTER*10 VARCLUS(NVARCLUS),VARFLOW(NVARFLOW),VARSURF(NVARSURF)
214 INTEGER NVARNG,NVARSG,NVARCG,NVARTG,NVARPG,NVARRG,NVARURG
215 INTEGER NVARNSG,NVARSPG
216 INTEGER NVARING,NVARRWG,NVARRBG,NVARACG,NVARSEG,NVARJOG
217 INTEGER NVARABG,NVARMG4,NVARMVG,NVARPAG,NVARFXG,NVARFXMG
218 INTEGER NVARF1G,NVARFRG,NVARGAUG,NVARFLOWG,NVARMVENTG,NVARSURFG
219 INTEGER NVARSLIPG,NVARRETG
220 parameter(nvarng = 8,nvarsg = 39767,nvarcg = 675,nvartg = 1)
221 parameter(nvarpg = 1,nvarrg = 1,nvarurg = 1)
222 parameter(nvarnsg = 1,nvarspg = 2)
223 parameter(nvaring = 6,nvarrwg = 3,nvarrbg = 6,nvarfxg = 1)
224 parameter(nvarfxmg =1)
225 parameter(nvaracg = 2,nvarseg = 7,nvarjog = 3)
226 parameter(nvarabg = 1,nvarmg4 = 1,nvarmvg = 62,nvarmventg=62)
227 parameter(nvarpag = 1)
228 parameter(nvarf1g = 3,nvarfrg = 4,nvargaug = 1,nvarclusg=2,nvarflowg=1,
229 . nvarsurfg = 1,nvarslipg=1,nvarretg = 1)
230 parameter(nvarsens = 1,nvarsensg = 1)
231 parameter(nvarchecksum = 1,nvarchecksumg = 1)
232 CHARACTER*10 VARNG(NVARNG)
233 CHARACTER*10 VARCG(NVARCG),VARTG(NVARTG)
234 CHARACTER*10 VARPG(NVARPG),VARRG(NVARRG),VARURG(NVARURG)
235 CHARACTER*10 VARNSG(NVARNSG),VARSPG(NVARSPG)
236 CHARACTER*10 VARING(NVARING),VARRWG(NVARRWG),VARRBG(NVARRBG)
237 CHARACTER*10 VARACG(NVARACG),VARSEG(NVARSEG),VARJOG(NVARJOG)
238 CHARACTER*10 VARABG(NVARABG),VARMVG4(NVARMG4),VARMVG(NVARMVG),
239 . VARMVENTG(NVARMVENTG),VARSENS(NVARSENS),VARSENSG(NVARSENSG),
240 . VARCHECKSUM(NVARCHECKSUM),VARCHECKSUMG(NVARCHECKSUMG)
241 CHARACTER*10 VARPAG(NVARPAG),VARFXG(NVARFXG),VARFXMG(NVARFXMG)
242 CHARACTER*10 VARF1G(NVARF1G),VARFRG(NVARFRG),VARGAUG(NVARGAUG)
243 CHARACTER*10 VARCLUSG(NVARCLUSG),VARFLOWG(NVARFLOWG)
244 CHARACTER*10 VARSURFG(NVARSURFG),VARSLIPG(NVARSLIPG),VARRETG(NVARRETG)
245 INTEGER,DIMENSION(:,:),ALLOCATABLE :: IVARSG
246 CHARACTER(LEN=10),DIMENSION(:),ALLOCATABLE:: VARSG
247C-----------------------------------------------
248 INTEGER IVARNG(18,NVARNG)
249 INTEGER IVARCG(18,NVARCG),IVARTG(18,NVARTG)
250 INTEGER IVARPG(18,NVARPG),IVARRG(18,NVARRG),IVARURG(18,NVARURG)
251 INTEGER IVARNSG(18,NVARNSG),IVARSPG(18,NVARSPG),IVNS2R(18,NVARNSG)
252 INTEGER IVARING(18,NVARING),IVARRWG(18,NVARRWG)
253 INTEGER IVARRBG(18,NVARRBG),IVARACG(18,NVARACG)
254 INTEGER IVARSEG(18,NVARSEG),IVARJOG(18,NVARJOG)
255 INTEGER IVARABG(18,NVARABG),IVARMG4(18,NVARMG4),
256 . IVARMVG(18,NVARMVG),IVARMVENTG(18,NVARMVENTG)
257 INTEGER IVARPAG(18,NVARPAG),IVARFXG(18,NVARFXG),IVARSENSG(18,NVARSENSG)
258 INTEGER IVARFXMG(18,NVARFXMG)
259 INTEGER IVARF1G(18,NVARF1G),IVARFRG(18,NVARFRG),
260 . IVARGAUG(18,NVARGAUG),IVARCLUSG(18,NVARCLUSG),
261 . IVARFLOWG(18,NVARFLOWG),IVARSURFG(18,NVARSURFG),
262 . IVARSLIPG(18,NVARSLIPG),IVARRETG(18,NVARRETG),
263 . IVARCHECKSUMG(18,NVARCHECKSUMG)
264 INTEGER NSMA
265 INTEGER IUN,INODACC,INODFRA
266 INTEGER NNOD, NSHEL, NTRUS, NBEAM, NSPRG, NCHTG, NCHINT, NCHJNT,
267 . CHNOD, CHELT, CNT1, CNT2, IVSG, IREC_0, OUT_VENT, NG, NEL,
268 . NFT, NUMPART, I15ATH, I15A, I15B, I15C, I15D, I15E, I15F,
269 . I15G, I15H, I15I, I15J, I15K, I15L, ITY
270 my_real
271 . bid
272 INTEGER IVCG
273 CHARACTER CHJJ*3,CHKK*2,CHII*2,CHJKK*4
274 CHARACTER CHI*1,CHJ*3,CHK*1,CHJS*1
275 CHARACTER*10 VARC2(574),VARC1(208),VARCG1(167),VARCG2(407), VAR ,
276 . VAR_TMP,VARC3(500),VARCG3(100),
277 . VARCGPINCH(1),VARCPINCH(7),VARCNLOC(2)
278 CHARACTER*80 STR, MES
279 LOGICAL :: IS_AVAILABLE
280 INTEGER, DIMENSION(:), POINTER :: IPARTEL_P
281 CHARACTER THFILE*4
282
283 CHARACTER(LEN=100),DIMENSION(:),ALLOCATABLE:: VARN1_TITLE,VARN1A_TITLE,VARN2_TITLE,VARNPINCH_TITLE,
284 . VARP_TITLE,VARR_TITLE,VART_TITLE,
285 . VARS1_TITLE,VARS2_TITLE,VARS3_TITLE,VARS4_TITLE,
286 . VARS5_TITLE,VARS6_TITLE,VARS7_TITLE,VARS8_TITLE,VARS9_TITLE,VARS10_TITLE,
287 . VARC_TITLE,
288 . VARSNLOC_TITLE,VARNS_TITLE,VARSPH_TITLE,
289 . VARIN_TITLE,VARRW_TITLE,VARRB_TITLE,
290 . VARMV_TITLE,VARSE_TITLE,VARAC_TITLE,
291 . VARJO_TITLE,VARMVENT_TITLE,
292 . VARPA_TITLE,VARFX_TITLE,
293 . VARGAU_TITLE,VARFR_TITLE,
294 . VARSLIP_TITLE,VARRET_TITLE,
295 . VARCLUS_TITLE,VARFLOW_TITLE,
296 . VARSURF_TITLE,VARSENS_TITLE
297
298 INTEGER :: LEN_TITLE
299 INTEGER :: LEN_CHECKSUM
300 CHARACTER(LEN=NCHARLINE):: CHECKSUM_TITLE
301 CHARACTER(LEN=64):: DECK_CHECKSUM
302C-----------------------------------------------
303C E x t e r n a l F u n c t i o n s
304C-----------------------------------------------
305 INTEGER USR2SYS,ULIST2S,LISTCNT
306C-----------------------------------------------
307C D a t a
308C-----------------------------------------------
309 DATA mess/'TH GROUP DEFINITION '/
310 DATA iun/1/
311
312C-----------------------------------------------
313C TH TITLES
314C DON T MISS TO ADD TH TITLES FOR NEW TH OPTIONS
315C-----------------------------------------------
316
317 nvarn = nvarn1 + nvarn1a + nvarn2 + nvarnpinch
318 ALLOCATE(varn(nvarn))
319
320
321 ALLOCATE(varn1_title(nvarn1),varn1a_title(nvarn1a),varn2_title(nvarn2),varnpinch_title(nvarnpinch),
322 . varp_title(nvarp),varr_title(nvarr),vart_title(nvart),
323 . vars1_title(nvars1),vars2_title(nvars2),vars3_title(nvars3),vars4_title(nvars4),
324 . vars5_title(nvars5),vars6_title(nvars6),vars7_title(nvars7),vars8_title(nvars8),vars9_title(nvars9),
325 . varc_title(nvarc),vars10_title(nvars10),
326 . varsnloc_title(nvarsnloc),varns_title(nvarns),varsph_title(nvarsph),
327 . varin_title(nvarin),varrw_title(nvarrw),varrb_title(nvarrb),
328 . varmv_title(nvarmv),varse_title(nvarse),varac_title(nvarac),
329 . varjo_title(nvarjo),varmvent_title(nvarmvent),
330 . varpa_title(nvarpa),varfx_title(nvarfx),
331 . vargau_title(nvargau),varfr_title(nvarfr),
332 . varslip_title(nvarslip),varret_title(nvarret),
333 . varclus_title(nvarclus),varflow_title(nvarflow),
334 . varsurf_title(nvarsurf),varsens_title(1))
335
336
337 CALL th_titles(
338 1 nvarn1 ,nvarn1a ,nvarn2 ,nvarnpinch ,nvars1 ,
339 1 nvars2 ,nvars3 ,nvars4 ,nvars5 ,nvars6 ,
340 1 nvars7 ,nvars8 ,nvars9 ,nvars10 ,nvarsnloc ,
341 1 nvarp ,nvarr ,nvart ,nvarns ,nvarsph ,
342 2 nvarin ,nvarrw ,nvarrb ,nvarmv ,nvarse ,
343 3 nvarac ,nvarjo ,nvarmvent ,nvarpa ,nvarfx ,
344 4 nvargau ,nvarfr ,nvarslip ,nvarret ,nvarclus ,
345 5 nvarflow ,nvarsurf ,nvarc ,nvarsens ,
346 a varn1_title ,varn1a_title ,varn2_title ,
347 a varnpinch_title,varp_title ,varr_title ,vart_title ,
348 6 vars1_title ,vars2_title ,vars3_title ,vars4_title ,vars5_title ,
349 6 vars6_title ,vars7_title ,vars8_title ,vars9_title ,varsnloc_title ,
350 6 varc_title ,vars10_title ,
351 7 varns_title ,varsph_title ,varin_title ,
352 8 varrw_title ,varrb_title ,varmv_title ,varse_title ,varac_title ,
353 9 varjo_title ,varmvent_title,varpa_title ,varfx_title ,vargau_title ,
354 a varfr_title ,varslip_title ,varret_title ,varclus_title ,varflow_title ,
355 b varsurf_title ,varsens_title )
356
357C-----------------------------------------------
358C TH VARIABLES
359C-----------------------------------------------
360 ! Allocate temporary Arrays
361 ! Character Arrays - ue ALLOCATEkm
362 ALLOCATE(vars6(97200))
363 vars6(1:97200) = ' '
364 ALLOCATE(vars7(97200))
365 vars7(1:97200) = ' '
366 ALLOCATE(varsg2(19782))
367 varsg2(1:19782) = ' '
368 ALLOCATE(varsg3(19782))
369 varsg3(1:19782) = ' '
370 ALLOCATE(vars(nvars))
371 vars(1:nvars) = ' '
372 ALLOCATE(varc(nvarc))
373 varc(1:nvarc) = ' '
374 ALLOCATE(varsg(nvarsg))
375 varsg(1:nvarsg) = ' '
376
377 ! INTEGER Arrays - use MY_ALLOC
378 CALL my_alloc(ivarsg,18,nvarsg)
379C nodes
380 DATA varn1/'DX ','DY ','DZ ',
381 . 'VX ','VY ','VZ ',
382 . 'AX ','AY ','AZ ',
383 . 'VRX ','VRY ','VRZ ',
384 . 'ARX ','ARY ','ARZ ',
385 . 'X ','Y ','Z ',
386 . 'TEMP '/
387 DATA varn2/'REACX ','REACY ','REACZ ',
388 . 'REACXX ','REACYY ','REACZZ ',
389 . 'DRX ','DRY ','DRZ ',
390 . 'PEXT '/
391 DATA varnpinch/
392 . 'APINCHX ','APINCHY ','APINCHZ ',
393 . 'VPINCHX ','VPINCHY ','VPINCHZ ',
394 . 'DPINCHX ','DPINCHY ','DPINCHZ '/
395C solids
396 DATA vars1/
397 .'OFF ','SX ','SY ','SZ ','SXY ',
398 .'SYZ ','SXZ ','IE ','DENS ','BULK ',
399 .'VOL ','PLAS ','TEMP ','PLSR ','DAM1 ',
400 .'DAM2 ','DAM3 ','DAM4 ','DAMA ','SA1 ',
401 .'SA2 ','SA3 ','CR ','CAP ','K0 ',
402 .'RK ','TD ','EFIB ','ISTA ','VPLA ',
403 .'BFRAC ','WPLA ','SFIB ','DAM5 ','LSX ',
404 .'LSY ','LSZ ','LSXY ','LSYZ ','LSXZ ',
405 .'SX1 ','SX2 ','SX3 ','SX4 ','SX5 ','SX6 ',
406 .'SX7 ','SX8 ','SY1 ','SY2 ','SY3 ','SY4 ',
407 .'SY5 ','SY6 ','SY7 ','SY8 ','SZ1 ','SZ2 ',
408 .'SZ3 ','SZ4 ','SZ5 ','SZ6 ','SZ7 ','SZ8 ',
409 .'SXY1 ','SXY2 ','SXY3 ','SXY4 ','SXY5 ','SXY6 ',
410 .'SXY7 ','SXY8 ','SYZ1 ','SYZ2 ','SYZ3 ','SYZ4 ',
411 .'SYZ5 ','SYZ6 ','SYZ7 ','SYZ8 ','SXZ1 ','SXZ2 ',
412 .'SXZ3 ','SXZ4 ','SXZ5 ','SXZ6 ','SXZ7 ','SXZ8 ',
413 .'LSX1 ','LSX2 ','LSX3 ','LSX4 ','LSX5 ','LSX6 ',
414 .'LSX7 ','LSX8 ','LSY1 ','LSY2 ','LSY3 ','LSY4 ',
415 .'LSY5 ','LSY6 ','LSY7 ','LSY8 ','LSZ1 ','LSZ2 ',
416 .'LSZ3 ','LSZ4 ','LSZ5 ','LSZ6 ','LSZ7 ','LSZ8 ',
417 .'LSXY1 ','LSXY2 ','LSXY3 ','LSXY4 ','LSXY5 ','LSXY6 ',
418 .'LSXY7 ','LSXY8 ','LSYZ1 ','LSYZ2 ','LSYZ3 ','LSYZ4 ',
419 .'LSYZ5 ','LSYZ6 ','LSYZ7 ','LSYZ8 ','LSXZ1 ','LSXZ2 ',
420 .'LSXZ3 ','LSXZ4 ','LSXZ5 ','LSXZ6 ','LSXZ7 ','LSXZ8 ',
421 .'USR1 ','USR2 ','USR3 ','USR4 ','USR5 ','USR6 ',
422 .'USR7 ','USR8 ','USR9 ','USR10 ','USR11 ','USR12 ',
423 .'USR13 ','USR14 ','USR15 ','USR16 ','USR17 ','USR18 ',
424 .'USR19 ','USR20 ','USR21 ','USR22 ','USR23 ','USR24 ',
425 .'USR25 ','USR26 ','USR27 ','USR28 ','USR29 ','USR30 ',
426 .'USR31 ','USR32 ','USR33 ','USR34 ','USR35 ','USR36 ',
427 .'USR37 ','USR38 ','USR39 ','USR40 ','USR41 ','USR42 ',
428 .'USR43 ','USR44 ','USR45 ','USR46 ','USR47 ','USR48 ',
429 .'USR49 ','USR50 ','USR51 ','USR52 ','USR53 ','USR54 ',
430 .'USR55 ','USR56 ','USR57 ','USR58 ','USR59 ','USR60 '/
431C
432 DATA vars2/
433 .'SX111 ','SY111 ','SZ111 ','SXY111 ','SYZ111 ','SXZ111 ',
434 .'EP111 ','SX211 ','SY211 ','SZ211 ','SXY211 ','SYZ211 ',
435 .'SXZ211 ','EP211 ','SX311 ','SY311 ','SZ311 ','SXY311 ',
436 .'SYZ311 ','SXZ311 ','EP311 ','SX121 ','SY121 ','SZ121 ',
437 .'SXY121 ','SYZ121 ','SXZ121 ','EP121 ','SX221 ','SY221 ',
438 .'SZ221 ','SXY221 ','SYZ221 ','SXZ221 ','EP221 ','SX321 ',
439 .'SY321 ','SZ321 ','SXY321 ','SYZ321 ','SXZ321 ','EP321 ',
440 .'SX131 ','SY131 ','SZ131 ','SXY131 ','SYZ131 ','SXZ131 ',
441 .'EP131 ','SX231 ','SY231 ','SZ231 ','SXY231 ','SYZ331 ',
442 .'SXZ231 ','EP231 ','SX331 ','SY331 ','SZ331 ','SXY331 ',
443 .'SYZ331 ','SXZ331 ','EP331 ','SX141 ','SY141 ','SZ141 ',
444 .'SXY141 ','SYZ141 ','SXZ141 ','EP141 ','SX241 ','SY241 ',
445 .'SZ241 ','SXY241 ','SYZ241 ','SXZ241 ','EP241 ','SX341 ',
446 .'SY341 ','SZ341 ','SXY341 ','SYZ341 ','SXZ341 ','EP341 ',
447 .'SX151 ','SY151 ','SZ151 ','SXY151 ','SYZ151 ','SXZ151 ',
448 .'EP151 ','SX251 ','SY251 ','SZ251 ','SXY251 ','SYZ251 ',
449 .'SXZ251 ','EP251 ','SX351 ','SY351 ','SZ351 ','SXY351 ',
450 .'SYZ351 ','SXZ351 ','EP351 ','SX161 ','SY161 ','SZ161 ',
451 .'SXY161 ','SYZ161 ','SXZ161 ','EP161 ','sx261 ','sy261 ',
452 .'sz261 ','sxy261 ','syz261 ','sxz261 ','ep261 ','sx361 ',
453 .'sy361 ','sz361 ','sxy361 ','syz361 ','sxz361 ','ep361 ',
454 .'sx171 ','sy171 ','sz171 ','sxy171 ','syz171 ','sxz171 ',
455 .'ep171 ','sx271 ','sy271 ','sz271 ','sxy271 ','syz271 ',
456 .'sxz271 ','ep271 ','sx371 ','sy371 ','sz371 ','sxy371 ',
457 .'syz371 ','sxz371 ','ep371 ','sx181 ','sy181 ','sz181 ',
458 .'sxy181 ','syz181 ','sxz181 ','ep181 ','sx281 ','sy281 ',
459 .'sz281 ','sxy281 ','syz281 ','sxz281 ','ep281 ','sx381 ',
460 .'sy381 ','sz381 ','sxy381 ','syz381 ','sxz381 ','ep381 ',
461 .'sx191 ','sy191 ','sz191 ','sxy191 ','syz191 ','sxz191 ',
462 .'ep191 ','sx291 ','sy291 ','sz291 ','sxy291 ','syz291 ',
463 .'sxz291 ','ep291 ','sx391 ','sy391 ','sz391 ','sxy391 ',
464 .'syz391 ','sxz391 ','ep391 ','sx112 ','sy112 ','sz112 ',
465 .'sxy112 ','syz112 ','sxz112 ','ep112 ','sx212 ','sy212 ',
466 .'sz212 ','sxy212 ','syz212 ','sxz212 ','ep212 ','sx312 ',
467 .'sy312 ','sz112 ','sxy112 ','syz312 ','sxz112 ','ep312 ',
468 .'sx122 ','sy122 ','sz122 ','sxy122 ','syz122 ','sxz122 ',
469 .'ep122 ','sx222 ','sy222 ','sz222 ','sxy222 ','syz222 ',
470 .'sxz222 ','ep222 ','sx322 ','sy322 ','sz322 ','sxy322 ',
471 .'syz322 ','sxz322 ','ep322 ','sx132 ','sy132 ','sz132 ',
472 .'sxy132 ','syz132 ','sxz132 ','ep132 ','sx232 ','sy232 ',
473 .'sz232 ','sxy232 ','syz332 ','sxz232 ','ep232 ','sx332 ',
474 .'sy332 ','sz332 ','sxy332 ','syz332 ','sxz332 ','ep332 ',
475 .'sx142 ','sy142 ','sz142 ','sxy142 ','syz142 ','sxz142 ',
476 .'ep142 ','sx242 ','SY242 ','SZ242 ','SXY242 ','SYZ242 ',
477 .'SXZ242 ','EP242 ','SX342 ','SY342 ','SZ342 ','SXY342 ',
478 .'SYZ342 ','SXZ342 ','EP342 ','SX152 ','SY152 ','SZ152 ',
479 .'SXY152 ','SYZ152 ','SXZ152 ','EP152 ','SX252 ','SY252 ',
480 .'SZ252 ','SXY252 ','SYZ252 ','SXZ252 ','EP252 ','SX352 ',
481 .'SY352 ','SZ352 ','SXY352 ','SYZ352 ','SXZ352 ','EP352 ',
482 .'SX162 ','SY162 ','SZ162 ','SXY162 ','SYZ162 ','SXZ162 ',
483 .'EP162 ','SX262 ','SY262 ','SZ262 ','SXY262 ','SYZ262 ',
484 .'SXZ262 ','EP262 ','SX362 ','SY362 ','SZ362 ','SXY362 ',
485 .'SYZ362 ','SXZ362 ','EP362 ','SX172 ','SY172 ','SZ172 ',
486 .'SXY172 ','SYZ172 ','SXZ172 ','EP172 ','SX272 ','SY272 ',
487 .'SZ272 ','SXY272 ','SYZ272 ','SXZ272 ','EP272 ','SX372 ',
488 .'SY372 ','SZ372 ','SXY372 ','SYZ372 ','SXZ372 ','EP372 ',
489 .'SX182 ','SY182 ','SZ182 ','SXY182 ','SYZ182 ','SXZ182 ',
490 .'EP182 ','SX282 ','SY282 ','SZ282 ','SXY282 ','SYZ282 ',
491 .'SXZ282 ','EP282 ','SX382 ','SY382 ','SZ382 ','SXY382 ',
492 .'SYZ382 ','SXZ382 ','EP382 ','SX192 ','SY192 ','SZ192 ',
493 .'SXY192 ','SYZ192 ','SXZ192 ','EP192 ','SX292 ','SY292 ',
494 .'SZ292 ','SXY292 ','SYZ292 ','SXZ292 ','EP292 ','SX392 ',
495 .'SY392 ','SZ392 ','SXY392 ','SYZ392 ','SXZ392 ','EP392 ',
496 .'SX113 ','SY113 ','SZ113 ','SXY113 ','SYZ113 ','SXZ113 ',
497 .'EP113 ','SX213 ','SY213 ','SZ213 ','SXY213 ','SYZ213 ',
498 .'SXZ213 ','EP213 ','SX313 ','SY313 ','SZ113 ','SXY113 ',
499 .'SYZ313 ','SXZ113 ','EP313 ','SX123 ','SY123 ','SZ123 ',
500 .'SXY123 ','SYZ123 ','SXZ123 ','EP123 ','SX223 ','SY223 ',
501 .'SZ223 ','SXY223 ','SYZ223 ','SXZ223 ','EP223 ','SX323 ',
502 .'SY323 ','SZ323 ','SXY323 ','SYZ323 ','SXZ323 ','EP323 ',
503 .'SX133 ','SY133 ','SZ133 ','SXY133 ','SYZ133 ','SXZ133 ',
504 .'EP133 ','SX233 ','SY233 ','SZ233 ','SXY233 ','SYZ333 ',
505 .'SXZ233 ','EP233 ','SX333 ','SY333 ','SZ333 ','SXY333 ',
506 .'SYZ333 ','SXZ333 ','EP333 ','SX143 ','SY143 ','SZ143 ',
507 .'SXY143 ','SYZ143 ','SXZ143 ','EP143 ','SX243 ','SY243 ',
508 .'SZ243 ','SXY243 ','SYZ243 ','SXZ243 ','EP243 ','SX343 ',
509 .'SY343 ','SZ343 ','SXY343 ','SYZ343 ','SXZ343 ','EP343 ',
510 .'SX153 ','SY153 ','SZ153 ','SXY153 ','SYZ153 ','SXZ153 ',
511 .'EP153 ','SX253 ','SY253 ','SZ253 ','SXY253 ','SYZ253 ',
512 .'SXZ253 ','EP253 ','SX353 ','SY353 ','SZ353 ','SXY353 ',
513 .'SYZ353 ','SXZ353 ','EP353 ','SX163 ','SY163 ','SZ163 ',
514 .'SXY163 ','SYZ163 ','SXZ163 ','EP163 ','SX263 ','SY263 ',
515 .'SZ263 ','SXY263 ','SYZ263 ','SXZ263 ','EP263 ','SX363 ',
516 .'SY363 ','SZ363 ','SXY363 ','SYZ363 ','SXZ363 ','EP363 ',
517 .'SX173 ','SY173 ','SZ173 ','SXY173 ','SYZ173 ','SXZ173 ',
518 .'EP173 ','SX273 ','SY273 ','SZ273 ','SXY273 ','SYZ273 ',
519 .'SXZ273 ','EP273 ','SX373 ','SY373 ','SZ373 ','SXY373 ',
520 .'SYZ373 ','SXZ373 ','EP373 ','SX183 ','SY183 ','SZ183 ',
521 .'SXY183 ','SYZ183 ','SXZ183 ','EP183 ','SX283 ','SY283 ',
522 .'SZ283 ','SXY283 ','SYZ283 ','SXZ283 ','EP283 ','SX383 ',
523 .'SY383 ','SZ383 ','SXY383 ','SYZ383 ','SXZ383 ','EP383 ',
524 .'SX193 ','SY193 ','SZ193 ','SXY193 ','SYZ193 ','SXZ193 ',
525 .'EP193 ','SX293 ','SY293 ','SZ293 ','SXY293 ','SYZ293 ',
526 .'SXZ293 ','EP293 ','SX393 ','SY393 ','SZ393 ','SXY393 ',
527 .'SYZ393 ','SXZ393 ','EP393 '/
528 DATA vars3/
529 .'SX1U1 ','SY1U1 ','SZ1U1 ','SXY1U1 ','SYZ1U1 ','SXZ1U1 ',
530 .'EP1U1 ','SX2U1 ','SY2U1 ','SZ2U1 ','SXY2U1 ','SYZ2U1 ',
531 .'SXZ2U1 ','EP2U1 ','SX3U1 ','SY3U1 ','SZ3U1 ','SXY3U1 ',
532 .'SYZ3U1 ','SXZ3U1 ','EP3U1 ','SX1U2 ','SY1U2 ','SZ1U2 ',
533 .'SXY1U2 ','SYZ1U2 ','SXZ1U2 ','EP1U2 ','SX2U2 ','SY2U2 ',
534 .'SZ2U2 ','SXY2U2 ','SYZ2U2 ','SXZ2U2 ','EP2U2 ','SX3U2 ',
535 .'SY3U2 ','SZ3U2 ','SXY3U2 ','SYZ3U2 ','SXZ3U2 ','EP3U2 ',
536 .'SX1U3 ','SY1U3 ','SZ1U3 ','SXY1U3 ','SYZ1U3 ','SXZ1U3 ',
537 .'EP1U3 ','SX2U3 ','SY2U3 ','SZ2U3 ','SXY2U3 ','SYZ3U3 ',
538 .'SXZ2U3 ','EP2U3 ','SX3U3 ','SY3U3 ','SZ3U3 ','SXY3U3 ',
539 .'SYZ3U3 ','SXZ3U3 ','EP3U3 ','SX1D1 ','SY1D1 ','SZ1D1 ',
540 .'SXY1D1 ','SYZ1D1 ','SXZ1D1 ','EP1D1 ','SX2D1 ','SY2U1 ',
541 .'sz2d1 ','sxy2d1 ','syz2d1 ','sxz2d1 ','ep2d1 ','sx3d1 ',
542 .'sy3d1 ','sz3d1 ','sxy3d1 ','syz3d1 ','sxz3d1 ','ep3d1 ',
543 .'sx1d2 ','sy1d2 ','sz1d2 ','sxy1d2 ','syz1d2 ','sxz1d2 ',
544 .'ep1d2 ','sx2d2 ','sy2d2 ','sz2d2 ','sxy2d2 ','syz2d2 ',
545 .'sxz2d2 ','ep2d2 ','sx3d2 ','sy3d2 ','sz3d2 ','sxy3d2 ',
546 .'syz3d2 ','sxz3d2 ','ep3d2 ','sx1d3 ','sy1d3 ','sz1d3 ',
547 .'sxy1d3 ','syz1d3 ','sxz1d3 ','ep1d3 ','sx2d3 ','sy2d3 ',
548 .'sz2d3 ','sxy2d3 ','syz3d3 ','sxz2d3 ','ep2d3 ','sx3d3 ',
549 .'sy3d3 ','sz3d3 ','sxy3d3 ','syz3d3 ','sxz3d3 ','ep3d3 ',
550 .'usr1_111','usr2_111','usr3_111','usr4_111','usr5_111','usr6_111',
551 .'usr7_111','usr8_111','usr9_111','usr1_211','usr2_211','usr3_211',
552 .'usr4_211','usr5_211','usr6_211','usr7_211','usr8_211','usr9_211',
553 .'usr1_311','usr3_311','usr3_311','usr4_311','usr5_311','usr6_311',
554 .'usr7_311','usr8_311','usr9_311','usr1_121','usr2_121','usr3_121',
555 .'usr4_121','usr5_121','usr6_121','usr7_121','usr8_121','usr9_121',
556 .'usr1_221','usr2_221','usr3_221','usr4_221','usr5_221','usr6_221',
557 .'usr7_221','usr8_221','usr9_221','usr1_321','usr3_321','usr3_321',
558 .'usr4_321','usr5_321','usr6_321','usr7_321','usr8_321','usr9_321',
559 .'usr1_131','usr2_131','usr3_131','usr4_131','usr5_131','usr6_131',
560 .'usr7_131','usr8_131','usr9_131','usr1_231','usr2_231','usr3_231',
561 .'usr4_231','usr5_231','usr6_231','usr7_231','usr8_231','usr9_231',
562 .'usr1_331','usr3_331','usr3_331','usr4_331','usr5_331','usr6_331',
563 .'usr7_331','usr8_331','usr9_331','usr1_141','usr2_141','usr3_141',
564 .'usr4_141','usr5_141','usr6_141','usr7_141','usr8_141','usr9_141',
565 .'usr1_241','usr2_241','usr3_241','usr4_241','usr5_241','usr6_241',
566 .'usr7_241','usr8_241','usr9_241','usr1_341','usr3_341','usr3_341',
567 .'usr4_341','usr5_341','usr6_341','usr7_341','usr8_341','usr9_341',
568 .'usr1_151','usr2_151','usr3_151','usr4_151','usr5_151','usr6_151',
569 .'usr7_151','usr8_151','usr9_151','usr1_251','usr2_251','usr3_251',
570 .'USR4_251','USR5_251','USR6_251','USR7_251','USR8_251','USR9_251',
571 .'USR1_351','USR3_351','USR3_351','USR4_351','USR5_351','USR6_351',
572 .'USR7_351','USR8_351','USR9_351','USR1_161','USR2_161','USR3_161',
573 .'USR4_161','USR5_161','USR6_161','USR7_161','USR8_161','USR9_161',
574 .'USR1_261','USR2_261','USR3_261','USR4_261','USR5_261','USR6_261',
575 .'USR7_261','USR8_261','USR9_261','USR1_361','USR3_361','USR3_361',
576 .'USR4_361','USR5_361','USR6_361','USR7_361','USR8_361','USR9_361',
577 .'USR1_171','USR2_171','USR3_171','USR4_171','USR5_171','USR6_171',
578 .'USR7_171','USR8_171','USR9_171','USR1_271','USR2_271','USR3_271',
579 .'USR4_271','USR5_271','USR6_271','USR7_271','USR8_271','USR9_271',
580 .'USR1_371','USR3_371','USR3_371','USR4_371','USR5_371','USR6_371',
581 .'USR7_371','USR8_371','USR9_371','USR1_181','USR2_181','USR3_181',
582 .'USR4_181','USR5_171','USR6_171','USR7_181','USR8_181','USR9_181',
583 .'USR1_281','USR2_281','USR3_281','USR4_281','USR5_271','USR6_271',
584 .'USR7_281','USR8_281','USR9_281','USR1_381','USR3_381','USR3_381',
585 .'USR4_381','USR5_371','USR6_371','USR7_381','USR8_381','USR9_381',
586 .'USR1_191','USR2_191','USR3_191','USR4_191','USR5_191','USR6_191',
587 .'USR7_191','USR8_191','USR9_191','USR1_291','USR2_291','USR3_291',
588 .'USR4_291','USR5_291','USR6_291','USR7_291','USR8_291','USR9_291',
589 .'USR1_391','USR3_391','USR3_391','USR4_391','USR5_391','USR6_391',
590 .'USR7_391','USR8_391','USR9_391'/
591 DATA vars4/
592 .'USR1_112','USR2_112','USR3_112','USR4_112','USR5_112','USR6_112',
593 .'USR7_112','USR8_112','USR9_112',
594 .'USR1_212','USR2_212','USR3_212','USR4_212','USR5_212','USR6_212',
595 .'USR7_212','USR8_212','USR9_212','USR1_312','USR3_312','USR3_312',
596 .'USR4_312','USR5_312','USR6_312','USR7_312','USR8_312','USR9_312',
597 .'USR1_122','USR2_122','USR3_122','USR4_122','USR5_122','USR6_122',
598 .'USR7_122','USR8_122','USR9_122','USR1_222','USR2_222','USR3_222',
599 .'USR4_222','USR5_222','USR6_222','USR7_222','USR8_222','USR9_222',
600 .'USR1_322','USR3_322','USR3_322','USR4_322','USR5_322','USR6_322',
601 .'USR7_322','USR8_322','USR9_322','USR1_132','USR2_132','USR3_132',
602 .'USR4_132','USR5_132','USR6_132','USR7_132','USR8_132','USR9_132',
603 .'USR1_232','USR2_232','USR3_232','USR4_232','USR5_232','USR6_232',
604 .'USR7_232','USR8_232','USR9_232','USR1_332','USR3_332','USR3_332',
605 .'USR4_332','USR5_332','USR6_332','USR7_332','USR8_332','USR9_332',
606 .'USR1_142','USR2_142','USR3_142','USR4_142','USR5_142','USR6_142',
607 .'usr7_142','usr8_142','usr9_142','usr1_242','usr2_242','usr3_242',
608 .'usr4_242','usr5_242','usr6_242','usr7_242','usr8_242','usr9_242',
609 .'usr1_342','usr3_342','usr3_342','usr4_342','usr5_342','usr6_342',
610 .'usr7_342','usr8_342','usr9_342','usr1_152','usr2_152','usr3_152',
611 .'usr4_152','usr5_152','usr6_152','usr7_152','usr8_152','usr9_152',
612 .'usr1_252','usr2_252','usr3_252','usr4_252','usr5_252','usr6_252',
613 .'usr7_252','usr8_252','usr9_252','usr1_352','usr3_352','usr3_352',
614 .'usr4_352','usr5_352','usr6_352','usr7_352','usr8_352','usr9_352',
615 .'usr1_162','usr2_162','usr3_162','usr4_162','usr5_162','usr6_162',
616 .'usr7_162','usr8_162','usr9_162','usr1_262','usr2_262','usr3_262',
617 .'usr4_262','usr5_262','usr6_262','usr7_262','usr8_262','usr9_262',
618 .'usr1_362','usr3_362','usr3_362','usr4_362','usr5_362','usr6_362',
619 .'usr7_362','usr8_362','usr9_362','usr1_172','usr2_172','usr3_172',
620 .'usr4_172','usr5_172','usr6_172','usr7_172','usr8_172','usr9_172',
621 .'usr1_272','usr2_272','usr3_272','usr4_272','usr5_272','usr6_272',
622 .'usr7_272','usr8_272','usr9_272','usr1_372','usr3_372','usr3_372',
623 .'usr4_372','usr5_372','usr6_372','usr7_372','usr8_372','USR9_372',
624 .'USR1_182','USR2_182','USR3_182','USR4_182','USR5_182','USR6_182',
625 .'USR7_182','USR8_182','USR9_182','USR1_282','USR2_282','USR3_282',
626 .'USR4_282','USR5_282','USR6_282','USR7_282','USR8_282','USR9_282',
627 .'USR1_382','USR3_382','USR3_382','USR4_382','USR5_382','USR6_382',
628 .'USR7_382','USR8_382','USR9_382','USR1_192','USR2_192','USR3_192',
629 .'USR4_192','USR5_192','USR6_192','USR7_192','usr8_192','usr9_192',
630 .'usr1_292','usr2_292','usr3_292','usr4_292','usr5_292','usr6_292',
631 .'usr7_292','usr8_292','usr9_292','usr1_392','usr3_392','usr3_392',
632 .'usr4_392','usr5_392','usr6_392','usr7_392','usr8_392','usr9_392',
633 .'usr1_113','usr2_113','usr3_113','usr4_113','usr5_113','usr6_113',
634 .'usr7_113','usr8_113','usr9_113','usr1_213','usr2_213','usr3_213',
635 .'usr4_213','usr5_213','usr6_213','usr7_213','usr8_213','usr9_213',
636 .'usr1_313','usr3_313','usr3_313','usr4_313','usr5_313','usr6_313',
637 .'usr7_313','usr8_313','usr9_313','usr1_123','usr2_123','usr3_123',
638 .'usr4_123','usr5_123','usr6_123','usr7_123','usr8_123','usr9_123',
639 .'usr1_223','usr2_223','usr3_223','usr4_223','usr5_223','usr6_223',
640 .'usr7_223','usr8_223','usr9_223','usr1_323','usr3_323','usr3_323',
641 .'usr4_323','usr5_323','usr6_323','usr7_323','usr8_323','usr9_323',
642 .'usr1_133','usr2_133','usr3_133','usr4_133','usr5_133','usr6_133',
643 .'usr7_133','usr8_133','usr9_133','usr1_233','usr2_233','usr3_233',
644 .'usr4_233','usr5_233','usr6_233','usr7_233','usr8_233','usr9_233',
645 .'usr1_333','usr3_333','usr3_333','usr4_333','usr5_333','usr6_333',
646 .'usr7_333','usr8_333','usr9_333','usr1_143','usr2_143','usr3_143',
647 .'usr4_143','usr5_143','usr6_143','usr7_143','usr8_143','usr9_143',
648 .'usr1_243','usr2_243','usr3_243','usr4_243','usr5_243','usr6_243',
649 .'usr7_243','usr8_243','usr9_243','usr1_343','usr3_343','usr3_343',
650 .'usr4_343','usr5_343','usr6_343','usr7_343','usr8_343','usr9_343',
651 .'usr1_153','usr2_153','usr3_153','usr4_153','usr5_153','usr6_153',
652 .'usr7_153','usr8_153','usr9_153','usr1_253','usr2_253','usr3_253',
653 .'usr4_253','usr5_253','usr6_253','usr7_253','usr8_253','usr9_253',
654 .'usr1_353','USR3_353','USR3_353','USR4_353','USR5_353','USR6_353',
655 .'USR7_353','USR8_353','USR9_353','USR1_163','USR2_163','USR3_163',
656 .'USR4_163','USR5_163','USR6_163','USR7_163','USR8_163','USR9_163',
657 .'USR1_263','USR2_263','USR3_263','USR4_263','USR5_263','USR6_263',
658 .'USR7_263','USR8_263','USR9_263','USR1_363','USR3_363','USR3_363',
659 .'USR4_363','USR5_363','USR6_363','USR7_363','USR8_363','USR9_363',
660 .'USR1_173','USR2_173','USR3_173','USR4_173','USR5_173','USR6_173',
661 .'USR7_173','USR8_173','USR9_173','USR1_273','USR2_273','USR3_273',
662 .'USR4_273','USR5_273','USR6_273','USR7_273','USR8_273','USR9_273',
663 .'USR1_373','USR3_373','USR3_373','USR4_373','USR5_373','USR6_373',
664 .'USR7_373','USR8_373','USR9_373','USR1_183','USR2_183','USR3_183',
665 .'USR4_183','USR5_183','USR6_183','USR7_183','USR8_183','USR9_183',
666 .'USR1_283','USR2_283','USR3_283','USR4_283','USR5_283','USR6_283',
667 .'USR7_283','USR8_283','USR9_283','USR1_383','USR3_383','USR3_383',
668 .'USR4_383','USR5_383','USR6_383','USR7_383','USR8_383','USR9_383',
669 .'USR1_193','USR2_193','USR3_193','USR4_193','USR5_193','USR6_193',
670 .'USR7_193','USR8_193','USR9_193','USR1_293','USR2_293','USR3_293',
671 .'USR4_293','USR5_293','USR6_293','USR7_293','USR8_293','USR9_293',
672 .'USR1_393','USR3_393','USR3_393','USR4_393','USR5_393','USR6_393',
673 .'USR7_393','USR8_393','USR9_393',
674 .'EPSXX ','EPSYY ','EPSZZ ','EPSXY ','EPSZX ','EPSYZ '/
675 DATA vars5/
676 .'SX_MAT1 ','SY_MAT1 ','SZ_MAT1 ','SXY_MAT1','SYZ_MAT1','SXZ_MAT1',
677 .'IE_MAT1 ','DEN_MAT1','VOL_MAT1','PLA_MAT1','TEM_MAT1',
678 .'SX_MAT2 ','SY_MAT2 ','SZ_MAT2 ','SXY_MAT2','SYZ_MAT2','SXZ_MAT2',
679 .'IE_MAT2 ','DEN_MAT2','VOL_MAT2','PLA_MAT2','TEM_MAT2'/
680C Definition of VARS6 : strain for /PROP/TYPE22
681C=======================================================================
682 DO j=1,200
683 IF (j <= 9) THEN
684 WRITE(chj,'(I1.1)')j
685 ELSEIF (j <= 99) THEN
686 WRITE(chj,'(I2.2)')j
687 ELSE
688 WRITE(chj,'(I3.3)')j
689 ENDIF
690 DO i=1,9
691 WRITE(chi,'(I1.1)')i
692 DO k=1,9
693 WRITE(chk,'(I1.1)')k
694 DO l=1,6
695 IF (l==1) var_tmp(1:5)='E11_'//chi
696 IF (l==2) var_tmp(1:5)='E12_'//chi
697 IF (l==3) var_tmp(1:5)='E13_'//chi
698 IF (l==4) var_tmp(1:5)='E22_'//chi
699 IF (l==5) var_tmp(1:5)='E23_'//chi
700 IF (l==6) var_tmp(1:5)='E33_'//chi
701 var_tmp=var_tmp(1:5)//chk//'_'//chj
702 vars6((j-1)*9*9*6+((i-1)*9+k-1)*6+l)=var_tmp
703 ENDDO
704 ENDDO
705 ENDDO
706 ENDDO
707C Definition of VARS7 : stress for /PROP/TYPE22
708 DO j=1,200
709 IF (j <= 9) THEN
710 WRITE(chj,'(I1.1)')j
711 ELSEIF (j <= 99) THEN
712 WRITE(chj,'(I2.2)')j
713 ELSE
714 WRITE(chj,'(I3.3)')j
715 ENDIF
716 DO i=1,9
717 WRITE(chi,'(I1.1)')i
718 DO k=1,9
719 WRITE(chk,'(I1.1)')k
720 DO l=1,6
721 IF (l==1) var_tmp(1:5)='S11_'//chi
722 IF (l==2) var_tmp(1:5)='S12_'//chi
723 IF (l==3) var_tmp(1:5)='S13_'//chi
724 IF (l==4) var_tmp(1:5)='S22_'//chi
725 IF (l==5) var_tmp(1:5)='S23_'//chi
726 IF (l==6) var_tmp(1:5)='S33_'//chi
727 var_tmp=var_tmp(1:5)//chk//'_'//chj
728 vars7((j-1)*9*9*6+((i-1)*9+k-1)*6+l)=var_tmp
729 ENDDO
730 ENDDO
731 ENDDO
732 ENDDO
733C Definition of VARS8 : strain matrix output
734
735 vars8(1) = 'LEPSX '
736 vars8(2) = 'LEPSY '
737 vars8(3) = 'LEPSZ '
738 vars8(4) = 'LEPSXY '
739 vars8(5) = 'LEPSYZ '
740 vars8(6) = 'LEPSZX '
741 vars8(7) = 'EPSX1 '
742 vars8(8) = 'EPSX2 '
743 vars8(9) = 'EPSX3 '
744 vars8(10) = 'EPSX4 '
745 vars8(11) = 'EPSY1 '
746 vars8(12) = 'EPSY2 '
747 vars8(13) = 'EPSY3 '
748 vars8(14) = 'EPSY4 '
749 vars8(15) = 'EPSZ1 '
750 vars8(16) = 'EPSZ2 '
751 vars8(17) = 'EPSZ3 '
752 vars8(18) = 'EPSZ4 '
753 vars8(19) = 'EPSXY1 '
754 vars8(20) = 'EPSXY2 '
755 vars8(21) = 'EPSXY3 '
756 vars8(22) = 'EPSXY4 '
757 vars8(23) = 'EPSYZ1 '
758 vars8(24) = 'EPSYZ2 '
759 vars8(25) = 'EPSYZ3 '
760 vars8(26) = 'EPSYZ4 '
761 vars8(27) = 'EPSZX1 '
762 vars8(28) = 'EPSZX2 '
763 vars8(29) = 'EPSZX3 '
764 vars8(30) = 'EPSZX4 '
765
766 DO k=1,3
767 DO j=1,9
768 DO i=1,3
769 WRITE(chi,'(I1.1)')i
770 WRITE(chjs,'(I1.1)')j
771 WRITE(chk,'(I1.1)')k
772 var_tmp = 'EPSX'
773 vars8(30+(k-1)*9*6*3+(j-1)*3*6+(i-1)*6+1)=var_tmp(1:4)//chi//chjs//chk
774 var_tmp = 'EPSY'
775 vars8(30+(k-1)*9*6*3+(j-1)*3*6+(i-1)*6+2)=var_tmp(1:4)//chi//chjs//chk
776 var_tmp = 'EPSZ'
777 vars8(30+(k-1)*9*6*3+(j-1)*3*6+(i-1)*6+3)=var_tmp(1:4)//chi//chjs//chk
778 var_tmp = 'EPSXY'
779 vars8(30+(k-1)*9*6*3+(j-1)*3*6+(i-1)*6+4)=var_tmp(1:5)//chi//chjs//chk
780 var_tmp = 'EPSYZ'
781 vars8(30+(k-1)*9*6*3+(j-1)*3*6+(i-1)*6+5)=var_tmp(1:5)//chi//chjs//chk
782 var_tmp = 'EPSZX'
783 vars8(30+(k-1)*9*6*3+(j-1)*3*6+(i-1)*6+6)=var_tmp(1:5)//chi//chjs//chk
784
785 ENDDO
786 ENDDO
787 ENDDO
788
789 DATA vars9/
790 . 'VX ','VY ','VZ ','SSP ','MACH ',
791 . 'VK '/
792
793 DATA varsnloc/
794 . 'NL_PLAS ','NL_PLSR '/
795
796 DATA vars10/
797 . 'VSTRAIN ' /
798
799
800C shells
801c VARC = VARC1 + VARC2 + VARC3
802 DATA varc1/
803 . 'F1 ','F2 ','F12 ','Q1 ','Q2 ',
804 . 'M1 ','M2 ','M12 ','IEM ','IEB ',
805 . 'OFF ','THIC ','EMIN ','EMAX ','E1 ',
806 . 'E2 ','E12 ','SH1 ','SH2 ','K1 ',
807 .'K2 ','K12 ','EPSD ',
808 .'USR1 ','USR2 ','USR3 ','USR4 ','USR5 ','USR6 ',
809 .'USR7 ','USR8 ','USR9 ','USR10 ','USR11 ','USR12 ',
810 .'USR13 ','USR14 ','USR15 ','USR16 ','USR17 ','USR18 ',
811 .'USR19 ','USR20 ','USR21 ','USR22 ','USR23 ','USR24 ',
812 .'USR25 ','USR26 ','USR27 ','USR28 ','USR29 ','USR30 ',
813 .'USR31 ','USR32 ','USR33 ','USR34 ','USR35 ','USR36 ',
814 .'USR37 ','USR38 ','USR39 ','USR40 ','USR41 ','USR42 ',
815 .'USR43 ','USR44 ','USR45 ','USR46 ','USR47 ','USR48 ',
816 .'USR49 ','USR50 ','USR51 ','USR52 ','USR53 ','USR54 ',
817 .'USR55 ','USR56 ','USR57 ','USR58 ','USR59 ','USR60 ',
818 .'USR1_1 ','USR1_2 ','USR1_3 ','USR1_4 ','USR1_5 ','USR2_1 ',
819 .'USR2_2 ','USR2_3 ','USR2_4 ','USR2_5 ','USR3_1 ','USR3_2 ',
820 .'USR3_3 ','USR3_4 ','USR3_5 ','USR4_1 ','USR4_2 ','USR4_3 ',
821 .'USR4_4 ','USR4_5 ','USR5_1 ','USR5_2 ','usr5_3 ','usr5_4 ',
822 .'usr5_5 ','usr6_1 ','usr6_2 ','usr6_3 ','usr6_4 ','usr6_5 ',
823 .'usr7_1 ','usr7_2 ','usr7_3 ','usr7_4 ','usr7_5 ','usr8_1 ',
824 .'usr8_2 ','usr8_3 ','USR8_4 ','USR8_5 ','USR9_1 ','USR9_2 ',
825 .'USR9_3 ','USR9_4 ','USR9_5 ','USR10_1 ','USR10_2 ','USR10_3 ',
826 .'USR10_4 ','USR10_5 ','USR11_1 ','USR11_2 ','USR11_3 ','USR11_4 ',
827 .'USR11_5 ','USR12_1 ','USR12_2 ','USR12_3 ','USR12_4 ','USR12_5 ',
828 .'USR13_1 ','USR13_2 ','USR13_3 ','USR13_4 ','USR13_5 ','USR14_1 ',
829 .'USR14_2 ','USR14_3 ','USR14_4 ','USR14_5 ','USR15_1 ','USR15_2 ',
830 .'USR15_3 ','USR15_4 ','USR15_5 ','USR16_1 ','USR16_2 ','USR16_3 ',
831 .'USR16_4 ','USR16_5 ','USR17_1 ','USR17_2 ','USR17_3 ','USR17_4 ',
832 .'USR17_5 ','USR18_1 ','USR18_2 ','USR18_3 ','USR18_4 ','USR18_5 ',
833 .'USR19_1 ','USR19_2 ','USR19_3 ','USR19_4 ','USR19_5 ','USR20_1 ',
834 .'USR20_2 ','USR20_3 ','USR20_4 ','USR20_5 ','SX_1 ','SY_1 ',
835 .'SXY_1 ','SYZ_1 ','SZX_1 ','SX_2 ','SY_2 ','SXY_2 ',
836 .'SYZ_2 ','SZX_2 ','SX_3 ','SY_3 ','SXY_3 ','SYZ_3 ',
837 .'SZX_3 ','SX_4 ','SY_4 ','SXY_4 ','SYZ_4 ','SZX_4 ',
838 .'SX_5 ','SY_5 ','SXY_5 ','SYZ_5 ','SZX_5 '/
839 DATA varc2/
840 .'SX_6 ','SY_6 ','SXY_6 ','SYZ_6 ','SZX_6 ',
841 .'SX_7 ','SY_7 ','SXY_7 ','SYZ_7 ','SZX_7 ','SX_8 ',
842 .'SY_8 ','SXY_8 ','SYZ_8 ','SZX_8 ','SX_9 ','SY_9 ',
843 .'SXY_9 ','SYZ_9 ','SZX_9 ','SX_10 ','SY_10 ','SXY_10 ',
844 .'SYZ_10 ','SZX_10 ','SX_11 ','SY_11 ','SXY_11 ','SYZ_11 ',
845 .'SZX_11 ','SX_12 ','SY_12 ','SXY_12 ','SYZ_12 ','SZX_12 ',
846 .'SX_13 ','SY_13 ','SXY_13 ','SYZ_13 ','SZX_13 ','SX_14 ',
847 .'SY_14 ','SXY_14 ','SYZ_14 ','SZX_14 ','SX_15 ','SY_15 ',
848 .'SXY_15 ','SYZ_15 ','SZX_15 ','SX_16 ','SY_16 ','SXY_16 ',
849 .'SYZ_16 ','SZX_16 ','SX_17 ','SY_17 ','SXY_17 ','SYZ_17 ',
850 .'SZX_17 ','SX_18 ','SY_18 ','SXY_18 ','SYZ_18 ','SZX_18 ',
851 .'SX_19 ','SY_19 ','SXY_19 ','SYZ_19 ','SZX_19 ','SX_20 ',
852 .'SY_20 ','SXY_20 ','SYZ_20 ','SZX_20 ','SX_21 ','SY_21 ',
853 .'SXY_21 ','SYZ_21 ','SZX_21 ','SX_22 ','SY_22 ','SXY_22 ',
854 .'SYZ_22 ','SZX_22 ','SX_23 ','SY_23 ','SXY_23 ','SYZ_23 ',
855 .'SZX_23 ','SX_24 ','SY_24 ','SXY_24 ','SYZ_24 ','SZX_24 ',
856 .'SX_25 ','SY_25 ','SXY_25 ','SYZ_25 ','SZX_25 ','SX_26 ',
857 .'SY_26 ','SXY_26 ','SYZ_26 ','SZX_26 ','SX_27 ','SY_27 ',
858 .'SXY_27 ','SYZ_27 ','SZX_27 ','SX_28 ','SY_28 ','SXY_28 ',
859 .'SYZ_28 ','SZX_28 ','SX_29 ','SY_29 ','SXY_29 ','SYZ_29 ',
860 .'SZX_29 ','SX_30 ','SY_30 ','SXY_30 ','SYZ_30 ','SZX_30 ',
861 .'SX_31 ','SY_31 ','SXY_31 ','SYZ_31 ','SZX_31 ','SX_32 ',
862 .'SY_32 ','SXY_32 ','SYZ_32 ','SZX_32 ','SX_33 ','SY_33 ',
863 .'SXY_33 ','SYZ_33 ','SZX_33 ','SX_34 ','SY_34 ','SXY_34 ',
864 .'SYZ_34 ','SZX_34 ','SX_35 ','SY_35 ','SXY_35 ','SYZ_35 ',
865 .'SZX_35 ','SX_36 ','SY_36 ','SXY_36 ','SYZ_36 ','SZX_36 ',
866 .'SX_37 ','SY_37 ','SXY_37 ','SYZ_37 ','SZX_37 ','SX_38 ',
867 .'SY_38 ','SXY_38 ','SYZ_38 ','SZX_38 ','SX_39 ','SY_39 ',
868 .'SXY_39 ','SYZ_39 ','SZX_39 ','SX_40 ','SY_40 ','SXY_40 ',
869 .'SYZ_40 ','SZX_40 ','SX_41 ','SY_41 ','SXY_41 ','SYZ_41 ',
870 .'SZX_41 ','SX_42 ','SY_42 ','SXY_42 ','SYZ_42 ','SZX_42 ',
871 .'SX_43 ','SY_43 ','SXY_43 ','SYZ_43 ','SZX_43 ','SX_44 ',
872 .'SY_44 ','SXY_44 ','SYZ_44 ','SZX_44 ','SX_45 ','SY_45 ',
873 .'SXY_45 ','SYZ_45 ','SZX_45 ','SX_46 ','SY_46 ','SXY_46 ',
874 .'SYZ_46 ','SZX_46 ','SX_47 ','SY_47 ','SXY_47 ','SYZ_47 ',
875 .'SZX_47 ','SX_48 ','SY_48 ','SXY_48 ','SYZ_48 ','SZX_48 ',
876 .'SX_49 ','SY_49 ','SXY_49 ','SYZ_49 ','SZX_49 ','SX_50 ',
877 .'SY_50 ','SXY_50 ','SYZ_50 ','SZX_50 ','SX_51 ','SY_51 ',
878 .'SXY_51 ','SYZ_51 ','SZX_51 ','SX_52 ','SY_52 ','SXY_52 ',
879 .'SYZ_52 ','SZX_52 ','SX_53 ','SY_53 ','SXY_53 ','SYZ_53 ',
880 .'SZX_53 ','SX_54 ','SY_54 ','SXY_54 ','SYZ_54 ','SZX_54 ',
881 .'SX_55 ','SY_55 ','SXY_55 ','SYZ_55 ','SZX_55 ','SX_56 ',
882 .'SY_56 ','SXY_56 ','SYZ_56 ','SZX_56 ','SX_57 ','SY_57 ',
883 .'SXY_57 ','SYZ_57 ','SZX_57 ','SX_58 ','SY_58 ','SXY_58 ',
884 .'SYZ_58 ','SZX_58 ','SX_59 ','SY_59 ','SXY_59 ','SYZ_59 ',
885 .'SZX_59 ','SX_60 ','SY_60 ','SXY_60 ','SYZ_60 ','SZX_60 ',
886 .'SX_61 ','SY_61 ','SXY_61 ','SYZ_61 ','SZX_61 ','SX_62 ',
887 .'SY_62 ','SXY_62 ','SYZ_62 ','SZX_62 ','SX_63 ','SY_63 ',
888 .'SXY_63 ','SYZ_63 ','SZX_63 ','SX_64 ','SY_64 ','SXY_64 ',
889 .'SYZ_64 ','SZX_64 ','SX_65 ','SY_65 ','SXY_65 ','SYZ_65 ',
890 .'SZX_65 ','SX_66 ','SY_66 ','SXY_66 ','SYZ_66 ','SZX_66 ',
891 .'SX_67 ','SY_67 ','SXY_67 ','SYZ_67 ','SZX_67 ','SX_68 ',
892 .'SY_68 ','SXY_68 ','SYZ_68 ','SZX_68 ','SX_69 ','SY_69 ',
893 .'SXY_69 ','SYZ_69 ','SZX_69 ','SX_70 ','SY_70 ','SXY_70 ',
894 .'SYZ_70 ','SZX_70 ','SX_71 ','SY_71 ','SXY_71 ','SYZ_71 ',
895 .'SZX_71 ','SX_72 ','SY_72 ','SXY_72 ','SYZ_72 ','SZX_72 ',
896 .'SX_73 ','SY_73 ','SXY_73 ','SYZ_73 ','SZX_73 ','SX_74 ',
897 .'SY_74 ','SXY_74 ','SYZ_74 ','SZX_74 ','SX_75 ','SY_75 ',
898 .'SXY_75 ','SYZ_75 ','SZX_75 ','SX_76 ','SY_76 ','SXY_76 ',
899 .'SYZ_76 ','SZX_76 ','SX_77 ','SY_77 ','SXY_77 ','SYZ_77 ',
900 .'SZX_77 ','SX_78 ','SY_78 ','SXY_78 ','SYZ_78 ','SZX_78 ',
901 .'SX_79 ','SY_79 ','SXY_79 ','SYZ_79 ','SZX_79 ','SX_80 ',
902 .'SY_80 ','SXY_80 ','SYZ_80 ','SZX_80 ','SX_81 ','SY_81 ',
903 .'SXY_81 ','SYZ_81 ','SZX_81 ','SX_82 ','SY_82 ','SXY_82 ',
904 .'SYZ_82 ','SZX_82 ','SX_83 ','SY_83 ','SXY_83 ','SYZ_83 ',
905 .'SZX_83 ','SX_84 ','SY_84 ','SXY_84 ','SYZ_84 ','SZX_84 ',
906 .'SX_85 ','SY_85 ','SXY_85 ','SYZ_85 ','SZX_85 ','SX_86 ',
907 .'SY_86 ','SXY_86 ','SYZ_86 ','SZX_86 ','SX_87 ','SY_87 ',
908 .'SXY_87 ','SYZ_87 ','SZX_87 ','SX_88 ','SY_88 ','SXY_88 ',
909 .'SYZ_88 ','SZX_88 ','SX_89 ','SY_89 ','SXY_89 ','SYZ_89 ',
910 .'SZX_89 ','SX_90 ','SY_90 ','SXY_90 ','SYZ_90 ','SZX_90 ',
911 .'SX_91 ','SY_91 ','SXY_91 ','SYZ_91 ','SZX_91 ','SX_92 ',
912 .'SY_92 ','SXY_92 ','SYZ_92 ','SZX_92 ','SX_93 ','SY_93 ',
913 .'SXY_93 ','SYZ_93 ','SZX_93 ','SX_94 ','SY_94 ','SXY_94 ',
914 .'SYZ_94 ','SZX_94 ','SX_95 ','SY_95 ','SXY_95 ','SYZ_95 ',
915 .'SZX_95 ','SX_96 ','SY_96 ','SXY_96 ','SYZ_96 ','SZX_96 ',
916 .'SX_97 ','SY_97 ','SXY_97 ','SYZ_97 ','SZX_97 ','SX_98 ',
917 .'SY_98 ','SXY_98 ','SYZ_98 ','SZX_98 ','SX_99 ','SY_99 ',
918 .'SXY_99 ','SYZ_99 ','SZX_99 ','NFAIL ','PFAIL ','FAIL_D1 ',
919 .'FAIL_D2 ','FAIL_EN ','WPLAY01 ','WPLAY02 ','WPLAY03 ','WPLAY04 ',
920 .'WPLAY05 ','WPLAY06 ','WPLAY07 ','WPLAY08 ','WPLAY09 ','WPLAY10 ',
921 .'WPLAY11 ','WPLAY12 ','WPLAY13 ','WPLAY14 ','WPLAY15 ','WPLAY16 ',
922 .'WPLAY17 ','WPLAY18 ','WPLAY19 ','WPLAY20 ','WPLAY21 ','WPLAY22 ',
923 .'WPLAY23 ','WPLAY24 ','WPLAY25 ','WPLAY26 ','WPLAY27 ','WPLAY28 ',
924 .'WPLAY29 ','WPLAY30 ','WPLAY31 ','WPLAY32 ','WPLAY33 ','WPLAY34 ',
925 .'WPLAY35 ','WPLAY36 ','WPLAY37 ','WPLAY38 ','WPLAY39 ','WPLAY40 ',
926 .'WPLAY41 ','WPLAY42 ','WPLAY43 ','WPLAY44 ','WPLAY45 ','WPLAY46 ',
927 .'WPLAY47 ','WPLAY48 ','WPLAY49 ','WPLAY50 ','WPLAY51 ','WPLAY52 ',
928 .'WPLAY53 ','WPLAY54 ','WPLAY55 ','WPLAY56 ','WPLAY57 ','WPLAY58 ',
929 .'WPLAY59 ','WPLAY60 ','WPLAY61 ','WPLAY62 ','WPLAY63 ','WPLAY64 ',
930 .'WPLAY65 ','WPLAY66 ','WPLAY67 ','WPLAY68 ','WPLAY69 ','WPLAY70 ',
931 .'WPLAY71 ','WPLAY72 ','WPLAY73 ','WPLAY74 ','WPLAY75 ','WPLAY76 ',
932 .'WPLAY77 ','WPLAY78 ','WPLAY79 ','WPLAY80 ','WPLAY81 ','WPLAY82 ',
933 .'WPLAY83 ','WPLAY84 ','WPLAY85 ','WPLAY86 ','WPLAY87 ','WPLAY88 ',
934 .'WPLAY89 ','WPLAY90 ','WPLAY91 ','WPLAY92 ','WPLAY93 ','WPLAY94 ',
935 .'WPLAY95 ','WPLAY96 ','WPLAY97 ','WPLAY98 ','WPLAY99 '/
936 DATA varc3/
937 .'SVX_1 ','SVY_1 ','SVXY_1 ','SVYZ_1 ','SVZX_1 ','SVX_2 ',
938 .'SVY_2 ','SVXY_2 ','SVYZ_2 ','SVZX_2 ','SVX_3 ','SVY_3 ',
939 .'SVXY_3 ','SVYZ_3 ','SVZX_3 ','SVX_4 ','SVY_4 ','SVXY_4 ',
940 .'SVYZ_4 ','SVZX_4 ','SVX_5 ','SVY_5 ','SVXY_5 ','SVYZ_5 ',
941 .'SVZX_5 ','SVX_6 ','SVY_6 ','SVXY_6 ','SVYZ_6 ','SVZX_6 ',
942 .'SVX_7 ','SVY_7 ','SVXY_7 ','SVYZ_7 ','SVZX_7 ','SVX_8 ',
943 .'SVY_8 ','SVXY_8 ','SVYZ_8 ','SVZX_8 ','SVX_9 ','SVY_9 ',
944 .'SVXY_9 ','SVYZ_9 ','SVZX_9 ','SVX_10 ','SVY_10 ','SVXY_10 ',
945 .'SVYZ_10 ','SVZX_10 ','SVX_11 ','SVY_11 ','SVXY_11 ','SVYZ_11 ',
946 .'SVZX_11 ','SVX_12 ','SVY_12 ','SVXY_12 ','SVYZ_12 ','SVZX_12 ',
947 .'SVX_13 ','SVY_13 ','SVXY_13 ','SVYZ_13 ','SVZX_13 ','SVX_14 ',
948 .'SVY_14 ','SVXY_14 ','SVYZ_14 ','SVZX_14 ','SVX_15 ','SVY_15 ',
949 .'SVXY_15 ','SVYZ_15 ','SVZX_15 ','SVX_16 ','SVY_16 ','SVXY_16 ',
950 .'SVYZ_16 ','SVZX_16 ','SVX_17 ','SVY_17 ','SVXY_17 ','SVYZ_17 ',
951 .'SVZX_17 ','SVX_18 ','SVY_18 ','SVXY_18 ','SVYZ_18 ','SVZX_18 ',
952 .'SVX_19 ','SVY_19 ','SVXY_19 ','SVYZ_19 ','SVZX_19 ','SVX_20 ',
953 .'SVY_20 ','SVXY_20 ','SVYZ_20 ','SVZX_20 ','SVX_21 ','SVY_21 ',
954 .'SVXY_21 ','SVYZ_21 ','SVZX_21 ','SVX_22 ','SVY_22 ','SVXY_22 ',
955 .'SVYZ_22 ','SVZX_22 ','SVX_23 ','SVY_23 ','SVXY_23 ','SVYZ_23 ',
956 .'SVZX_23 ','SVX_24 ','SVY_24 ','SVXY_24 ','SVYZ_24 ','SVZX_24 ',
957 .'SVX_25 ','SVY_25 ','SVXY_25 ','SVYZ_25 ','SVZX_25 ','SVX_26 ',
958 .'SVY_26 ','SVXY_26 ','SVYZ_26 ','SVZX_26 ','SVX_27 ','SVY_27 ',
959 .'SVXY_27 ','SVYZ_27 ','SVZX_27 ','SVX_28 ','SVY_28 ','SVXY_28 ',
960 .'SVYZ_28 ','SVZX_28 ','SVX_29 ','SVY_29 ','SVXY_29 ','SVYZ_29 ',
961 .'SVZX_29 ','SVX_30 ','SVY_30 ','SVXY_30 ','SVYZ_30 ','svzx_30 ',
962 .'svx_31 ','svy_31 ','svxy_31 ','svyz_31 ','svzx_31 ','svx_32 ',
963 .'svy_32 ','svxy_32 ','svyz_32 ','SVZX_32 ','SVX_33 ','SVY_33 ',
964 .'SVXY_33 ','SVYZ_33 ','SVZX_33 ','SVX_34 ','SVY_34 ','SVXY_34 ',
965 .'SVYZ_34 ','SVZX_34 ','SVX_35 ','SVY_35 ','SVXY_35 ','SVYZ_35 ',
966 .'SVZX_35 ','SVX_36 ','SVY_36 ','SVXY_36 ','SVYZ_36 ','SVZX_36 ',
967 .'SVX_37 ','SVY_37 ','SVXY_37 ','SVYZ_37 ','SVZX_37 ','SVX_38 ',
968 .'SVY_38 ','SVXY_38 ','SVYZ_38 ','SVZX_38 ','SVX_39 ','SVY_39 ',
969 .'SVXY_39 ','SVYZ_39 ','SVZX_39 ','SVX_40 ','SVY_40 ','svxy_40 ',
970 .'svyz_40 ','svzx_40 ','svx_41 ','svy_41 ','svxy_41 ','svyz_41 ',
971 .'svzx_41 ','svx_42 ','svy_42 ','svxy_42 ','svyz_42 ','svzx_42 ',
972 .'svx_43 ','svy_43 ','svxy_43 ','svyz_43 ','svzx_43 ','svx_44 ',
973 .'svy_44 ','svxy_44 ','svyz_44 ','svzx_44 ','svx_45 ','svy_45 ',
974 .'svxy_45 ','svyz_45 ','svzx_45 ','svx_46 ','svy_46 ','svxy_46 ',
975 .'SVYZ_46 ','SVZX_46 ','SVX_47 ','SVY_47 ','SVXY_47 ','SVYZ_47 ',
976 .'SVZX_47 ','SVX_48 ','SVY_48 ','SVXY_48 ','SVYZ_48 ','SVZX_48 ',
977 .'SVX_49 ','SVY_49 ','SVXY_49 ','SVYZ_49 ','SVZX_49 ','SVX_50 ',
978 .'SVY_50 ','SVXY_50 ','SVYZ_50 ','SVZX_50 ','SVX_51 ','SVY_51 ',
979 .'SVXY_51 ','SVYZ_51 ','SVZX_51 ','SVX_52 ','SVY_52 ','SVXY_52 ',
980 .'SVYZ_52 ','SVZX_52 ','SVX_53 ','SVY_53 ','SVXY_53 ','SVYZ_53 ',
981 .'SVZX_53 ','SVX_54 ','SVY_54 ','SVXY_54 ','SVYZ_54 ','SVZX_54 ',
982 .'SVX_55 ','SVY_55 ','SVXY_55 ','SVYZ_55 ','SVZX_55 ','SVX_56 ',
983 .'SVY_56 ','SVXY_56 ','SVYZ_56 ','SVZX_56 ','SVX_57 ','SVY_57 ',
984 .'SVXY_57 ','SVYZ_57 ','SVZX_57 ','SVX_58 ','SVY_58 ','SVXY_58 ',
985 .'SVYZ_58 ','SVZX_58 ','SVX_59 ','SVY_59 ','SVXY_59 ','SVYZ_59 ',
986 .'SVZX_59 ','SVX_60 ','SVY_60 ','SVXY_60 ','SVYZ_60 ','SVZX_60 ',
987 .'SVX_61 ','SVY_61 ','SVXY_61 ','SVYZ_61 ','SVZX_61 ','SVX_62 ',
988 .'SVY_62 ','SVXY_62 ','SVYZ_62 ','SVZX_62 ','SVX_63 ','SVY_63 ',
989 .'SVXY_63 ','SVYZ_63 ','SVZX_63 ','SVX_64 ','SVY_64 ','SVXY_64 ',
990 .'SVYZ_64 ','SVZX_64 ','SVX_65 ','SVY_65 ','SVXY_65 ','SVYZ_65 ',
991 .'SVZX_65 ','SVX_66 ','SVY_66 ','SVXY_66 ','SVYZ_66 ','SVZX_66 ',
992 .'SVX_67 ','SVY_67 ','SVXY_67 ','SVYZ_67 ','SVZX_67 ','SVX_68 ',
993 .'SVY_68 ','SVXY_68 ','SVYZ_68 ','SVZX_68 ','SVX_69 ','SVY_69 ',
994 .'SVXY_69 ','SVYZ_69 ','SVZX_69 ','SVX_70 ','SVY_70 ','SVXY_70 ',
995 .'SVYZ_70 ','SVZX_70 ','SVX_71 ','SVY_71 ','SVXY_71 ','SVYZ_71 ',
996 .'SVZX_71 ','SVX_72 ','SVY_72 ','SVXY_72 ','SVYZ_72 ','SVZX_72 ',
997 .'SVX_73 ','SVY_73 ','SVXY_73 ','SVYZ_73 ','SVZX_73 ','SVX_74 ',
998 .'SVY_74 ','SVXY_74 ','SVYZ_74 ','SVZX_74 ','SVX_75 ','SVY_75 ',
999 .'SVXY_75 ','SVYZ_75 ','SVZX_75 ','SVX_76 ','SVY_76 ','SVXY_76 ',
1000 .'SVYZ_76 ','SVZX_76 ','SVX_77 ','SVY_77 ','SVXY_77 ','SVYZ_77 ',
1001 .'SVZX_77 ','SVX_78 ','SVY_78 ','SVXY_78 ','SVYZ_78 ','SVZX_78 ',
1002 .'SVX_79 ','SVY_79 ','SVXY_79 ','SVYZ_79 ','SVZX_79 ','SVX_80 ',
1003 .'SVY_80 ','SVXY_80 ','SVYZ_80 ','SVZX_80 ','SVX_81 ','SVY_81 ',
1004 .'SVXY_81 ','SVYZ_81 ','SVZX_81 ','SVX_82 ','SVY_82 ','SVXY_82 ',
1005 .'SVYZ_82 ','SVZX_82 ','SVX_83 ','SVY_83 ','SVXY_83 ','SVYZ_83 ',
1006 .'SVZX_83 ','SVX_84 ','SVY_84 ','SVXY_84 ','SVYZ_84 ','SVZX_84 ',
1007 .'SVX_85 ','SVY_85 ','SVXY_85 ','SVYZ_85 ','SVZX_85 ','SVX_86 ',
1008 .'SVY_86 ','SVXY_86 ','SVYZ_86 ','SVZX_86 ','SVX_87 ','SVY_87 ',
1009 .'SVXY_87 ','SVYZ_87 ','SVZX_87 ','SVX_88 ','SVY_88 ','SVXY_88 ',
1010 .'SVYZ_88 ','SVZX_88 ','SVX_89 ','SVY_89 ','SVXY_89 ','SVYZ_89 ',
1011 .'SVZX_89 ','SVX_90 ','SVY_90 ','SVXY_90 ','SVYZ_90 ','SVZX_90 ',
1012 .'SVX_91 ','SVY_91 ','SVXY_91 ','SVYZ_91 ','SVZX_91 ','SVX_92 ',
1013 .'SVY_92 ','SVXY_92 ','SVYZ_92 ','SVZX_92 ','SVX_93 ','SVY_93 ',
1014 .'SVXY_93 ','SVYZ_93 ','SVZX_93 ','SVX_94 ','SVY_94 ','SVXY_94 ',
1015 .'SVYZ_94 ','SVZX_94 ','SVX_95 ','SVY_95 ','SVXY_95 ','SVYZ_95 ',
1016 .'SVZX_95 ','SVX_96 ','SVY_96 ','SVXY_96 ','SVYZ_96 ','SVZX_96 ',
1017 .'SVX_97 ','SVY_97 ','SVXY_97 ','SVYZ_97 ','SVZX_97 ','SVX_98 ',
1018 .'SVY_98 ','SVXY_98 ','SVYZ_98 ','SVZX_98 ','SVX_99 ','SVY_99 ',
1019 .'SVXY_99 ','SVYZ_99 ','SVZX_99 ','SVX_100 ','SVY_100 ','SVXY_100',
1020 .'SVYZ_100','SVZX_100'/
1021 DATA varcpinch/
1022 .'EPINCHXZ','EPINCHYZ','EPINCHZZ','FPINCH ','MPINCHXZ','MPINCHYZ',
1023 .'THK '/
1024 DATA varcnloc/
1025 .'NL_PLAS ','NL_EPSD '/
1026C trusses
1027 DATA vart/
1028 . 'OFF ','F ','IE ','A ','L ',
1029 . 'PLAS '/
1030C beams
1031 DATA varp/
1032 . 'OFF ','F1 ','F2 ','F3 ','M1 ',
1033 . 'M2 ','M3 ','IE ','SX ','SXY ',
1034 . 'SZX ','SX_1 ','SXY_1 ','SZX_1 ','SX_2 ',
1035 . 'SXY_2 ','SZX_2 ','SX_3 ','SXY_3 ','SZX_3 ',
1036 . 'SX_4 ','SXY_4 ','SZX_4 ','SX_5 ','SXY_5 ',
1037 . 'SZX_5 ','SX_6 ','SXY_6 ','SZX_6 ','SX_7 ',
1038 . 'SXY_7 ','SZX_7 ','SX_8 ','SXY_8 ','SZX_8 ',
1039 . 'SX_9 ','SXY_9 ','SZX_9 ','SX_10 ','SXY_10 ',
1040 . 'SZX_10 ','SX_11 ','SXY_11 ','SZX_11 ','SX_12 ',
1041 . 'SXY_12 ','SZX_12 ','SX_13 ','SXY_13 ','SZX_13 ',
1042 . 'SX_14 ','SXY_14 ','SZX_14 ','SX_15 ','SXY_15 ',
1043 . 'SZX_15 ','SX_16 ','SXY_16 ','SZX_16 ','SX_17 ',
1044 . 'sxy_17 ','szx_17 ','sx_18 ','sxy_18 ','szx_18 ',
1045 . 'sx_19 ','sxy_19 ','szx_19 ','sx_20 ','sxy_20 ',
1046 . 'szx_20 ','sx_21 ','sxy_21 ','szx_21 ','sx_22 ',
1047 . 'sxy_22 ','szx_22 ','sx_23 ','sxy_23 ','szx_23 ',
1048 . 'sx_24 ','sxy_24 ','szx_24 ','sx_25 ','sxy_25 ',
1049 . 'szx_25 ','sx_26 ','sxy_26 ','szx_26 ','sx_27 ',
1050 . 'sxy_27 ','szx_27 ','SX_28 ','SXY_28 ','SZX_28 ',
1051 . 'SX_29 ','SXY_29 ','SZX_29 ','SX_30 ','SXY_30 ',
1052 . 'SZX_30 ','SX_31 ','SXY_31 ','SZX_31 ','SX_32 ',
1053 . 'SXY_32 ','SZX_32 ','SX_33 ','SXY_33 ','SZX_33 ',
1054 . 'SX_34 ','SXY_34 ','SZX_34 ','SX_35 ','SXY_35 ',
1055 . 'SZX_35 ','SX_36 ','SXY_36 ','SZX_36 ','SX_37 ',
1056 . 'SXY_37 ','SZX_37 ','SX_38 ','SXY_38 ','SZX_38 ',
1057 . 'SX_39 ','SXY_39 ','SZX_39 ','SX_40 ','SXY_40 ',
1058 . 'SZX_40 ','SX_41 ','SXY_41 ','SZX_41 ','SX_42 ',
1059 . 'SXY_42 ','SZX_42 ','SX_43 ','SXY_43 ','SZX_43 ',
1060 . 'SX_44 ','SXY_44 ','SZX_44 ','SX_45 ','SXY_45 ',
1061 . 'SZX_45 ','SX_46 ','SXY_46 ','SZX_46 ','SX_47 ',
1062 . 'SXY_47 ','SZX_47 ','SX_48 ','SXY_48 ','SZX_48 ',
1063 . 'SX_49 ','SXY_49 ','SZX_49 ','SX_50 ','SXY_50 ',
1064 . 'SZX_50 ','SX_51 ','SXY_51 ','SZX_51 ','SX_52 ',
1065 . 'SXY_52 ','SZX_52 ','SX_53 ','SXY_53 ','SZX_53 ',
1066 . 'SX_54 ','SXY_54 ','SZX_54 ','SX_55 ','SXY_55 ',
1067 . 'SZX_55 ','SX_56 ','SXY_56 ','SZX_56 ','SX_57 ',
1068 . 'SXY_57 ','SZX_57 ','SX_58 ','SXY_58 ','SZX_58 ',
1069 . 'SX_59 ','SXY_59 ','SZX_59 ','SX_60 ','SXY_60 ',
1070 . 'SZX_60 ','SX_61 ','SXY_61 ','SZX_61 ','SX_62 ',
1071 . 'SXY_62 ','SZX_62 ','SX_63 ','SXY_63 ','SZX_63 ',
1072 . 'sx_64 ','sxy_64 ','szx_64 ','sx_65 ','sxy_65 ',
1073 . 'szx_65 ','sx_66 ','sxy_66 ','szx_66 ','sx_67 ',
1074 . 'sxy_67 ','szx_67 ','sx_68 ','sxy_68 ','szx_68 ',
1075 . 'sx_69 ','sxy_69 ','szx_69 ','sx_70 ','sxy_70 ',
1076 . 'szx_70 ','sx_71 ','sxy_71 ','szx_71 ','sx_72 ',
1077 . 'sxy_72 ','szx_72 ','sx_73 ','SXY_73 ','SZX_73 ',
1078 . 'SX_74 ','SXY_74 ','SZX_74 ','SX_75 ','SXY_75 ',
1079 . 'SZX_75 ','SX_76 ','SXY_76 ','SZX_76 ','SX_77 ',
1080 . 'SXY_77 ','SZX_77 ','SX_78 ','SXY_78 ','SZX_78 ',
1081 . 'SX_79 ','SXY_79 ','SZX_79 ','SX_80 ','SXY_80 ',
1082 . 'SZX_80 ','SX_81 ','SXY_81 ','SZX_81 ','EPSP ',
1083 . 'EPSP_1 ','EPSP_2 ','EPSP_3 ','EPSP_4 ','EPSP_5 ',
1084 . 'EPSP_6 ','EPSP_7 ','EPSP_8 ','EPSP_9 ','EPSP_10 ',
1085 . 'EPSP_11 ','EPSP_12 ','EPSP_13 ','EPSP_14 ','EPSP_15 ',
1086 . 'EPSP_16 ','EPSP_17 ','EPSP_18 ','EPSP_19 ','EPSP_20 ',
1087 . 'EPSP_21 ','EPSP_22 ','EPSP_23 ','EPSP_24 ','EPSP_25 ',
1088 . 'EPSP_26 ','EPSP_27 ','EPSP_28 ','EPSP_29 ','EPSP_30 ',
1089 . 'epsp_31 ','epsp_32 ','epsp_33 ','epsp_34 ','epsp_35 ',
1090 . 'epsp_36 ','epsp_37 ','epsp_38 ','epsp_39 ','epsp_40 ',
1091 . 'epsp_41 ','epsp_42 ','epsp_43 ','epsp_44 ','epsp_45 ',
1092 . 'epsp_46 ','epsp_47 ','epsp_48 ','epsp_49 ','epsp_50 ',
1093 . 'EPSP_51 ','EPSP_52 ','EPSP_53 ','EPSP_54 ','EPSP_55 ',
1094 . 'EPSP_56 ','EPSP_57 ','EPSP_58 ','EPSP_59 ','EPSP_60 ',
1095 . 'EPSP_61 ','EPSP_62 ','EPSP_63 ','EPSP_64 ','EPSP_65 ',
1096 . 'EPSP_66 ','EPSP_67 ','EPSP_68 ','EPSP_69 ','EPSP_70 ',
1097 . 'EPSP_71 ','EPSP_72 ','EPSP_73 ','EPSP_74 ','EPSP_75 ',
1098 . 'EPSP_76 ','EPSP_77 ','epsp_78 ','epsp_79 ','epsp_80 ',
1099 . 'epsp_81 ','epsd '/
1100C springs
1101 DATA VARR/
1102 . 'off ','fx ','fy ','fz ','mx ',
1103 . 'my ','mz ','lx ','ly ','lz ',
1104 . 'rx ','ry ','RZ ','IE ','F1 ',
1105 . 'F2 ','F1G_X ','F1G_Y ','F1G_Z ','F2G_X ',
1106 . 'F2G_Y ','F2G_Z ','F1L_X ','F1L_Y ','F1L_Z ',
1107 . 'F2L_X ','F2L_Y ','F2L_Z ','M1G_X ','M1G_Y ',
1108 . 'M1G_Z ','M2G_X ','M2G_Y ','M2G_Z ','M1L_X ',
1109 . 'M1L_Y ','M1L_Z ','M2L_X ','M2L_Y ','M2L_Z ',
1110 . 'D1G_X ','D1G_Y ','D1G_Z ','D2G_X ','D2G_Y ',
1111 . 'D2G_Z ','D1L_X ','D1L_Y ','D1L_Z ','D2L_X ',
1112 . 'D2L_Y ','D2L_Z ','R1G_X ','R1G_Y ','R1G_Z ',
1113 . 'R2G_X ','R2G_Y ','R2G_Z ','R1L_X ','R1L_Y ',
1114 . 'R1L_Z ','R2L_X ','R2L_Y ','R2L_Z ','LENGTH ',
1115 . 'FAIL '/
1116
1117 DATA varns/
1118 . 'OFF ','FX ','LX ','IE '/
1119
1120 DATA varsph/
1121 . 'OFF ','SX ','SY ','SZ ','SXY ',
1122 . 'SYZ ','SXZ ','IE ','DENS ','WFVIS ',
1123 . 'VOL ','PLAS ','TEMP ','PLSR ','DAM1 ',
1124 . 'DAM2 ','DAM3 ','DAM4 ','DAMA ','SA1 ',
1125 . 'SA2 ','SA3 ','CR ','CAP ','K0 ',
1126 . 'RK ','TD ','EFIB ','ISTA ','VPLA ',
1127 . 'BFRAC ','WPLA ','SFIB ','DAM5 ','LSX ',
1128 . 'LSY ','LSZ ','LSXY ','LSYZ ','LSXZ ',
1129 . 'DIAMETER'/
1130 DATA varur/
1131 . 'OFF ','EINT ','EINR ','XN1 ','XN2 ',
1132 . 'XN3 ','DX1 ','DX2 ','DX3 ','XM1 ',
1133 . 'XM2 ','XM3 '/
1134C interfaces
1135 DATA varin/
1136 . 'FNX ','FNY ','FNZ ','FTX ','FTY ',
1137 . 'FTZ ','SFW ','|FNX| ','|FNY| ','|FNZ| ',
1138 . '||FN|| ','|FX| ','|FY| ','|FZ| ','||F|| ',
1139 . 'PVOL ','PSURF ','PMED ','DELTAP ','VOL ',
1140 . 'SURF ','MX ','MY ','MZ ','QFRIC ',
1141 . 'CE_ELAST','CE_FRIC ','CE_DAMP ','CAREA '/
1142C rigid wall
1143 DATA varrw/
1144 . 'FNX ','FNY ','FNZ ','FTX ','FTY ',
1145 . 'FTZ '/
1146C rigid body
1147 DATA varrb/
1148 . 'FX ','FY ','FZ ','MX ','MY ',
1149 . 'MZ ','RX ','RY ','RZ ',
1150 . 'FXI ','FYI ','FZI ',
1151 . 'MXI ','MYI ','MZI '/
1152C flexible body
1153 DATA varfx/
1154 . 'IE ','KE ','EFW ','DE '/
1155C flexible body mode
1156 DATA varfxm/
1157 . 'D ','V ','A '/
1158C accel
1159 DATA varac/
1160 . 'AX ','AY ','AZ ','WX ','WY ',
1161 . 'WZ '/
1162C section
1163 DATA varse/
1164 . 'FNX ','FNY ','FNZ ','FTX ','FTY ',
1165 . 'FTZ ','M1 ','M2 ','M3 ','WORK ',
1166 . 'DFX ','DFY ','DFZ ','DF2 ','WORKR ',
1167 . 'DMX ','DMY ','DMZ ','DM2 ','KIN ',
1168 . 'KINR ','DMVX ','DMVY ','DMVZ ','DKIN ',
1169 . 'DMVRX ','DMVRY ','DMVRZ ','DKINR ','TFEXT ',
1170 . 'MX ','MY ','MZ ','F1 ','F2 ',
1171 . 'F3 ','CX ','CY ','CZ '/
1172C joint cyl
1173 DATA varjo/
1174 . 'FX ','FY ','FZ ','MX ','MY ',
1175 . 'MZ '/
1176C airbag
1177 DATA varab/
1178 . 'MASS ','VOL ','P ','RO ','IE ',
1179 . 'OE ','OU '/
1180C monitored volume
1181 DATA varmv4/
1182 . 'MASS ','VOL ','P ','A ','T ',
1183 . 'AO ','UO ','AC ','UC '/
1184 DATA varmv/
1185 . 'MASS ','VOL ','P ','A ','T ',
1186 . 'AO ','UO ','AC ','UC ','CP ',
1187 . 'CV ','GAMA ','DTBAG ','NFV ','MASS-IN ',
1188 . 'ENTHA-IN','ENER-INT','WORK ','UPCRIT '/
1189C
1190C Caution : all variables added or defined in VARMV/ must be added also in VARMVENT/
1191C
1192 DATA varmvent/
1193 . 'MASS ','VOL ','P ','A ','T ',
1194 . 'AO ','UO ','AC ','UC ','CP ',
1195 . 'CV ','GAMA ','DTBAG ','NFV ','MASS-IN ',
1196 . 'ENTHA-IN','ENER-INT','WORK ','UPCRIT ','UNUSED ',
1197 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1198 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1199 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1200 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1201 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1202 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1203 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1204 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1205 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1206 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1207 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1208 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1209 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1210 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1211 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1212 . 'UNUSED ','UNUSED ','unused ','unused ','unused ',
1213 . 'aout1 ','bout1 ','uout1 ','mout1 ','hout1 ',
1214 . 'aout2 ','bout2 ','uout2 ','mout2 ','hout2 ',
1215 . 'aout3 ','bout3 ','uout3 ','mout3 ','hout3 ',
1216 . 'aout4 ','bout4 ','uout4 ','mout4 ','hout4 ',
1217 . 'aout5 ','bout5 ','uout5 ','mout5 ','hout5 ',
1218 . 'aout6 ','bout6 ','uout6 ','mout6 ','hout6 ',
1219 . 'aout7 ','bout7 ','uout7 ','MOUT7 ','HOUT7 ',
1220 . 'AOUT8 ','BOUT8 ','UOUT8 ','MOUT8 ','HOUT8 ',
1221 . 'AOUT9 ','BOUT9 ','UOUT9 ','MOUT9 ','HOUT9 ',
1222 . 'AOUT10 ','BOUT10 ','UOUT10 ','MOUT10 ','HOUT10 '/
1223C part
1224 DATA varpa/
1225 . 'IE ','KE ','XMOM ','YMOM ','ZMOM ',
1226 . 'MASS ','HE ','TURBKE ','XCG ','YCG ',
1227 . 'ZCG ','XXMOM ','yymom ','zzmom ','ixx ',
1228 . 'iyy ','izz ','ixy ','iyz ','izx ',
1229 . 'rie ','kerb ','rkerb ','rke ','eroded ',
1230 . 'unused ','unused ','heat ','vx ','vy ',
1231 . 'vz ','pw '/
1232C rivets
1233 DATA VARRIV/
1234 . 'off ','f1x ','f1y ','f1z ','f2x ',
1235 . 'f2y ','f2z ','fx ','fy ','fz '/
1236C frames
1237 DATA VARF1/
1238 . 'ox ','oy ','oz ','r11 ','r12 ',
1239 . 'r13 ','r21 ','r22 ','r23 ','r31 ',
1240 . 'r32 ','r33 ','vx ','vy ','vz ',
1241 . 'vrx ','vry ','VRZ '/
1242 DATA varfr/
1243 . 'OX ','OY ','OZ ','R11 ','R12 ',
1244 . 'R13 ','R21 ','R22 ','R23 ','R31 ',
1245 . 'R32 ','R33 ','VX ','VY ','VZ ',
1246 . 'VRX ','VRY ','VRZ ','AX ','AY ',
1247 . 'AZ ','ARX ','ARY ','ARZ '/
1248 DATA vargau/
1249 . 'P ','IE ','DENS ','T ','PA ',
1250 . 'PAX ','PAY ','PAZ '/
1251C cluster element
1252 DATA varclus/
1253 . 'FX ','FY ','FZ ','MX ','MY ',
1254 . 'MZ ','FS ','FN ','MS ','MN ',
1255 . 'FAIL '/
1256 DATA varflow/
1257 . 'FLOW '/
1258 DATA varsurf/
1259 . 'AREA ','MASSFLOW','VELOCITY','P ','A ','MASS '/
1260C
1261 DATA varslip/
1262 . 'RINGSLIP','FN ','F1 ','F2 ','THETA ','GAMMA '/
1263C
1264 DATA varret/
1265 . 'SLIP','FN ','LOCK '/
1266c
1267 DATA varsens/
1268 . 'STATUS '/
1269c
1270 DATA varchecksum/
1271 . 'CHECKSUM '/
1272C-----------------------------------------------
1273C TH GENERIC VARIABLES
1274C-----------------------------------------------
1275 DATA varng/'DEF ','D ','V ',
1276 . 'A ','VR ','AR ','XYZ ',
1277 . 'PINCH '/
1278 DATA ivarng
1279 . / 1, 2, 3, 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1280 . 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1281 . 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1282 . 7, 8, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1283 . 10,11,12, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1284 . 13,14,15, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1285 . 16,17,18, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1286 . 630,631,632,633,634,635,636,637,638, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1287 DATA varsg1/
1288 .'DEF ','STRESS ','LOCSTRS ','SXGP ','SYGP ','SZGP ',
1289 .'SXYGP ','SYZGP ','SXZGP ','LSXGP ','LSYGP ','LSZGP ',
1290 .'LSXYGP ','LSYZGP ','LSXZGP ','USERA ','USERB ','USERC ',
1291 .'USERD ',
1292 .'SIG111 ','SIG211 ','SIG311 ','SIG121 ','SIG221 ','SIG321 ',
1293 .'SIG131 ','SIG231 ','SIG331 ','SIG141 ','SIG241 ','SIG341 ',
1294 .'SIG151 ','SIG251 ','SIG351 ','SIG161 ','SIG261 ','SIG361 ',
1295 .'SIG171 ','SIG271 ','SIG371 ','SIG181 ','SIG281 ','SIG381 ',
1296 .'SIG191 ','SIG291 ','SIG391 ',
1297 .'SIG112 ','SIG212 ','SIG312 ','SIG122 ','SIG222 ','SIG322 ',
1298 .'SIG132 ','SIG232 ','SIG332 ','SIG142 ','SIG242 ','SIG342 ',
1299 .'SIG152 ','SIG252 ','SIG352 ','SIG162 ','SIG262 ','SIG362 ',
1300 .'SIG172 ','SIG272 ','SIG372 ','SIG182 ','SIG282 ','SIG382 ',
1301 .'SIG192 ','SIG292 ','SIG392 ',
1302 .'SIG113 ','SIG213 ','SIG313 ','SIG123 ','SIG223 ','SIG323 ',
1303 .'SIG133 ','SIG233 ','SIG333 ','SIG143 ','SIG243 ','SIG343 ',
1304 .'SIG153 ','SIG253 ','SIG353 ','SIG163 ','SIG263 ','SIG363 ',
1305 .'SIG173 ','SIG273 ','SIG373 ','SIG183 ','SIG283 ','SIG383 ',
1306 .'SIG193 ','SIG293 ','SIG393 ',
1307 .'SIG1U1 ','SIG2U1 ','SIG3U1 ','SIG1U2 ','SIG2U2 ','SIG3U2 ',
1308 .'SIG1U3 ','SIG2U3 ','SIG3U3 ',
1309 .'SIG1D1 ','SIG2D1 ','SIG3D1 ','SIG1D2 ','SIG2D2 ','SIG3D2 ',
1310 .'SIG1D3 ','SIG2D3 ','SIG3D3 ','USER111 ','USER211 ','USER311 ',
1311 .'USER121 ','USER221 ','USER321 ','USER131 ','USER231 ','USER331 ',
1312 .'USER141 ','USER241 ','USER341 ','USER151 ','USER251 ','USER351 ',
1313 .'USER161 ','USER261 ','USER361 ','USER171 ','USER271 ','USER371 ',
1314 .'USER181 ','USER281 ','USER381 ','USER191 ','USER291 ','USER391 ',
1315 .'USER112 ','USER212 ','USER312 ','USER122 ','USER222 ','USER322 ',
1316 .'USER132 ','USER232 ','USER332 ','USER142 ','USER242 ','USER342 ',
1317 .'USER152 ','USER252 ','USER352 ','USER162 ','USER262 ','USER362 ',
1318 .'USER172 ','USER272 ','USER372 ','USER182 ','USER282 ','USER382 ',
1319 .'USER192 ','USER292 ','USER392 ','USER113 ','USER213 ','USER313 ',
1320 .'USER123 ','USER223 ','USER323 ','USER133 ','USER233 ','USER333 ',
1321 .'USER143 ','USER243 ','USER343 ','USER153 ','USER253 ','USER353 ',
1322 .'USER163 ','USER263 ','USER363 ','USER173 ','USER273 ','USER373 ',
1323 .'USER183 ','USER283 ','USER383 ','USER193 ','USER293 ','USER393 ',
1324 .'STRAIN ','DEF_MAT1','DEF_MAT2'/
1325C
1326c DATA VARSG2
1327C Definition of VARSG2 : strain for /PROP/TYPE22
1328 DO j=1,200
1329 IF (j <= 9) THEN
1330 WRITE(chj,'(I1.1)')j
1331 ELSEIF (j <= 99) THEN
1332 WRITE(chj,'(I2.2)')j
1333 ELSE
1334 WRITE(chj,'(I3.3)')j
1335 ENDIF
1336 DO i=1,9
1337 WRITE(chi,'(I1.1)')i
1338 DO k=1,9
1339 WRITE(chk,'(I1.1)')k
1340 var_tmp(1:5)='EPS_'//chi
1341 var_tmp=var_tmp(1:5)//chk//'_'//chj
1342 varsg2((j-1)*99+(i-1)*9+k)=var_tmp
1343 ENDDO
1344 ENDDO
1345 ENDDO
1346c DATA VARSG3
1347C Definition of VARSG3 : stress for /PROP/TYPE22
1348 DO j=1,200
1349 IF (j <= 9) THEN
1350 WRITE(chj,'(I1.1)')j
1351 ELSEIF (j <= 99) THEN
1352 WRITE(chj,'(I2.2)')j
1353 ELSE
1354 WRITE(chj,'(I3.3)')j
1355 ENDIF
1356 DO i=1,9
1357 WRITE(chi,'(I1.1)')i
1358 DO k=1,9
1359 WRITE(chk,'(I1.1)')k
1360 var_tmp(1:5)='SIG_'//chi
1361 var_tmp=var_tmp(1:5)//chk//'_'//chj
1362 varsg3((j-1)*9*9+(i-1)*9+k)=var_tmp
1363 ENDDO
1364 ENDDO
1365 ENDDO
1366c DATA VARSG4
1367C Definition of VARSG4
1368C
1369 DATA varsg4/
1370 .'LOCSTRN '/
1371
1372c VARCG = VARCG1 + VARCG2+ VARCG3
1373 DATA varcg1/
1374 . 'DEF ','STRESS ','STRAIN ','PLAS ','USERA ',
1375 . 'USERB ','USERC ','USERD ','USER1 ','USER2 ',
1376 . 'USER3 ','USER4 ','USER5 ','USER6 ','USER7 ',
1377 . 'USER8 ','USER9 ','USER10 ','USER11 ','USER12 ',
1378 . 'USER13 ','USER14 ','USER15 ','USER16 ','USER17 ',
1379 . 'USER18 ','USER19 ','USER20 ',
1380 . 'SIG_1 ','SIG_2 ','SIG_3 ','SIG_4 ','SIG_5 ',
1381 . 'SIG_6 ','SIG_7 ','SIG_8 ','SIG_9 ','SIG_10 ',
1382 . 'SIG_11 ','SIG_12 ','SIG_13 ','SIG_14 ','SIG_15 ',
1383 . 'SIG_16 ','SIG_17 ','SIG_18 ','SIG_19 ','SIG_10 ',
1384 . 'SIG_21 ','SIG_22 ','SIG_23 ','SIG_24 ','SIG_25 ',
1385 . 'SIG_26 ','SIG_27 ','SIG_28 ','SIG_29 ','SIG_30 ',
1386 . 'SIG_31 ','SIG_32 ','SIG_33 ','SIG_34 ','SIG_35 ',
1387 . 'SIG_36 ','SIG_37 ','SIG_38 ','SIG_39 ','SIG_40 ',
1388 . 'SIG_41 ','SIG_42 ','SIG_43 ','SIG_44 ','SIG_45 ',
1389 . 'SIG_46 ','SIG_47 ','SIG_48 ','SIG_49 ','SIG_50 ',
1390 . 'SIG_51 ','SIG_52 ','SIG_53 ','SIG_54 ','SIG_55 ',
1391 . 'SIG_56 ','SIG_57 ','SIG_58 ','SIG_59 ','SIG_60 ',
1392 . 'SIG_61 ','SIG_62 ','SIG_63 ','SIG_64 ','SIG_65 ',
1393 . 'SIG_66 ','SIG_67 ','SIG_68 ','SIG_69 ','SIG_70 ',
1394 . 'SIG_71 ','SIG_72 ','SIG_73 ','SIG_74 ','SIG_75 ',
1395 . 'SIG_76 ','SIG_77 ','SIG_78 ','SIG_79 ','SIG_80 ',
1396 . 'SIG_81 ','SIG_82 ','SIG_83 ','SIG_84 ','SIG_85 ',
1397 . 'SIG_86 ','SIG_87 ','SIG_88 ','SIG_89 ','SIG_90 ',
1398 . 'SIG_91 ','SIG_92 ','SIG_93 ','SIG_94 ','SIG_95 ',
1399 . 'SIG_96 ','SIG_97 ','SIG_98 ','SIG_99 ',
1400 . 'USER21 ','USER22 ',
1401 . 'USER23 ','USER24 ','USER25 ','USER26 ','USER27 ',
1402 . 'USER28 ','USER29 ','USER30 ','USER31 ','USER32 ',
1403 . 'USER33 ','USER34 ','USER35 ','USER36 ','USER37 ',
1404 . 'USER38 ','USER39 ','USER40 ','USER41 ','USER42 ',
1405 . 'USER43 ','USER44 ','USER45 ','USER46 ','USER47 ',
1406 . 'USER48 ','USER49 ','USER50 ','USER51 ','USER52 ',
1407 . 'USER53 ','USER54 ','USER55 ','USER56 ','USER57 ',
1408 . 'USER58 ','USER59 ','USER60 '/
1409 DATA varcg2/
1410 . 'USLAY011','USLAY012','USLAY013','USLAY014','USLAY021',
1411 . 'USLAY022','USLAY023','USLAY024','USLAY031','USLAY032',
1412 . 'USLAY033','USLAY034','USLAY041','USLAY042','USLAY043',
1413 . 'USLAY044','USLAY051','USLAY012','USLAY013','USLAY014',
1414 . 'USLAY061','USLAY062','USLAY063','USLAY064','USLAY071',
1415 . 'USLAY072','USLAY073','USLAY074','USLAY081','USLAY082',
1416 . 'USLAY083','USLAY084','USLAY091','USLAY092','USLAY093',
1417 . 'USLAY094','USLAY101','USLAY102','USLAY103','USLAY104',
1418 . 'USLAY111','USLAY112','USLAY113','USLAY114','USLAY121',
1419 . 'USLAY122','USLAY123','USLAY124','USLAY131','USLAY132',
1420 . 'USLAY133','USLAY134','USLAY141','USLAY142','USLAY143',
1421 . 'USLAY144','USLAY151','USLAY152','USLAY153','USLAY154',
1422 . 'USLAY161','USLAY162','USLAY163','USLAY164','USLAY171',
1423 . 'USLAY172','USLAY173','USLAY174','USLAY181','USLAY182',
1424 . 'USLAY183','USLAY184','USLAY191','USLAY192','USLAY193',
1425 . 'USLAY194','USLAY201','USLAY202','USLAY203','USLAY204',
1426 . 'USLAY211','USLAY212','USLAY213','USLAY214','USLAY221',
1427 . 'USLAY222','USLAY223','USLAY224','USLAY231','USLAY232',
1428 . 'USLAY233','USLAY234','USLAY241','USLAY242','USLAY243',
1429 . 'USLAY244','USLAY251','USLAY252','USLAY253','USLAY254',
1430 . 'USLAY261','USLAY262','USLAY263','USLAY264','USLAY271',
1431 . 'USLAY272','USLAY273','USLAY274','USLAY281','USLAY282',
1432 . 'USLAY283','uslay284','uslay291','uslay292','uslay293',
1433 . 'uslay294','uslay301','uslay302','uslay303','uslay304',
1434 . 'uslay311','uslay312','uslay313','uslay314','uslay321',
1435 . 'uslay322','uslay323','uslay324','uslay331','uslay332',
1436 . 'uslay333','uslay334','uslay341','uslay342','uslay343',
1437 . 'uslay344','uslay351','uslay352','uslay353','uslay354',
1438 . 'uslay361','uslay362','uslay363','uslay364','uslay371',
1439 . 'uslay372','uslay373','uslay374','uslay381','uslay382',
1440 . 'uslay383','uslay384','uslay391','uslay392','uslay393',
1441 . 'uslay394','uslay401','uslay402','uslay403','uslay404',
1442 . 'uslay411','uslay412','uslay413','uslay414','uslay421',
1443 . 'uslay422','uslay423','uslay424','uslay431','uslay432',
1444 . 'uslay433','uslay434','uslay441','uslay442','uslay443',
1445 . 'uslay444','uslay451','uslay452','uslay453','uslay454',
1446 . 'uslay461','uslay462','uslay463','uslay464','uslay471',
1447 . 'uslay472','uslay473','uslay474','uslay481','uslay482',
1448 . 'uslay483','uslay484','uslay491','uslay492','uslay493',
1449 . 'uslay494','uslay501','uslay502','uslay503','uslay504',
1450 . 'uslay511','uslay512','uslay513','uslay514','uslay521',
1451 . 'uslay522','uslay523','uslay524','uslay531','uslay532',
1452 . 'uslay533','uslay534','uslay541','uslay542','uslay543',
1453 . 'uslay544','uslay551','uslay552','uslay553','uslay554',
1454 . 'uslay561','uslay562','uslay563','uslay564','uslay571',
1455 . 'uslay572','uslay573','uslay574','uslay581','uslay582',
1456 . 'uslay583','uslay584','uslay591','uslay592','uslay593',
1457 . 'uslay594','uslay601','uslay602','uslay603','uslay604',
1458 . 'uslay611','uslay612','uslay613','uslay614','uslay621',
1459 . 'uslay622','uslay623','uslay624','uslay631','uslay632',
1460 . 'uslay633','uslay634','uslay641','uslay642','uslay643',
1461 . 'uslay644','uslay651','uslay652','uslay653','uslay654',
1462 . 'uslay661','uslay662','uslay663','uslay664','uslay671',
1463 . 'uslay672','uslay673','uslay674','uslay681','uslay682',
1464 . 'uslay683','uslay684','uslay691','uslay692','uslay693',
1465 . 'uslay694','uslay711','uslay712','uslay713','uslay714',
1466 . 'uslay701','uslay702','uslay703','uslay704','uslay721',
1467 . 'uslay722','uslay723','uslay724','uslay731','uslay732',
1468 . 'uslay733','uslay734','uslay741','uslay742','uslay743',
1469 . 'uslay744','uslay751','uslay752','uslay753','uslay754',
1470 . 'uslay761','uslay762','uslay763','uslay764','uslay771',
1471 . 'uslay772','uslay773','uslay774','uslay781','uslay782',
1472 . 'uslay783','uslay784','uslay791','uslay792','uslay793',
1473 . 'uslay794','uslay801','uslay802','uslay803','uslay804',
1474 . 'uslay811','uslay812','uslay813','uslay814','uslay821',
1475 . 'uslay822','uslay823','uslay824','uslay831','uslay832',
1476 . 'uslay833','uslay834','uslay841','uslay842','uslay843',
1477 . 'uslay844','uslay851','uslay852','uslay853','uslay854',
1478 . 'uslay861','uslay862','uslay863','uslay864','uslay871',
1479 . 'uslay872','uslay873','uslay874','uslay881','uslay882',
1480 . 'uslay883','uslay884','uslay891','uslay892','uslay893',
1481 . 'uslay894','uslay901','uslay902','uslay903','uslay904',
1482 . 'uslay911','uslay912','uslay913','uslay914','uslay921',
1483 . 'uslay922','uslay923','uslay924','uslay931','uslay932',
1484 . 'uslay933','uslay934','uslay941','uslay942','uslay943',
1485 . 'uslay944','uslay951','uslay952','uslay953','uslay954',
1486 . 'uslay961','uslay962','uslay963','uslay964','uslay971',
1487 . 'uslay972','uslay973','uslay974','uslay981','uslay982',
1488 . 'uslay983','uslay984','uslay991','uslay992','uslay993',
1489 . 'uslay994','failure ','wpl01_10','wpl11_20','wpl21_30',
1490 . 'wpl31_40','wpl41_50','wpl51_60','wpl61_70','wpl71_80',
1491 . 'wpl81_90','wpl91_99'/
1492 DATA VARCG3/
1493 . 'sigv_1 ','sigv_2 ','sigv_3 ','sigv_4 ','sigv_5 ',
1494 . 'sigv_6 ','sigv_7 ','sigv_8 ','sigv_9 ','sigv_10 ',
1495 . 'sigv_11 ','sigv_12 ','sigv_13 ','sigv_14 ','sigv_15 ',
1496 . 'sigv_16 ','sigv_17 ','sigv_18 ','sigv_19 ','sigv_10 ',
1497 . 'sigv_21 ','sigv_22 ','sigv_23 ','sigv_24 ','sigv_25 ',
1498 . 'sigv_26 ','sigv_27 ','sigv_28 ','sigv_29 ','sigv_30 ',
1499 . 'sigv_31 ','sigv_32 ','sigv_33 ','sigv_34 ','sigv_35 ',
1500 . 'sigv_36 ','sigv_37 ','sigv_38 ','sigv_39 ','sigv_40 ',
1501 . 'sigv_41 ','sigv_42 ','sigv_43 ','sigv_44 ','sigv_45 ',
1502 . 'sigv_46 ','sigv_47 ','sigv_48 ','sigv_49 ','sigv_50 ',
1503 . 'sigv_51 ','sigv_52 ','sigv_53 ','sigv_54 ','sigv_55 ',
1504 . 'sigv_56 ','sigv_57 ','sigv_58 ','sigv_59 ','sigv_60 ',
1505 . 'sigv_61 ','sigv_62 ','sigv_63 ','sigv_64 ','sigv_65 ',
1506 . 'sigv_66 ','sigv_67 ','sigv_68 ','sigv_69 ','sigv_70 ',
1507 . 'sigv_71 ','sigv_72 ','sigv_73 ','sigv_74 ','sigv_75 ',
1508 . 'sigv_76 ','sigv_77 ','sigv_78 ','sigv_79 ','sigv_80 ',
1509 . 'sigv_81 ','sigv_82 ','sigv_83 ','sigv_84 ','sigv_85 ',
1510 . 'sigv_86 ','sigv_87 ','sigv_88 ','sigv_89 ','sigv_90 ',
1511 . 'sigv_91 ','sigv_92 ','sigv_93 ','sigv_94 ','sigv_95 ',
1512 . 'sigv_96 ','sigv_97 ','sigv_98 ','sigv_99 ','sigv_100'/
1513 DATA VARCGPINCH/
1514 . 'pinch '/
1515 DATA VARTG/'def '/
1516 DATA IVARTG
1517 . / 1, 2, 3, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1518 DATA VARPG/'def '/
1519 DATA IVARPG
1520 . / 1, 2, 3, 4, 5, 6, 7, 8, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1521 DATA VARRG/'def '/
1522 DATA IVARRG
1523 . / 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,65, 0, 0, 0/
1524 DATA VARURG/'def '/
1525 DATA IVARURG
1526 . / 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12, 0, 0, 0, 0, 0, 0/
1527 DATA VARNSG/'def '/
1528 DATA IVARNSG
1529 . / 1, 2, 3, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1530 DATA IVNS2R
1531 . / 1, 2, 8,14, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1532 DATA VARSPG/
1533 . 'def ','stress '/
1534 DATA IVARSPG
1535 . / 1, 2, 3, 4, 5, 6, 7, 8, 9,12,13,41, 0, 0, 0, 0, 0, 0,
1536 . 2, 3, 4, 5, 6, 7, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1537C
1538 DATA VARING/'def ','fn ','ft ','sfw ',
1539 . '|fn| ','|f| '/
1540 DATA IVARING
1541 . / 1, 2, 3, 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1542 . 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1543 . 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1544 . 7, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1545 . 8, 9,10,11, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1546 . 12,13,14,15, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1547 DATA VARRWG/'def ','fn ','ft '/
1548 DATA IVARRWG
1549 . / 1, 2, 3, 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1550 . 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1551 . 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1552 DATA VARRBG/'def ','f ','m ','r ',
1553 . 'fi ','mi '/
1554 DATA IVARRBG
1555 . / 1, 2, 3, 4, 5, 6, 7, 8, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1556 . 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1557 . 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1558 . 7, 8, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1559 . 10,11,12, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1560 . 13,14,15, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1561 DATA VARFXG/'def '/
1562 DATA IVARFXG
1563 . / 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1564 DATA VARFXMG/'def '/
1565 DATA IVARFXMG
1566 . / 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1567 DATA VARACG/'def ','w '/
1568 DATA IVARACG
1569 . / 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1570 . 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1571 DATA VARSEG/'def ','fn ','ft ','m ','center ',
1572 . 'global ','local '/
1573 DATA IVARSEG
1574 . / 1, 2, 3, 4, 5, 6, 7, 8, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1575 . 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1576 . 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1577 . 7, 8, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1578 . 37,38,39, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1579 . 1, 2, 3, 4, 5, 6,31,32,33, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1580 . 7, 8, 9,34,35,36, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1581 DATA VARJOG/'def ','f ','m '/
1582 DATA IVARJOG
1583 . / 1, 2, 3, 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1584 . 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1585 . 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1586 DATA VARABG/'def '/
1587 DATA IVARABG
1588 . / 1, 2, 3, 4, 5, 6, 7,13,14, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1589 DATA VARMVG4/'def '/,VARMVG/'def ','gas ',
1590 . 'out1 ','out2 ','out3 ','out4 ','out5 ',
1591 . 'out6 ','out7 ','out8 ','out9 ','out10 ',
1592 . 'ao1 ','bo1 ','uo1 ','mo1 ','ho1 ',
1593 . 'ao2 ','bo2 ','uo2 ','mo2 ','ho2 ',
1594 . 'ao3 ','bo3 ','uo3 ','mo3 ','ho3 ',
1595 . 'ao4 ','bo4 ','uo4 ','mo4 ','ho4 ',
1596 . 'ao5 ','bo5 ','uo5 ','mo5 ','ho5 ',
1597 . 'ao6 ','bo6 ','uo6 ','mo6 ','ho6 ',
1598 . 'ao7 ','bo7 ','uo7 ','mo7 ','ho7 ',
1599 . 'ao8 ','bo8 ','uo8 ','mo8 ','ho8 ',
1600 . 'ao9 ','bo9 ','uo9 ','mo9 ','ho9 ',
1601 . 'ao10 ','bo10 ','uo10 ','mo10 ','ho10 '/
1602 DATA VARMVENTG/'def ','gas ',
1603 . 'out1 ','out2 ','out3 ','out4 ','out5 ',
1604 . 'out6 ','out7 ','out8 ','out9 ','out10 ',
1605 . 'ao1 ','bo1 ','uo1 ','mo1 ','ho1 ',
1606 . 'ao2 ','bo2 ','uo2 ','mo2 ','ho2 ',
1607 . 'ao3 ','bo3 ','uo3 ','mo3 ','ho3 ',
1608 . 'ao4 ','bo4 ','uo4 ','mo4 ','ho4 ',
1609 . 'ao5 ','bo5 ','uo5 ','mo5 ','ho5 ',
1610 . 'ao6 ','bo6 ','uo6 ','mo6 ','ho6 ',
1611 . 'ao7 ','bo7 ','uo7 ','mo7 ','ho7 ',
1612 . 'ao8 ','bo8 ','uo8 ','mo8 ','ho8 ',
1613 . 'ao9 ','bo9 ','uo9 ','mo9 ','ho9 ',
1614 . 'ao10 ','bo10 ','uo10 ','mo10 ','ho10 '/
1615 DATA IVARMG4
1616 . / 1, 2, 3, 4, 5, 6, 7,13,14, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1617 DATA IVARMVG
1618 . / 1, 2, 3, 4, 5, 6, 7,13,14, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1619 . 10,11,12, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1620 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1621 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1622 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1623 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1624 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1625 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1626 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1627 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1628 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1629 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1630 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1631 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1632 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1633 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1634 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1635 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1636 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1637 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1638 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1639 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1640 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1641 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1642 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1643 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1644 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1645 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1646 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1647 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1648 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1649 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1650 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1651 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1652 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1653 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1654 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1655 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1656 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1657 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1658 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1659 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1660 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1661 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1662 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1663 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1664 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1665 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1666 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1667 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1668 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1669 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1670 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1671 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1672 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1673 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1674 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1675 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1676 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1677 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1678 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1679 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1680
1681
1682 DATA IVARMVENTG
1683 . / 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1684 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1685 . 101,102,103,104,105, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1686 . 106,107,108,109,110, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1687 . 111,112,113,114,115, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1688 . 116,117,118,119,120, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1689 . 121,122,123,124,125, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1690 . 126,127,128,129,130, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1691 . 131,132,133,134,135, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1692 . 136,137,138,139,140, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1693 . 141,142,143,144,145, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1694 . 146,147,148,149,150, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1695 . 101, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1696 . 102, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1697 . 103, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1698 . 104, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1699 . 105, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1700 . 106, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1701 . 107, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1702 . 108, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1703 . 109, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1704 . 110, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1705 . 111, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1706 . 112, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1707 . 113, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1708 . 114, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1709 . 115, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1710 . 116, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1711 . 117, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1712 . 118, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1713 . 119, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1714 . 120, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1715 . 121, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1716 . 122, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1717 . 123, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1718 . 124, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1719 . 125, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1720 . 126, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1721 . 127, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1722 . 128, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1723 . 129, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1724 . 130, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1725 . 131, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1726 . 132, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1727 . 133, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1728 . 134, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1729 . 135, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1730 . 136, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1731 . 137, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1732 . 138, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1733 . 139, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1734 . 140, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1735 . 141, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1736 . 142, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1737 . 143, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1738 . 144, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1739 . 145, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1740 . 146, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1741 . 147, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1742 . 148, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1743 . 149, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1744 . 150, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1745 DATA VARPAG/'def '/
1746 DATA IVARPAG
1747 . / 1, 2, 3, 4, 5, 6, 7,28, 25, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1748 DATA VARRIVG/'def '/
1749 DATA IVARRIV
1750 . / 1, 2, 3, 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1751C frames
1752 DATA VARF1G/'def ','o+r ','v+vr '/
1753 DATA IVARF1G
1754 . / 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,15,16,17,18,
1755 . 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12, 0, 0, 0, 0, 0, 0,
1756 . 13,14,15,16,17,18, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1757 DATA VARFRG/'def ','o+r ','v+vr ','a+ar '/
1758 DATA IVARFRG
1759 . / 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,15,16,17,18,
1760 . 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12, 0, 0, 0, 0, 0, 0,
1761 . 13,14,15,16,17,18, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1762 . 19,20,21,22,23,24, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1763 DATA VARGAUG/'def '/
1764 DATA IVARGAUG
1765 . / 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1766 DATA VARCLUSG/'def ','floc '/
1767 DATA IVARCLUSG
1768 . / 1, 2, 3, 4, 5, 6,11, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1769 . 7, 8, 9,10, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1770 DATA VARFLOWG/'def '/
1771 DATA IVARFLOWG
1772 . / 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1773 DATA VARSURFG/'def '/
1774 DATA IVARSURFG
1775 . / 1, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1776 DATA VARSLIPG/'def '/
1777 DATA IVARSLIPG
1778 . / 1, 2, 3, 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1779 DATA VARRETG/'def '/
1780 DATA IVARRETG
1781 . / 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1782 DATA VARSENSG/
1783 . 'def '/
1784 DATA IVARSENSG
1785 . / 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1786 DATA VARCHECKSUMG/
1787 . 'def '/
1788 DATA IVARCHECKSUMG
1789 . / 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1790C=======================================================================
1791 TAGP=0
1792 TAGS=0
1793 IBID = 0
1794 BID = 0
1795 NSNE = 0
1796C-----------------------------------------------
1797C Brick
1798C ------------------------------------------
1799C (VARS1, VARS2, VARS3, VARS4, VARS5, VARS6, VARS7 )----> VARS
1800 DO I=1,NVARS1
1801 VARS(I) = VARS1(I)
1802 ENDDO
1803C
1804 DO I=1,NVARS2
1805 VARS(NVARS1+ I) = VARS2(I)
1806 ENDDO
1807C
1808 DO I=1,NVARS3
1809 VARS(NVARS1+ NVARS2 + I) = VARS3(I)
1810 ENDDO
1811 DO I=1,NVARS4
1812 VARS(NVARS1+ NVARS2 + NVARS3 + I) = VARS4(I)
1813 ENDDO
1814 DO I=1,NVARS5
1815 VARS(1624 + I) = VARS5(I)
1816 ENDDO
1817 DO I=1,NVARS6
1818 VARS(1646 + I) = VARS6(I)
1819 ENDDO
1820 DO I=1,NVARS7
1821 VARS(98846 + I) = VARS7(I)
1822 ENDDO
1823 DO I=1,NVARS8
1824 VARS(239030 + I) = VARS8(I)
1825 ENDDO
1826 DO I=1,NVARS9
1827 VARS(239546 + I) = VARS9(I)
1828 ENDDO
1829 DO I=1,NVARSNLOC
1830 VARS(239552 + I) = VARSNLOC(I)
1831 ENDDO
1832 DO I=1,NVARS10
1833 VARS(239554 + I) = VARS10(I)
1834 ENDDO
1835C (VARSG1, VARSG2, VARSG3 )----> VARSG
1836 DO I=1,202
1837 VARSG(I) = VARSG1(I)
1838 ENDDO
1839C
1840 DO I=1,16200
1841 VARSG(202 + I) = VARSG2(I)
1842 ENDDO
1843C
1844 DO I=1,16200
1845 VARSG(16402 + I) = VARSG3(I)
1846 ENDDO
1847c
1848 DO I=1,1
1849 VARSG(39766 + I) = VARSG4(I)
1850 ENDDO
1851c
1852 DO I=1,18
1853 DO J = 1,NVARSG
1854 IVARSG(I,J)=0
1855 ENDDO
1856 ENDDO
1857C Option 1
1858 DO J=1,9
1859 IVARSG(J,1)=J
1860 ENDDO
1861 IVARSG(10,1)= 12
1862 IVARSG(11,1)= 13
1863C Option 2+3
1864 IVSG = 34
1865 DO J=1,6
1866 IVARSG(J,2)= J+1
1867 IVARSG(J,3)= IVSG + J
1868 ENDDO
1869 IVSG = IVSG + 6
1870C Option 4-->15 IVSG = 40
1871 DO I=4,15
1872 DO J=1,8
1873 IVARSG(J,I)= IVSG + J
1874 ENDDO
1875 IVSG = IVSG + 8
1876 ENDDO
1877C Option 16--->19 IVSG = 136
1878 DO I=16,18
1879 DO J=1,18
1880 IVARSG(J,I)= IVSG + J
1881 ENDDO
1882 IVSG= IVSG + 18
1883 ENDDO
1884C Option 19 IVSG = 190
1885 DO J=1,6
1886 IVARSG(J,19)= IVSG + J
1887 ENDDO
1888 IVSG = IVSG + 6
1889C options 20--->100 IVSG=196
1890 DO I=20,100
1891 DO J=1,7
1892 IVARSG(J,I)= IVSG + J
1893 ENDDO
1894 IVSG = IVSG + 7
1895 ENDDO
1896C Options 101 ---> 118 IVSG = 763
1897 DO I=101,118
1898 DO J=1,7
1899 IVARSG(J,I)= IVSG + J
1900 ENDDO
1901 IVSG = IVSG + 7
1902 ENDDO
1903C Options 118 ---> 199 IVSG = 889
1904 DO I=119,199
1905 DO J=1,9
1906 IVARSG(J,I)= IVSG + J
1907 ENDDO
1908 IVSG = IVSG + 9
1909 ENDDO
1910 DO J=1,6
1911 IVARSG(J,200)= IVSG + J
1912 ENDDO
1913 IVSG = IVSG + 6
1914C Options 201 ---> 202 IVSG = 1624
1915 DO I=201,202
1916 DO J=1,11
1917 IVARSG(J,I)= IVSG + J
1918 ENDDO
1919 IVSG = IVSG + 11
1920 ENDDO
1921C Options 203 ---> 19984 IVSG = 1646
1922 DO I=203,19984
1923 DO J=1,6
1924 IVARSG(J,I)= IVSG + J
1925 ENDDO
1926 IVSG = IVSG + 6
1927 ENDDO
1928C Options 203 ---> 19984 IVSG = 120338
1929 DO I=19985,39766
1930 DO J=1,6
1931 IVARSG(J,I)= IVSG + J
1932 ENDDO
1933 IVSG = IVSG + 6
1934 ENDDO
1935C Locstrn option
1936 DO J=1,6
1937 IVARSG(J,39767)= 239030 + J
1938 ENDDO
1939C
1940 DO I=1,208
1941 VARC(I) = VARC1(I)
1942 ENDDO
1943 DO I=1,470
1944 VARC(208 + I) = VARC2(I)
1945 ENDDO
1946Ccc USRII_KK (II=1,20, KK=6,99
1947 DO II=1,20
1948 WRITE(CHII,'(i2.2)')II
1949 DO KK= 6,99
1950 WRITE(CHKK,'(i2.2)')KK
1951 IF(CHII(1:1) == '0')THEN
1952 VAR = 'usr'//CHII(2:2)//'_'
1953 IF(CHKK(1:1) == '0')THEN
1954 VAR = VAR(1:5)//CHKK(2:2)
1955 ELSE
1956 VAR = VAR(1:5)//CHKK
1957 ENDIF
1958 ELSE
1959 VAR ='usr'//CHII//'_'
1960 IF(CHKK(1:1) == '0')THEN
1961 VAR = VAR(1:6)//CHKK(2:2)
1962 ELSE
1963 VAR = VAR(1:6)//CHKK
1964 ENDIF
1965 ENDIF
1966 VARC(678 + (II -1)*94 + KK - 5) = VAR
1967 ENDDO
1968 ENDDO
1969Ccc USRII_KK(II=21,60, KK=1,99)
1970 DO II=21,60
1971 WRITE(CHII,'(i2.2)')II
1972 DO KK= 1,99
1973 WRITE(CHKK,'(i2.2)')KK
1974 IF(CHII(1:1) == '0')THEN
1975 VAR = 'usr'//CHII(2:2)//'_'
1976 IF(CHKK(1:1) == '0')THEN
1977 VAR = VAR(1:5)//CHKK(2:2)
1978 ELSE
1979 VAR = VAR(1:5)//CHKK
1980 ENDIF
1981 ELSE
1982 VAR ='usr'//CHII//'_'
1983 IF(CHKK(1:1) == '0')THEN
1984 VAR = VAR(1:6)//CHKK(2:2)
1985 ELSE
1986 VAR = VAR(1:6)//CHKK
1987 ENDIF
1988 ENDIF
1989 VARC(2558 + (II - 21)*99 + KK ) = VAR
1990 ENDDO
1991 ENDDO
1992C USII_jKK (II=01, 60, j=1,4, KK=01,99)
1993 DO KK = 1, 99
1994 WRITE(CHKK,'(i2.2)')KK
1995 DO J = 1,4
1996 JJ = J*100 + KK
1997 WRITE(CHJJ,'(i3.3)')JJ
1998 DO II = 1,60
1999 WRITE(CHII,'(i2.2)')II
2000 VAR ='us'//CHII//'_'//CHJJ
2001 VARC(6518 + II + (KK - 1)*240 + (J-1)*60)=VAR
2002 ENDDO
2003 ENDDO
2004 ENDDO
2005 DO I = 1,104
2006 VARC(30278 + I ) = VARC2(470 + I)
2007 ENDDO
2008C
2009C viscous stress for shell 5*100
2010CC USII_jKK (II=01, 60, j=1,4, KK=01,99)
2011 DO I = 1,500
2012 VARC(30382 + I ) = VARC3(I)
2013 ENDDO
2014C NVC_old = 30882
2015C stress, strain, damage for interply
2016C SZILYjKK (j=1,4, KK=001,200)
2017C SYZILYjKK (j=1,4, KK=001,200)
2018C SXZILYjKK (j=1,4, KK=001,200)
2019C EZILYjKK (j=1,4, KK=001,200)
2020C EYZILYjKK (j=1,4, KK=001,200)
2021C EXZILYjKK (j=1,4, KK=001,200)
2022C SZILY_KK (KK=001,200)
2023C SYZILY_KK ( KK=001,200)
2024C SXZILY_KK ( KK=001,200)
2025C EZILY_KK (KK=001,200)
2026C EYZILY_KK ( KK=001,200)
2027C EXZILY_KK ( KK=001,200)
2028C DAM1ILYKK ( KK=001,200)
2029C DAM2ILYKK (KK=001,200)
2030C DAM3ILYKK ( KK=001,200)
2031C DAM4ILYKK ( KK=001,200)
2032C EINTILKK
2033 DO KK = 1, 199
2034C 3 stress + 3 strain + 3 average stress + 3 average strain + 4 max damage
2035C + EINT
2036 II = 30882 + (KK -1)*(6*4 + 11)
2037 DO J = 1,4
2038 JJ = J*1000 + KK
2039 WRITE(CHJKK ,'(i4.4)')JJ
2040 DO L=1,6
2041 IF (L==1) VAR= 'szily'//CHJKK
2042 IF (L==2) VAR ='syzily'//CHJKK
2043 IF (L==3) VAR ='sxzily'//CHJKK
2044C
2045 IF (L==4) VAR= 'ezily'//CHJKK
2046 IF (L==5) VAR ='eyzily'//CHJKK
2047 IF (L==6) VAR ='exzily'//CHJKK
2048 VARC(II + L ) = VAR
2049 ENDDO
2050 II = II + 6
2051 ENDDO
2052 WRITE(CHJJ,'(i3.3)')KK
2053 DO L = 1,6
2054 IF (L==1) VAR = 'szily'//'_'//CHJJ
2055 IF (L==2) VAR = 'syzily'//'_'//CHJJ
2056 IF (L==3) VAR = 'sxzily'//'_'//CHJJ
2057 IF (L==4) VAR = 'ezily'//'_'//CHJJ
2058 IF (L==5) VAR = 'eyzily'//'_'//CHJJ
2059 IF (L==6) VAR = 'exzily'//'_'//CHJJ
2060 VARC(II + L ) = VAR
2061 ENDDO
2062 II =II +6
2063 DO L = 1,4
2064 IF (L==1) VAR = 'damily'//CHJJ
2065 IF (L==2) VAR = 'dam2ily'//CHJJ
2066 IF (L==3) VAR = 'dam3ily'//CHJJ
2067 IF (L==4) VAR = 'dam4ily'//CHJJ
2068 VARC(II + L ) = VAR
2069 ENDDO
2070 II = II + 4
2071 VAR = 'eintily'//CHJJ
2072 VARC(II + 1 ) = VAR
2073 II = II + 1
2074 ENDDO
2075C---pinching
2076 DO I=1,7
2077 VARC(37847+I) = VARCPINCH(I)
2078 ENDDO
2079C---non-local
2080 DO I=1,2
2081 VARC(37854+I) = VARCNLOC(I)
2082 ENDDO
2083C
2084C NVC = 30882 + 199*34 = 37847
2085Cccc array dimension of VARC =
2086C dimension_varc =83 + 99*60 + 5*99 + 4*99*60 + 104 + 500 ! 30882
2087C + (4*6 + 11)*199 = 37847
2088C + (4*6 + 11)*199 = 37847
2089C-----------------------------------------------------
2090C GENERALIZED KEYWORDS
2091C-----------------------------------------------------
2092 DO I =1,167
2093 VARCG(I)= VARCG1(I)
2094 ENDDO
2095 DO I = 1,407
2096 VARCG(167 + I) = VARCG2(I)
2097 ENDDO
2098 DO I = 1,100
2099 VARCG(574 + I) = VARCG3(I)
2100 ENDDO
2101 VARCG(675) = VARCGPINCH(1)
2102C
2103 DO I=1,18
2104 DO J = 1, NVARCG
2105 IVARCG(I,J)=0
2106 ENDDO
2107 ENDDO
2108C keyword 1
2109 DO I=1,3
2110 IVARCG(I,1) = I
2111 IVARCG(I+3,1) = 5 + I
2112 IVARCG(I+6,1) = 8 + I
2113 ENDDO
2114 IVARCG(10,1) = 13
2115 IVARCG(11,1) = 14
2116C keyword 2 + 3
2117 DO I = 1,8
2118 IVARCG(I,2) = I
2119 IVARCG(I,3) = 14 + I
2120 ENDDO
2121C keyword 4
2122 IVARCG(1,4) = 13
2123 IVARCG(2,4) = 14
2124C keyword 5,6,7
2125 DO I =1,18
2126 IVARCG(I,5)= 23 + I
2127 IVARCG(I,6)= 41 + I
2128 IVARCG(I,7)= 59 + I
2129 ENDDO
2130C keyword 8
2131 DO I =1,6
2132 IVARCG(I,8)= 77 + I
2133 ENDDO
2134 IVCG = 83
2135C keyword 9 -----> 28
2136 II = 678
2137 DO J = 9, 28
2138 DO I = 1,5
2139 IVARCG(I,J) = IVCG + I
2140 ENDDO
2141 IVCG = IVCG + 5
2142 DO I=6,18
2143 IVARCG(I,J) = II + I - 5
2144 ENDDO
2145 II = II + 94
2146 ENDDO
2147C keyword 29 -----> 127
2148 DO J = 29, 127
2149 DO I = 1,5
2150 IVARCG(I,J)= IVCG + I
2151 ENDDO
2152 IVCG = IVCG + 5
2153 ENDDO
2154C c ivcg = 183 + 5*99 + 20*94
2155 IVCG = IVCG + 20*94
2156C keyword 128 - 167
2157 DO J=128,167
2158 DO I=1,18
2159 IVARCG(I,J) = IVCG + I
2160 ENDDO
2161 IVCG = IVCG + 99
2162 ENDDO
2163C keyword 167 -----> 563
2164 DO J = 168, 563
2165 DO I = 1,18
2166 IVARCG(I,J)= IVCG + I
2167 ENDDO
2168 IVCG = IVCG + 60
2169 ENDDO
2170C keyword 564
2171 DO I = 1,5
2172 IVARCG(I,564)= IVCG + I
2173 ENDDO
2174 IVCG = IVCG + 5
2175c keyword 565 -----> 573
2176 DO J = 565,573
2177 DO I = 1,10
2178 IVARCG(I,J)= IVCG + 1
2179 IVCG = IVCG + 1
2180 ENDDO
2181 ENDDO
2182c keyword 574
2183 DO I = 1,9
2184 IVARCG(I,574)= IVCG + 1
2185 IVCG = IVCG + 1
2186 ENDDO
2187C keyword de 575 -----> 674
2188 DO J = 575, 674
2189 DO I = 1,5
2190 IVARCG(I,J)= IVCG + I
2191 ENDDO
2192 IVCG = IVCG + 5
2193 ENDDO
2194C pinching
2195 DO I = 1,7
2196 IVARCG(I,675)= 37847 + I
2197 ENDDO
2198C
2199C to ctrole the value of ivcg
2200C icvg = 30882
2201C end . ivcg = ivcg + 4*99*60 + 5 + 99 + 100
2202C--------------fin --------------------
2203C--------------fin --------------------
2204C ---- node + displacement by ply for shell
2205 DO I=1,NVARN1
2206 VARN(I) = VARN1(I)
2207 ENDDO
2208 DO J= 1,200
2209 WRITE(CHJJ,'(i3.3)')J
2210 VAR ='dx_'//CHJJ
2211 VARN(NVARN1 + (J-1)*3 + 1 )= VAR
2212 VARN1A((J-1)*3 + 1 )= VAR
2213 VAR ='dy_'//CHJJ
2214 VARN(NVARN1 + (J-1)*3 + 2) = VAR
2215 VARN1A((J-1)*3 + 2 )= VAR
2216 VAR ='dz_'//CHJJ
2217 VARN(NVARN1 + (J-1)*3 + 3) = VAR
2218 VARN1A((J-1)*3 + 3 )= VAR
2219 ENDDO
2220 DO I=1,NVARN2
2221 VARN(NVARN1+NVARN1A+I) = VARN2(I)
2222 ENDDO
2223 DO I=1,NVARNPINCH
2224 VARN(NVARN1+NVARN1A+NVARN2+I) = VARNPINCH(I)
2225 ENDDO
2226 NNOD = 0
2227 CHNOD = 78
2228 NSHEL = 0
2229 NTRUS = 0
2230 NBEAM = 0
2231 NSPRG = 0
2232 NCHTG = 0
2233 NCHINT = 0
2234 NCHJNT = 0
2235C-----------------------------------------------
2236 NSNOD=0
2237 NSELS=0
2238 NSELQ=0
2239 NSELC=0
2240 NSELT=0
2241 NSELP=0
2242 NSELR=0
2243 NSELTG=0
2244 NSELU=0
2245 NSMAT=0
2246 NSRBY=0
2247C
2248 NSFXBY=0
2249C
2250 NSMOD=0
2251C for parts and subsets:
2252 NTHWA=100
2253 NSELSP=0
2254 NSELRT=0
2255C
2256 OK=0
2257 IS_AVAILABLE = .FALSE.
2258 IGS=0
2259 IAD=1
2260 HM_NTHGRP=0
2261c
2262 ! Choose TH file type
2263 IF (IFLAG == 0) THEN
2264 THFILE = '/th'
2265 ELSEIF (IFLAG == 1) THEN
2266 THFILE = '/ath'
2267 ELSEIF (IFLAG == 2) THEN
2268 THFILE = '/bth'
2269 ELSEIF (IFLAG == 3) THEN
2270 THFILE = '/cth'
2271 ELSEIF (IFLAG == 4) THEN
2272 THFILE = '/dth'
2273 ELSEIF (IFLAG == 5) THEN
2274 THFILE = '/eth'
2275 ELSEIF (IFLAG == 6) THEN
2276 THFILE = '/fth'
2277 ELSEIF (IFLAG == 7) THEN
2278 THFILE = '/gth'
2279 ELSEIF (IFLAG == 8) THEN
2280 THFILE = '/hth'
2281 ELSEIF (IFLAG == 9) THEN
2282 THFILE = '/ith'
2283 ENDIF
2284C
2285C
2286C-----------------------------------------------
2287! ************************** !
2288! /TH read with hm reader !
2289! ************************** !
2290C-----------------------------------------------
2291 !-------------------------------------------
2292 ! /TH/ACCEL
2293 !-------------------------------------------
2294 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/accel' ,NTHACCEL )
2295 IF (NTHACCEL > 0) THEN
2296 CALL HM_OPTION_START(TRIM(THFILE)//'/accel')
2297 HM_NTHGRP = HM_NTHGRP + NTHACCEL
2298 DO I = 1, NTHACCEL
2299 TITR = ''
2300 CALL HM_OPTION_READ_KEY(LSUBMODEL,
2301 . OPTION_ID = ID,
2302 . OPTION_TITR = TITR)
2303C
2304 IGS = IGS+1
2305 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2306 ITHGRP(1,IGS) = ID
2307 NVALL = 3
2308 KEY = 'accel '
2309 CALL HM_READ_THGRKI(
2310 1 108 ,KEY ,PTR_NOPT_ACC ,
2311 2 IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2312 3 NVARAC ,VARAC ,NACCELM ,VARACG ,
2313 4 NVARACG ,IVARACG ,NSNE ,NVALL,ITHVAR ,
2314 5 FLAGABF ,NVARABF ,NOM_OPT,IGS ,LSUBMODEL,LITHBUFMX)
2315 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2316 ENDDO
2317 ENDIF
2318 !-------------------------------------------
2319 ! /TH/INTER
2320 !-------------------------------------------
2321 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/inter' ,NTHINTER )
2322 IF (NTHINTER > 0) THEN
2323 CALL HM_OPTION_START(TRIM(THFILE)//'/inter')
2324 HM_NTHGRP = HM_NTHGRP + NTHINTER
2325 DO I = 1, NTHINTER
2326 TITR = ''
2327 CALL HM_OPTION_READ_KEY(LSUBMODEL,
2328 . OPTION_ID = ID,
2329 . OPTION_TITR = TITR)
2330C
2331 IGS = IGS+1
2332 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2333 ITHGRP(1,IGS) = ID
2334 CHELT = 258
2335 CNT1 = 0
2336 NVALL = 7
2337 KEY = 'inter '
2338 CALL HM_READ_THGRKI(
2339 1 101 ,KEY ,PTR_NOPT_INTER,
2340 2 IAD ,IFI ,ITHGRP(1,IGS) ,ITHBUF ,
2341 3 NVARIN ,VARIN ,NINTER+NINTSUB,VARING ,NVARING ,
2342 4 IVARING ,NSNE ,NVALL ,ITHVAR,FLAGABF ,NVARABF ,
2343 5 NOM_OPT ,IGS ,LSUBMODEL, LITHBUFMX)
2344 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2345 DO JJ = ITHGRP(7,IGS),ITHGRP(5,IGS)
2346 IF(ITHBUF(JJ) ==29) INTERFACES%PARAMETERS%INTCAREA =1
2347 ENDDO
2348 IF(INTERFACES%PARAMETERS%INTCAREA > 0 ) THEN
2349 DO JJ = ITHGRP(5,IGS),ITHGRP(5,IGS)+ITHGRP(4,IGS)-1
2350 NIN = ITHBUF(JJ)
2351 IF(NIN <= NINTER) THEN
2352.OR. IF(IPARI(7,NIN)==24 IPARI(7,NIN)==25) IPARI(99,NIN) = 1
2353 ENDIF
2354 ENDDO
2355 ENDIF
2356 ENDDO
2357 ENDIF
2358 !-------------------------------------------
2359 ! /TH/RWALL
2360 !-------------------------------------------
2361 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/rwall' ,NTHRWALL )
2362 IF (NTHRWALL > 0) THEN
2363 CALL HM_OPTION_START(TRIM(THFILE)//'/rwall')
2364 HM_NTHGRP = HM_NTHGRP + NTHRWALL
2365 DO I = 1, NTHRWALL
2366 TITR = ''
2367 CALL HM_OPTION_READ_KEY(LSUBMODEL,
2368 . OPTION_ID = ID,
2369 . OPTION_TITR = TITR)
2370C
2371 IGS = IGS+1
2372 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2373 ITHGRP(1,IGS) = ID
2374 NVALL = 6
2375 KEY = 'rwall '
2376 CALL HM_READ_THGRKI(
2377 1 102 ,KEY ,PTR_NOPT_RWALL,
2378 2 IAD ,IFI ,ITHGRP(1,IGS) ,ITHBUF ,
2379 3 NVARRW ,VARRW ,NRWALL ,VARRWG ,NVARRWG ,
2380 4 IVARRWG ,NSNE ,NVALL,ITHVAR ,FLAGABF ,NVARABF ,
2381 5 NOM_OPT ,IGS ,LSUBMODEL, LITHBUFMX)
2382 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2383 ENDDO
2384 ENDIF
2385 !-------------------------------------------
2386 ! /TH/SECTIO
2387 !-------------------------------------------
2388 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/sectio' ,NTHSECTIO)
2389 IF (NTHSECTIO > 0) THEN
2390 CALL HM_OPTION_START(TRIM(THFILE)//'/sectio')
2391 HM_NTHGRP = HM_NTHGRP + NTHSECTIO
2392 DO I = 1, NTHSECTIO
2393 TITR = ''
2394 CALL HM_OPTION_READ_KEY(LSUBMODEL,
2395 . OPTION_ID = ID,
2396 . OPTION_TITR = TITR)
2397C
2398 IGS = IGS+1
2399 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2400 ITHGRP(1,IGS) = ID
2401 NVALL = 13
2402 KEY = 'sectio '
2403 CALL HM_READ_THGRKI(
2404 1 104 ,KEY ,PTR_NOPT_SECT,
2405 2 IAD ,IFI ,ITHGRP(1,IGS) ,ITHBUF ,
2406 3 NVARSE ,VARSE ,NSECT ,VARSEG ,NVARSEG ,
2407 4 IVARSEG ,NSNE ,NVALL ,ITHVAR ,FLAGABF ,NVARABF ,
2408 5 NOM_OPT ,IGS ,LSUBMODEL, LITHBUFMX)
2409 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2410 ENDDO
2411 ENDIF
2412 !-------------------------------------------
2413 ! /TH/CLUSTER
2414 !-------------------------------------------
2415 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/cluster' ,NTHCLUS)
2416 IF (NTHCLUS > 0) THEN
2417 CALL HM_OPTION_START(TRIM(THFILE)//'/cluster')
2418 HM_NTHGRP = HM_NTHGRP + NTHCLUS
2419 DO I = 1, NTHCLUS
2420 TITR = ''
2421 CALL HM_OPTION_READ_KEY(LSUBMODEL,
2422 . OPTION_ID = ID,
2423 . OPTION_TITR = TITR)
2424C
2425 IGS = IGS+1
2426 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2427 ITHGRP(1,IGS) = ID
2428 NVALL = 11
2429 ITYP = 114
2430 KEY = 'cluster '
2431 CALL HM_READ_THGRKI(
2432 1 ITYP ,KEY ,PTR_NOPT_CLUS,
2433 3 IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2434 4 NVARCLUS ,VARCLUS ,NCLUSTER ,VARCLUSG ,NVARCLUSG ,
2435 5 IVARCLUSG,NSNE ,NVALL ,ITHVAR,FLAGABF ,NVARABF ,
2436 6 NOM_OPT ,IGS ,LSUBMODEL, LITHBUFMX )
2437 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2438 ENDDO
2439 ENDIF
2440 !-------------------------------------------
2441 ! /TH/BEAM
2442 !-------------------------------------------
2443 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/beam' ,NTHBEAM)
2444 IF (NTHBEAM > 0) THEN
2445 CALL HM_OPTION_START(TRIM(THFILE)//'/beam')
2446 HM_NTHGRP = HM_NTHGRP + NTHBEAM
2447 DO I = 1, NTHBEAM
2448 TITR = ''
2449 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
2450 IGS = IGS+1
2451 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2452 ITHGRP(1,IGS) = ID
2453 NVALL = 11
2454 ITYP = 5
2455 KEY = 'beam '
2456 CALL HM_READ_THGRNE(
2457 1 ITYP ,KEY ,ITAB ,ITABM1,IXP ,
2458 3 NIXP ,IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2459 4 NVARP ,VARP ,NUMELP,VARPG ,NVARPG ,
2460 5 IVARPG,NSELP,NVALL,IMERGE,ITHVAR ,ITHERM_FE,
2461 6 FLAGABF,NVARABF,RFI,LSUBMODEL,MAP_TABLES%IBEAMM, NUMELP)
2462
2463 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2464 ENDDO
2465 ENDIF
2466 !-------------------------------------------
2467 ! /TH/TRUS
2468 !-------------------------------------------
2469 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/truss' ,NTHTRUS)
2470 IF (NTHTRUS > 0) THEN
2471 CALL HM_OPTION_START(TRIM(THFILE)//'/truss')
2472 HM_NTHGRP = HM_NTHGRP + NTHTRUS
2473 DO I = 1, NTHTRUS
2474 TITR = ''
2475 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
2476 IGS = IGS+1
2477 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2478 ITHGRP(1,IGS) = ID
2479 NVALL = 6
2480 ITYP = 4
2481 KEY = 'truss '
2482 CALL HM_READ_THGRNE(
2483 1 ITYP ,KEY ,ITAB ,ITABM1,IXT ,
2484 3 NIXT ,IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2485 4 NVART ,VART ,NUMELT,VARTG ,NVARTG ,
2486 5 IVARTG,NSELT,NVALL,IMERGE,ITHVAR ,ITHERM_FE,
2487 6 FLAGABF,NVARABF,RFI,LSUBMODEL,MAP_TABLES%ITRUSSM, NUMELT)
2488 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2489 ENDDO
2490 ENDIF
2491 !-------------------------------------------
2492 ! /TH/SPRI
2493 !-------------------------------------------
2494 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/spring' ,NTHSPRING)
2495 IF (NTHSPRING > 0) THEN
2496 CALL HM_OPTION_START(TRIM(THFILE)//'/spring')
2497 HM_NTHGRP = HM_NTHGRP + NTHSPRING
2498 DO I = 1, NTHSPRING
2499 TITR = ''
2500 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
2501 IGS = IGS+1
2502 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2503 ITHGRP(1,IGS) = ID
2504 NVALL = 16
2505 ITYP = 6
2506 KEY = 'spring '
2507 CALL HM_READ_THGRNE(
2508 1 ITYP ,KEY ,ITAB ,ITABM1,IXR ,
2509 3 NIXR ,IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2510 4 NVARR ,VARR ,NUMELR,VARRG ,NVARRG ,
2511 5 IVARRG,NSELR,NVALL,IMERGE,ITHVAR ,ITHERM_FE,
2512 6 FLAGABF,NVARABF,RFI,LSUBMODEL,MAP_TABLES%ISPRINGM, NUMELR)
2513 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2514 ENDDO
2515 ENDIF
2516 !-------------------------------------------
2517 ! /TH/BRIC
2518 !-------------------------------------------
2519 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/bric' ,NTHBRIC)
2520 IF (NTHBRIC > 0) THEN
2521 CALL HM_OPTION_START(TRIM(THFILE)//'/bric')
2522 hm_nthgrp = hm_nthgrp + nthbric
2523 DO i = 1, nthbric
2524 titr = ''
2525 CALL hm_option_read_key(lsubmodel, option_id=id, option_titr=titr)
2526 igs = igs+1
2527 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2528 ithgrp(1,igs) = id
2529 nvall = 40
2530 ityp = 1
2531 key = 'BRICK '
2532 CALL hm_read_thgrne(
2533 1 ityp ,key ,itab ,itabm1,ixs ,
2534 3 nixs ,iad ,ifi ,ithgrp(1,igs),ithbuf ,
2535 4 nvars ,vars ,numels,varsg ,nvarsg ,
2536 5 ivarsg,nsels,nvall,imerge,ithvar ,itherm_fe,
2537 6 flagabf,nvarabf,rfi,lsubmodel,map_tables%ISOLM, numels)
2538 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2539 ENDDO
2540 ENDIF
2541 !-------------------------------------------
2542 ! /TH/NODE
2543 !-------------------------------------------
2544 CALL hm_option_count(trim(thfile)//'/NODE' ,nthnode)
2545 IF (nthnode > 0) THEN
2546 CALL hm_option_start(trim(thfile)//'/NODE')
2547 hm_nthgrp = hm_nthgrp + nthnode
2548 DO i = 1, nthnode
2549 titr = ''
2550 CALL hm_option_read_key(lsubmodel, option_id=id, option_titr=titr)
2551 igs = igs+1
2552 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2553 ithgrp(1,igs) = id
2554 nvall = 9
2555 ityp = 0
2556 key = 'NODE '
2557 CALL hm_read_thgrne(
2558 1 ityp ,key ,itab ,itabm1 ,ixs ,
2559 3 nixs ,iad ,ifi ,ithgrp(1,igs),ithbuf ,
2560 4 nvarn ,varn ,0 ,varng ,nvarng ,
2561 5 ivarng,nsnod ,nvall ,imerge ,ithvar ,itherm_fe,
2562 6 flagabf,nvarabf,rfi,lsubmodel,ibid,ibid)
2563C decode les numeros de skews.
2564 nne =ithgrp(4,igs)
2565 iad2=ithgrp(5,igs)+nne
2566 DO k=1,nne
2567 isk=ithbuf(iad2)
2568 DO 100 j=1,(numskw+1)+nsubmod
2569 IF(isk == iskwn(4,j)) THEN
2570 ithbuf(iad2)=j-1
2571 GO TO 110
2572 ENDIF
2573 100 CONTINUE
2574C decode les numeros de frames.
2575 DO 101 j=(numskw+1)+min(iun,nspcond)*numsph+1+nsubmod,
2576 . (numskw+1)+min(iun,nspcond)*numsph+nsubmod+
2577 . (numfram+1)
2578 IF(isk == iskwn(4,j)) THEN
2579 ithbuf(iad2)=j-1
2580 GO TO 110
2581 ENDIF
2582 101 CONTINUE
2583 CALL ancmsg(msgid=434,
2584 . msgtype=msgerror,
2585 . anmode=aninfo_blind_1,
2586 . i1=id,
2587 . c1=titr,
2588 . i2=isk)
2589 110 CONTINUE
2590 iad2=iad2+1
2591 ENDDO
2592 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2593 ENDDO
2594 ENDIF
2595 !-------------------------------------------
2596 ! /TH/SHEL
2597 !-------------------------------------------
2598 CALL hm_option_count(trim(thfile)//'/SHEL' ,nthshel)
2599 IF (nthshel > 0) THEN
2600 CALL hm_option_start(trim(thfile)//'/SHEL')
2601 hm_nthgrp = hm_nthgrp + nthshel
2602 DO i = 1, nthshel
2603 titr = ''
2604 CALL hm_option_read_key(lsubmodel, option_id=id, option_titr=titr)
2605 igs = igs+1
2606 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2607 ithgrp(1,igs) = id
2608 nvall = 22
2609 ityp = 3
2610 key = 'SHELL '
2611 CALL hm_read_thgrne(
2612 1 ityp ,key ,itab ,itabm1,ixc ,
2613 3 nixc ,iad ,ifi ,ithgrp(1,igs),ithbuf ,
2614 4 nvarc ,varc ,numelc,varcg ,nvarcg ,
2615 5 ivarcg,nselc,nvall,imerge,ithvar ,itherm_fe,
2616 6 flagabf,nvarabf,rfi,lsubmodel,map_tables%ISH4NM, numelc)
2617C decode les numeros de skews.
2618 nne =ithgrp(4,igs)
2619 iad2=ithgrp(5,igs)+3*nne
2620 DO k=1,nne
2621 isk=ithbuf(iad2)
2622 DO 200 j=1,(numskw+1)+nsubmod
2623 IF(isk == iskwn(4,j)) THEN
2624 ithbuf(iad2)=j-1
2625 GO TO 210
2626 ENDIF
2627 200 CONTINUE
2628C frame numbers user -> internal
2629 DO 201 j=(numskw+1)+min(iun,nspcond)*numsph+1+nsubmod,
2630 . (numskw+1)+min(iun,nspcond)*numsph+nsubmod+
2631 . (numfram+1)
2632 IF(isk == iskwn(4,j)) THEN
2633 ithbuf(iad2)=j-1
2634 GO TO 210
2635 ENDIF
2636 201 CONTINUE
2637 CALL ancmsg(msgid=434,
2638 . msgtype=msgerror,
2639 . anmode=aninfo_blind_1,
2640 . i1=id,
2641 . c1=titr,
2642 . i2=isk)
2643 210 CONTINUE
2644 iad2=iad2+1
2645 ENDDO
2646 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2647 ENDDO
2648 ENDIF
2649 !-------------------------------------------
2650 ! /TH/SH3N
2651 !-------------------------------------------
2652 CALL hm_option_count(trim(thfile)//'/SH3N' ,nthsh3n)
2653 IF (nthsh3n > 0) THEN
2654 CALL hm_option_start(trim(thfile)//'/SH3N')
2655 hm_nthgrp = hm_nthgrp + nthsh3n
2656 DO i = 1, nthsh3n
2657 titr = ''
2658 CALL hm_option_read_key(lsubmodel, option_id=id, option_titr=titr)
2659 igs = igs+1
2660 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2661 ithgrp(1,igs) = id
2662 nvall = 22
2663 ityp = 7
2664 key = 'SH3N '
2665 CALL hm_read_thgrne(
2666 1 ityp ,key ,itab ,itabm1,ixtg ,
2667 3 nixtg ,iad ,ifi ,ithgrp(1,igs),ithbuf ,
2668 4 nvarc ,varc ,numeltg,varcg ,nvarcg ,
2669 5 ivarcg,nseltg,nvall,imerge,ithvar ,itherm_fe,
2670 6 flagabf,nvarabf,rfi,lsubmodel,map_tables%ISH3NM, numeltg)
2671C skew numbers user -> internal
2672 nne =ithgrp(4,igs)
2673 iad2=ithgrp(5,igs)+3*nne
2674 DO k=1,nne
2675 isk=ithbuf(iad2)
2676 DO 300 j=1,(numskw+1)+nsubmod
2677 IF(isk == iskwn(4,j)) THEN
2678 ithbuf(iad2)=j-1
2679 GO TO 310
2680 ENDIF
2681 300 CONTINUE
2682C frame numbers user -> internal
2683 DO 301 j=(numskw+1)+min(iun,nspcond)*numsph+1+nsubmod,
2684 . (numskw+1)+min(iun,nspcond)*numsph+nsubmod+
2685 . (numfram+1)
2686 IF(isk == iskwn(4,j)) THEN
2687 ithbuf(iad2)=j-1
2688 GO TO 310
2689 ENDIF
2690 301 CONTINUE
2691 CALL ancmsg(msgid=434,
2692 . msgtype=msgerror,
2693 . anmode=aninfo_blind_1,
2694 . i1=id,
2695 . c1=titr,
2696 . i2=isk)
2697 310 CONTINUE
2698 iad2=iad2+1
2699 ENDDO
2700 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2701 ENDDO
2702 ENDIF
2703 !-------------------------------------------
2704 ! /TH/RBODY
2705 !-------------------------------------------
2706 CALL hm_option_count(trim(thfile)//'/RBODY' ,nthrbody)
2707 IF (nthrbody > 0) THEN
2708 CALL hm_option_start(trim(thfile)//'/RBODY')
2709 hm_nthgrp = hm_nthgrp + nthrbody
2710 DO i = 1, nthrbody
2711 titr = ''
2712 CALL hm_option_read_key(lsubmodel, option_id=id, option_titr=titr)
2713 igs = igs+1
2714 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2715 ithgrp(1,igs) = id
2716 nvall = 15
2717 ityp = 103
2718 key = 'RBODY '
2719 CALL hm_read_thgrki_rbody(
2720 1 ityp ,key ,0,
2721 3 iad ,ifi ,ithgrp(1,igs),ithbuf ,
2722 4 nvarrb ,varrb ,nrbody ,varrbg ,nvarrbg,
2723 5 ivarrbg,nsrby ,nvall ,ithvar ,flagabf,nvarabf,
2724 6 nom_opt,igs ,npby ,lsubmodel )
2725 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2726 ENDDO
2727 ENDIF
2728 !-------------------------------------------
2729 ! /TH/MONVOL
2730 !-------------------------------------------
2731 CALL hm_option_count(trim(thfile)//'/MONV' ,nthmonvol)
2732 IF (nthmonvol > 0) THEN
2733 CALL hm_option_start(trim(thfile)//'/MONV')
2734 hm_nthgrp = hm_nthgrp + nthmonvol
2735 DO i = 1, nthmonvol
2736 titr = ''
2737 CALL hm_option_read_key(lsubmodel, option_id=id, option_titr=titr)
2738 igs = igs+1
2739 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2740 ithgrp(1,igs) = id
2741 nvall = 9
2742 ityp = 107
2743 key = 'MONVOL '
2744 CALL hm_read_thgrki(
2745 1 ityp ,key ,ptr_nopt_monv,
2746 3 iad ,ifi ,ithgrp(1,igs),ithbuf ,
2747 4 nvarmv ,varmv ,nvolu ,varmvg,nvarmvg,
2748 5 ivarmvg,nsne,nvall,ithvar,flagabf,nvarabf,
2749 6 nom_opt,igs ,lsubmodel,lithbufmx )
2750 nvall = 9
2751 igs=igs+1
2752 titr1(1:5)='VENT '
2753 DO j=1,nchartitle-5
2754 titr1(j+5:j+5) = titr(j:j)
2755 ENDDO
2756 CALL fretitl(titr1,ithgrp(nithgr-ltitr+1,igs),ltitr)
2757 ithgrp(1,igs)=id_limit%TH
2758 id_limit%TH = id_limit%TH + 1
2759 CALL hm_thgrki_vent(
2760 1 ityp ,key ,ptr_nopt_monv,
2761 3 iad ,ifi ,ithgrp(1,igs),ithbuf ,
2762 4 nvarmvent ,varmvent ,nvolu ,varmventg, nvarmventg,
2763 5 ivarmventg,nsne ,nvall ,ithvar , flagabf ,nvarabf,
2764 6 nom_opt ,igs ,t_monvol ,nvarmvent, lsubmodel)
2765 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2766 ENDDO
2767 ENDIF
2768 !-------------------------------------------
2769 ! /TH/PART
2770 !-------------------------------------------
2771 !HM_NTHPART: total number of /TH/PART in 0000.rad
2772 CALL hm_option_count(trim(thfile)//'/PART' ,hm_nthpart)
2773 IF (hm_nthpart > 0) THEN
2774 CALL hm_option_start(trim(thfile)//'/PART')
2775 hm_nthgrp = hm_nthgrp + hm_nthpart
2776 numthpart = npart + nthpart
2777 DO i = 1, hm_nthpart
2778 CALL hm_read_thgrpa(
2779 . ipart ,lipart1 ,ithbuf ,iad ,ifi ,
2780 . varpa ,nvarpa ,varpag ,nvarpag ,numthpart ,
2781 . ivarpag ,pathid ,tagp ,iparth ,nparth ,
2782 . nvparth ,lsubmodel)
2783 ENDDO
2784 ENDIF
2785 !-------------------------------------------
2786 ! /TH/SUBS
2787 !-------------------------------------------
2788 !HM_NTHSUBS: total number of /TH/SUBS in 0000.rad
2789 CALL hm_option_count(trim(thfile)//'/SUBS' ,hm_nthsubs)
2790 IF (hm_nthsubs > 0) THEN
2791 CALL hm_option_start(trim(thfile)//'/SUBS')
2792 hm_nthgrp = hm_nthgrp + hm_nthsubs
2793 ! HM_NTHGRP: Number of /TH/ Group in 0000.rad
2794 DO i = 1, hm_nthsubs
2795 ! TAGS: Total number of SUBS for the TIME history output
2796 CALL hm_read_thgrpa_sub(
2797 . iad ,ifi ,ithbuf ,nvarpa ,varpa ,
2798 . varpag ,nvarpag ,ivarpag ,pathid ,
2799 . suthid ,tags ,subset ,ithflag ,lsubmodel)
2800 ENDDO
2801 ENDIF
2802 !-------------------------------------------
2803 ! /TH/QUAD
2804 !-------------------------------------------
2805 th_strain = 0 ! Time history strain flag for QUADS -> TH output STRAIN TENSOR
2806 CALL hm_option_count(trim(thfile)//'/QUAD' ,hm_nthquad )
2807 IF (hm_nthquad > 0) THEN
2808 CALL hm_option_start(trim(thfile)//'/QUAD')
2809 hm_nthgrp = hm_nthgrp + hm_nthquad
2810 DO i = 1, hm_nthquad
2811 titr = ''
2812 CALL hm_option_read_key(lsubmodel, option_id = id, option_titr = titr)
2813 igs = igs+1
2814 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2815 ithgrp(1,igs) = id
2816 nvall = 40
2817 key = 'QUAD '
2818 CALL hm_read_thgrne(
2819 1 2 ,key ,itab ,itabm1,ixq ,
2820 3 nixq ,iad ,ifi ,ithgrp(1,igs),ithbuf ,
2821 4 nvars ,vars ,numelq,varsg ,nvarsg ,
2822 5 ivarsg,nselq,nvall,imerge,ithvar ,itherm_fe,
2823 6 flagabf,nvarabf,rfi,lsubmodel,map_tables%IQUADM, numelq)
2824 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2825 ENDDO
2826 ENDIF
2827 !-------------------------------------------
2828 ! /TH/SPHCEL
2829 !-------------------------------------------
2830 CALL hm_option_count(trim(thfile)//'/SPHCEL' ,hm_nthsphcel )
2831 IF (hm_nthsphcel > 0) THEN
2832 CALL hm_option_start(trim(thfile)//'/SPHCEL')
2833 hm_nthgrp = hm_nthgrp + hm_nthsphcel
2834 DO i = 1, hm_nthsphcel
2835 titr = ''
2836 CALL hm_option_read_key(lsubmodel, option_id = id, option_titr = titr)
2837 igs = igs+1
2838 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2839 ithgrp(1,igs) = id
2840 nvall = 41
2841 key = 'SPHCEL '
2842 CALL hm_read_thgrne(
2843 1 51 ,key ,itab ,itabm1,kxsp ,
2844 2 nisp ,iad ,ifi ,ithgrp(1,igs),ithbuf ,
2845 4 nvarsph,varsph,numsph,varspg,nvarspg ,
2846 5 ivarspg,nselsp,nvall,imerge,ithvar ,itherm_fe,
2847 6 flagabf,nvarabf,rfi,lsubmodel,map_tables%ISPHM,numsph)
2848 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2849 ENDDO
2850 ENDIF
2851 !-------------------------------------------
2852 ! /TH/CYL_JO
2853 !-------------------------------------------
2854 CALL hm_option_count(trim(thfile)//'/CYL_JO' ,hm_nthcyljo)
2855 IF (hm_nthcyljo > 0) THEN
2856 CALL hm_option_start(trim(thfile)//'/CYL_JO')
2857 hm_nthgrp = hm_nthgrp + hm_nthcyljo
2858 DO i = 1, hm_nthcyljo
2859 titr = ''
2860 CALL hm_option_read_key(lsubmodel, option_id = id, option_titr = titr)
2861 igs = igs+1
2862 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2863 ithgrp(1,igs) = id
2864 nvall = 6
2865 key = 'CYL_JO '
2866 CALL hm_read_thgrki(
2867 1 105 ,key ,ptr_nopt_joint,
2868 2 iad ,ifi ,ithgrp(1,igs) ,ithbuf ,
2869 3 nvarjo ,varjo ,njoint,varjog,nvarjog,
2870 4 ivarjog ,nsne ,nvall ,ithvar ,flagabf ,nvarabf ,
2871 5 nom_opt ,igs ,lsubmodel, lithbufmx)
2872 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2873 ENDDO
2874 ENDIF
2875 !-------------------------------------------
2876 ! /TH/FXBODY
2877 !-------------------------------------------
2878 CALL hm_option_count(trim(thfile)//'/FXBODY' ,hm_nthfxbody)
2879 IF (hm_nthfxbody > 0) THEN
2880 CALL hm_option_start(trim(thfile)//'/FXBODY')
2881 hm_nthgrp = hm_nthgrp + hm_nthfxbody
2882 DO i = 1, hm_nthfxbody
2883 titr = ''
2884 CALL hm_option_read_key(lsubmodel, option_id = id, option_titr = titr)
2885 igs = igs+1
2886 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2887 ithgrp(1,igs) = id
2888 nvall = 0
2889 key = 'FXBODY '
2890 CALL hm_read_thgrki(
2891 1 111 ,key ,ptr_nopt_fxby,
2892 2 iad ,ifi ,ithgrp(1,igs) ,ithbuf ,
2893 3 nvarfx ,varfx ,nfxbody ,varfxg ,nvarfxg,
2894 4 ivarfxg ,nsne ,nvall ,ithvar ,flagabf ,nvarabf ,
2895 5 nom_opt ,igs ,lsubmodel, lithbufmx )
2896 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2897 ENDDO
2898 ENDIF
2899 !-------------------------------------------
2900 ! /TH/FRAME
2901 !-------------------------------------------
2902 CALL hm_option_count(trim(thfile)//'/FRAME' ,hm_nthframe)
2903 IF (hm_nthframe > 0) THEN
2904 CALL hm_option_start(trim(thfile)//'/FRAME')
2905 hm_nthgrp = hm_nthgrp + hm_nthframe
2906 DO i = 1, hm_nthframe
2907 titr = ''
2908 CALL hm_option_read_key(lsubmodel, option_id = id, option_titr = titr)
2909 igs = igs+1
2910 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2911 ithgrp(1,igs) = id
2912 nvall = 24
2913 key = 'FRAME '
2914 CALL hm_read_thgrki(
2915 1 110 ,key ,ptr_nopt_skw+numskw+1,
2916 2 iad ,ifi ,ithgrp(1,igs) ,ithbuf ,
2917 3 nvarfr ,varfr ,numfram+1 ,varfrg, nvarfrg ,
2918 4 ivarfrg ,nsne ,nvall ,ithvar ,flagabf ,nvarabf ,
2919 5 nom_opt ,igs ,lsubmodel, lithbufmx)
2920 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2921 ENDDO
2922 ENDIF
2923 !-------------------------------------------
2924 ! /TH/GAUGE
2925 !-------------------------------------------
2926 !HM_NTHGAUGE: total number of /TH/GAUGE in 0000.rad
2927 CALL hm_option_count(trim(thfile)//'/GAUGE' ,hm_nthgauge)
2928 IF (hm_nthgauge > 0) THEN
2929 CALL hm_option_start(trim(thfile)//'/GAUGE')
2930 ! HM_NTHGRP: Number of /TH/ Group in 0000.rad
2931 hm_nthgrp = hm_nthgrp + hm_nthgauge
2932 DO i = 1, hm_nthgauge
2933 titr = ''
2934 CALL hm_option_read_key(lsubmodel, option_id = id, option_titr = titr)
2935 igs = igs+1
2936 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2937 ithgrp(1,igs) = id
2938 nvall = 3
2939 key = 'GAUGE '
2940 CALL hm_read_thgrki(
2941 1 113 ,key ,ptr_nopt_gau,
2942 2 iad ,ifi ,ithgrp(1,igs) ,ithbuf ,
2943 3 nvargau ,vargau ,nbgauge ,vargaug,nvargaug ,
2944 4 ivargaug ,nsne ,nvall ,ithvar ,flagabf ,nvarabf ,
2945 5 nom_opt ,igs ,lsubmodel, lithbufmx)
2946 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2947 ENDDO
2948 ENDIF
2949 !-------------------------------------------
2950 ! /TH/SPH_FLOW
2951 !-------------------------------------------
2952 CALL hm_option_count(trim(thfile)//'/SPH_FLOW' ,hm_nthsphflow)
2953 IF (hm_nthsphflow > 0) THEN
2954 CALL hm_option_start(trim(thfile)//'/SPH_FLOW')
2955 hm_nthgrp = hm_nthgrp + hm_nthsphflow
2956 DO i = 1, hm_nthsphflow
2957 titr = ''
2958 CALL hm_option_read_key(lsubmodel, option_id = id, option_titr = titr)
2959 igs = igs+1
2960 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2961 ithgrp(1,igs) = id
2962 nvall = 1
2963 key = 'SPH_FLOW '
2964 CALL hm_read_thgrki(
2965 1 115 ,key ,ptr_nopt_sphio,
2966 2 iad ,ifi ,ithgrp(1,igs) ,ithbuf ,
2967 3 nvarflow,varflow,nsphio,varflowg,nvarflowg ,
2968 5 ivarflowg,nsne ,nvall ,ithvar ,flagabf ,nvarabf ,
2969 5 nom_opt ,igs ,lsubmodel,lithbufmx)
2970C error message if type of inout is not correct
2971 DO j=1,ithgrp(4,igs)
2972 inout = ithbuf(ithgrp(5,igs)+j-1)
2973 IF (inout>0) THEN
2974 IF (isphio(1,inout) == 1) THEN
2975 CALL ancmsg(msgid=1111,
2976 . msgtype=msgerror,
2977 . anmode=aninfo_blind_1,
2978 . i1=ithgrp(1,igs),
2979 . c1=titr,
2980 . i2=nom_opt(1,ptr_nopt_sphio+inout))
2981 ENDIF
2982 ENDIF
2983 END DO
2984 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2985 ENDDO
2986 ENDIF
2987 !-------------------------------------------
2988 ! /TH/NSTRAND
2989 !-------------------------------------------
2990 CALL hm_option_count(trim(thfile)//'/NSTRAND' ,hm_nthnstrand )
2991 IF (hm_nthnstrand > 0) THEN
2992 CALL hm_option_start(trim(thfile)//'/NSTRAND')
2993 hm_nthgrp = hm_nthgrp + hm_nthnstrand
2994 DO i = 1, hm_nthnstrand
2995 titr = ''
2996 CALL hm_option_read_key(lsubmodel, option_id=id, option_titr=titr)
2997 igs = igs+1
2998 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2999 ithgrp(1,igs) = id
3000 nvall = 4
3001 ityp = 100
3002 key = 'NSTRAND '
3003 CALL hm_read_thgrns(
3004 1 100 ,key ,itab ,itabm1,kxx ,
3005 3 ixx ,iad ,ifi ,ithgrp(1,igs),ithbuf ,
3006 4 nvarns ,varns ,varnsg ,nvarnsg ,ivarnsg,
3007 5 nselr ,ivns2r,nvall ,id,titr,
3008 6 ithvar,flagabf,nvarabf,lsubmodel)
3009 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
3010 ENDDO
3011 ENDIF
3012 !-------------------------------------------
3013 ! /TH/SURF
3014 !-------------------------------------------
3015 CALL hm_option_count(trim(thfile)//'/SURF' ,hm_nthsurf )
3016 IF (hm_nthsurf> 0) THEN
3017 CALL hm_option_start(trim(thfile)//'/SURF')
3018 hm_nthgrp = hm_nthgrp + hm_nthsurf
3019 DO i = 1, hm_nthsurf
3020 titr = ''
3021 CALL hm_option_read_key(lsubmodel, option_id=id, option_titr=titr)
3022 igs = igs+1
3023 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
3024 ithgrp(1,igs) = id
3025 nvall = 1
3026 ityp = 116
3027 key = 'SURF '
3028 CALL hm_read_thgrsurf(
3029 1 ityp ,key ,
3030 3 iad ,ifi ,ithgrp(1,igs),ithbuf ,
3031 4 nvarsurf,varsurf,nsurf,varsurfg,nvarsurfg ,
3032 5 ivarsurfg,nsne,nvall,ithvar,flagabf,nvarabf,
3033 6 igrsurf,igs,lsubmodel)
3034 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
3035 ENDDO
3036 ENDIF
3037 !-------------------------------------------
3038 ! /TH/TRIA
3039 !-------------------------------------------
3040 CALL hm_option_count(trim(thfile)//'/TRIA' ,hm_nthtria )
3041 IF (hm_nthtria > 0) THEN
3042 CALL hm_option_start(trim(thfile)//'/TRIA')
3043 hm_nthgrp = hm_nthgrp + hm_nthtria
3044 DO i = 1, hm_nthtria
3045 titr = ''
3046 CALL hm_option_read_key(lsubmodel, option_id=id, option_titr=titr)
3047 igs = igs+1
3048 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
3049 ithgrp(1,igs) = id
3050 nvall = 40
3051 ityp = 117
3052 key = 'TRIA '
3053 CALL hm_read_thgrne(
3054 1 ityp ,key ,itab ,itabm1,ixtg ,
3055 3 nixtg ,iad ,ifi ,ithgrp(1,igs),ithbuf ,
3056 4 nvars ,vars ,numeltg,varsg ,nvarsg ,
3057 5 ivarsg,nseltg,nvall,imerge,ithvar ,itherm_fe,
3058 6 flagabf,nvarabf,rfi,lsubmodel,map_tables%ITRIAM, numeltg)
3059 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
3060 ENDDO
3061 ENDIF
3062 !-------------------------------------------
3063 ! /TH/SLIPRING
3064 !-------------------------------------------
3065 CALL hm_option_count(trim(thfile)//'/SLIPRING' ,hm_nthslipring)
3066 IF (hm_nthslipring > 0) THEN
3067 CALL hm_option_start(trim(thfile)//'/SLIPRING')
3068 hm_nthgrp = hm_nthgrp + hm_nthslipring
3069 DO i = 1, hm_nthslipring
3070 titr = ''
3071 CALL hm_option_read_key(lsubmodel, option_id = id, option_titr = titr)
3072 igs = igs+1
3073 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
3074 ithgrp(1,igs) = id
3075 nvall = 6
3076 ityp = 118
3077 key = 'SLIPRING '
3078 CALL hm_read_thgrki(
3079 1 ityp ,key ,ptr_nopt_slipring,
3080 2 iad ,ifi ,ithgrp(1,igs) ,ithbuf ,
3081 3 nvarslip ,varslip ,nslipring ,varslipg,nvarslipg ,
3082 4 ivarslipg,nsne ,nvall ,ithvar ,flagabf ,nvarabf ,
3083 5 nom_opt ,igs ,lsubmodel, lithbufmx)
3084 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
3085 ENDDO
3086 ENDIF
3087 !-------------------------------------------
3088 ! /TH/RETRACTOR
3089 !-------------------------------------------
3090 CALL hm_option_count(trim(thfile)//'/RETRACTOR' ,hm_nthretractor)
3091 IF (hm_nthretractor > 0) THEN
3092 CALL hm_option_start(trim(thfile)//'/RETRACTOR')
3093 hm_nthgrp = hm_nthgrp + hm_nthretractor
3094 DO i = 1, hm_nthretractor
3095 titr = ''
3096 CALL hm_option_read_key(lsubmodel, option_id = id, option_titr = titr)
3097 igs = igs+1
3098 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
3099 ithgrp(1,igs) = id
3100 nvall = 3
3101 ityp = 119
3102 key = 'RETRACTOR '
3103 CALL hm_read_thgrki(
3104 1 ityp ,key ,ptr_nopt_retractor,
3105 2 iad ,ifi ,ithgrp(1,igs) ,ithbuf ,
3106 3 nvarret ,varret ,nretractor ,varretg,nvarretg ,
3107 4 ivarretg ,nsne ,nvall ,ithvar ,flagabf ,nvarabf ,
3108 5 nom_opt ,igs ,lsubmodel, lithbufmx)
3109 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
3110 ENDDO
3111 ENDIF
3112 !-------------------------------------------
3113 ! /TH/SENSOR
3114 !-------------------------------------------
3115 CALL hm_option_count(trim(thfile)//'/SENSOR' ,hm_nthsens )
3116 IF (hm_nthsens > 0) THEN
3117 CALL hm_option_start(trim(thfile)//'/SENSOR')
3118 hm_nthgrp = hm_nthgrp + hm_nthsens
3119 DO i = 1, hm_nthsens
3120 titr = ''
3121 CALL hm_option_read_key(lsubmodel, option_id=id, option_titr=titr)
3122 igs = igs+1
3123 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
3124 ithgrp(1,igs) = id
3125 nvall = 1
3126 ityp = 120
3127 key = 'SENSOR '
3128 CALL hm_read_thgrsens(sensors,
3129 . ityp ,key ,igs ,lithbufmx ,ithbuf ,
3130 . iad ,ifi ,ithgrp(1,igs),ithvar,nvall ,
3131 . nvarsens ,nvarsensg,varsens ,varsensg ,ivarsensg ,
3132 . nsne ,nvarabf ,lsubmodel )
3133c
3134 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
3135 ENDDO
3136 ENDIF
3137 !-------------------------------------------
3138 ! /TH/CHECKSUM ( activated automaticaly if /CHECKSUM is used )
3139 !-------------------------------------------
3140 hm_nthchecksum=checksum%checksum_count
3141 IF (hm_nthchecksum > 0) THEN
3142 hm_nthgrp = hm_nthgrp + 1
3143 igs = igs+1
3144 checksum_title = ''
3145 deck_checksum = ''
3146
3147 titr = 'CHECKSUM '
3148 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
3149
3150
3151 ithgrp(1,igs) = 100000
3152 nvall = 1
3153 ityp = 121
3154 key = 'CHECKSUM '
3155 CALL hm_read_thchecksum(
3156 1 121 ,key ,iad ,ifi ,nsne ,
3157 2 ithgrp(1:,igs) ,ithbuf ,nvarchecksum ,varchecksum ,lithbufmx,
3158 3 checksum ,nsubdom ,nithgr ,ltitr ,ipri )
3159c
3160 nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
3161 ENDIF
3162C
3163 DO n=1,npart
3164 IF(ipart(8,n)>0) nsmat=nsmat+1
3165 ENDDO
3166C=======================================================================
3167C Tag groups with shares in th part and subset not def
3168C=======================================================================
3169 DO ng=1,ngroup
3170 nel=iparg(2,ng)
3171 nft=iparg(3,ng)
3172 ity=iparg(5,ng)
3173 IF (ity == 1) THEN
3174 ipartel_p => iparts ! 3D solids
3175 ELSEIF (ity == 2) THEN
3176 ipartel_p => ipartq ! Quads
3177 ELSEIF (ity == 3) THEN
3178 ipartel_p => ipartc ! Shells
3179 ELSEIF (ity == 4) THEN
3180 ipartel_p => ipartt ! Truss
3181 ELSEIF (ity == 5) THEN
3182 ipartel_p => ipartp ! Beams
3183 ELSEIF (ity == 6) THEN
3184 ipartel_p => ipartr ! Springs
3185 ELSEIF (ity == 7 .OR. ity == 8) THEN
3186 ipartel_p => ipartg ! SH3N or Tria
3187 ELSEIF (ity == 51) THEN
3188 ipartel_p => ipartsp ! Sph
3189 ELSEIF (ity == 100) THEN
3190 ipartel_p => ipartx ! Xelems
3191 ELSEIF (ity == 101) THEN
3192 ipartel_p => ipartig3d ! Ig3d
3193 ENDIF
3194 DO j=nft+1,nft+nel
3195 IF(iparg(80,ng) == 1) EXIT
3196 numpart = ipartel_p(j)
3197 IF(iparth(nvparth,numpart)>0) THEN
3198 nvar =iparth(nvparth,numpart)
3199 iad =iparth(nvparth+1,numpart)
3200 DO n=iad,iad+nvar-1
3201 k=ithbuf(n)
3202 IF(k>7.AND.k/=28) THEN
3203 iparg(80,ng) = 1
3204 EXIT
3205 ENDIF
3206 ENDDO
3207 ENDIF
3208 ENDDO
3209 ENDDO
3210C=======================================================================
3211C print group of PART and SUBSET
3212C=======================================================================
3213 CALL thprin(
3214 1 'PART ',ipart ,lipart1 ,4 ,
3215 2 nvparth ,ithbuf ,varpa ,npart + nthpart )
3216 CALL thprin_sub(
3217 1 'SUBSET ',subset ,ithflag ,1 ,
3218 2 nvsubth ,ithbuf ,varpa ,nsubs )
3219C
3220C TH GROUPS PART AND SUBSET ARE NO LONGER USED
3221C
3222 nthgrp2=igs
3223C------------------------------------------------
3224C Search for shares or subsets in several th
3225C------------------------------------------------
3226 DO i=1,tagp
3227 DO j=1, tagp
3228 IF((i < j).AND.(pathid(i) == pathid(j)))THEN
3229 CALL ancmsg(msgid=553,
3230 . msgtype=msgerror,
3231 . anmode=aninfo_blind_1,
3232 . c1='PART',
3233 . i1=pathid(i))
3234 ENDIF
3235 ENDDO
3236 ENDDO
3237 DO i=1,tags
3238 DO j=1, tags
3239 IF((i < j).AND.(suthid(i) == suthid(j)))THEN
3240 CALL ancmsg(msgid=553,
3241 . msgtype=msgerror,
3242 . anmode=aninfo_blind_1,
3243 . c1='SUBSET',
3244 . i1=suthid(i))
3245 ENDIF
3246 ENDDO
3247 ENDDO
3248C-------------------------------------
3249C Recherche des ID doubles
3250C-------------------------------------
3251 CALL udouble(ithgrp,nithgr,nthgrp2,mess,0,bid)
3252C
3253C-------------------------------------
3254C Write th.nms1 file
3255C-------------------------------------
3256
3257 IF(dump_thnms1_file == 1) THEN
3258 CALL write_thnms1(
3259 1 nvarn1 ,nvarn1a ,nvarn2 ,nvarnpinch ,nvars1 ,
3260 2 nvars2 ,nvars3 ,nvars4 ,nvars5 ,nvars6 ,
3261 3 nvars7 ,nvars8 ,nvars9 ,nvars10 ,nvarsnloc ,
3262 4 nvarp ,nvarr ,nvart ,nvarns ,nvarsph ,
3263 5 nvarin ,nvarrw ,nvarrb ,nvarse ,
3264 5 nvarac ,nvarjo ,nvarmvent ,nvarpa ,nvarfx ,
3265 6 nvargau ,nvarfr ,nvarslip ,nvarret ,nvarclus ,
3266 7 nvarflow ,nvarsurf ,nvarc ,nvarsens ,
3267 8 varn1_title ,varn1a_title ,varn2_title ,
3268 9 varnpinch_title,varp_title ,varr_title ,vart_title ,
3269 a vars1_title ,vars2_title ,vars3_title ,vars4_title ,vars5_title ,
3270 b vars6_title ,vars7_title ,vars8_title ,vars9_title ,varsnloc_title ,
3271 c varc_title ,vars10_title ,
3272 d varns_title ,varsph_title ,varin_title ,
3273 e varrw_title ,varrb_title ,varse_title ,varac_title ,
3274 f varjo_title ,varmvent_title,varpa_title ,varfx_title ,vargau_title ,
3275 g varfr_title ,varslip_title ,varret_title ,varclus_title ,varflow_title ,
3276 h varsurf_title ,varsens_title ,
3277 i varn1 ,varn1a ,varn2 ,varnpinch ,
3278 j varp ,varr ,vart ,vars1 ,vars2 ,
3279 k vars3 ,vars4 ,vars5 ,vars6 ,vars7 ,
3280 l vars8 ,vars9 ,vars10 ,varsnloc ,
3281 m varc ,
3282 n varns ,varsph ,varin ,
3283 o varrw ,varrb ,varse ,varac ,
3284 p varjo ,varmvent ,varpa ,varfx ,vargau ,
3285 q varfr ,varslip ,varret ,varclus ,varflow ,
3286 r varsurf ,varsens )
3287 ENDIF
3288C-------------------------------------
3289
3290
3291
3292 ! Deallocate Variables
3293 DEALLOCATE(vars6)
3294 DEALLOCATE(vars7)
3295 DEALLOCATE(varsg2)
3296 DEALLOCATE(varsg3)
3297 DEALLOCATE(vars)
3298 DEALLOCATE(varc)
3299 DEALLOCATE(varsg)
3300 DEALLOCATE(ivarsg)
3301
3302 DEALLOCATE(varn1_title,varn1a_title,varn2_title,varnpinch_title,
3303 . varp_title,varr_title,vart_title,
3304 . vars1_title,vars2_title,vars3_title,vars4_title,
3305 . vars5_title,vars6_title,vars7_title,vars8_title,vars9_title,
3306 . varc_title,vars10_title,
3307 . varns_title,varsph_title,
3308 . varin_title,varrw_title,varrb_title,
3309 . varmv_title,varse_title,varac_title,
3310 . varjo_title,varmvent_title,
3311 . varpa_title,varfx_title,
3312 . vargau_title,varfr_title,
3313 . varslip_title,varret_title,
3314 . varclus_title,varflow_title,
3315 . varsurf_title,varsens_title)
3316
3317 DEALLOCATE(varn)
3318
3319
3320 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
hm_option_count
subroutine hm_option_count(entity_type, hm_option_number)
Definition hm_option_count.F:105
hm_option_start
subroutine hm_option_start(entity_type)
Definition hm_option_start.F:253
hm_read_thgrki
subroutine hm_read_thgrki(ityp, key, inopt1, iad, ifi, ithgrp, ithbuf, nv, vare, num, varg, nvg, ivarg, nsne, nv0, ithvar, flagabf, nvarabf, nom_opt, igs, lsubmodel, lithbufmx)
Definition hm_read_thgrki.F:48
hm_read_thgrki_rbody
subroutine hm_read_thgrki_rbody(ityp, key, inopt1, iad, ifi, ithgrp, ithbuf, nv, vare, num, varg, nvg, ivarg, nsne, nv0, ithvar, flagabf, nvarabf, nom_opt, igs, npby, lsubmodel)
Definition hm_read_thgrki_rbody.F:48
hm_read_thgrne
subroutine hm_read_thgrne(ityp, key, itab, itabm1, ix, nix, iad, ifi, ithgrp, ithbuf, nv, vare, numel, varg, nvg, ivarg, nsne, nv0, imerge, ithvar, itherm_fe, flagabf, nvarabf, rfi, lsubmodel, map, mapsize)
Definition hm_read_thgrne.F:56
hm_read_thgrns
subroutine hm_read_thgrns(ityp, key, itab, itabm1, kxx, ixx, iad, ifi, ithgrp, ithbuf, nv, vare, varg, nvg, ivarg, nsne, ivns2r, nv0, id, titr, ithvar, flagabf, nvarabf, lsubmodel)
Definition hm_read_thgrns.F:48
hm_read_thgrpa
subroutine hm_read_thgrpa(ipasu, npasu, ithbuf, iad, ifi, varpa, nvarpa, varg, nvarg, numthpart, ivarpag, pathid, tagp, iparth, nparth, nvparth, lsubmodel)
Definition hm_read_thgrpa.F:44
hm_read_thgrpa_sub
subroutine hm_read_thgrpa_sub(iad, ifi, ithbuf, nvarpa, varpa, varg, nvarg, ivarpag, pathid, suthid, tags, subset, ithflag, lsubmodel)
Definition hm_read_thgrpa.F:196
hm_read_thgrsens
subroutine hm_read_thgrsens(sensors, ityp, key, igs, lithbufmx, ithbuf, iad, ifi, ithgrp, ithvar, nvall, nvare, nvarg, vare, varg, ivarg, nsne, nvarabf, lsubmodel)
Definition hm_read_thgrsens.F:47
hm_read_thgrsurf
subroutine hm_read_thgrsurf(ityp, key, iad, ifi, ithgrp, ithbuf, nv, vare, num, varg, nvg, ivarg, nsne, nv0, ithvar, flagabf, nvarabf, igrsurf, igs, lsubmodel)
Definition hm_read_thgrsurf.F:48
hm_thgrki_vent
subroutine hm_thgrki_vent(ityp, key, inopt1, iad, ifi, ithgrp, ithbuf, nv, vare, num, varg, nvg, ivarg, nsne, nv0, ithvar, flagabf, nvarabf, nom_opt, igs, t_monvol, nvarmvent, lsubmodel)
Definition hm_thgrki_vent.F:58
hm_option_read_mod::hm_option_read_key
Definition hm_option_read_mod.F:414
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
id
initmumps id
Definition multiplerhs_example.m:4
groupdef_mod
Definition groupdef_mod.F:668
hm_option_read_mod
Definition hm_option_read_mod.F:412
mapping_option_mod
Definition dichotomy_mod.F:36
message_mod
Definition message_mod.F:1257
monvol_struct_mod
Definition monvol_struct_mod.F:66
names_and_titles_mod
Definition names_and_titles_mod.F:1002
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1008
names_and_titles_mod::ncharline
integer, parameter ncharline
Definition names_and_titles_mod.F:1006
submodel_mod
Definition submodel_mod.F:717
submodel_mod::nsubmod
integer nsubmod
Definition submodel_mod.F:725
nvar
integer function nvar(text)
Definition nvar.F:32
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
fretitl
subroutine fretitl(titr, iasc, l)
Definition freform.F:615
udouble
subroutine udouble(list, ilist, nlist, mess, ir, rlist)
Definition sysfus.F:573
thprin_sub
subroutine thprin_sub(key, subset, ithflag, nid, nvps, ithbuf, vare, num)
Definition thprin.F:113
thprin
subroutine thprin(key, ipasu, npasu, nid, nvps, ithbuf, vare, num)
Definition thprin.F:37