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
95C-----------------------------------------------
96C I m p l i c i t T y p e s
97C-----------------------------------------------
98#include "implicit_f.inc"
99C-----------------------------------------------
100C C o m m o n B l o c k s
101C-----------------------------------------------
102#include "scr23_c.inc"
103#include "scr17_c.inc"
104#include "com10_c.inc"
105#include "com01_c.inc"
106#include "com04_c.inc"
107#include "param_c.inc"
108#include "sphcom.inc"
109#include "fxbcom.inc"
110C-----------------------------------------------
111C D u m m y A r g u m e n t s
112C-----------------------------------------------
113 INTEGER,INTENT(IN) :: NSUBDOM
114 INTEGER,INTENT(IN) :: IPRI
115 INTEGER,INTENT(IN) :: LITHBUFMX
116 INTEGER,INTENT(IN) :: ITHERM_FE
117 INTEGER ITAB(NUMNOD),ITABM1(*), IFLAG,
118 . IXS(NIXS,NUMELS),IXQ(NIXQ,NUMELQ),IXC(NIXC,NUMELC),IXT(NIXT,NUMELT),
119 . IXP(NIXP,NUMELP),IXR(NIXR,NUMELR),IXTG(NIXTG,NUMELTG),KXX(NIXX,*),
120 . IXX(*),IPART(LIPART1,*),ITHGRP(NITHGR,*),ITHBUF(LITHBUFMX),
121 . IMERGE(*),
122 . IFI,NTHWA,KXSP(NISP,*),IXRI(4,*),
123 . IFRAME(LISKN,NUMFRAM+1), NTHGRP2,PATHID(*),SUTHID(*),
124 . FXBIPM(NBIPM,*),ISKWN(LISKN,*),
125 . NPARTH,IPARTH(NPARTH,*),NVPARTH,ITYP,
126 . NVSUBTH,ITHVAR(*),FLAGABF,NVARABF,
127 . RFI,ISPHIO(NISPHIO,*),ITHFLAG,IPARG(NPARG,*)
128 INTEGER NOM_OPT(LNOPT1,*),PTR_NOPT_FXBY,PTR_NOPT_INTER,
129 . PTR_NOPT_RWALL,PTR_NOPT_SECT,PTR_NOPT_JOINT,
130 . PTR_NOPT_MONV,PTR_NOPT_ACC,PTR_NOPT_SKW,PTR_NOPT_GAU,
131 . PTR_NOPT_CLUS,PTR_NOPT_SPHIO,PTR_NOPT_SURF,PTR_NOPT_RETRACTOR,
132 . PTR_NOPT_SLIPRING
133 INTEGER NPBY(*)
134 INTEGER, TARGET :: IPARTS(NUMELS),IPARTQ(NUMELQ),IPARTC(NUMELC),
135 . IPARTT(NUMELT),IPARTP(NUMELP),IPARTR(NUMELR),IPARTG(NUMELTG),
136 . IPARTX(NUMELX),IPARTSP(NUMSPH),IPARTIG3D(NUMELIG3D)
137 INTEGER,INTENT(INOUT) :: IPARI(NPARI,NINTER)
138 TYPE(MONVOL_STRUCT_), DIMENSION(NVOLU), INTENT(IN) :: T_MONVOL
139 TYPE (SUBMODEL_DATA) ,DIMENSION(NSUBMOD) ,INTENT(IN) :: LSUBMODEL
140 TYPE(MAPPING_STRUCT_),INTENT(IN) :: MAP_TABLES
141 TYPE (INTERFACES_) ,INTENT(INOUT):: INTERFACES
142 TYPE (SENSORS_) ,INTENT(IN) :: SENSORS
143C-----------------------------------------------
144 TYPE (SUBSET_) , DIMENSION(NSUBS) :: SUBSET
145 TYPE (SURF_) , DIMENSION(NSURF) :: IGRSURF
146 INTEGER, INTENT(IN):: DUMP_THNMS1_FILE
147 TYPE(checksum_option_) , INTENT(IN) :: checksum
148C-----------------------------------------------
149C L o c a l V a r i a b l e s
150C-----------------------------------------------
151 INTEGER J,JJ,II, I,ISU,ID,NOSYS,J10(10),NTOT,KK,IER,
152 . OK,IGS,IGRS,NSU,K,L,CONT,IAD0,IADV,NTRI,NL,
153 . IFITMP,IADFIN,NUMEL,NV,NVAR,M,N,IAD2,IAD,NSNE,NNE,ISK,
154 . NVALL,TAGP,TAGS,IFIX_TMP,IBID,NIN
155 CHARACTER(LEN=NCHARTITLE) :: TITR,TITR1
156 CHARACTER MESS*40
157 CHARACTER(LEN=NCHARLINE) :: KEY,COPT
158 INTEGER NVARN,NVARN1,NVARN1A,NVARN2,NVARNPINCH,NVARS,NVARC,NVART,NVARP,NVARR,NVARUR
159 INTEGER NVARS1,NVARS2,NVARS3,NVARS4,NVARS5,NVARS6,NVARS7,NVARS8,NVARS9,NVARS10,NVARSNLOC
160 INTEGER NVARNS,NVARSPH
161 INTEGER NVARIN,NVARRW,NVARRB,NVARAC,NVARSE,NVARJO,NVARFX,NVARFXM
162 INTEGER NVARAB,NVARMV4,NVARMV,NVARPA,NVARMVENT,NVARSENS,NVARSENSG
163 INTEGER NVARCHECKSUM,NVARCHECKSUMG
164 INTEGER NVARF1,NVARFR,NVARGAU,NVARCLUS,NVARCLUSG,NVARFLOW,NVARSURF
165 INTEGER NVARRIV,NVARRIVG,NVARSLIP,NVARRET,NSELRT,INOUT
166 INTEGER HM_NTHGRP,NTHACCEL,NTHINTER,NTHRWALL,NTHSECTIO,NTHCLUS
167 INTEGER NTHBEAM,NTHTRUS,NTHBRIC,NTHNODE,NTHSHEL,NTHSH3N,NTHSPRING,NTHRBODY
168 INTEGER NTHMONVOL,HM_NTHPART ,NUMTHPART, HM_NTHSUBS, HM_NTHSURF
169 INTEGER HM_NTHSPHCEL, HM_NTHQUAD, HM_NTHSPHFLOW
170 INTEGER HM_NTHGAUGE, HM_NTHFXBODY, HM_NTHFRAME, HM_NTHCYLJO, HM_NTHNSTRAND
171 INTEGER HM_NTHTRIA,HM_NTHSLIPRING,HM_NTHRETRACTOR,HM_NTHSENS,HM_NTHCHECKSUM
172 parameter(nvarriv = 10,nvarrivg = 1)
173 INTEGER IVARRIV(18,NVARRIVG)
174 CHARACTER*10 VARRIV(NVARRIV),VARRIVG(NVARRIVG)
175C-----------------------------------------------
176 parameter(nvarn1 = 19,nvarn1a = 600,nvarn2 = 10 ,nvarnpinch = 9)
177 parameter(nvars1 = 196,nvars2 = 567,nvars3 = 369,nvars4 = 492,nvars5 = 22)
178 parameter(nvars6 = 97200,nvars7 = 97200,nvars8 = 516,nvars9 = 6,nvars10 = 1,nvarsnloc = 2)
179 parameter(nvars =239555 ,nvarc = 37856,nvart = 6)
180 parameter(nvarp = 337,nvarr = 66,nvarur = 12)
181 parameter(nvarns = 4,nvarsph = 41)
182 parameter(nvarin = 29,nvarrw = 6,nvarrb =15,nvarfx =4)
183 parameter(nvarfxm = 3)
184 parameter(nvarac = 6,nvarse =39,nvarjo = 6)
185 parameter(nvarab = 7,nvarmv4= 9,nvarmv = 19,nvarmvent = 150)
186 parameter(nvarpa = 32)
187 parameter(nvarf1 = 18,nvarfr = 24,nvargau = 8,nvarclus=11,nvarflow=1,nvarsurf = 6,nvarslip=6,nvarret=3)
188 CHARACTER*10, DIMENSION(:), ALLOCATABLE :: VARN
189 CHARACTER*10 :: VART(NVART)
190 CHARACTER*10 VARP(NVARP),VARR(NVARR),VARUR(NVARUR)
191 CHARACTER(LEN=10),DIMENSION(:),ALLOCATABLE::VARS
192 CHARACTER(LEN=10),DIMENSION(:),ALLOCATABLE::VARC
193 CHARACTER*10 VARS1(NVARS1),VARS2(NVARS2), VARS3(NVARS3), VARS4(NVARS4),
194 . VARS5(NVARS5),VARS8(NVARS8),VARS9(NVARS9),VARS10(NVARS10),
195 . VARN1(NVARN1),VARN1A(NVARN1A),
196 . VARN2(NVARN2),
197 . VARNPINCH(NVARNPINCH),
198 . VARSG1(202),VARSG4(1),
199 . VARSNLOC(NVARSNLOC)
200 CHARACTER(LEN=10),DIMENSION(:),ALLOCATABLE::VARS6
201 CHARACTER(LEN=10),DIMENSION(:),ALLOCATABLE::VARS7
202 CHARACTER(LEN=10),DIMENSION(:),ALLOCATABLE::VARSG2
203 CHARACTER(LEN=10),DIMENSION(:),ALLOCATABLE::VARSG3
204 CHARACTER*10 VARNS(NVARNS),VARSPH(NVARSPH)
205 CHARACTER*10 VARIN(NVARIN),VARRW(NVARRW),VARRB(NVARRB)
206 CHARACTER*10 VARAC(NVARAC),VARSE(NVARSE),VARJO(NVARJO)
207 CHARACTER*10 VARAB(NVARAB),VARMV4(NVARMV4),VARMV(NVARMV),
208 . VARMVENT(NVARMVENT)
209 CHARACTER*10 VARPA(NVARPA),VARFX(NVARFX),VARFXM(NVARFXM)
210 CHARACTER*10 VARF1(NVARF1),VARFR(NVARFR),VARGAU(NVARGAU),
211 . VARSLIP(NVARSLIP),VARRET(NVARRET)
212 CHARACTER*10 VARCLUS(NVARCLUS),VARFLOW(NVARFLOW),VARSURF(NVARSURF)
213 INTEGER NVARNG,NVARSG,NVARCG,NVARTG,NVARPG,NVARRG,NVARURG
214 INTEGER NVARNSG,NVARSPG
215 INTEGER NVARING,NVARRWG,NVARRBG,NVARACG,NVARSEG,NVARJOG
216 INTEGER NVARABG,NVARMG4,NVARMVG,NVARPAG,NVARFXG,NVARFXMG
217 INTEGER NVARF1G,NVARFRG,NVARGAUG,NVARFLOWG,NVARMVENTG,NVARSURFG
218 INTEGER NVARSLIPG,NVARRETG
219 parameter(nvarng = 8,nvarsg = 39767,nvarcg = 675,nvartg = 1)
220 parameter(nvarpg = 1,nvarrg = 1,nvarurg = 1)
221 parameter(nvarnsg = 1,nvarspg = 2)
222 parameter(nvaring = 6,nvarrwg = 3,nvarrbg = 6,nvarfxg = 1)
223 parameter(nvarfxmg =1)
224 parameter(nvaracg = 2,nvarseg = 7,nvarjog = 3)
225 parameter(nvarabg = 1,nvarmg4 = 1,nvarmvg = 62,nvarmventg=62)
226 parameter(nvarpag = 1)
227 parameter(nvarf1g = 3,nvarfrg = 4,nvargaug = 1,nvarclusg=2,nvarflowg=1,
228 . nvarsurfg = 1,nvarslipg=1,nvarretg = 1)
229 parameter(nvarsens = 1,nvarsensg = 1)
230 parameter(nvarchecksum = 1,nvarchecksumg = 1)
231 CHARACTER*10 VARNG(NVARNG)
232 CHARACTER*10 VARCG(NVARCG),VARTG(NVARTG)
233 CHARACTER*10 VARPG(NVARPG),VARRG(NVARRG),VARURG(NVARURG)
234 CHARACTER*10 VARNSG(NVARNSG),VARSPG(NVARSPG)
235 CHARACTER*10 VARING(NVARING),VARRWG(NVARRWG),VARRBG(NVARRBG)
236 CHARACTER*10 VARACG(NVARACG),VARSEG(NVARSEG),VARJOG(NVARJOG)
237 CHARACTER*10 VARABG(NVARABG),VARMVG4(NVARMG4),VARMVG(NVARMVG),
238 . VARMVENTG(NVARMVENTG),VARSENS(NVARSENS),VARSENSG(NVARSENSG),
239 . VARCHECKSUM(NVARCHECKSUM),VARCHECKSUMG(NVARCHECKSUMG)
240 CHARACTER*10 VARPAG(NVARPAG),VARFXG(NVARFXG),VARFXMG(NVARFXMG)
241 CHARACTER*10 VARF1G(NVARF1G),VARFRG(NVARFRG),VARGAUG(NVARGAUG)
242 CHARACTER*10 VARCLUSG(NVARCLUSG),VARFLOWG(NVARFLOWG)
243 CHARACTER*10 VARSURFG(NVARSURFG),VARSLIPG(NVARSLIPG),VARRETG(NVARRETG)
244 INTEGER,DIMENSION(:,:),ALLOCATABLE :: IVARSG
245 CHARACTER(LEN=10),DIMENSION(:),ALLOCATABLE:: VARSG
246C-----------------------------------------------
247 INTEGER IVARNG(18,NVARNG)
248 INTEGER IVARCG(18,NVARCG),IVARTG(18,NVARTG)
249 INTEGER IVARPG(18,NVARPG),IVARRG(18,NVARRG),IVARURG(18,NVARURG)
250 INTEGER IVARNSG(18,NVARNSG),IVARSPG(18,NVARSPG),IVNS2R(18,NVARNSG)
251 INTEGER IVARING(18,NVARING),IVARRWG(18,NVARRWG)
252 INTEGER IVARRBG(18,NVARRBG),IVARACG(18,NVARACG)
253 INTEGER IVARSEG(18,NVARSEG),IVARJOG(18,NVARJOG)
254 INTEGER IVARABG(18,NVARABG),IVARMG4(18,NVARMG4),
255 . IVARMVG(18,NVARMVG),IVARMVENTG(18,NVARMVENTG)
256 INTEGER IVARPAG(18,NVARPAG),IVARFXG(18,NVARFXG),IVARSENSG(18,NVARSENSG)
257 INTEGER IVARFXMG(18,NVARFXMG)
258 INTEGER IVARF1G(18,NVARF1G),IVARFRG(18,NVARFRG),
259 . IVARGAUG(18,NVARGAUG),IVARCLUSG(18,NVARCLUSG),
260 . IVARFLOWG(18,NVARFLOWG),IVARSURFG(18,NVARSURFG),
261 . IVARSLIPG(18,NVARSLIPG),IVARRETG(18,NVARRETG),
262 . IVARCHECKSUMG(18,NVARCHECKSUMG)
263 INTEGER NSMA
264 INTEGER IUN,INODACC,INODFRA
265 INTEGER NNOD, NSHEL, NTRUS, NBEAM, NSPRG, NCHTG, NCHINT, NCHJNT,
266 . CHNOD, CHELT, CNT1, CNT2, IVSG, IREC_0, OUT_VENT, NG, NEL,
267 . NFT, NUMPART, I15ATH, I15A, I15B, I15C, I15D, I15E, I15F,
268 . I15G, I15H, I15I, I15J, I15K, I15L, ITY
269 my_real
270 . bid
271 INTEGER IVCG
272 CHARACTER CHJJ*3,CHKK*2,CHII*2,CHJKK*4
273 CHARACTER CHI*1,CHJ*3,CHK*1,CHJS*1
274 CHARACTER*10 VARC2(574),VARC1(208),VARCG1(167),VARCG2(407), VAR ,
275 . VAR_TMP,VARC3(500),VARCG3(100),
276 . VARCGPINCH(1),VARCPINCH(7),VARCNLOC(2)
277 CHARACTER*80 STR, MES
278 LOGICAL :: IS_AVAILABLE
279 INTEGER, DIMENSION(:), POINTER :: IPARTEL_P
280 CHARACTER THFILE*4
281
282 CHARACTER(LEN=100),DIMENSION(:),ALLOCATABLE:: VARN1_TITLE,VARN1A_TITLE,VARN2_TITLE,VARNPINCH_TITLE,
283 . VARP_TITLE,VARR_TITLE,VART_TITLE,
284 . VARS1_TITLE,VARS2_TITLE,VARS3_TITLE,VARS4_TITLE,
285 . VARS5_TITLE,VARS6_TITLE,VARS7_TITLE,VARS8_TITLE,VARS9_TITLE,VARS10_TITLE,
286 . VARC_TITLE,
287 . VARSNLOC_TITLE,VARNS_TITLE,VARSPH_TITLE,
288 . VARIN_TITLE,VARRW_TITLE,VARRB_TITLE,
289 . VARMV_TITLE,VARSE_TITLE,VARAC_TITLE,
290 . VARJO_TITLE,VARMVENT_TITLE,
291 . VARPA_TITLE,VARFX_TITLE,
292 . VARGAU_TITLE,VARFR_TITLE,
293 . VARSLIP_TITLE,VARRET_TITLE,
294 . VARCLUS_TITLE,VARFLOW_TITLE,
295 . VARSURF_TITLE,VARSENS_TITLE
296
297 INTEGER :: LEN_TITLE
298 INTEGER :: LEN_CHECKSUM
299 CHARACTER(LEN=NCHARLINE):: CHECKSUM_TITLE
300 CHARACTER(LEN=64):: DECK_CHECKSUM
301C-----------------------------------------------
302C E x t e r n a l F u n c t i o n s
303C-----------------------------------------------
304 INTEGER USR2SYS,ULIST2S,LISTCNT
305C-----------------------------------------------
306C D a t a
307C-----------------------------------------------
308 DATA mess/'TH GROUP DEFINITION '/
309 DATA iun/1/
310
311C-----------------------------------------------
312C TH TITLES
313C DON T MISS TO ADD TH TITLES FOR NEW TH OPTIONS
314C-----------------------------------------------
315
316 nvarn = nvarn1 + nvarn1a + nvarn2 + nvarnpinch
317 ALLOCATE(varn(nvarn))
318
319
320 ALLOCATE(varn1_title(nvarn1),varn1a_title(nvarn1a),varn2_title(nvarn2),varnpinch_title(nvarnpinch),
321 . varp_title(nvarp),varr_title(nvarr),vart_title(nvart),
322 . vars1_title(nvars1),vars2_title(nvars2),vars3_title(nvars3),vars4_title(nvars4),
323 . vars5_title(nvars5),vars6_title(nvars6),vars7_title(nvars7),vars8_title(nvars8),vars9_title(nvars9),
324 . varc_title(nvarc),vars10_title(nvars10),
325 . varsnloc_title(nvarsnloc),varns_title(nvarns),varsph_title(nvarsph),
326 . varin_title(nvarin),varrw_title(nvarrw),varrb_title(nvarrb),
327 . varmv_title(nvarmv),varse_title(nvarse),varac_title(nvarac),
328 . varjo_title(nvarjo),varmvent_title(nvarmvent),
329 . varpa_title(nvarpa),varfx_title(nvarfx),
330 . vargau_title(nvargau),varfr_title(nvarfr),
331 . varslip_title(nvarslip),varret_title(nvarret),
332 . varclus_title(nvarclus),varflow_title(nvarflow),
333 . varsurf_title(nvarsurf),varsens_title(1))
334
335
336 CALL th_titles(
337 1 nvarn1 ,nvarn1a ,nvarn2 ,nvarnpinch ,nvars1 ,
338 1 nvars2 ,nvars3 ,nvars4 ,nvars5 ,nvars6 ,
339 1 nvars7 ,nvars8 ,nvars9 ,nvars10 ,nvarsnloc ,
340 1 nvarp ,nvarr ,nvart ,nvarns ,nvarsph ,
341 2 nvarin ,nvarrw ,nvarrb ,nvarmv ,nvarse ,
342 3 nvarac ,nvarjo ,nvarmvent ,nvarpa ,nvarfx ,
343 4 nvargau ,nvarfr ,nvarslip ,nvarret ,nvarclus ,
344 5 nvarflow ,nvarsurf ,nvarc ,nvarsens ,
345 a varn1_title ,varn1a_title ,varn2_title ,
346 a varnpinch_title,varp_title ,varr_title ,vart_title ,
347 6 vars1_title ,vars2_title ,vars3_title ,vars4_title ,vars5_title ,
348 6 vars6_title ,vars7_title ,vars8_title ,vars9_title ,varsnloc_title ,
349 6 varc_title ,vars10_title ,
350 7 varns_title ,varsph_title ,varin_title ,
351 8 varrw_title ,varrb_title ,varmv_title ,varse_title ,varac_title ,
352 9 varjo_title ,varmvent_title,varpa_title ,varfx_title ,vargau_title ,
353 a varfr_title ,varslip_title ,varret_title ,varclus_title ,varflow_title ,
354 b varsurf_title ,varsens_title )
355
356C-----------------------------------------------
357C TH VARIABLES
358C-----------------------------------------------
359 ! Allocate temporary Arrays
360 ! Character Arrays - ue ALLOCATEkm
361 ALLOCATE(vars6(97200))
362 vars6(1:97200) = ' '
363 ALLOCATE(vars7(97200))
364 vars7(1:97200) = ' '
365 ALLOCATE(varsg2(19782))
366 varsg2(1:19782) = ' '
367 ALLOCATE(varsg3(19782))
368 varsg3(1:19782) = ' '
369 ALLOCATE(vars(nvars))
370 vars(1:nvars) = ' '
371 ALLOCATE(varc(nvarc))
372 varc(1:nvarc) = ' '
373 ALLOCATE(varsg(nvarsg))
374 varsg(1:nvarsg) = ' '
375
376 ! INTEGER Arrays - use MY_ALLOC
377 CALL my_alloc(ivarsg,18,nvarsg)
378C nodes
379 DATA varn1/'DX ','DY ','DZ ',
380 . 'VX ','VY ','VZ ',
381 . 'AX ','AY ','AZ ',
382 . 'VRX ','VRY ','VRZ ',
383 . 'ARX ','ARY ','ARZ ',
384 . 'X ','Y ','Z ',
385 . 'TEMP '/
386 DATA varn2/'REACX ','REACY ','REACZ ',
387 . 'REACXX ','REACYY ','REACZZ ',
388 . 'DRX ','DRY ','DRZ ',
389 . 'PEXT '/
390 DATA varnpinch/
391 . 'APINCHX ','APINCHY ','apinchz ',
392 . 'vpinchx ','vpinchy ','vpinchz ',
393 . 'dpinchx ','dpinchy ','dpinchz '/
394C solids
395 DATA VARS1/
396 .'off ','sx ','sy ','sz ','sxy ',
397 .'syz ','sxz ','IE ','DENS ','BULK ',
398 .'VOL ','PLAS ','TEMP ','PLSR ','DAM1 ',
399 .'DAM2 ','DAM3 ','DAM4 ','DAMA ','SA1 ',
400 .'SA2 ','SA3 ','CR ','CAP ','K0 ',
401 .'RK ','TD ','EFIB ','ISTA ','VPLA ',
402 .'BFRAC ','WPLA ','SFIB ','DAM5 ','LSX ',
403 .'LSY ','LSZ ','LSXY ','LSYZ ','LSXZ ',
404 .'SX1 ','SX2 ','SX3 ','SX4 ','SX5 ','SX6 ',
405 .'SX7 ','SX8 ','SY1 ','SY2 ','SY3 ','SY4 ',
406 .'SY5 ','SY6 ','SY7 ','SY8 ','SZ1 ','SZ2 ',
407 .'SZ3 ','SZ4 ','SZ5 ','SZ6 ','SZ7 ','SZ8 ',
408 .'SXY1 ','SXY2 ','SXY3 ','SXY4 ','SXY5 ','SXY6 ',
409 .'SXY7 ','SXY8 ','SYZ1 ','SYZ2 ','SYZ3 ','SYZ4 ',
410 .'SYZ5 ','SYZ6 ','SYZ7 ','SYZ8 ','SXZ1 ','SXZ2 ',
411 .'SXZ3 ','SXZ4 ','SXZ5 ','SXZ6 ','SXZ7 ','SXZ8 ',
412 .'LSX1 ','lsx2 ','lsx3 ','lsx4 ','lsx5 ','lsx6 ',
413 .'lsx7 ','lsx8 ','lsy1 ','lsy2 ','LSY3 ','LSY4 ',
414 .'LSY5 ','LSY6 ','LSY7 ','LSY8 ','LSZ1 ','LSZ2 ',
415 .'LSZ3 ','LSZ4 ','LSZ5 ','LSZ6 ','LSZ7 ','LSZ8 ',
416 .'LSXY1 ','LSXY2 ','LSXY3 ','LSXY4 ','LSXY5 ','LSXY6 ',
417 .'LSXY7 ','LSXY8 ','LSYZ1 ','LSYZ2 ','LSYZ3 ','lsyz4 ',
418 .'lsyz5 ','lsyz6 ','lsyz7 ','lsyz8 ','lsxz1 ','lsxz2 ',
419 .'lsxz3 ','lsxz4 ','lsxz5 ','lsxz6 ','lsxz7 ','lsxz8 ',
420 .'usr1 ','usr2 ','usr3 ','usr4 ','usr5 ','usr6 ',
421 .'usr7 ','usr8 ','usr9 ','usr10 ','usr11 ','usr12 ',
422 .'usr13 ','usr14 ','usr15 ','usr16 ','usr17 ','usr18 ',
423 .'usr19 ','usr20 ','usr21 ','usr22 ','usr23 ','usr24 ',
424 .'usr25 ','usr26 ','usr27 ','usr28 ','usr29 ','usr30 ',
425 .'usr31 ','usr32 ','usr33 ','usr34 ','usr35 ','usr36 ',
426 .'usr37 ','usr38 ','usr39 ','usr40 ','usr41 ','usr42 ',
427 .'usr43 ','usr44 ','usr45 ','usr46 ','usr47 ','usr48 ',
428 .'usr49 ','usr50 ','usr51 ','usr52 ','usr53 ','usr54 ',
429 .'usr55 ','usr56 ','usr57 ','usr58 ','usr59 ','usr60 '/
430C
431 DATA VARS2/
432 .'sx111 ','sy111 ','sz111 ','sxy111 ','syz111 ','sxz111 ',
433 .'ep111 ','sx211 ','sy211 ','sz211 ','sxy211 ','syz211 ',
434 .'sxz211 ','ep211 ','sx311 ','sy311 ','sz311 ','sxy311 ',
435 .'syz311 ','sxz311 ','ep311 ','sx121 ','sy121 ','sz121 ',
436 .'sxy121 ','syz121 ','sxz121 ','ep121 ','sx221 ','sy221 ',
437 .'sz221 ','sxy221 ','syz221 ','sxz221 ','ep221 ','sx321 ',
438 .'sy321 ','sz321 ','sxy321 ','syz321 ','sxz321 ','ep321 ',
439 .'sx131 ','sy131 ','sz131 ','sxy131 ','syz131 ','sxz131 ',
440 .'ep131 ','sx231 ','sy231 ','sz231 ','sxy231 ','syz331 ',
441 .'sxz231 ','ep231 ','sx331 ','sy331 ','sz331 ','sxy331 ',
442 .'syz331 ','sxz331 ','ep331 ','sx141 ','sy141 ','sz141 ',
443 .'sxy141 ','syz141 ','sxz141 ','EP141 ','SX241 ','SY241 ',
444 .'SZ241 ','SXY241 ','SYZ241 ','SXZ241 ','EP241 ','SX341 ',
445 .'SY341 ','SZ341 ','SXY341 ','SYZ341 ','SXZ341 ','EP341 ',
446 .'SX151 ','SY151 ','SZ151 ','SXY151 ','SYZ151 ','SXZ151 ',
447 .'EP151 ','SX251 ','SY251 ','SZ251 ','SXY251 ','SYZ251 ',
448 .'SXZ251 ','EP251 ','SX351 ','SY351 ','SZ351 ','SXY351 ',
449 .'SYZ351 ','SXZ351 ','EP351 ','SX161 ','SY161 ','SZ161 ',
450 .'SXY161 ','SYZ161 ','SXZ161 ','EP161 ','SX261 ','SY261 ',
451 .'SZ261 ','SXY261 ','SYZ261 ','SXZ261 ','EP261 ','SX361 ',
452 .'SY361 ','SZ361 ','SXY361 ','SYZ361 ','SXZ361 ','EP361 ',
453 .'SX171 ','SY171 ','SZ171 ','SXY171 ','SYZ171 ','SXZ171 ',
454 .'EP171 ','SX271 ','SY271 ','SZ271 ','SXY271 ','SYZ271 ',
455 .'SXZ271 ','EP271 ','SX371 ','SY371 ','SZ371 ','SXY371 ',
456 .'SYZ371 ','SXZ371 ','EP371 ','SX181 ','SY181 ','SZ181 ',
457 .'SXY181 ','SYZ181 ','SXZ181 ','EP181 ','SX281 ','SY281 ',
458 .'SZ281 ','SXY281 ','SYZ281 ','SXZ281 ','EP281 ','SX381 ',
459 .'SY381 ','SZ381 ','SXY381 ','SYZ381 ','SXZ381 ','EP381 ',
460 .'SX191 ','SY191 ','SZ191 ','SXY191 ','SYZ191 ','SXZ191 ',
461 .'EP191 ','SX291 ','SY291 ','SZ291 ','SXY291 ','SYZ291 ',
462 .'SXZ291 ','EP291 ','SX391 ','SY391 ','SZ391 ','SXY391 ',
463 .'SYZ391 ','SXZ391 ','EP391 ','SX112 ','SY112 ','SZ112 ',
464 .'SXY112 ','SYZ112 ','SXZ112 ','ep112 ','sx212 ','sy212 ',
465 .'sz212 ','sxy212 ','syz212 ','SXZ212 ','EP212 ','SX312 ',
466 .'SY312 ','SZ112 ','SXY112 ','SYZ312 ','SXZ112 ','EP312 ',
467 .'SX122 ','SY122 ','SZ122 ','SXY122 ','SYZ122 ','SXZ122 ',
468 .'EP122 ','SX222 ','SY222 ','SZ222 ','SXY222 ','SYZ222 ',
469 .'SXZ222 ','EP222 ','SX322 ','SY322 ','SZ322 ','SXY322 ',
470 .'SYZ322 ','SXZ322 ','EP322 ','SX132 ','SY132 ','SZ132 ',
471 .'SXY132 ','SYZ132 ','SXZ132 ','EP132 ','SX232 ','SY232 ',
472 .'SZ232 ','SXY232 ','SYZ332 ','SXZ232 ','EP232 ','SX332 ',
473 .'SY332 ','SZ332 ','SXY332 ','SYZ332 ','SXZ332 ','EP332 ',
474 .'SX142 ','SY142 ','SZ142 ','SXY142 ','SYZ142 ','SXZ142 ',
475 .'EP142 ','SX242 ','SY242 ','SZ242 ','SXY242 ','SYZ242 ',
476 .'SXZ242 ','EP242 ','SX342 ','SY342 ','SZ342 ','SXY342 ',
477 .'SYZ342 ','SXZ342 ','EP342 ','SX152 ','SY152 ','SZ152 ',
478 .'SXY152 ','SYZ152 ','SXZ152 ','EP152 ','SX252 ','SY252 ',
479 .'SZ252 ','SXY252 ','SYZ252 ','SXZ252 ','EP252 ','SX352 ',
480 .'SY352 ','SZ352 ','SXY352 ','SYZ352 ','SXZ352 ','EP352 ',
481 .'SX162 ','SY162 ','SZ162 ','SXY162 ','SYZ162 ','SXZ162 ',
482 .'EP162 ','SX262 ','SY262 ','SZ262 ','SXY262 ','SYZ262 ',
483 .'SXZ262 ','EP262 ','SX362 ','SY362 ','SZ362 ','SXY362 ',
484 .'SYZ362 ','SXZ362 ','EP362 ','SX172 ','SY172 ','SZ172 ',
485 .'SXY172 ','SYZ172 ','SXZ172 ','EP172 ','SX272 ','SY272 ',
486 .'SZ272 ','SXY272 ','SYZ272 ','SXZ272 ','EP272 ','SX372 ',
487 .'SY372 ','SZ372 ','SXY372 ','SYZ372 ','SXZ372 ','EP372 ',
488 .'sx182 ','sy182 ','sz182 ','sxy182 ','syz182 ','sxz182 ',
489 .'ep182 ','sx282 ','sy282 ','sz282 ','sxy282 ','syz282 ',
490 .'sxz282 ','ep282 ','sx382 ','sy382 ','sz382 ','SXY382 ',
491 .'SYZ382 ','SXZ382 ','EP382 ','SX192 ','SY192 ','SZ192 ',
492 .'SXY192 ','SYZ192 ','SXZ192 ','EP192 ','SX292 ','SY292 ',
493 .'SZ292 ','SXY292 ','SYZ292 ','SXZ292 ','EP292 ','SX392 ',
494 .'SY392 ','SZ392 ','SXY392 ','SYZ392 ','SXZ392 ','EP392 ',
495 .'SX113 ','SY113 ','SZ113 ','SXY113 ','SYZ113 ','SXZ113 ',
496 .'EP113 ','SX213 ','SY213 ','SZ213 ','SXY213 ','SYZ213 ',
497 .'SXZ213 ','EP213 ','SX313 ','SY313 ','SZ113 ','SXY113 ',
498 .'SYZ313 ','SXZ113 ','EP313 ','SX123 ','SY123 ','SZ123 ',
499 .'SXY123 ','SYZ123 ','SXZ123 ','EP123 ','SX223 ','SY223 ',
500 .'SZ223 ','SXY223 ','SYZ223 ','SXZ223 ','EP223 ','SX323 ',
501 .'SY323 ','SZ323 ','SXY323 ','SYZ323 ','SXZ323 ','EP323 ',
502 .'SX133 ','SY133 ','SZ133 ','SXY133 ','SYZ133 ','SXZ133 ',
503 .'EP133 ','SX233 ','SY233 ','SZ233 ','SXY233 ','SYZ333 ',
504 .'SXZ233 ','EP233 ','SX333 ','SY333 ','SZ333 ','SXY333 ',
505 .'SYZ333 ','SXZ333 ','EP333 ','SX143 ','SY143 ','SZ143 ',
506 .'SXY143 ','SYZ143 ','SXZ143 ','EP143 ','SX243 ','SY243 ',
507 .'SZ243 ','SXY243 ','SYZ243 ','SXZ243 ','EP243 ','SX343 ',
508 .'SY343 ','SZ343 ','SXY343 ','SYZ343 ','SXZ343 ','EP343 ',
509 .'SX153 ','SY153 ','SZ153 ','SXY153 ','SYZ153 ','SXZ153 ',
510 .'EP153 ','SX253 ','SY253 ','sz253 ','sxy253 ','syz253 ',
511 .'sxz253 ','ep253 ','sx353 ','sy353 ','sz353 ','sxy353 ',
512 .'syz353 ','sxz353 ','ep353 ','sx163 ','sy163 ','sz163 ',
513 .'sxy163 ','syz163 ','sxz163 ','ep163 ','sx263 ','sy263 ',
514 .'sz263 ','sxy263 ','syz263 ','sxz263 ','ep263 ','sx363 ',
515 .'sy363 ','sz363 ','sxy363 ','syz363 ','sxz363 ','ep363 ',
516 .'sx173 ','sy173 ','sz173 ','sxy173 ','syz173 ','sxz173 ',
517 .'ep173 ','sx273 ','sy273 ','sz273 ','sxy273 ','syz273 ',
518 .'sxz273 ','ep273 ','sx373 ','sy373 ','sz373 ','sxy373 ',
519 .'syz373 ','sxz373 ','ep373 ','sx183 ','sy183 ','sz183 ',
520 .'sxy183 ','syz183 ','sxz183 ','ep183 ','sx283 ','sy283 ',
521 .'sz283 ','sxy283 ','syz283 ','sxz283 ','ep283 ','sx383 ',
522 .'sy383 ','sz383 ','sxy383 ','SYZ383 ','SXZ383 ','EP383 ',
523 .'SX193 ','SY193 ','SZ193 ','SXY193 ','SYZ193 ','SXZ193 ',
524 .'EP193 ','SX293 ','SY293 ','SZ293 ','SXY293 ','SYZ293 ',
525 .'SXZ293 ','EP293 ','SX393 ','SY393 ','SZ393 ','SXY393 ',
526 .'SYZ393 ','SXZ393 ','EP393 '/
527 DATA vars3/
528 .'SX1U1 ','SY1U1 ','SZ1U1 ','SXY1U1 ','SYZ1U1 ','SXZ1U1 ',
529 .'EP1U1 ','SX2U1 ','SY2U1 ','SZ2U1 ','SXY2U1 ','SYZ2U1 ',
530 .'SXZ2U1 ','EP2U1 ','SX3U1 ','SY3U1 ','SZ3U1 ','SXY3U1 ',
531 .'SYZ3U1 ','SXZ3U1 ','EP3U1 ','SX1U2 ','SY1U2 ','SZ1U2 ',
532 .'SXY1U2 ','SYZ1U2 ','SXZ1U2 ','EP1U2 ','SX2U2 ','SY2U2 ',
533 .'SZ2U2 ','SXY2U2 ','SYZ2U2 ','SXZ2U2 ','EP2U2 ','SX3U2 ',
534 .'SY3U2 ','SZ3U2 ','SXY3U2 ','SYZ3U2 ','SXZ3U2 ','EP3U2 ',
535 .'SX1U3 ','SY1U3 ','SZ1U3 ','SXY1U3 ','SYZ1U3 ','SXZ1U3 ',
536 .'EP1U3 ','SX2U3 ','SY2U3 ','SZ2U3 ','sxy2u3 ','syz3u3 ',
537 .'sxz2u3 ','ep2u3 ','sx3u3 ','sy3u3 ','sz3u3 ','sxy3u3 ',
538 .'syz3u3 ','sxz3u3 ','ep3u3 ','sx1d1 ','sy1d1 ','sz1d1 ',
539 .'sxy1d1 ','syz1d1 ','sxz1d1 ','ep1d1 ','sx2d1 ','sy2u1 ',
540 .'sz2d1 ','sxy2d1 ','syz2d1 ','sxz2d1 ','ep2d1 ','sx3d1 ',
541 .'sy3d1 ','sz3d1 ','sxy3d1 ','syz3d1 ','sxz3d1 ','ep3d1 ',
542 .'sx1d2 ','sy1d2 ','sz1d2 ','sxy1d2 ','syz1d2 ','sxz1d2 ',
543 .'ep1d2 ','sx2d2 ','sy2d2 ','sz2d2 ','sxy2d2 ','syz2d2 ',
544 .'sxz2d2 ','ep2d2 ','sx3d2 ','sy3d2 ','sz3d2 ','sxy3d2 ',
545 .'syz3d2 ','sxz3d2 ','ep3d2 ','sx1d3 ','sy1d3 ','sz1d3 ',
546 .'sxy1d3 ','syz1d3 ','sxz1d3 ','ep1d3 ','sx2d3 ','sy2d3 ',
547 .'sz2d3 ','sxy2d3 ','syz3d3 ','sxz2d3 ','ep2d3 ','sx3d3 ',
548 .'sy3d3 ','sz3d3 ','sxy3d3 ','syz3d3 ','sxz3d3 ','ep3d3 ',
549 .'usr1_111','usr2_111','usr3_111','usr4_111','usr5_111','usr6_111',
550 .'usr7_111','usr8_111','usr9_111','usr1_211','usr2_211','usr3_211',
551 .'usr4_211','usr5_211','usr6_211','usr7_211','usr8_211','usr9_211',
552 .'usr1_311','usr3_311','usr3_311','usr4_311','usr5_311','usr6_311',
553 .'usr7_311','usr8_311','usr9_311','usr1_121','usr2_121','USR3_121',
554 .'USR4_121','USR5_121','USR6_121','USR7_121','USR8_121','USR9_121',
555 .'USR1_221','USR2_221','USR3_221','USR4_221','USR5_221','USR6_221',
556 .'USR7_221','USR8_221','USR9_221','USR1_321','USR3_321','USR3_321',
557 .'USR4_321','USR5_321','USR6_321','USR7_321','USR8_321','USR9_321',
558 .'USR1_131','USR2_131','USR3_131','USR4_131','USR5_131','USR6_131',
559 .'USR7_131','USR8_131','USR9_131','USR1_231','USR2_231','USR3_231',
560 .'USR4_231','USR5_231','USR6_231','USR7_231','USR8_231','USR9_231',
561 .'USR1_331','USR3_331','USR3_331','USR4_331','USR5_331','USR6_331',
562 .'USR7_331','USR8_331','USR9_331','USR1_141','USR2_141','USR3_141',
563 .'USR4_141','USR5_141','USR6_141','USR7_141','USR8_141','USR9_141',
564 .'USR1_241','USR2_241','USR3_241','USR4_241','USR5_241','USR6_241',
565 .'USR7_241','USR8_241','USR9_241','USR1_341','USR3_341','USR3_341',
566 .'USR4_341','USR5_341','USR6_341','USR7_341','USR8_341','USR9_341',
567 .'USR1_151','USR2_151','USR3_151','USR4_151','USR5_151','USR6_151',
568 .'USR7_151','USR8_151','USR9_151','USR1_251','USR2_251','USR3_251',
569 .'USR4_251','USR5_251','USR6_251','USR7_251','USR8_251','USR9_251',
570 .'USR1_351','USR3_351','USR3_351','USR4_351','USR5_351','USR6_351',
571 .'USR7_351','USR8_351','USR9_351','USR1_161','USR2_161','USR3_161',
572 .'USR4_161','USR5_161','USR6_161','USR7_161','USR8_161','USR9_161',
573 .'USR1_261','USR2_261','USR3_261','USR4_261','USR5_261','USR6_261',
574 .'USR7_261','USR8_261','USR9_261','USR1_361','USR3_361','USR3_361',
575 .'USR4_361','USR5_361','USR6_361','USR7_361','USR8_361','USR9_361',
576 .'USR1_171','USR2_171','USR3_171','USR4_171','USR5_171','USR6_171',
577 .'USR7_171','USR8_171','USR9_171','USR1_271','USR2_271','USR3_271',
578 .'USR4_271','USR5_271','USR6_271','USR7_271','USR8_271','USR9_271',
579 .'USR1_371','USR3_371','USR3_371','USR4_371','USR5_371','USR6_371',
580 .'USR7_371','USR8_371','USR9_371','USR1_181','USR2_181','USR3_181',
581 .'USR4_181','USR5_171','USR6_171','USR7_181','USR8_181','USR9_181',
582 .'USR1_281','USR2_281','USR3_281','USR4_281','USR5_271','usr6_271',
583 .'usr7_281','usr8_281','usr9_281','usr1_381','usr3_381','usr3_381',
584 .'usr4_381','usr5_371','usr6_371','usr7_381','usr8_381','usr9_381',
585 .'usr1_191','usr2_191','usr3_191','usr4_191','usr5_191','usr6_191',
586 .'usr7_191','usr8_191','usr9_191','USR1_291','USR2_291','USR3_291',
587 .'USR4_291','USR5_291','USR6_291','USR7_291','USR8_291','USR9_291',
588 .'USR1_391','USR3_391','USR3_391','USR4_391','USR5_391','USR6_391',
589 .'USR7_391','USR8_391','USR9_391'/
590 DATA vars4/
591 .'USR1_112','USR2_112','USR3_112','USR4_112','USR5_112','USR6_112',
592 .'USR7_112','USR8_112','USR9_112',
593 .'USR1_212','usr2_212','usr3_212','usr4_212','usr5_212','usr6_212',
594 .'usr7_212','usr8_212','usr9_212','usr1_312','usr3_312','usr3_312',
595 .'usr4_312','usr5_312','usr6_312','usr7_312','usr8_312','USR9_312',
596 .'USR1_122','USR2_122','USR3_122','USR4_122','USR5_122','USR6_122',
597 .'USR7_122','USR8_122','USR9_122','USR1_222','USR2_222','USR3_222',
598 .'USR4_222','USR5_222','USR6_222','USR7_222','USR8_222','USR9_222',
599 .'USR1_322','USR3_322','USR3_322','USR4_322','USR5_322','USR6_322',
600 .'USR7_322','USR8_322','USR9_322','USR1_132','USR2_132','USR3_132',
601 .'USR4_132','USR5_132','USR6_132','USR7_132','USR8_132','USR9_132',
602 .'USR1_232','USR2_232','USR3_232','USR4_232','USR5_232','USR6_232',
603 .'USR7_232','USR8_232','USR9_232','USR1_332','USR3_332','usr3_332',
604 .'usr4_332','usr5_332','usr6_332','usr7_332','usr8_332','usr9_332',
605 .'usr1_142','usr2_142','usr3_142','usr4_142','usr5_142','usr6_142',
606 .'usr7_142','usr8_142','usr9_142','usr1_242','usr2_242','usr3_242',
607 .'usr4_242','usr5_242','usr6_242','usr7_242','usr8_242','usr9_242',
608 .'usr1_342','usr3_342','usr3_342','usr4_342','usr5_342','usr6_342',
609 .'usr7_342','usr8_342','usr9_342','usr1_152','usr2_152','usr3_152',
610 .'usr4_152','usr5_152','usr6_152','usr7_152','usr8_152','usr9_152',
611 .'usr1_252','usr2_252','usr3_252','usr4_252','usr5_252','usr6_252',
612 .'usr7_252','usr8_252','usr9_252','usr1_352','usr3_352','usr3_352',
613 .'usr4_352','usr5_352','usr6_352','usr7_352','usr8_352','usr9_352',
614 .'usr1_162','usr2_162','usr3_162','usr4_162','usr5_162','usr6_162',
615 .'usr7_162','usr8_162','usr9_162','usr1_262','usr2_262','usr3_262',
616 .'usr4_262','usr5_262','USR6_262','USR7_262','USR8_262','USR9_262',
617 .'USR1_362','USR3_362','USR3_362','USR4_362','USR5_362','USR6_362',
618 .'USR7_362','USR8_362','USR9_362','USR1_172','USR2_172','USR3_172',
619 .'USR4_172','USR5_172','USR6_172','USR7_172','USR8_172','USR9_172',
620 .'USR1_272','USR2_272','USR3_272','USR4_272','USR5_272','USR6_272',
621 .'USR7_272','USR8_272','USR9_272','USR1_372','USR3_372','USR3_372',
622 .'USR4_372','USR5_372','USR6_372','USR7_372','USR8_372','USR9_372',
623 .'USR1_182','USR2_182','USR3_182','USR4_182','USR5_182','USR6_182',
624 .'USR7_182','USR8_182','USR9_182','USR1_282','USR2_282','USR3_282',
625 .'usr4_282','usr5_282','usr6_282','usr7_282','usr8_282','usr9_282',
626 .'usr1_382','usr3_382','usr3_382','usr4_382','usr5_382','usr6_382',
627 .'usr7_382','usr8_382','usr9_382','usr1_192','usr2_192','usr3_192',
628 .'usr4_192','usr5_192','usr6_192','usr7_192','usr8_192','usr9_192',
629 .'usr1_292','usr2_292','usr3_292','usr4_292','usr5_292','usr6_292',
630 .'usr7_292','usr8_292','usr9_292','usr1_392','usr3_392','usr3_392',
631 .'usr4_392','usr5_392','usr6_392','usr7_392','usr8_392','usr9_392',
632 .'usr1_113','usr2_113','usr3_113','usr4_113','usr5_113','usr6_113',
633 .'usr7_113','usr8_113','usr9_113','usr1_213','usr2_213','usr3_213',
634 .'usr4_213','usr5_213','usr6_213','usr7_213','usr8_213','usr9_213',
635 .'usr1_313','usr3_313','usr3_313','usr4_313','usr5_313','usr6_313',
636 .'usr7_313','usr8_313','usr9_313','usr1_123','usr2_123','usr3_123',
637 .'usr4_123','usr5_123','usr6_123','usr7_123','usr8_123','usr9_123',
638 .'usr1_223','usr2_223','usr3_223','usr4_223','usr5_223','usr6_223',
639 .'usr7_223','usr8_223','usr9_223','usr1_323','usr3_323','usr3_323',
640 .'usr4_323','usr5_323','usr6_323','usr7_323','usr8_323','usr9_323',
641 .'usr1_133','usr2_133','usr3_133','usr4_133','usr5_133','usr6_133',
642 .'usr7_133','usr8_133','usr9_133','usr1_233','usr2_233','usr3_233',
643 .'usr4_233','usr5_233','usr6_233','usr7_233','usr8_233','usr9_233',
644 .'usr1_333','usr3_333','usr3_333','usr4_333','usr5_333','usr6_333',
645 .'usr7_333','usr8_333','usr9_333','usr1_143','usr2_143','usr3_143',
646 .'usr4_143','usr5_143','usr6_143','usr7_143','usr8_143','usr9_143',
647 .'usr1_243','usr2_243','usr3_243','usr4_243','usr5_243','USR6_243',
648 .'USR7_243','USR8_243','USR9_243','USR1_343','USR3_343','USR3_343',
649 .'USR4_343','USR5_343','USR6_343','USR7_343','USR8_343','USR9_343',
650 .'USR1_153','USR2_153','USR3_153','USR4_153','USR5_153','USR6_153',
651 .'USR7_153','USR8_153','USR9_153','USR1_253','USR2_253','USR3_253',
652 .'USR4_253','USR5_253','USR6_253','USR7_253','USR8_253','USR9_253',
653 .'USR1_353','USR3_353','USR3_353','USR4_353','USR5_353','USR6_353',
654 .'USR7_353','USR8_353','USR9_353','USR1_163','USR2_163','USR3_163',
655 .'USR4_163','USR5_163','USR6_163','USR7_163','USR8_163','USR9_163',
656 .'USR1_263','USR2_263','USR3_263','USR4_263','USR5_263','USR6_263',
657 .'USR7_263','USR8_263','USR9_263','USR1_363','USR3_363','USR3_363',
658 .'USR4_363','USR5_363','USR6_363','USR7_363','USR8_363','USR9_363',
659 .'USR1_173','USR2_173','USR3_173','USR4_173','USR5_173','USR6_173',
660 .'USR7_173','USR8_173','USR9_173','USR1_273','USR2_273','USR3_273',
661 .'USR4_273','USR5_273','USR6_273','USR7_273','USR8_273','USR9_273',
662 .'USR1_373','USR3_373','USR3_373','USR4_373','USR5_373','USR6_373',
663 .'USR7_373','USR8_373','USR9_373','USR1_183','USR2_183','USR3_183',
664 .'USR4_183','USR5_183','USR6_183','USR7_183','USR8_183','USR9_183',
665 .'USR1_283','USR2_283','USR3_283','USR4_283','USR5_283','USR6_283',
666 .'USR7_283','USR8_283','USR9_283','USR1_383','USR3_383','USR3_383',
667 .'USR4_383','USR5_383','USR6_383','USR7_383','USR8_383','USR9_383',
668 .'USR1_193','USR2_193','USR3_193','USR4_193','USR5_193','USR6_193',
669 .'USR7_193','USR8_193','USR9_193','USR1_293','USR2_293','USR3_293',
670 .'USR4_293','USR5_293','USR6_293','USR7_293','USR8_293','USR9_293',
671 .'USR1_393','USR3_393','USR3_393','USR4_393','USR5_393','USR6_393',
672 .'USR7_393','USR8_393','USR9_393',
673 .'EPSXX ','EPSYY ','EPSZZ ','EPSXY ','EPSZX ','EPSYZ '/
674 DATA vars5/
675 .'SX_MAT1 ','SY_MAT1 ','SZ_MAT1 ','SXY_MAT1','SYZ_MAT1','SXZ_MAT1',
676 .'IE_MAT1 ','DEN_MAT1','VOL_MAT1','PLA_MAT1','TEM_MAT1',
677 .'SX_MAT2 ','SY_MAT2 ','SZ_MAT2 ','SXY_MAT2','SYZ_MAT2','SXZ_MAT2',
678 .'IE_MAT2 ','DEN_MAT2','VOL_MAT2','PLA_MAT2','TEM_MAT2'/
679C Definition of VARS6 : strain for /PROP/TYPE22
680C=======================================================================
681 DO j=1,200
682 IF (j <= 9) THEN
683 WRITE(chj,'(I1.1)')j
684 ELSEIF (j <= 99) THEN
685 WRITE(chj,'(I2.2)')j
686 ELSE
687 WRITE(chj,'(I3.3)')j
688 ENDIF
689 DO i=1,9
690 WRITE(chi,'(I1.1)')i
691 DO k=1,9
692 WRITE(chk,'(I1.1)')k
693 DO l=1,6
694 IF (l==1) var_tmp(1:5)='E11_'//chi
695 IF (l==2) var_tmp(1:5)='E12_'//chi
696 IF (l==3) var_tmp(1:5)='E13_'//chi
697 IF (l==4) var_tmp(1:5)='E22_'//chi
698 IF (l==5) var_tmp(1:5)='E23_'//chi
699 IF (l==6) var_tmp(1:5)='E33_'//chi
700 var_tmp=var_tmp(1:5)//chk//'_'//chj
701 vars6((j-1)*9*9*6+((i-1)*9+k-1)*6+l)=var_tmp
702 ENDDO
703 ENDDO
704 ENDDO
705 ENDDO
706C Definition of VARS7 : stress for /PROP/TYPE22
707 DO j=1,200
708 IF (j <= 9) THEN
709 WRITE(chj,'(I1.1)')j
710 ELSEIF (j <= 99) THEN
711 WRITE(chj,'(I2.2)')j
712 ELSE
713 WRITE(chj,'(I3.3)')j
714 ENDIF
715 DO i=1,9
716 WRITE(chi,'(I1.1)')i
717 DO k=1,9
718 WRITE(chk,'(I1.1)')k
719 DO l=1,6
720 IF (l==1) var_tmp(1:5)='S11_'//chi
721 IF (l==2) var_tmp(1:5)='S12_'//chi
722 IF (l==3) var_tmp(1:5)='S13_'//chi
723 IF (l==4) var_tmp(1:5)='S22_'//chi
724 IF (l==5) var_tmp(1:5)='S23_'//chi
725 IF (l==6) var_tmp(1:5)='S33_'//chi
726 var_tmp=var_tmp(1:5)//chk//'_'//chj
727 vars7((j-1)*9*9*6+((i-1)*9+k-1)*6+l)=var_tmp
728 ENDDO
729 ENDDO
730 ENDDO
731 ENDDO
732C Definition of VARS8 : strain matrix output
733
734 vars8(1) = 'LEPSX '
735 vars8(2) = 'LEPSY '
736 vars8(3) = 'LEPSZ '
737 vars8(4) = 'LEPSXY '
738 vars8(5) = 'LEPSYZ '
739 vars8(6) = 'LEPSZX '
740 vars8(7) = 'EPSX1 '
741 vars8(8) = 'EPSX2 '
742 vars8(9) = 'EPSX3 '
743 vars8(10) = 'EPSX4 '
744 vars8(11) = 'EPSY1 '
745 vars8(12) = 'EPSY2 '
746 vars8(13) = 'EPSY3 '
747 vars8(14) = 'EPSY4 '
748 vars8(15) = 'EPSZ1 '
749 vars8(16) = 'EPSZ2 '
750 vars8(17) = 'EPSZ3 '
751 vars8(18) = 'EPSZ4 '
752 vars8(19) = 'EPSXY1 '
753 vars8(20) = 'EPSXY2 '
754 vars8(21) = 'EPSXY3 '
755 vars8(22) = 'EPSXY4 '
756 vars8(23) = 'EPSYZ1 '
757 vars8(24) = 'EPSYZ2 '
758 vars8(25) = 'EPSYZ3 '
759 vars8(26) = 'EPSYZ4 '
760 vars8(27) = 'EPSZX1 '
761 vars8(28) = 'EPSZX2 '
762 vars8(29) = 'EPSZX3 '
763 vars8(30) = 'EPSZX4 '
764
765 DO k=1,3
766 DO j=1,9
767 DO i=1,3
768 WRITE(chi,'(I1.1)')i
769 WRITE(chjs,'(I1.1)')j
770 WRITE(chk,'(I1.1)')k
771 var_tmp = 'EPSX'
772 vars8(30+(k-1)*9*6*3+(j-1)*3*6+(i-1)*6+1)=var_tmp(1:4)//chi//chjs//chk
773 var_tmp = 'EPSY'
774 vars8(30+(k-1)*9*6*3+(j-1)*3*6+(i-1)*6+2)=var_tmp(1:4)//chi//chjs//chk
775 var_tmp = 'EPSZ'
776 vars8(30+(k-1)*9*6*3+(j-1)*3*6+(i-1)*6+3)=var_tmp(1:4)//chi//chjs//chk
777 var_tmp = 'EPSXY'
778 vars8(30+(k-1)*9*6*3+(j-1)*3*6+(i-1)*6+4)=var_tmp(1:5)//chi//chjs//chk
779 var_tmp = 'EPSYZ'
780 vars8(30+(k-1)*9*6*3+(j-1)*3*6+(i-1)*6+5)=var_tmp(1:5)//chi//chjs//chk
781 var_tmp = 'EPSZX'
782 vars8(30+(k-1)*9*6*3+(j-1)*3*6+(i-1)*6+6)=var_tmp(1:5)//chi//chjs//chk
783
784 ENDDO
785 ENDDO
786 ENDDO
787
788 DATA vars9/
789 . 'VX ','VY ','VZ ','SSP ','MACH ',
790 . 'VK '/
791
792 DATA varsnloc/
793 . 'NL_PLAS ','NL_PLSR '/
794
795 DATA vars10/
796 . 'VSTRAIN ' /
797
798
799C shells
800c VARC = VARC1 + VARC2 + VARC3
801 DATA varc1/
802 . 'F1 ','F2 ','F12 ','Q1 ','Q2 ',
803 . 'M1 ','M2 ','M12 ','IEM ','IEB ',
804 . 'OFF ','THIC ','EMIN ','EMAX ','E1 ',
805 . 'E2 ','E12 ','SH1 ','SH2 ','K1 ',
806 .'K2 ','K12 ','EPSD ',
807 .'USR1 ','USR2 ','USR3 ','USR4 ','USR5 ','USR6 ',
808 .'USR7 ','USR8 ','USR9 ','USR10 ','USR11 ','USR12 ',
809 .'USR13 ','USR14 ','USR15 ','USR16 ','USR17 ','USR18 ',
810 .'USR19 ','USR20 ','USR21 ','USR22 ','USR23 ','USR24 ',
811 .'USR25 ','USR26 ','USR27 ','USR28 ','USR29 ','USR30 ',
812 .'USR31 ','USR32 ','USR33 ','USR34 ','USR35 ','USR36 ',
813 .'USR37 ','USR38 ','USR39 ','USR40 ','USR41 ','USR42 ',
814 .'USR43 ','USR44 ','USR45 ','USR46 ','USR47 ','USR48 ',
815 .'USR49 ','USR50 ','USR51 ','USR52 ','USR53 ','USR54 ',
816 .'USR55 ','USR56 ','USR57 ','USR58 ','USR59 ','USR60 ',
817 .'USR1_1 ','USR1_2 ','USR1_3 ','USR1_4 ','USR1_5 ','USR2_1 ',
818 .'USR2_2 ','USR2_3 ','USR2_4 ','USR2_5 ','USR3_1 ','USR3_2 ',
819 .'USR3_3 ','USR3_4 ','USR3_5 ','USR4_1 ','USR4_2 ','USR4_3 ',
820 .'USR4_4 ','USR4_5 ','USR5_1 ','USR5_2 ','USR5_3 ','USR5_4 ',
821 .'USR5_5 ','USR6_1 ','USR6_2 ','USR6_3 ','USR6_4 ','USR6_5 ',
822 .'USR7_1 ','USR7_2 ','USR7_3 ','USR7_4 ','USR7_5 ','USR8_1 ',
823 .'USR8_2 ','USR8_3 ','USR8_4 ','USR8_5 ','USR9_1 ','USR9_2 ',
824 .'USR9_3 ','USR9_4 ','USR9_5 ','USR10_1 ','USR10_2 ','USR10_3 ',
825 .'USR10_4 ','USR10_5 ','USR11_1 ','USR11_2 ','USR11_3 ','USR11_4 ',
826 .'USR11_5 ','USR12_1 ','USR12_2 ','USR12_3 ','USR12_4 ','USR12_5 ',
827 .'USR13_1 ','USR13_2 ','USR13_3 ','USR13_4 ','USR13_5 ','USR14_1 ',
828 .'USR14_2 ','USR14_3 ','USR14_4 ','USR14_5 ','USR15_1 ','USR15_2 ',
829 .'USR15_3 ','USR15_4 ','USR15_5 ','USR16_1 ','USR16_2 ','USR16_3 ',
830 .'USR16_4 ','USR16_5 ','USR17_1 ','USR17_2 ','USR17_3 ','USR17_4 ',
831 .'USR17_5 ','USR18_1 ','USR18_2 ','USR18_3 ','USR18_4 ','USR18_5 ',
832 .'USR19_1 ','USR19_2 ','USR19_3 ','USR19_4 ','USR19_5 ','USR20_1 ',
833 .'USR20_2 ','USR20_3 ','USR20_4 ','USR20_5 ','SX_1 ','SY_1 ',
834 .'SXY_1 ','SYZ_1 ','SZX_1 ','SX_2 ','SY_2 ','SXY_2 ',
835 .'SYZ_2 ','SZX_2 ','SX_3 ','SY_3 ','SXY_3 ','SYZ_3 ',
836 .'SZX_3 ','SX_4 ','SY_4 ','SXY_4 ','SYZ_4 ','SZX_4 ',
837 .'SX_5 ','SY_5 ','SXY_5 ','SYZ_5 ','SZX_5 '/
838 DATA varc2/
839 .'SX_6 ','SY_6 ','SXY_6 ','SYZ_6 ','SZX_6 ',
840 .'SX_7 ','SY_7 ','SXY_7 ','SYZ_7 ','SZX_7 ','SX_8 ',
841 .'SY_8 ','SXY_8 ','SYZ_8 ','SZX_8 ','SX_9 ','SY_9 ',
842 .'SXY_9 ','SYZ_9 ','SZX_9 ','SX_10 ','SY_10 ','SXY_10 ',
843 .'SYZ_10 ','SZX_10 ','SX_11 ','SY_11 ','SXY_11 ','SYZ_11 ',
844 .'SZX_11 ','SX_12 ','SY_12 ','SXY_12 ','SYZ_12 ','SZX_12 ',
845 .'SX_13 ','SY_13 ','SXY_13 ','SYZ_13 ','SZX_13 ','SX_14 ',
846 .'SY_14 ','SXY_14 ','SYZ_14 ','SZX_14 ','SX_15 ','SY_15 ',
847 .'SXY_15 ','SYZ_15 ','SZX_15 ','SX_16 ','SY_16 ','SXY_16 ',
848 .'SYZ_16 ','SZX_16 ','SX_17 ','SY_17 ','SXY_17 ','SYZ_17 ',
849 .'SZX_17 ','SX_18 ','SY_18 ','SXY_18 ','SYZ_18 ','SZX_18 ',
850 .'SX_19 ','SY_19 ','SXY_19 ','SYZ_19 ','SZX_19 ','SX_20 ',
851 .'SY_20 ','SXY_20 ','SYZ_20 ','SZX_20 ','SX_21 ','SY_21 ',
852 .'SXY_21 ','SYZ_21 ','SZX_21 ','SX_22 ','SY_22 ','SXY_22 ',
853 .'SYZ_22 ','SZX_22 ','SX_23 ','SY_23 ','SXY_23 ','SYZ_23 ',
854 .'SZX_23 ','SX_24 ','SY_24 ','SXY_24 ','SYZ_24 ','SZX_24 ',
855 .'SX_25 ','SY_25 ','SXY_25 ','SYZ_25 ','SZX_25 ','SX_26 ',
856 .'SY_26 ','SXY_26 ','SYZ_26 ','SZX_26 ','SX_27 ','SY_27 ',
857 .'SXY_27 ','SYZ_27 ','SZX_27 ','SX_28 ','SY_28 ','SXY_28 ',
858 .'SYZ_28 ','SZX_28 ','SX_29 ','SY_29 ','SXY_29 ','SYZ_29 ',
859 .'SZX_29 ','SX_30 ','SY_30 ','SXY_30 ','SYZ_30 ','SZX_30 ',
860 .'SX_31 ','SY_31 ','SXY_31 ','SYZ_31 ','SZX_31 ','SX_32 ',
861 .'SY_32 ','SXY_32 ','SYZ_32 ','SZX_32 ','SX_33 ','SY_33 ',
862 .'SXY_33 ','SYZ_33 ','SZX_33 ','SX_34 ','SY_34 ','SXY_34 ',
863 .'SYZ_34 ','SZX_34 ','SX_35 ','SY_35 ','SXY_35 ','SYZ_35 ',
864 .'SZX_35 ','SX_36 ','SY_36 ','SXY_36 ','SYZ_36 ','SZX_36 ',
865 .'SX_37 ','SY_37 ','SXY_37 ','SYZ_37 ','SZX_37 ','SX_38 ',
866 .'SY_38 ','SXY_38 ','SYZ_38 ','SZX_38 ','SX_39 ','SY_39 ',
867 .'SXY_39 ','SYZ_39 ','SZX_39 ','SX_40 ','SY_40 ','SXY_40 ',
868 .'SYZ_40 ','SZX_40 ','SX_41 ','SY_41 ','SXY_41 ','SYZ_41 ',
869 .'SZX_41 ','SX_42 ','SY_42 ','SXY_42 ','SYZ_42 ','SZX_42 ',
870 .'SX_43 ','SY_43 ','SXY_43 ','SYZ_43 ','SZX_43 ','SX_44 ',
871 .'SY_44 ','SXY_44 ','SYZ_44 ','SZX_44 ','SX_45 ','SY_45 ',
872 .'SXY_45 ','SYZ_45 ','SZX_45 ','SX_46 ','SY_46 ','SXY_46 ',
873 .'SYZ_46 ','SZX_46 ','SX_47 ','SY_47 ','SXY_47 ','SYZ_47 ',
874 .'SZX_47 ','SX_48 ','SY_48 ','SXY_48 ','SYZ_48 ','SZX_48 ',
875 .'SX_49 ','SY_49 ','SXY_49 ','SYZ_49 ','SZX_49 ','SX_50 ',
876 .'SY_50 ','SXY_50 ','SYZ_50 ','SZX_50 ','SX_51 ','SY_51 ',
877 .'SXY_51 ','SYZ_51 ','SZX_51 ','SX_52 ','SY_52 ','SXY_52 ',
878 .'SYZ_52 ','SZX_52 ','SX_53 ','SY_53 ','SXY_53 ','SYZ_53 ',
879 .'SZX_53 ','SX_54 ','SY_54 ','SXY_54 ','SYZ_54 ','SZX_54 ',
880 .'SX_55 ','SY_55 ','SXY_55 ','SYZ_55 ','SZX_55 ','SX_56 ',
881 .'SY_56 ','SXY_56 ','SYZ_56 ','SZX_56 ','SX_57 ','SY_57 ',
882 .'SXY_57 ','SYZ_57 ','SZX_57 ','SX_58 ','SY_58 ','SXY_58 ',
883 .'SYZ_58 ','SZX_58 ','SX_59 ','SY_59 ','SXY_59 ','SYZ_59 ',
884 .'SZX_59 ','SX_60 ','SY_60 ','SXY_60 ','SYZ_60 ','SZX_60 ',
885 .'SX_61 ','SY_61 ','SXY_61 ','SYZ_61 ','SZX_61 ','SX_62 ',
886 .'SY_62 ','SXY_62 ','SYZ_62 ','SZX_62 ','SX_63 ','SY_63 ',
887 .'SXY_63 ','SYZ_63 ','SZX_63 ','SX_64 ','SY_64 ','SXY_64 ',
888 .'SYZ_64 ','SZX_64 ','SX_65 ','SY_65 ','SXY_65 ','SYZ_65 ',
889 .'SZX_65 ','SX_66 ','SY_66 ','SXY_66 ','SYZ_66 ','SZX_66 ',
890 .'SX_67 ','SY_67 ','SXY_67 ','SYZ_67 ','SZX_67 ','SX_68 ',
891 .'SY_68 ','SXY_68 ','SYZ_68 ','SZX_68 ','SX_69 ','SY_69 ',
892 .'SXY_69 ','SYZ_69 ','SZX_69 ','SX_70 ','SY_70 ','SXY_70 ',
893 .'SYZ_70 ','SZX_70 ','SX_71 ','SY_71 ','SXY_71 ','SYZ_71 ',
894 .'SZX_71 ','SX_72 ','SY_72 ','SXY_72 ','SYZ_72 ','SZX_72 ',
895 .'SX_73 ','SY_73 ','SXY_73 ','SYZ_73 ','SZX_73 ','SX_74 ',
896 .'SY_74 ','SXY_74 ','SYZ_74 ','SZX_74 ','SX_75 ','SY_75 ',
897 .'SXY_75 ','SYZ_75 ','SZX_75 ','SX_76 ','SY_76 ','SXY_76 ',
898 .'SYZ_76 ','SZX_76 ','SX_77 ','SY_77 ','SXY_77 ','SYZ_77 ',
899 .'SZX_77 ','SX_78 ','SY_78 ','SXY_78 ','SYZ_78 ','SZX_78 ',
900 .'SX_79 ','SY_79 ','SXY_79 ','SYZ_79 ','SZX_79 ','SX_80 ',
901 .'SY_80 ','SXY_80 ','SYZ_80 ','SZX_80 ','SX_81 ','SY_81 ',
902 .'SXY_81 ','SYZ_81 ','SZX_81 ','SX_82 ','SY_82 ','SXY_82 ',
903 .'SYZ_82 ','SZX_82 ','SX_83 ','SY_83 ','SXY_83 ','SYZ_83 ',
904 .'SZX_83 ','SX_84 ','SY_84 ','SXY_84 ','SYZ_84 ','SZX_84 ',
905 .'SX_85 ','SY_85 ','SXY_85 ','SYZ_85 ','SZX_85 ','SX_86 ',
906 .'SY_86 ','SXY_86 ','SYZ_86 ','SZX_86 ','SX_87 ','SY_87 ',
907 .'SXY_87 ','SYZ_87 ','SZX_87 ','SX_88 ','SY_88 ','SXY_88 ',
908 .'SYZ_88 ','SZX_88 ','SX_89 ','SY_89 ','SXY_89 ','SYZ_89 ',
909 .'SZX_89 ','SX_90 ','SY_90 ','SXY_90 ','SYZ_90 ','SZX_90 ',
910 .'SX_91 ','SY_91 ','SXY_91 ','SYZ_91 ','SZX_91 ','SX_92 ',
911 .'SY_92 ','SXY_92 ','SYZ_92 ','SZX_92 ','SX_93 ','SY_93 ',
912 .'SXY_93 ','SYZ_93 ','SZX_93 ','SX_94 ','SY_94 ','SXY_94 ',
913 .'SYZ_94 ','SZX_94 ','SX_95 ','SY_95 ','SXY_95 ','SYZ_95 ',
914 .'SZX_95 ','SX_96 ','SY_96 ','SXY_96 ','SYZ_96 ','SZX_96 ',
915 .'SX_97 ','SY_97 ','SXY_97 ','SYZ_97 ','SZX_97 ','SX_98 ',
916 .'SY_98 ','SXY_98 ','SYZ_98 ','SZX_98 ','SX_99 ','SY_99 ',
917 .'SXY_99 ','SYZ_99 ','SZX_99 ','NFAIL ','PFAIL ','FAIL_D1 ',
918 .'FAIL_D2 ','FAIL_EN ','WPLAY01 ','WPLAY02 ','WPLAY03 ','WPLAY04 ',
919 .'WPLAY05 ','WPLAY06 ','WPLAY07 ','WPLAY08 ','WPLAY09 ','WPLAY10 ',
920 .'WPLAY11 ','WPLAY12 ','WPLAY13 ','WPLAY14 ','WPLAY15 ','WPLAY16 ',
921 .'WPLAY17 ','WPLAY18 ','WPLAY19 ','WPLAY20 ','WPLAY21 ','WPLAY22 ',
922 .'WPLAY23 ','WPLAY24 ','WPLAY25 ','WPLAY26 ','WPLAY27 ','WPLAY28 ',
923 .'WPLAY29 ','WPLAY30 ','WPLAY31 ','WPLAY32 ','WPLAY33 ','WPLAY34 ',
924 .'WPLAY35 ','WPLAY36 ','WPLAY37 ','WPLAY38 ','WPLAY39 ','WPLAY40 ',
925 .'WPLAY41 ','WPLAY42 ','WPLAY43 ','WPLAY44 ','WPLAY45 ','WPLAY46 ',
926 .'WPLAY47 ','WPLAY48 ','WPLAY49 ','WPLAY50 ','WPLAY51 ','WPLAY52 ',
927 .'WPLAY53 ','WPLAY54 ','WPLAY55 ','WPLAY56 ','WPLAY57 ','WPLAY58 ',
928 .'WPLAY59 ','WPLAY60 ','WPLAY61 ','WPLAY62 ','WPLAY63 ','WPLAY64 ',
929 .'WPLAY65 ','WPLAY66 ','WPLAY67 ','WPLAY68 ','WPLAY69 ','WPLAY70 ',
930 .'WPLAY71 ','WPLAY72 ','WPLAY73 ','WPLAY74 ','WPLAY75 ','WPLAY76 ',
931 .'WPLAY77 ','WPLAY78 ','WPLAY79 ','WPLAY80 ','WPLAY81 ','WPLAY82 ',
932 .'WPLAY83 ','WPLAY84 ','WPLAY85 ','WPLAY86 ','WPLAY87 ','WPLAY88 ',
933 .'WPLAY89 ','WPLAY90 ','WPLAY91 ','WPLAY92 ','WPLAY93 ','WPLAY94 ',
934 .'WPLAY95 ','WPLAY96 ','WPLAY97 ','WPLAY98 ','WPLAY99 '/
935 DATA varc3/
936 .'SVX_1 ','SVY_1 ','SVXY_1 ','SVYZ_1 ','SVZX_1 ','SVX_2 ',
937 .'SVY_2 ','SVXY_2 ','SVYZ_2 ','SVZX_2 ','SVX_3 ','SVY_3 ',
938 .'SVXY_3 ','SVYZ_3 ','SVZX_3 ','SVX_4 ','SVY_4 ','SVXY_4 ',
939 .'SVYZ_4 ','SVZX_4 ','SVX_5 ','SVY_5 ','SVXY_5 ','SVYZ_5 ',
940 .'SVZX_5 ','SVX_6 ','SVY_6 ','SVXY_6 ','SVYZ_6 ','SVZX_6 ',
941 .'SVX_7 ','SVY_7 ','SVXY_7 ','SVYZ_7 ','SVZX_7 ','SVX_8 ',
942 .'SVY_8 ','SVXY_8 ','SVYZ_8 ','SVZX_8 ','SVX_9 ','SVY_9 ',
943 .'SVXY_9 ','SVYZ_9 ','SVZX_9 ','SVX_10 ','SVY_10 ','SVXY_10 ',
944 .'SVYZ_10 ','SVZX_10 ','SVX_11 ','SVY_11 ','SVXY_11 ','SVYZ_11 ',
945 .'SVZX_11 ','SVX_12 ','SVY_12 ','SVXY_12 ','SVYZ_12 ','SVZX_12 ',
946 .'SVX_13 ','SVY_13 ','SVXY_13 ','SVYZ_13 ','SVZX_13 ','SVX_14 ',
947 .'SVY_14 ','SVXY_14 ','SVYZ_14 ','SVZX_14 ','SVX_15 ','SVY_15 ',
948 .'SVXY_15 ','SVYZ_15 ','SVZX_15 ','SVX_16 ','SVY_16 ','SVXY_16 ',
949 .'SVYZ_16 ','SVZX_16 ','SVX_17 ','SVY_17 ','SVXY_17 ','SVYZ_17 ',
950 .'SVZX_17 ','SVX_18 ','SVY_18 ','SVXY_18 ','SVYZ_18 ','SVZX_18 ',
951 .'SVX_19 ','SVY_19 ','SVXY_19 ','SVYZ_19 ','SVZX_19 ','SVX_20 ',
952 .'SVY_20 ','SVXY_20 ','SVYZ_20 ','SVZX_20 ','SVX_21 ','SVY_21 ',
953 .'SVXY_21 ','SVYZ_21 ','SVZX_21 ','SVX_22 ','SVY_22 ','SVXY_22 ',
954 .'SVYZ_22 ','SVZX_22 ','SVX_23 ','SVY_23 ','SVXY_23 ','SVYZ_23 ',
955 .'SVZX_23 ','SVX_24 ','SVY_24 ','SVXY_24 ','SVYZ_24 ','SVZX_24 ',
956 .'SVX_25 ','SVY_25 ','SVXY_25 ','SVYZ_25 ','SVZX_25 ','SVX_26 ',
957 .'SVY_26 ','SVXY_26 ','SVYZ_26 ','SVZX_26 ','SVX_27 ','SVY_27 ',
958 .'SVXY_27 ','SVYZ_27 ','SVZX_27 ','SVX_28 ','SVY_28 ','SVXY_28 ',
959 .'SVYZ_28 ','SVZX_28 ','SVX_29 ','SVY_29 ','SVXY_29 ','SVYZ_29 ',
960 .'SVZX_29 ','SVX_30 ','SVY_30 ','SVXY_30 ','SVYZ_30 ','SVZX_30 ',
961 .'SVX_31 ','SVY_31 ','SVXY_31 ','SVYZ_31 ','SVZX_31 ','SVX_32 ',
962 .'SVY_32 ','SVXY_32 ','SVYZ_32 ','SVZX_32 ','SVX_33 ','SVY_33 ',
963 .'SVXY_33 ','SVYZ_33 ','SVZX_33 ','SVX_34 ','SVY_34 ','SVXY_34 ',
964 .'SVYZ_34 ','SVZX_34 ','SVX_35 ','SVY_35 ','SVXY_35 ','SVYZ_35 ',
965 .'SVZX_35 ','SVX_36 ','SVY_36 ','SVXY_36 ','SVYZ_36 ','SVZX_36 ',
966 .'SVX_37 ','SVY_37 ','SVXY_37 ','SVYZ_37 ','SVZX_37 ','SVX_38 ',
967 .'SVY_38 ','SVXY_38 ','SVYZ_38 ','SVZX_38 ','SVX_39 ','SVY_39 ',
968 .'SVXY_39 ','SVYZ_39 ','SVZX_39 ','SVX_40 ','SVY_40 ','SVXY_40 ',
969 .'SVYZ_40 ','SVZX_40 ','SVX_41 ','SVY_41 ','SVXY_41 ','SVYZ_41 ',
970 .'SVZX_41 ','SVX_42 ','SVY_42 ','SVXY_42 ','SVYZ_42 ','SVZX_42 ',
971 .'SVX_43 ','SVY_43 ','SVXY_43 ','SVYZ_43 ','SVZX_43 ','SVX_44 ',
972 .'SVY_44 ','SVXY_44 ','SVYZ_44 ','SVZX_44 ','SVX_45 ','SVY_45 ',
973 .'SVXY_45 ','SVYZ_45 ','SVZX_45 ','SVX_46 ','SVY_46 ','SVXY_46 ',
974 .'SVYZ_46 ','SVZX_46 ','SVX_47 ','SVY_47 ','SVXY_47 ','SVYZ_47 ',
975 .'SVZX_47 ','SVX_48 ','SVY_48 ','SVXY_48 ','SVYZ_48 ','SVZX_48 ',
976 .'SVX_49 ','SVY_49 ','SVXY_49 ','SVYZ_49 ','SVZX_49 ','SVX_50 ',
977 .'SVY_50 ','SVXY_50 ','SVYZ_50 ','SVZX_50 ','SVX_51 ','SVY_51 ',
978 .'SVXY_51 ','SVYZ_51 ','SVZX_51 ','SVX_52 ','SVY_52 ','SVXY_52 ',
979 .'SVYZ_52 ','SVZX_52 ','SVX_53 ','SVY_53 ','SVXY_53 ','SVYZ_53 ',
980 .'SVZX_53 ','SVX_54 ','SVY_54 ','SVXY_54 ','SVYZ_54 ','SVZX_54 ',
981 .'SVX_55 ','SVY_55 ','SVXY_55 ','SVYZ_55 ','SVZX_55 ','SVX_56 ',
982 .'SVY_56 ','SVXY_56 ','SVYZ_56 ','SVZX_56 ','SVX_57 ','SVY_57 ',
983 .'SVXY_57 ','SVYZ_57 ','SVZX_57 ','SVX_58 ','SVY_58 ','SVXY_58 ',
984 .'SVYZ_58 ','SVZX_58 ','SVX_59 ','SVY_59 ','SVXY_59 ','SVYZ_59 ',
985 .'SVZX_59 ','SVX_60 ','SVY_60 ','SVXY_60 ','SVYZ_60 ','SVZX_60 ',
986 .'SVX_61 ','SVY_61 ','SVXY_61 ','SVYZ_61 ','SVZX_61 ','SVX_62 ',
987 .'SVY_62 ','SVXY_62 ','SVYZ_62 ','SVZX_62 ','SVX_63 ','SVY_63 ',
988 .'SVXY_63 ','SVYZ_63 ','SVZX_63 ','SVX_64 ','SVY_64 ','SVXY_64 ',
989 .'SVYZ_64 ','SVZX_64 ','SVX_65 ','SVY_65 ','SVXY_65 ','SVYZ_65 ',
990 .'SVZX_65 ','SVX_66 ','SVY_66 ','SVXY_66 ','SVYZ_66 ','SVZX_66 ',
991 .'SVX_67 ','SVY_67 ','SVXY_67 ','SVYZ_67 ','SVZX_67 ','SVX_68 ',
992 .'SVY_68 ','SVXY_68 ','SVYZ_68 ','SVZX_68 ','SVX_69 ','SVY_69 ',
993 .'SVXY_69 ','SVYZ_69 ','SVZX_69 ','SVX_70 ','SVY_70 ','SVXY_70 ',
994 .'SVYZ_70 ','SVZX_70 ','SVX_71 ','SVY_71 ','SVXY_71 ','SVYZ_71 ',
995 .'SVZX_71 ','SVX_72 ','SVY_72 ','SVXY_72 ','SVYZ_72 ','SVZX_72 ',
996 .'SVX_73 ','SVY_73 ','SVXY_73 ','SVYZ_73 ','SVZX_73 ','SVX_74 ',
997 .'SVY_74 ','SVXY_74 ','SVYZ_74 ','SVZX_74 ','SVX_75 ','SVY_75 ',
998 .'SVXY_75 ','SVYZ_75 ','SVZX_75 ','SVX_76 ','SVY_76 ','SVXY_76 ',
999 .'SVYZ_76 ','SVZX_76 ','SVX_77 ','SVY_77 ','SVXY_77 ','SVYZ_77 ',
1000 .'SVZX_77 ','SVX_78 ','SVY_78 ','SVXY_78 ','SVYZ_78 ','SVZX_78 ',
1001 .'SVX_79 ','SVY_79 ','SVXY_79 ','SVYZ_79 ','SVZX_79 ','SVX_80 ',
1002 .'SVY_80 ','SVXY_80 ','SVYZ_80 ','SVZX_80 ','SVX_81 ','SVY_81 ',
1003 .'SVXY_81 ','SVYZ_81 ','SVZX_81 ','SVX_82 ','SVY_82 ','SVXY_82 ',
1004 .'SVYZ_82 ','SVZX_82 ','SVX_83 ','SVY_83 ','SVXY_83 ','SVYZ_83 ',
1005 .'SVZX_83 ','SVX_84 ','SVY_84 ','SVXY_84 ','SVYZ_84 ','SVZX_84 ',
1006 .'SVX_85 ','SVY_85 ','SVXY_85 ','SVYZ_85 ','SVZX_85 ','SVX_86 ',
1007 .'SVY_86 ','SVXY_86 ','SVYZ_86 ','SVZX_86 ','SVX_87 ','SVY_87 ',
1008 .'SVXY_87 ','SVYZ_87 ','SVZX_87 ','SVX_88 ','SVY_88 ','SVXY_88 ',
1009 .'SVYZ_88 ','SVZX_88 ','SVX_89 ','SVY_89 ','SVXY_89 ','SVYZ_89 ',
1010 .'SVZX_89 ','SVX_90 ','SVY_90 ','SVXY_90 ','SVYZ_90 ','SVZX_90 ',
1011 .'SVX_91 ','SVY_91 ','SVXY_91 ','SVYZ_91 ','SVZX_91 ','SVX_92 ',
1012 .'SVY_92 ','SVXY_92 ','SVYZ_92 ','SVZX_92 ','svx_93 ','svy_93 ',
1013 .'svxy_93 ','svyz_93 ','svzx_93 ','svx_94 ','svy_94 ','svxy_94 ',
1014 .'svyz_94 ','svzx_94 ','svx_95 ','svy_95 ','svxy_95 ','svyz_95 ',
1015 .'svzx_95 ','svx_96 ','svy_96 ','svxy_96 ','svyz_96 ','svzx_96 ',
1016 .'svx_97 ','svy_97 ','svxy_97 ','svyz_97 ','svzx_97 ','svx_98 ',
1017 .'svy_98 ','svxy_98 ','svyz_98 ','svzx_98 ','svx_99 ','svy_99 ',
1018 .'svxy_99 ','svyz_99 ','svzx_99 ','svx_100 ','svy_100 ','svxy_100',
1019 .'svyz_100','svzx_100'/
1020 DATA VARCPINCH/
1021 .'epinchxz','epinchyz','epinchzz','fpinch ','mpinchxz','mpinchyz',
1022 .'thk '/
1023 DATA VARCNLOC/
1024 .'nl_plas ','nl_epsd '/
1025C trusses
1026 DATA VART/
1027 . 'off ','f ','ie ','a ','l ',
1028 . 'plas '/
1029C beams
1030 DATA VARP/
1031 . 'off ','f1 ','f2 ','f3 ','m1 ',
1032 . 'm2 ','m3 ','ie ','sx ','sxy ',
1033 . 'szx ','sx_1 ','sxy_1 ','szx_1 ','sx_2 ',
1034 . 'sxy_2 ','szx_2 ','sx_3 ','sxy_3 ','szx_3 ',
1035 . 'sx_4 ','sxy_4 ','szx_4 ','sx_5 ','sxy_5 ',
1036 . 'szx_5 ','sx_6 ','sxy_6 ','SZX_6 ','SX_7 ',
1037 . 'SXY_7 ','SZX_7 ','SX_8 ','SXY_8 ','SZX_8 ',
1038 . 'SX_9 ','SXY_9 ','SZX_9 ','SX_10 ','SXY_10 ',
1039 . 'SZX_10 ','SX_11 ','SXY_11 ','SZX_11 ','SX_12 ',
1040 . 'SXY_12 ','SZX_12 ','SX_13 ','SXY_13 ','SZX_13 ',
1041 . 'SX_14 ','SXY_14 ','SZX_14 ','SX_15 ','SXY_15 ',
1042 . 'SZX_15 ','SX_16 ','SXY_16 ','SZX_16 ','SX_17 ',
1043 . 'SXY_17 ','SZX_17 ','SX_18 ','SXY_18 ','SZX_18 ',
1044 . 'SX_19 ','SXY_19 ','SZX_19 ','SX_20 ','SXY_20 ',
1045 . 'SZX_20 ','SX_21 ','SXY_21 ','SZX_21 ','SX_22 ',
1046 . 'SXY_22 ','SZX_22 ','SX_23 ','SXY_23 ','SZX_23 ',
1047 . 'SX_24 ','SXY_24 ','SZX_24 ','SX_25 ','SXY_25 ',
1048 . 'SZX_25 ','SX_26 ','SXY_26 ','SZX_26 ','SX_27 ',
1049 . 'SXY_27 ','SZX_27 ','SX_28 ','SXY_28 ','SZX_28 ',
1050 . 'SX_29 ','SXY_29 ','SZX_29 ','SX_30 ','SXY_30 ',
1051 . 'SZX_30 ','SX_31 ','SXY_31 ','SZX_31 ','SX_32 ',
1052 . 'SXY_32 ','SZX_32 ','SX_33 ','SXY_33 ','SZX_33 ',
1053 . 'SX_34 ','SXY_34 ','SZX_34 ','SX_35 ','SXY_35 ',
1054 . 'SZX_35 ','SX_36 ','SXY_36 ','SZX_36 ','SX_37 ',
1055 . 'SXY_37 ','SZX_37 ','SX_38 ','SXY_38 ','SZX_38 ',
1056 . 'SX_39 ','SXY_39 ','SZX_39 ','SX_40 ','SXY_40 ',
1057 . 'SZX_40 ','SX_41 ','SXY_41 ','SZX_41 ','SX_42 ',
1058 . 'SXY_42 ','SZX_42 ','SX_43 ','SXY_43 ','SZX_43 ',
1059 . 'SX_44 ','SXY_44 ','SZX_44 ','SX_45 ','SXY_45 ',
1060 . 'SZX_45 ','SX_46 ','SXY_46 ','SZX_46 ','SX_47 ',
1061 . 'SXY_47 ','SZX_47 ','SX_48 ','SXY_48 ','SZX_48 ',
1062 . 'SX_49 ','sxy_49 ','szx_49 ','sx_50 ','sxy_50 ',
1063 . 'szx_50 ','sx_51 ','sxy_51 ','SZX_51 ','SX_52 ',
1064 . 'SXY_52 ','SZX_52 ','SX_53 ','SXY_53 ','SZX_53 ',
1065 . 'SX_54 ','SXY_54 ','SZX_54 ','SX_55 ','SXY_55 ',
1066 . 'SZX_55 ','SX_56 ','SXY_56 ','SZX_56 ','SX_57 ',
1067 . 'SXY_57 ','SZX_57 ','SX_58 ','SXY_58 ','SZX_58 ',
1068 . 'SX_59 ','SXY_59 ','SZX_59 ','SX_60 ','SXY_60 ',
1069 . 'SZX_60 ','SX_61 ','SXY_61 ','SZX_61 ','SX_62 ',
1070 . 'SXY_62 ','SZX_62 ','SX_63 ','SXY_63 ','SZX_63 ',
1071 . 'SX_64 ','SXY_64 ','SZX_64 ','SX_65 ','SXY_65 ',
1072 . 'SZX_65 ','SX_66 ','SXY_66 ','SZX_66 ','SX_67 ',
1073 . 'SXY_67 ','SZX_67 ','sx_68 ','sxy_68 ','szx_68 ',
1074 . 'sx_69 ','sxy_69 ','szx_69 ','sx_70 ','sxy_70 ',
1075 . 'szx_70 ','sx_71 ','sxy_71 ','szx_71 ','SX_72 ',
1076 . 'SXY_72 ','SZX_72 ','SX_73 ','SXY_73 ','SZX_73 ',
1077 . 'SX_74 ','SXY_74 ','SZX_74 ','SX_75 ','SXY_75 ',
1078 . 'SZX_75 ','SX_76 ','SXY_76 ','SZX_76 ','SX_77 ',
1079 . 'SXY_77 ','SZX_77 ','SX_78 ','SXY_78 ','SZX_78 ',
1080 . 'SX_79 ','SXY_79 ','SZX_79 ','SX_80 ','SXY_80 ',
1081 . 'SZX_80 ','SX_81 ','SXY_81 ','SZX_81 ','EPSP ',
1082 . 'EPSP_1 ','EPSP_2 ','EPSP_3 ','EPSP_4 ','EPSP_5 ',
1083 . 'EPSP_6 ','EPSP_7 ','EPSP_8 ','EPSP_9 ','EPSP_10 ',
1084 . 'EPSP_11 ','EPSP_12 ','EPSP_13 ','EPSP_14 ','EPSP_15 ',
1085 . 'EPSP_16 ','EPSP_17 ','EPSP_18 ','EPSP_19 ','EPSP_20 ',
1086 . 'EPSP_21 ','EPSP_22 ','EPSP_23 ','EPSP_24 ','EPSP_25 ',
1087 . 'EPSP_26 ','EPSP_27 ','EPSP_28 ','EPSP_29 ','EPSP_30 ',
1088 . 'EPSP_31 ','EPSP_32 ','EPSP_33 ','EPSP_34 ','EPSP_35 ',
1089 . 'EPSP_36 ','EPSP_37 ','EPSP_38 ','EPSP_39 ','EPSP_40 ',
1090 . 'EPSP_41 ','epsp_42 ','epsp_43 ','epsp_44 ','epsp_45 ',
1091 . 'epsp_46 ','epsp_47 ','epsp_48 ','epsp_49 ','epsp_50 ',
1092 . 'epsp_51 ','epsp_52 ','EPSP_53 ','EPSP_54 ','EPSP_55 ',
1093 . 'EPSP_56 ','EPSP_57 ','EPSP_58 ','EPSP_59 ','EPSP_60 ',
1094 . 'EPSP_61 ','EPSP_62 ','EPSP_63 ','EPSP_64 ','EPSP_65 ',
1095 . 'EPSP_66 ','EPSP_67 ','EPSP_68 ','EPSP_69 ','EPSP_70 ',
1096 . 'EPSP_71 ','EPSP_72 ','EPSP_73 ','EPSP_74 ','EPSP_75 ',
1097 . 'EPSP_76 ','EPSP_77 ','EPSP_78 ','EPSP_79 ','EPSP_80 ',
1098 . 'EPSP_81 ','EPSD '/
1099C springs
1100 DATA varr/
1101 . 'OFF ','FX ','FY ','FZ ','MX ',
1102 . 'MY ','MZ ','LX ','LY ','LZ ',
1103 . 'rx ','ry ','rz ','ie ','f1 ',
1104 . 'f2 ','f1g_x ','f1g_y ','f1g_z ','f2g_x ',
1105 . 'f2g_y ','f2g_z ','f1l_x ','f1l_y ','f1l_z ',
1106 . 'f2l_x ','f2l_y ','f2l_z ','m1g_x ','m1g_y ',
1107 . 'm1g_z ','m2g_x ','m2g_y ','m2g_z ','m1l_x ',
1108 . 'M1L_Y ','M1L_Z ','M2L_X ','M2L_Y ','M2L_Z ',
1109 . 'D1G_X ','D1G_Y ','d1g_z ','d2g_x ','d2g_y ',
1110 . 'd2g_z ','d1l_x ','d1l_y ','d1l_z ','d2l_x ',
1111 . 'd2l_y ','d2l_z ','r1g_x ','R1G_Y ','R1G_Z ',
1112 . 'R2G_X ','R2G_Y ','R2G_Z ','R1L_X ','R1L_Y ',
1113 . 'R1L_Z ','R2L_X ','R2L_Y ','R2L_Z ','LENGTH ',
1114 . 'FAIL '/
1115
1116 DATA varns/
1117 . 'OFF ','FX ','LX ','IE '/
1118
1119 DATA varsph/
1120 . 'OFF ','SX ','SY ','SZ ','SXY ',
1121 . 'SYZ ','SXZ ','IE ','DENS ','WFVIS ',
1122 . 'VOL ','PLAS ','TEMP ','PLSR ','DAM1 ',
1123 . 'DAM2 ','DAM3 ','DAM4 ','DAMA ','SA1 ',
1124 . 'SA2 ','SA3 ','CR ','CAP ','K0 ',
1125 . 'RK ','TD ','EFIB ','ISTA ','VPLA ',
1126 . 'BFRAC ','WPLA ','SFIB ','DAM5 ','LSX ',
1127 . 'LSY ','LSZ ','LSXY ','LSYZ ','LSXZ ',
1128 . 'DIAMETER'/
1129 DATA varur/
1130 . 'OFF ','EINT ','EINR ','XN1 ','XN2 ',
1131 . 'XN3 ','DX1 ','DX2 ','DX3 ','XM1 ',
1132 . 'XM2 ','XM3 '/
1133C interfaces
1134 DATA varin/
1135 . 'FNX ','FNY ','FNZ ','FTX ','FTY ',
1136 . 'FTZ ','SFW ','|FNX| ','|FNY| ','|FNZ| ',
1137 . '||FN|| ','|FX| ','|FY| ','|FZ| ','||F|| ',
1138 . 'PVOL ','PSURF ','PMED ','DELTAP ','VOL ',
1139 . 'SURF ','MX ','MY ','MZ ','QFRIC ',
1140 . 'CE_ELAST','CE_FRIC ','CE_DAMP ','CAREA '/
1141C rigid wall
1142 DATA varrw/
1143 . 'FNX ','FNY ','FNZ ','FTX ','FTY ',
1144 . 'FTZ '/
1145C rigid body
1146 DATA varrb/
1147 . 'FX ','FY ','FZ ','MX ','MY ',
1148 . 'MZ ','RX ','RY ','RZ ',
1149 . 'FXI ','FYI ','FZI ',
1150 . 'MXI ','MYI ','MZI '/
1151C flexible body
1152 DATA varfx/
1153 . 'IE ','KE ','EFW ','DE '/
1154C flexible body mode
1155 DATA varfxm/
1156 . 'D ','V ','A '/
1157C accel
1158 DATA varac/
1159 . 'AX ','AY ','AZ ','WX ','WY ',
1160 . 'WZ '/
1161C section
1162 DATA varse/
1163 . 'FNX ','FNY ','FNZ ','FTX ','FTY ',
1164 . 'FTZ ','M1 ','M2 ','M3 ','WORK ',
1165 . 'DFX ','DFY ','DFZ ','DF2 ','WORKR ',
1166 . 'DMX ','DMY ','DMZ ','DM2 ','KIN ',
1167 . 'KINR ','DMVX ','DMVY ','DMVZ ','DKIN ',
1168 . 'DMVRX ','DMVRY ','DMVRZ ','DKINR ','TFEXT ',
1169 . 'MX ','MY ','MZ ','F1 ','F2 ',
1170 . 'F3 ','CX ','CY ','CZ '/
1171C joint cyl
1172 DATA varjo/
1173 . 'FX ','FY ','FZ ','MX ','MY ',
1174 . 'MZ '/
1175C airbag
1176 DATA varab/
1177 . 'MASS ','VOL ','P ','RO ','IE ',
1178 . 'OE ','OU '/
1179C monitored volume
1180 DATA varmv4/
1181 . 'MASS ','VOL ','P ','A ','T ',
1182 . 'AO ','UO ','AC ','UC '/
1183 DATA varmv/
1184 . 'MASS ','VOL ','P ','A ','T ',
1185 . 'AO ','UO ','AC ','UC ','CP ',
1186 . 'CV ','GAMA ','DTBAG ','NFV ','MASS-IN ',
1187 . 'ENTHA-IN','ENER-INT','WORK ','UPCRIT '/
1188 DATA varmvent/
1189 . 'MASS ','VOL ','P ','A ','T ',
1190 . 'AO ','UO ','AC ','UC ','CP ',
1191 . 'CV ','GAMA ','DTBAG ','NFV ','MASS-IN ',
1192 . 'ENTHA-IN','ENER-INT','WORK ','UPCRIT ','UNUSED ',
1193 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1194 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1195 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','UNUSED ',
1196 . 'UNUSED ','UNUSED ','UNUSED ','UNUSED ','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 . 'AOUT1 ','BOUT1 ','UOUT1 ','MOUT1 ','HOUT1 ',
1210 . 'AOUT2 ','BOUT2 ','UOUT2 ','MOUT2 ','HOUT2 ',
1211 . 'AOUT3 ','BOUT3 ','UOUT3 ','MOUT3 ','HOUT3 ',
1212 . 'AOUT4 ','BOUT4 ','UOUT4 ','MOUT4 ','HOUT4 ',
1213 . 'AOUT5 ','BOUT5 ','UOUT5 ','MOUT5 ','HOUT5 ',
1214 . 'AOUT6 ','BOUT6 ','UOUT6 ','MOUT6 ','HOUT6 ',
1215 . 'AOUT7 ','BOUT7 ','UOUT7 ','MOUT7 ','HOUT7 ',
1216 . 'AOUT8 ','BOUT8 ','UOUT8 ','MOUT8 ','HOUT8 ',
1217 . 'AOUT9 ','BOUT9 ','UOUT9 ','MOUT9 ','HOUT9 ',
1218 . 'AOUT10 ','BOUT10 ','UOUT10 ','MOUT10 ','HOUT10 '/
1219C part
1220 DATA varpa/
1221 . 'IE ','KE ','XMOM ','YMOM ','ZMOM ',
1222 . 'MASS ','HE ','TURBKE ','XCG ','YCG ',
1223 . 'ZCG ','XXMOM ','YYMOM ','ZZMOM ','IXX ',
1224 . 'IYY ','IZZ ','IXY ','IYZ ','IZX ',
1225 . 'RIE ','KERB ','RKERB ','RKE ','ERODED ',
1226 . 'UNUSED ','UNUSED ','HEAT ','VX ','VY ',
1227 . 'VZ ','PW '/
1228C rivets
1229 DATA varriv/
1230 . 'OFF ','F1X ','F1Y ','F1Z ','F2X ',
1231 . 'F2Y ','F2Z ','FX ','FY ','FZ '/
1232C frames
1233 DATA varf1/
1234 . 'OX ','OY ','OZ ','R11 ','R12 ',
1235 . 'R13 ','R21 ','R22 ','R23 ','R31 ',
1236 . 'R32 ','R33 ','VX ','VY ','VZ ',
1237 . 'VRX ','VRY ','VRZ '/
1238 DATA varfr/
1239 . 'OX ','OY ','OZ ','R11 ','R12 ',
1240 . 'R13 ','R21 ','R22 ','R23 ','R31 ',
1241 . 'R32 ','R33 ','VX ','VY ','VZ ',
1242 . 'VRX ','VRY ','VRZ ','AX ','AY ',
1243 . 'AZ ','ARX ','ARY ','ARZ '/
1244 DATA vargau/
1245 . 'P ','IE ','DENS ','T ','PA ',
1246 . 'PAX ','PAY ','PAZ '/
1247C element cluster
1248 DATA varclus/
1249 . 'FX ','FY ','FZ ','MX ','MY ',
1250 . 'MZ ','FS ','FN ','MS ','MN ',
1251 . 'FAIL '/
1252 DATA varflow/
1253 . 'FLOW '/
1254 DATA varsurf/
1255 . 'AREA ','MASSFLOW','VELOCITY','P ','A ','MASS '/
1256C
1257 DATA varslip/
1258 . 'RINGSLIP','FN ','F1 ','F2 ','THETA ','GAMMA '/
1259C
1260 DATA varret/
1261 . 'SLIP','FN ','LOCK '/
1262c
1263 DATA varsens/
1264 . 'STATUS '/
1265c
1266 DATA varchecksum/
1267 . 'CHECKSUM '/
1268C-----------------------------------------------
1269C TH GENERIC VARIABLES
1270C-----------------------------------------------
1271 DATA varng/'DEF ','D ','V ',
1272 . 'A ','VR ','AR ','XYZ ',
1273 . 'PINCH '/
1274 DATA ivarng
1275 . / 1, 2, 3, 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1276 . 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1277 . 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1278 . 7, 8, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1279 . 10,11,12, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1280 . 13,14,15, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1281 . 16,17,18, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1282 . 630,631,632,633,634,635,636,637,638, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1283 DATA varsg1/
1284 .'DEF ','STRESS ','LOCSTRS ','SXGP ','SYGP ','SZGP ',
1285 .'SXYGP ','SYZGP ','SXZGP ','LSXGP ','LSYGP ','LSZGP ',
1286 .'LSXYGP ','LSYZGP ','LSXZGP ','USERA ','USERB ','USERC ',
1287 .'USERD ',
1288 .'SIG111 ','SIG211 ','SIG311 ','SIG121 ','SIG221 ','SIG321 ',
1289 .'SIG131 ','SIG231 ','SIG331 ','SIG141 ','SIG241 ','SIG341 ',
1290 .'SIG151 ','SIG251 ','SIG351 ','SIG161 ','SIG261 ','SIG361 ',
1291 .'SIG171 ','SIG271 ','SIG371 ','SIG181 ','SIG281 ','SIG381 ',
1292 .'SIG191 ','SIG291 ','SIG391 ',
1293 .'SIG112 ','SIG212 ','SIG312 ','SIG122 ','SIG222 ','SIG322 ',
1294 .'SIG132 ','SIG232 ','SIG332 ','SIG142 ','SIG242 ','SIG342 ',
1295 .'SIG152 ','SIG252 ','SIG352 ','SIG162 ','SIG262 ','SIG362 ',
1296 .'SIG172 ','SIG272 ','SIG372 ','SIG182 ','SIG282 ','SIG382 ',
1297 .'SIG192 ','SIG292 ','SIG392 ',
1298 .'SIG113 ','SIG213 ','SIG313 ','SIG123 ','SIG223 ','SIG323 ',
1299 .'SIG133 ','SIG233 ','SIG333 ','SIG143 ','SIG243 ','SIG343 ',
1300 .'SIG153 ','SIG253 ','SIG353 ','SIG163 ','SIG263 ','SIG363 ',
1301 .'SIG173 ','SIG273 ','SIG373 ','SIG183 ','SIG283 ','SIG383 ',
1302 .'SIG193 ','SIG293 ','SIG393 ',
1303 .'SIG1U1 ','SIG2U1 ','SIG3U1 ','SIG1U2 ','SIG2U2 ','SIG3U2 ',
1304 .'SIG1U3 ','SIG2U3 ','SIG3U3 ',
1305 .'SIG1D1 ','SIG2D1 ','SIG3D1 ','SIG1D2 ','SIG2D2 ','SIG3D2 ',
1306 .'SIG1D3 ','SIG2D3 ','SIG3D3 ','USER111 ','USER211 ','USER311 ',
1307 .'USER121 ','USER221 ','USER321 ','USER131 ','USER231 ','USER331 ',
1308 .'USER141 ','USER241 ','USER341 ','USER151 ','USER251 ','USER351 ',
1309 .'USER161 ','USER261 ','USER361 ','USER171 ','USER271 ','USER371 ',
1310 .'USER181 ','USER281 ','USER381 ','USER191 ','USER291 ','USER391 ',
1311 .'USER112 ','USER212 ','USER312 ','USER122 ','USER222 ','USER322 ',
1312 .'USER132 ','USER232 ','USER332 ','USER142 ','USER242 ','USER342 ',
1313 .'USER152 ','USER252 ','USER352 ','USER162 ','USER262 ','USER362 ',
1314 .'user172 ','user272 ','user372 ','user182 ','user282 ','user382 ',
1315 .'user192 ','user292 ','user392 ','user113 ','user213 ','user313 ',
1316 .'user123 ','user223 ','user323 ','user133 ','user233 ','user333 ',
1317 .'user143 ','user243 ','user343 ','user153 ','user253 ','user353 ',
1318 .'user163 ','user263 ','user363 ','user173 ','user273 ','user373 ',
1319 .'user183 ','user283 ','user383 ','user193 ','user293 ','user393 ',
1320 .'strain ','def_mat1','def_mat2'/
1321C
1322c DATA VARSG2
1323C Definition of VARSG2 : strain for /PROP/TYPE22
1324 DO J=1,200
1325 IF (J <= 9) THEN
1326 WRITE(CHJ,'(i1.1)')J
1327 ELSEIF (J <= 99) THEN
1328 WRITE(CHJ,'(i2.2)')J
1329 ELSE
1330 WRITE(CHJ,'(i3.3)')J
1331 ENDIF
1332 DO I=1,9
1333 WRITE(CHI,'(i1.1)')I
1334 DO K=1,9
1335 WRITE(CHK,'(i1.1)')K
1336 VAR_TMP(1:5)='eps_'//CHI
1337 VAR_TMP=VAR_TMP(1:5)//CHK//'_'//CHJ
1338 VARSG2((J-1)*99+(I-1)*9+K)=VAR_TMP
1339 ENDDO
1340 ENDDO
1341 ENDDO
1342c DATA VARSG3
1343C Definition of VARSG3 : stress for /PROP/TYPE22
1344 DO J=1,200
1345 IF (J <= 9) THEN
1346 WRITE(CHJ,'(i1.1)')J
1347 ELSEIF (J <= 99) THEN
1348 WRITE(CHJ,'(i2.2)')J
1349 ELSE
1350 WRITE(CHJ,'(i3.3)')J
1351 ENDIF
1352 DO I=1,9
1353 WRITE(CHI,'(i1.1)')I
1354 DO K=1,9
1355 WRITE(CHK,'(i1.1)')K
1356 VAR_TMP(1:5)='sig_'//CHI
1357 VAR_TMP=VAR_TMP(1:5)//CHK//'_'//CHJ
1358 VARSG3((J-1)*9*9+(I-1)*9+K)=VAR_TMP
1359 ENDDO
1360 ENDDO
1361 ENDDO
1362c DATA VARSG4
1363C Definition of VARSG4
1364C
1365 DATA VARSG4/
1366 .'locstrn '/
1367
1368c VARCG = VARCG1 + VARCG2+ VARCG3
1369 DATA VARCG1/
1370 . 'def ','stress ','strain ','plas ','usera ',
1371 . 'userb ','userc ','userd ','user1 ','user2 ',
1372 . 'user3 ','user4 ','user5 ','user6 ','user7 ',
1373 . 'user8 ','user9 ','user10 ','user11 ','user12 ',
1374 . 'user13 ','user14 ','user15 ','user16 ','user17 ',
1375 . 'user18 ','user19 ','user20 ',
1376 . 'sig_1 ','sig_2 ','sig_3 ','sig_4 ','sig_5 ',
1377 . 'sig_6 ','sig_7 ','sig_8 ','sig_9 ','sig_10 ',
1378 . 'sig_11 ','sig_12 ','sig_13 ','sig_14 ','sig_15 ',
1379 . 'sig_16 ','sig_17 ','sig_18 ','sig_19 ','sig_10 ',
1380 . 'sig_21 ','sig_22 ','sig_23 ','sig_24 ','sig_25 ',
1381 . 'sig_26 ','sig_27 ','sig_28 ','sig_29 ','sig_30 ',
1382 . 'sig_31 ','sig_32 ','sig_33 ','sig_34 ','sig_35 ',
1383 . 'sig_36 ','sig_37 ','sig_38 ','sig_39 ','sig_40 ',
1384 . 'sig_41 ','sig_42 ','sig_43 ','sig_44 ','sig_45 ',
1385 . 'sig_46 ','sig_47 ','sig_48 ','sig_49 ','sig_50 ',
1386 . 'sig_51 ','sig_52 ','sig_53 ','sig_54 ','sig_55 ',
1387 . 'sig_56 ','sig_57 ','sig_58 ','sig_59 ','sig_60 ',
1388 . 'sig_61 ','sig_62 ','sig_63 ','sig_64 ','sig_65 ',
1389 . 'sig_66 ','sig_67 ','sig_68 ','sig_69 ','sig_70 ',
1390 . 'sig_71 ','sig_72 ','sig_73 ','sig_74 ','sig_75 ',
1391 . 'sig_76 ','sig_77 ','sig_78 ','sig_79 ','sig_80 ',
1392 . 'sig_81 ','sig_82 ','sig_83 ','sig_84 ','sig_85 ',
1393 . 'sig_86 ','sig_87 ','sig_88 ','sig_89 ','sig_90 ',
1394 . 'sig_91 ','sig_92 ','sig_93 ','sig_94 ','sig_95 ',
1395 . 'sig_96 ','sig_97 ','sig_98 ','sig_99 ',
1396 . 'user21 ','user22 ',
1397 . 'user23 ','user24 ','user25 ','user26 ','user27 ',
1398 . 'user28 ','user29 ','user30 ','user31 ','user32 ',
1399 . 'user33 ','user34 ','user35 ','user36 ','user37 ',
1400 . 'user38 ','user39 ','user40 ','user41 ','user42 ',
1401 . 'USER43 ','USER44 ','USER45 ','USER46 ','USER47 ',
1402 . 'USER48 ','USER49 ','USER50 ','USER51 ','USER52 ',
1403 . 'USER53 ','USER54 ','USER55 ','USER56 ','USER57 ',
1404 . 'USER58 ','USER59 ','USER60 '/
1405 DATA varcg2/
1406 . 'USLAY011','USLAY012','USLAY013','USLAY014','USLAY021',
1407 . 'USLAY022','USLAY023','USLAY024','USLAY031','USLAY032',
1408 . 'USLAY033','USLAY034','USLAY041','USLAY042','USLAY043',
1409 . 'USLAY044','USLAY051','USLAY012','USLAY013','USLAY014',
1410 . 'USLAY061','USLAY062','USLAY063','USLAY064','USLAY071',
1411 . 'USLAY072','USLAY073','USLAY074','USLAY081','USLAY082',
1412 . 'USLAY083','USLAY084','USLAY091','USLAY092','USLAY093',
1413 . 'USLAY094','USLAY101','USLAY102','USLAY103','USLAY104',
1414 . 'USLAY111','USLAY112','USLAY113','USLAY114','USLAY121',
1415 . 'USLAY122','USLAY123','USLAY124','USLAY131','USLAY132',
1416 . 'USLAY133','USLAY134','USLAY141','USLAY142','USLAY143',
1417 . 'USLAY144','USLAY151','USLAY152','USLAY153','USLAY154',
1418 . 'USLAY161','USLAY162','USLAY163','USLAY164','USLAY171',
1419 . 'USLAY172','USLAY173','USLAY174','USLAY181','USLAY182',
1420 . 'USLAY183','USLAY184','USLAY191','USLAY192','USLAY193',
1421 . 'USLAY194','USLAY201','USLAY202','USLAY203','USLAY204',
1422 . 'USLAY211','USLAY212','USLAY213','USLAY214','USLAY221',
1423 . 'USLAY222','USLAY223','USLAY224','USLAY231','USLAY232',
1424 . 'USLAY233','USLAY234','USLAY241','USLAY242','USLAY243',
1425 . 'USLAY244','USLAY251','USLAY252','USLAY253','USLAY254',
1426 . 'USLAY261','USLAY262','USLAY263','USLAY264','USLAY271',
1427 . 'USLAY272','USLAY273','USLAY274','USLAY281','USLAY282',
1428 . 'USLAY283','USLAY284','USLAY291','USLAY292','USLAY293',
1429 . 'USLAY294','USLAY301','USLAY302','USLAY303','USLAY304',
1430 . 'USLAY311','USLAY312','USLAY313','USLAY314','USLAY321',
1431 . 'USLAY322','USLAY323','USLAY324','USLAY331','USLAY332',
1432 . 'USLAY333','uslay334','uslay341','uslay342','uslay343',
1433 . 'uslay344','USLAY351','USLAY352','USLAY353','USLAY354',
1434 . 'USLAY361','USLAY362','USLAY363','USLAY364','USLAY371',
1435 . 'USLAY372','USLAY373','USLAY374','USLAY381','USLAY382',
1436 . 'USLAY383','USLAY384','USLAY391','USLAY392','USLAY393',
1437 . 'USLAY394','USLAY401','USLAY402','USLAY403','uslay404',
1438 . 'uslay411','uslay412','uslay413','uslay414','uslay421',
1439 . 'uslay422','uslay423','uslay424','uslay431','uslay432',
1440 . 'uslay433','uslay434','uslay441','uslay442','uslay443',
1441 . 'uslay444','uslay451','uslay452','USLAY453','USLAY454',
1442 . 'USLAY461','USLAY462','USLAY463','USLAY464','USLAY471',
1443 . 'USLAY472','USLAY473','USLAY474','USLAY481','USLAY482',
1444 . 'USLAY483','USLAY484','USLAY491','USLAY492','USLAY493',
1445 . 'USLAY494','USLAY501','uslay502','uslay503','uslay504',
1446 . 'uslay511','uslay512','uslay513','uslay514','uslay521',
1447 . 'uslay522','uslay523','uslay524','uslay531','uslay532',
1448 . 'uslay533','uslay534','uslay541','uslay542','uslay543',
1449 . 'uslay544','uslay551','uslay552','uslay553','uslay554',
1450 . 'uslay561','uslay562','uslay563','uslay564','uslay571',
1451 . 'uslay572','uslay573','uslay574','uslay581','uslay582',
1452 . 'uslay583','uslay584','uslay591','uslay592','uslay593',
1453 . 'uslay594','uslay601','uslay602','uslay603','uslay604',
1454 . 'uslay611','uslay612','uslay613','uslay614','uslay621',
1455 . 'uslay622','uslay623','uslay624','uslay631','uslay632',
1456 . 'uslay633','uslay634','uslay641','uslay642','uslay643',
1457 . 'uslay644','uslay651','uslay652','uslay653','uslay654',
1458 . 'uslay661','uslay662','uslay663','uslay664','uslay671',
1459 . 'uslay672','uslay673','uslay674','uslay681','uslay682',
1460 . 'uslay683','uslay684','uslay691','uslay692','uslay693',
1461 . 'uslay694','uslay711','uslay712','uslay713','uslay714',
1462 . 'uslay701','uslay702','uslay703','uslay704','uslay721',
1463 . 'uslay722','uslay723','uslay724','uslay731','uslay732',
1464 . 'uslay733','uslay734','uslay741','uslay742','uslay743',
1465 . 'uslay744','uslay751','uslay752','uslay753','uslay754',
1466 . 'uslay761','uslay762','uslay763','uslay764','uslay771',
1467 . 'uslay772','uslay773','uslay774','uslay781','uslay782',
1468 . 'uslay783','uslay784','uslay791','uslay792','uslay793',
1469 . 'uslay794','uslay801','uslay802','uslay803','uslay804',
1470 . 'uslay811','uslay812','uslay813','uslay814','uslay821',
1471 . 'uslay822','uslay823','uslay824','uslay831','uslay832',
1472 . 'uslay833','uslay834','uslay841','uslay842','uslay843',
1473 . 'uslay844','uslay851','uslay852','uslay853','uslay854',
1474 . 'uslay861','uslay862','uslay863','uslay864','uslay871',
1475 . 'uslay872','uslay873','uslay874','uslay881','uslay882',
1476 . 'uslay883','uslay884','uslay891','uslay892','uslay893',
1477 . 'uslay894','uslay901','uslay902','uslay903','uslay904',
1478 . 'uslay911','uslay912','uslay913','uslay914','uslay921',
1479 . 'uslay922','uslay923','uslay924','uslay931','uslay932',
1480 . 'uslay933','uslay934','uslay941','uslay942','uslay943',
1481 . 'uslay944','uslay951','uslay952','uslay953','uslay954',
1482 . 'uslay961','uslay962','uslay963','uslay964','uslay971',
1483 . 'uslay972','uslay973','uslay974','uslay981','USLAY982',
1484 . 'USLAY983','USLAY984','USLAY991','USLAY992','USLAY993',
1485 . 'USLAY994','FAILURE ','wpl01_10','wpl11_20','wpl21_30',
1486 . 'wpl31_40','wpl41_50','wpl51_60','wpl61_70','wpl71_80',
1487 . 'wpl81_90','wpl91_99'/
1488 DATA VARCG3/
1489 . 'sigv_1 ','sigv_2 ','sigv_3 ','sigv_4 ','sigv_5 ',
1490 . 'sigv_6 ','sigv_7 ','sigv_8 ','sigv_9 ','sigv_10 ',
1491 . 'sigv_11 ','sigv_12 ','sigv_13 ','sigv_14 ','sigv_15 ',
1492 . 'sigv_16 ','sigv_17 ','sigv_18 ','sigv_19 ','sigv_10 ',
1493 . 'sigv_21 ','sigv_22 ','sigv_23 ','sigv_24 ','sigv_25 ',
1494 . 'sigv_26 ','sigv_27 ','sigv_28 ','sigv_29 ','sigv_30 ',
1495 . 'sigv_31 ','sigv_32 ','sigv_33 ','sigv_34 ','sigv_35 ',
1496 . 'sigv_36 ','sigv_37 ','sigv_38 ','sigv_39 ','sigv_40 ',
1497 . 'sigv_41 ','sigv_42 ','sigv_43 ','sigv_44 ','sigv_45 ',
1498 . 'sigv_46 ','sigv_47 ','sigv_48 ','sigv_49 ','sigv_50 ',
1499 . 'sigv_51 ','sigv_52 ','sigv_53 ','sigv_54 ','sigv_55 ',
1500 . 'sigv_56 ','sigv_57 ','sigv_58 ','sigv_59 ','sigv_60 ',
1501 . 'sigv_61 ','sigv_62 ','sigv_63 ','sigv_64 ','sigv_65 ',
1502 . 'sigv_66 ','sigv_67 ','sigv_68 ','sigv_69 ','sigv_70 ',
1503 . 'sigv_71 ','sigv_72 ','sigv_73 ','sigv_74 ','sigv_75 ',
1504 . 'sigv_76 ','sigv_77 ','sigv_78 ','sigv_79 ','sigv_80 ',
1505 . 'sigv_81 ','sigv_82 ','sigv_83 ','sigv_84 ','sigv_85 ',
1506 . 'sigv_86 ','sigv_87 ','sigv_88 ','sigv_89 ','sigv_90 ',
1507 . 'sigv_91 ','sigv_92 ','sigv_93 ','sigv_94 ','sigv_95 ',
1508 . 'sigv_96 ','sigv_97 ','sigv_98 ','sigv_99 ','sigv_100'/
1509 DATA VARCGPINCH/
1510 . 'pinch '/
1511 DATA VARTG/'def '/
1512 DATA IVARTG
1513 . / 1, 2, 3, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1514 DATA VARPG/'def '/
1515 DATA IVARPG
1516 . / 1, 2, 3, 4, 5, 6, 7, 8, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1517 DATA VARRG/'def '/
1518 DATA IVARRG
1519 . / 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,65, 0, 0, 0/
1520 DATA VARURG/'def '/
1521 DATA IVARURG
1522 . / 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12, 0, 0, 0, 0, 0, 0/
1523 DATA VARNSG/'def '/
1524 DATA IVARNSG
1525 . / 1, 2, 3, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1526 DATA IVNS2R
1527 . / 1, 2, 8,14, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1528 DATA VARSPG/
1529 . 'def ','stress '/
1530 DATA IVARSPG
1531 . / 1, 2, 3, 4, 5, 6, 7, 8, 9,12,13,41, 0, 0, 0, 0, 0, 0,
1532 . 2, 3, 4, 5, 6, 7, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1533C
1534 DATA VARING/'def ','fn ','ft ','sfw ',
1535 . '|fn| ','|f| '/
1536 DATA IVARING
1537 . / 1, 2, 3, 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1538 . 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1539 . 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1540 . 7, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1541 . 8, 9,10,11, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1542 . 12,13,14,15, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1543 DATA VARRWG/'def ','fn ','ft '/
1544 DATA IVARRWG
1545 . / 1, 2, 3, 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1546 . 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1547 . 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1548 DATA VARRBG/'def ','f ','m ','r ',
1549 . 'fi ','mi '/
1550 DATA IVARRBG
1551 . / 1, 2, 3, 4, 5, 6, 7, 8, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1552 . 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1553 . 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1554 . 7, 8, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1555 . 10,11,12, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1556 . 13,14,15, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1557 DATA VARFXG/'def '/
1558 DATA IVARFXG
1559 . / 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1560 DATA VARFXMG/'def '/
1561 DATA IVARFXMG
1562 . / 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1563 DATA VARACG/'def ','w '/
1564 DATA IVARACG
1565 . / 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1566 . 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1567 DATA VARSEG/'def ','fn ','ft ','m ','center ',
1568 . 'global ','local '/
1569 DATA IVARSEG
1570 . / 1, 2, 3, 4, 5, 6, 7, 8, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1571 . 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1572 . 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1573 . 7, 8, 9, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1574 . 37,38,39, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1575 . 1, 2, 3, 4, 5, 6,31,32,33, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1576 . 7, 8, 9,34,35,36, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1577 DATA VARJOG/'def ','f ','m '/
1578 DATA IVARJOG
1579 . / 1, 2, 3, 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1580 . 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1581 . 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1582 DATA VARABG/'def '/
1583 DATA IVARABG
1584 . / 1, 2, 3, 4, 5, 6, 7,13,14, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1585 DATA VARMVG4/'def '/,VARMVG/'def ','gas ',
1586 . 'out1 ','out2 ','out3 ','out4 ','out5 ',
1587 . 'out6 ','out7 ','out8 ','out9 ','out10 ',
1588 . 'ao1 ','bo1 ','uo1 ','mo1 ','ho1 ',
1589 . 'ao2 ','bo2 ','uo2 ','mo2 ','ho2 ',
1590 . 'ao3 ','bo3 ','uo3 ','mo3 ','ho3 ',
1591 . 'ao4 ','bo4 ','uo4 ','mo4 ','ho4 ',
1592 . 'ao5 ','bo5 ','uo5 ','mo5 ','ho5 ',
1593 . 'ao6 ','bo6 ','uo6 ','mo6 ','ho6 ',
1594 . 'ao7 ','bo7 ','uo7 ','mo7 ','ho7 ',
1595 . 'ao8 ','bo8 ','uo8 ','mo8 ','ho8 ',
1596 . 'ao9 ','bo9 ','uo9 ','mo9 ','ho9 ',
1597 . 'ao10 ','bo10 ','uo10 ','mo10 ','ho10 '/
1598 DATA VARMVENTG/'def ','gas ',
1599 . 'out1 ','out2 ','out3 ','out4 ','out5 ',
1600 . 'out6 ','out7 ','out8 ','out9 ','out10 ',
1601 . 'ao1 ','bo1 ','uo1 ','mo1 ','ho1 ',
1602 . 'ao2 ','bo2 ','uo2 ','mo2 ','ho2 ',
1603 . 'ao3 ','bo3 ','uo3 ','mo3 ','ho3 ',
1604 . 'ao4 ','bo4 ','uo4 ','mo4 ','ho4 ',
1605 . 'ao5 ','bo5 ','uo5 ','mo5 ','ho5 ',
1606 . 'ao6 ','bo6 ','uo6 ','mo6 ','ho6 ',
1607 . 'ao7 ','bo7 ','uo7 ','mo7 ','ho7 ',
1608 . 'ao8 ','bo8 ','uo8 ','mo8 ','ho8 ',
1609 . 'ao9 ','bo9 ','uo9 ','mo9 ','ho9 ',
1610 . 'ao10 ','bo10 ','uo10 ','mo10 ','ho10 '/
1611 DATA IVARMG4
1612 . / 1, 2, 3, 4, 5, 6, 7,13,14, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1613 DATA IVARMVG
1614 . / 1, 2, 3, 4, 5, 6, 7,13,14, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1615 . 10,11,12, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1616 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1617 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1618 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1619 . 0, 0, 0, 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
1677
1678 DATA IVARMVENTG
1679 . / 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1680 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1681 . 101,102,103,104,105, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1682 . 106,107,108,109,110, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1683 . 111,112,113,114,115, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1684 . 116,117,118,119,120, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1685 . 121,122,123,124,125, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1686 . 126,127,128,129,130, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1687 . 131,132,133,134,135, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1688 . 136,137,138,139,140, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1689 . 141,142,143,144,145, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1690 . 146,147,148,149,150, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1691 . 101, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1692 . 102, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1693 . 103, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1694 . 104, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1695 . 105, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1696 . 106, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1697 . 107, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1698 . 108, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1699 . 109, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1700 . 110, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1701 . 111, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1702 . 112, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1703 . 113, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1704 . 114, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1705 . 115, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1706 . 116, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1707 . 117, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1708 . 118, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1709 . 119, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1710 . 120, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1711 . 121, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1712 . 122, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1713 . 123, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1714 . 124, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1715 . 125, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1716 . 126, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1717 . 127, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1718 . 128, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1719 . 129, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1720 . 130, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1721 . 131, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1722 . 132, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1723 . 133, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1724 . 134, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1725 . 135, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1726 . 136, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1727 . 137, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1728 . 138, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1729 . 139, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1730 . 140, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1731 . 141, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1732 . 142, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1733 . 143, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1734 . 144, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1735 . 145, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1736 . 146, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1737 . 147, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1738 . 148, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1739 . 149, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1740 . 150, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1741 DATA VARPAG/'def '/
1742 DATA IVARPAG
1743 . / 1, 2, 3, 4, 5, 6, 7,28, 25, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1744 DATA VARRIVG/'def '/
1745 DATA IVARRIV
1746 . / 1, 2, 3, 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1747C frames
1748 DATA VARF1G/'def ','o+r ','v+vr '/
1749 DATA IVARF1G
1750 . / 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,15,16,17,18,
1751 . 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12, 0, 0, 0, 0, 0, 0,
1752 . 13,14,15,16,17,18, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1753 DATA VARFRG/'def ','o+r ','v+vr ','a+ar '/
1754 DATA IVARFRG
1755 . / 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,15,16,17,18,
1756 . 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12, 0, 0, 0, 0, 0, 0,
1757 . 13,14,15,16,17,18, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1758 . 19,20,21,22,23,24, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1759 DATA VARGAUG/'def '/
1760 DATA IVARGAUG
1761 . / 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1762 DATA VARCLUSG/'def ','floc '/
1763 DATA IVARCLUSG
1764 . / 1, 2, 3, 4, 5, 6,11, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
1765 . 7, 8, 9,10, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1766 DATA VARFLOWG/'def '/
1767 DATA IVARFLOWG
1768 . / 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1769 DATA VARSURFG/'def '/
1770 DATA IVARSURFG
1771 . / 1, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1772 DATA VARSLIPG/'def '/
1773 DATA IVARSLIPG
1774 . / 1, 2, 3, 4, 5, 6, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1775 DATA VARRETG/'def '/
1776 DATA IVARRETG
1777 . / 1, 2, 3, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1778 DATA VARSENSG/
1779 . 'def '/
1780 DATA IVARSENSG
1781 . / 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1782 DATA VARCHECKSUMG/
1783 . 'def '/
1784 DATA IVARCHECKSUMG
1785 . / 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
1786C=======================================================================
1787 TAGP=0
1788 TAGS=0
1789 IBID = 0
1790 BID = 0
1791 NSNE = 0
1792C-----------------------------------------------
1793C Brick
1794C ------------------------------------------
1795C (VARS1, VARS2, VARS3, VARS4, VARS5, VARS6, VARS7 )----> VARS
1796 DO I=1,NVARS1
1797 VARS(I) = VARS1(I)
1798 ENDDO
1799C
1800 DO I=1,NVARS2
1801 VARS(NVARS1+ I) = VARS2(I)
1802 ENDDO
1803C
1804 DO I=1,NVARS3
1805 VARS(NVARS1+ NVARS2 + I) = VARS3(I)
1806 ENDDO
1807 DO I=1,NVARS4
1808 VARS(NVARS1+ NVARS2 + NVARS3 + I) = VARS4(I)
1809 ENDDO
1810 DO I=1,NVARS5
1811 VARS(1624 + I) = VARS5(I)
1812 ENDDO
1813 DO I=1,NVARS6
1814 VARS(1646 + I) = VARS6(I)
1815 ENDDO
1816 DO I=1,NVARS7
1817 VARS(98846 + I) = VARS7(I)
1818 ENDDO
1819 DO I=1,NVARS8
1820 VARS(239030 + I) = VARS8(I)
1821 ENDDO
1822 DO I=1,NVARS9
1823 VARS(239546 + I) = VARS9(I)
1824 ENDDO
1825 DO I=1,NVARSNLOC
1826 VARS(239552 + I) = VARSNLOC(I)
1827 ENDDO
1828 DO I=1,NVARS10
1829 VARS(239554 + I) = VARS10(I)
1830 ENDDO
1831C (VARSG1, VARSG2, VARSG3 )----> VARSG
1832 DO I=1,202
1833 VARSG(I) = VARSG1(I)
1834 ENDDO
1835C
1836 DO I=1,16200
1837 VARSG(202 + I) = VARSG2(I)
1838 ENDDO
1839C
1840 DO I=1,16200
1841 VARSG(16402 + I) = VARSG3(I)
1842 ENDDO
1843c
1844 DO I=1,1
1845 VARSG(39766 + I) = VARSG4(I)
1846 ENDDO
1847c
1848 DO I=1,18
1849 DO J = 1,NVARSG
1850 IVARSG(I,J)=0
1851 ENDDO
1852 ENDDO
1853C Option 1
1854 DO J=1,9
1855 IVARSG(J,1)=J
1856 ENDDO
1857 IVARSG(10,1)= 12
1858 IVARSG(11,1)= 13
1859C Option 2+3
1860 IVSG = 34
1861 DO J=1,6
1862 IVARSG(J,2)= J+1
1863 IVARSG(J,3)= IVSG + J
1864 ENDDO
1865 IVSG = IVSG + 6
1866C Option 4-->15 IVSG = 40
1867 DO I=4,15
1868 DO J=1,8
1869 IVARSG(J,I)= IVSG + J
1870 ENDDO
1871 IVSG = IVSG + 8
1872 ENDDO
1873C Option 16--->19 IVSG = 136
1874 DO I=16,18
1875 DO J=1,18
1876 IVARSG(J,I)= IVSG + J
1877 ENDDO
1878 IVSG= IVSG + 18
1879 ENDDO
1880C Option 19 IVSG = 190
1881 DO J=1,6
1882 IVARSG(J,19)= IVSG + J
1883 ENDDO
1884 IVSG = IVSG + 6
1885C options 20--->100 IVSG=196
1886 DO I=20,100
1887 DO J=1,7
1888 IVARSG(J,I)= IVSG + J
1889 ENDDO
1890 IVSG = IVSG + 7
1891 ENDDO
1892C Options 101--->118 IVSG= 763
1893 DO I=101,118
1894 DO J=1,7
1895 IVARSG(J,I)= IVSG + J
1896 ENDDO
1897 IVSG = IVSG + 7
1898 ENDDO
1899C Options 118--->199 IVSG= 889
1900 DO I=119,199
1901 DO J=1,9
1902 IVARSG(J,I)= IVSG + J
1903 ENDDO
1904 IVSG = IVSG + 9
1905 ENDDO
1906 DO J=1,6
1907 IVARSG(J,200)= IVSG + J
1908 ENDDO
1909 IVSG = IVSG + 6
1910C Options 201--->202 IVSG= 1624
1911 DO I=201,202
1912 DO J=1,11
1913 IVARSG(J,I)= IVSG + J
1914 ENDDO
1915 IVSG = IVSG + 11
1916 ENDDO
1917C Options 203--->19984 IVSG= 1646
1918 DO I=203,19984
1919 DO J=1,6
1920 IVARSG(J,I)= IVSG + J
1921 ENDDO
1922 IVSG = IVSG + 6
1923 ENDDO
1924C Options 203--->19984 IVSG= 120338
1925 DO I=19985,39766
1926 DO J=1,6
1927 IVARSG(J,I)= IVSG + J
1928 ENDDO
1929 IVSG = IVSG + 6
1930 ENDDO
1931C Option LOCSTRN
1932 DO J=1,6
1933 IVARSG(J,39767)= 239030 + J
1934 ENDDO
1935C
1936 DO I=1,208
1937 VARC(I) = VARC1(I)
1938 ENDDO
1939 DO I=1,470
1940 VARC(208 + I) = VARC2(I)
1941 ENDDO
1942Ccc USRII_KK (II=1,20, KK=6,99
1943 DO II=1,20
1944 WRITE(CHII,'(i2.2)')II
1945 DO KK= 6,99
1946 WRITE(CHKK,'(i2.2)')KK
1947 IF(CHII(1:1) == '0')THEN
1948 VAR = 'usr'//CHII(2:2)//'_'
1949 IF(CHKK(1:1) == '0')THEN
1950 VAR = VAR(1:5)//CHKK(2:2)
1951 ELSE
1952 VAR = VAR(1:5)//CHKK
1953 ENDIF
1954 ELSE
1955 VAR ='usr'//CHII//'_'
1956 IF(CHKK(1:1) == '0')THEN
1957 VAR = VAR(1:6)//CHKK(2:2)
1958 ELSE
1959 VAR = VAR(1:6)//CHKK
1960 ENDIF
1961 ENDIF
1962 VARC(678 + (II -1)*94 + KK - 5) = VAR
1963 ENDDO
1964 ENDDO
1965Ccc USRII_KK(II=21,60, KK=1,99)
1966 DO II=21,60
1967 WRITE(CHII,'(i2.2)')II
1968 DO KK= 1,99
1969 WRITE(CHKK,'(i2.2)')KK
1970 IF(CHII(1:1) == '0')THEN
1971 VAR = 'usr'//CHII(2:2)//'_'
1972 IF(CHKK(1:1) == '0')THEN
1973 VAR = VAR(1:5)//CHKK(2:2)
1974 ELSE
1975 VAR = VAR(1:5)//CHKK
1976 ENDIF
1977 ELSE
1978 VAR ='usr'//CHII//'_'
1979 IF(CHKK(1:1) == '0')THEN
1980 VAR = VAR(1:6)//CHKK(2:2)
1981 ELSE
1982 VAR = VAR(1:6)//CHKK
1983 ENDIF
1984 ENDIF
1985 VARC(2558 + (II - 21)*99 + KK ) = VAR
1986 ENDDO
1987 ENDDO
1988C USII_jKK (II=01, 60, j=1,4, KK=01,99)
1989 DO KK = 1, 99
1990 WRITE(CHKK,'(i2.2)')KK
1991 DO J = 1,4
1992 JJ = J*100 + KK
1993 WRITE(CHJJ,'(i3.3)')JJ
1994 DO II = 1,60
1995 WRITE(CHII,'(i2.2)')II
1996 VAR ='us'//CHII//'_'//CHJJ
1997 VARC(6518 + II + (KK - 1)*240 + (J-1)*60)=VAR
1998 ENDDO
1999 ENDDO
2000 ENDDO
2001 DO I = 1,104
2002 VARC(30278 + I ) = VARC2(470 + I)
2003 ENDDO
2004C
2005C viscous stress for shell 5*100
2006CC USII_jKK (II=01, 60, j=1,4, KK=01,99)
2007 DO I = 1,500
2008 VARC(30382 + I ) = VARC3(I)
2009 ENDDO
2010C NVC_old = 30882
2011C stress, strain, damage for interply
2012C SZILYjKK (j=1,4, KK=001,200)
2013C SYZILYjKK (j=1,4, KK=001,200)
2014C SXZILYjKK (j=1,4, KK=001,200)
2015C EZILYjKK (j=1,4, KK=001,200)
2016C EYZILYjKK (j=1,4, KK=001,200)
2017C EXZILYjKK (j=1,4, KK=001,200)
2018C SZILY_KK (KK=001,200)
2019C SYZILY_KK ( KK=001,200)
2020C SXZILY_KK ( KK=001,200)
2021C EZILY_KK (KK=001,200)
2022C EYZILY_KK ( KK=001,200)
2023C EXZILY_KK ( KK=001,200)
2024C DAM1ILYKK ( KK=001,200)
2025C DAM2ILYKK (KK=001,200)
2026C DAM3ILYKK ( KK=001,200)
2027C DAM4ILYKK ( KK=001,200)
2028C EINTILKK
2029 DO KK = 1, 199
2030C 3 stress + 3 strain + 3 average stress + 3 average strain + 4 max damage
2031C + EINT
2032 II = 30882 + (KK -1)*(6*4 + 11)
2033 DO J = 1,4
2034 JJ = J*1000 + KK
2035 WRITE(CHJKK ,'(i4.4)')JJ
2036 DO L=1,6
2037 IF (L==1) VAR= 'szily'//CHJKK
2038 IF (L==2) VAR ='syzily'//CHJKK
2039 IF (L==3) VAR ='sxzily'//CHJKK
2040C
2041 IF (L==4) VAR= 'ezily'//CHJKK
2042 IF (L==5) VAR ='eyzily'//CHJKK
2043 IF (L==6) VAR ='exzily'//CHJKK
2044 VARC(II + L ) = VAR
2045 ENDDO
2046 II = II + 6
2047 ENDDO
2048 WRITE(CHJJ,'(i3.3)')KK
2049 DO L = 1,6
2050 IF (L==1) VAR = 'szily'//'_'//CHJJ
2051 IF (L==2) VAR = 'syzily'//'_'//CHJJ
2052 IF (L==3) VAR = 'sxzily'//'_'//CHJJ
2053 IF (L==4) VAR = 'ezily'//'_'//CHJJ
2054 IF (L==5) VAR = 'eyzily'//'_'//CHJJ
2055 IF (L==6) VAR = 'exzily'//'_'//CHJJ
2056 VARC(II + L ) = VAR
2057 ENDDO
2058 II =II +6
2059 DO L = 1,4
2060 IF (L==1) VAR = 'damily'//CHJJ
2061 IF (L==2) VAR = 'dam2ily'//CHJJ
2062 IF (L==3) VAR = 'dam3ily'//CHJJ
2063 IF (L==4) VAR = 'dam4ily'//CHJJ
2064 VARC(II + L ) = VAR
2065 ENDDO
2066 II = II + 4
2067 VAR = 'eintily'//CHJJ
2068 VARC(II + 1 ) = VAR
2069 II = II + 1
2070 ENDDO
2071C---pinching
2072 DO I=1,7
2073 VARC(37847+I) = VARCPINCH(I)
2074 ENDDO
2075C---non-local
2076 DO I=1,2
2077 VARC(37854+I) = VARCNLOC(I)
2078 ENDDO
2079C
2080C NVC = 30882 + 199*34 = 37847
2081Cccc array dimension of VARC =
2082C dimension_varc =83 + 99*60 + 5*99 + 4*99*60 + 104 + 500 ! 30882
2083C + (4*6 + 11)*199 = 37847
2084C + (4*6 + 11)*199 = 37847
2085C-----------------------------------------------------
2086C GENERALIZED KEYWORDS
2087C-----------------------------------------------------
2088 DO I =1,167
2089 VARCG(I)= VARCG1(I)
2090 ENDDO
2091 DO I = 1,407
2092 VARCG(167 + I) = VARCG2(I)
2093 ENDDO
2094 DO I = 1,100
2095 VARCG(574 + I) = VARCG3(I)
2096 ENDDO
2097 VARCG(675) = VARCGPINCH(1)
2098C
2099 DO I=1,18
2100 DO J = 1, NVARCG
2101 IVARCG(I,J)=0
2102 ENDDO
2103 ENDDO
2104C keyword 1
2105 DO I=1,3
2106 IVARCG(I,1) = I
2107 IVARCG(I+3,1) = 5 + I
2108 IVARCG(I+6,1) = 8 + I
2109 ENDDO
2110 IVARCG(10,1) = 13
2111 IVARCG(11,1) = 14
2112C keyword 2 + 3
2113 DO I = 1,8
2114 IVARCG(I,2) = I
2115 IVARCG(I,3) = 14 + I
2116 ENDDO
2117C keyword 4
2118 IVARCG(1,4) = 13
2119 IVARCG(2,4) = 14
2120C keyword 5,6,7
2121 DO I =1,18
2122 IVARCG(I,5)= 23 + I
2123 IVARCG(I,6)= 41 + I
2124 IVARCG(I,7)= 59 + I
2125 ENDDO
2126C keyword 8
2127 DO I =1,6
2128 IVARCG(I,8)= 77 + I
2129 ENDDO
2130 IVCG = 83
2131C keyword 9 -----> 28
2132 II = 678
2133 DO J = 9, 28
2134 DO I = 1,5
2135 IVARCG(I,J) = IVCG + I
2136 ENDDO
2137 IVCG = IVCG + 5
2138 DO I=6,18
2139 IVARCG(I,J) = II + I - 5
2140 ENDDO
2141 II = II + 94
2142 ENDDO
2143C keyword 29 -----> 127
2144 DO J = 29, 127
2145 DO I = 1,5
2146 IVARCG(I,J)= IVCG + I
2147 ENDDO
2148 IVCG = IVCG + 5
2149 ENDDO
2150C c ivcg = 183 + 5*99 + 20*94
2151 IVCG = IVCG + 20*94
2152C keyword 128 - 167
2153 DO J=128,167
2154 DO I=1,18
2155 IVARCG(I,J) = IVCG + I
2156 ENDDO
2157 IVCG = IVCG + 99
2158 ENDDO
2159C keyword 167 -----> 563
2160 DO J = 168, 563
2161 DO I = 1,18
2162 IVARCG(I,J)= IVCG + I
2163 ENDDO
2164 IVCG = IVCG + 60
2165 ENDDO
2166C keyword 564
2167 DO I = 1,5
2168 IVARCG(I,564)= IVCG + I
2169 ENDDO
2170 IVCG = IVCG + 5
2171c keyword 565 -----> 573
2172 DO J = 565,573
2173 DO I = 1,10
2174 IVARCG(I,J)= IVCG + 1
2175 IVCG = IVCG + 1
2176 ENDDO
2177 ENDDO
2178c keyword 574
2179 DO I = 1,9
2180 IVARCG(I,574)= IVCG + 1
2181 IVCG = IVCG + 1
2182 ENDDO
2183C keyword de 575 -----> 674
2184 DO J = 575, 674
2185 DO I = 1,5
2186 IVARCG(I,J)= IVCG + I
2187 ENDDO
2188 IVCG = IVCG + 5
2189 ENDDO
2190C pinching
2191 DO I = 1,7
2192 IVARCG(I,675)= 37847 + I
2193 ENDDO
2194C
2195C to ctrole the value of ivcg
2196C icvg = 30882
2197C end . ivcg = ivcg + 4*99*60 + 5 + 99 + 100
2198C--------------fin --------------------
2199C--------------fin --------------------
2200C ---- node + displacement by ply for shell
2201 DO I=1,NVARN1
2202 VARN(I) = VARN1(I)
2203 ENDDO
2204 DO J= 1,200
2205 WRITE(CHJJ,'(i3.3)')J
2206 VAR ='dx_'//CHJJ
2207 VARN(NVARN1 + (J-1)*3 + 1 )= VAR
2208 VARN1A((J-1)*3 + 1 )= VAR
2209 VAR ='dy_'//CHJJ
2210 VARN(NVARN1 + (J-1)*3 + 2) = VAR
2211 VARN1A((J-1)*3 + 2 )= VAR
2212 VAR ='dz_'//CHJJ
2213 VARN(NVARN1 + (J-1)*3 + 3) = VAR
2214 VARN1A((J-1)*3 + 3 )= VAR
2215 ENDDO
2216 DO I=1,NVARN2
2217 VARN(NVARN1+NVARN1A+I) = VARN2(I)
2218 ENDDO
2219 DO I=1,NVARNPINCH
2220 VARN(NVARN1+NVARN1A+NVARN2+I) = VARNPINCH(I)
2221 ENDDO
2222 NNOD = 0
2223 CHNOD = 78
2224 NSHEL = 0
2225 NTRUS = 0
2226 NBEAM = 0
2227 NSPRG = 0
2228 NCHTG = 0
2229 NCHINT = 0
2230 NCHJNT = 0
2231C-----------------------------------------------
2232 NSNOD=0
2233 NSELS=0
2234 NSELQ=0
2235 NSELC=0
2236 NSELT=0
2237 NSELP=0
2238 NSELR=0
2239 NSELTG=0
2240 NSELU=0
2241 NSMAT=0
2242 NSRBY=0
2243C
2244 NSFXBY=0
2245C
2246 NSMOD=0
2247C for parts and subsets:
2248 NTHWA=100
2249 NSELSP=0
2250 NSELRT=0
2251C
2252 OK=0
2253 IS_AVAILABLE = .FALSE.
2254 IGS=0
2255 IAD=1
2256 HM_NTHGRP=0
2257c
2258 ! Choose TH file type
2259 IF (IFLAG == 0) THEN
2260 THFILE = '/th'
2261 ELSEIF (IFLAG == 1) THEN
2262 THFILE = '/ath'
2263 ELSEIF (IFLAG == 2) THEN
2264 THFILE = '/bth'
2265 ELSEIF (IFLAG == 3) THEN
2266 THFILE = '/cth'
2267 ELSEIF (IFLAG == 4) THEN
2268 THFILE = '/dth'
2269 ELSEIF (IFLAG == 5) THEN
2270 THFILE = '/eth'
2271 ELSEIF (IFLAG == 6) THEN
2272 THFILE = '/fth'
2273 ELSEIF (IFLAG == 7) THEN
2274 THFILE = '/gth'
2275 ELSEIF (IFLAG == 8) THEN
2276 THFILE = '/hth'
2277 ELSEIF (IFLAG == 9) THEN
2278 THFILE = '/ith'
2279 ENDIF
2280C
2281C
2282C-----------------------------------------------
2283! ************************** !
2284! /TH read with hm reader !
2285! ************************** !
2286C-----------------------------------------------
2287 !-------------------------------------------
2288 ! /TH/ACCEL
2289 !-------------------------------------------
2290 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/accel' ,NTHACCEL )
2291 IF (NTHACCEL > 0) THEN
2292 CALL HM_OPTION_START(TRIM(THFILE)//'/accel')
2293 HM_NTHGRP = HM_NTHGRP + NTHACCEL
2294 DO I = 1, NTHACCEL
2295 TITR = ''
2296 CALL HM_OPTION_READ_KEY(LSUBMODEL,
2297 . OPTION_ID = ID,
2298 . OPTION_TITR = TITR)
2299C
2300 IGS = IGS+1
2301 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2302 ITHGRP(1,IGS) = ID
2303 NVALL = 3
2304 KEY = 'accel '
2305 CALL HM_READ_THGRKI(
2306 1 108 ,KEY ,PTR_NOPT_ACC ,
2307 2 IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2308 3 NVARAC ,VARAC ,NACCELM ,VARACG ,
2309 4 NVARACG ,IVARACG ,NSNE ,NVALL,ITHVAR ,
2310 5 FLAGABF ,NVARABF ,NOM_OPT,IGS ,LSUBMODEL,LITHBUFMX)
2311 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2312 ENDDO
2313 ENDIF
2314 !-------------------------------------------
2315 ! /TH/INTER
2316 !-------------------------------------------
2317 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/inter' ,NTHINTER )
2318 IF (NTHINTER > 0) THEN
2319 CALL HM_OPTION_START(TRIM(THFILE)//'/inter')
2320 HM_NTHGRP = HM_NTHGRP + NTHINTER
2321 DO I = 1, NTHINTER
2322 TITR = ''
2323 CALL HM_OPTION_READ_KEY(LSUBMODEL,
2324 . OPTION_ID = ID,
2325 . OPTION_TITR = TITR)
2326C
2327 IGS = IGS+1
2328 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2329 ITHGRP(1,IGS) = ID
2330 CHELT = 258
2331 CNT1 = 0
2332 NVALL = 7
2333 KEY = 'inter '
2334 CALL HM_READ_THGRKI(
2335 1 101 ,KEY ,PTR_NOPT_INTER,
2336 2 IAD ,IFI ,ITHGRP(1,IGS) ,ITHBUF ,
2337 3 NVARIN ,VARIN ,NINTER+NINTSUB,VARING ,NVARING ,
2338 4 IVARING ,NSNE ,NVALL ,ITHVAR,FLAGABF ,NVARABF ,
2339 5 NOM_OPT ,IGS ,LSUBMODEL, LITHBUFMX)
2340 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2341 DO JJ = ITHGRP(7,IGS),ITHGRP(5,IGS)
2342 IF(ITHBUF(JJ) ==29) INTERFACES%PARAMETERS%INTCAREA =1
2343 ENDDO
2344 IF(INTERFACES%PARAMETERS%INTCAREA > 0 ) THEN
2345 DO JJ = ITHGRP(5,IGS),ITHGRP(5,IGS)+ITHGRP(4,IGS)-1
2346 NIN = ITHBUF(JJ)
2347 IF(NIN <= NINTER) THEN
2348.OR. IF(IPARI(7,NIN)==24 IPARI(7,NIN)==25) IPARI(99,NIN) = 1
2349 ENDIF
2350 ENDDO
2351 ENDIF
2352 ENDDO
2353 ENDIF
2354 !-------------------------------------------
2355 ! /TH/RWALL
2356 !-------------------------------------------
2357 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/rwall' ,NTHRWALL )
2358 IF (NTHRWALL > 0) THEN
2359 CALL HM_OPTION_START(TRIM(THFILE)//'/rwall')
2360 HM_NTHGRP = HM_NTHGRP + NTHRWALL
2361 DO I = 1, NTHRWALL
2362 TITR = ''
2363 CALL HM_OPTION_READ_KEY(LSUBMODEL,
2364 . OPTION_ID = ID,
2365 . OPTION_TITR = TITR)
2366C
2367 IGS = IGS+1
2368 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2369 ITHGRP(1,IGS) = ID
2370 NVALL = 6
2371 KEY = 'rwall '
2372 CALL HM_READ_THGRKI(
2373 1 102 ,KEY ,PTR_NOPT_RWALL,
2374 2 IAD ,IFI ,ITHGRP(1,IGS) ,ITHBUF ,
2375 3 NVARRW ,VARRW ,NRWALL ,VARRWG ,NVARRWG ,
2376 4 IVARRWG ,NSNE ,NVALL,ITHVAR ,FLAGABF ,NVARABF ,
2377 5 NOM_OPT ,IGS ,LSUBMODEL, LITHBUFMX)
2378 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2379 ENDDO
2380 ENDIF
2381 !-------------------------------------------
2382 ! /TH/SECTIO
2383 !-------------------------------------------
2384 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/sectio' ,NTHSECTIO)
2385 IF (NTHSECTIO > 0) THEN
2386 CALL HM_OPTION_START(TRIM(THFILE)//'/sectio')
2387 HM_NTHGRP = HM_NTHGRP + NTHSECTIO
2388 DO I = 1, NTHSECTIO
2389 TITR = ''
2390 CALL HM_OPTION_READ_KEY(LSUBMODEL,
2391 . OPTION_ID = ID,
2392 . OPTION_TITR = TITR)
2393C
2394 IGS = IGS+1
2395 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2396 ITHGRP(1,IGS) = ID
2397 NVALL = 13
2398 KEY = 'sectio '
2399 CALL HM_READ_THGRKI(
2400 1 104 ,KEY ,PTR_NOPT_SECT,
2401 2 IAD ,IFI ,ITHGRP(1,IGS) ,ITHBUF ,
2402 3 NVARSE ,VARSE ,NSECT ,VARSEG ,NVARSEG ,
2403 4 IVARSEG ,NSNE ,NVALL ,ITHVAR ,FLAGABF ,NVARABF ,
2404 5 NOM_OPT ,IGS ,LSUBMODEL, LITHBUFMX)
2405 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2406 ENDDO
2407 ENDIF
2408 !-------------------------------------------
2409 ! /TH/CLUSTER
2410 !-------------------------------------------
2411 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/cluster' ,NTHCLUS)
2412 IF (NTHCLUS > 0) THEN
2413 CALL HM_OPTION_START(TRIM(THFILE)//'/cluster')
2414 HM_NTHGRP = HM_NTHGRP + NTHCLUS
2415 DO I = 1, NTHCLUS
2416 TITR = ''
2417 CALL HM_OPTION_READ_KEY(LSUBMODEL,
2418 . OPTION_ID = ID,
2419 . OPTION_TITR = TITR)
2420C
2421 IGS = IGS+1
2422 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2423 ITHGRP(1,IGS) = ID
2424 NVALL = 11
2425 ITYP = 114
2426 KEY = 'cluster '
2427 CALL HM_READ_THGRKI(
2428 1 ITYP ,KEY ,PTR_NOPT_CLUS,
2429 3 IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2430 4 NVARCLUS ,VARCLUS ,NCLUSTER ,VARCLUSG ,NVARCLUSG ,
2431 5 IVARCLUSG,NSNE ,NVALL ,ITHVAR,FLAGABF ,NVARABF ,
2432 6 NOM_OPT ,IGS ,LSUBMODEL, LITHBUFMX )
2433 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2434 ENDDO
2435 ENDIF
2436 !-------------------------------------------
2437 ! /TH/BEAM
2438 !-------------------------------------------
2439 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/beam' ,NTHBEAM)
2440 IF (NTHBEAM > 0) THEN
2441 CALL HM_OPTION_START(TRIM(THFILE)//'/beam')
2442 HM_NTHGRP = HM_NTHGRP + NTHBEAM
2443 DO I = 1, NTHBEAM
2444 TITR = ''
2445 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
2446 IGS = IGS+1
2447 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2448 ITHGRP(1,IGS) = ID
2449 NVALL = 11
2450 ITYP = 5
2451 KEY = 'beam '
2452 CALL HM_READ_THGRNE(
2453 1 ITYP ,KEY ,ITAB ,ITABM1,IXP ,
2454 3 NIXP ,IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2455 4 NVARP ,VARP ,NUMELP,VARPG ,NVARPG ,
2456 5 IVARPG,NSELP,NVALL,IMERGE,ITHVAR ,ITHERM_FE,
2457 6 FLAGABF,NVARABF,RFI,LSUBMODEL,MAP_TABLES%IBEAMM, NUMELP)
2458
2459 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2460 ENDDO
2461 ENDIF
2462 !-------------------------------------------
2463 ! /TH/TRUS
2464 !-------------------------------------------
2465 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/truss' ,NTHTRUS)
2466 IF (NTHTRUS > 0) THEN
2467 CALL HM_OPTION_START(TRIM(THFILE)//'/truss')
2468 hm_nthgrp = hm_nthgrp + nthtrus
2469 DO i = 1, nthtrus
2470 titr = ''
2471 CALL hm_option_read_key(lsubmodel, option_id=id, option_titr=titr)
2472 igs = igs+1
2473 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2474 ithgrp(1,igs) = id
2475 nvall = 6
2476 ityp = 4
2477 key = 'truss '
2478 CALL HM_READ_THGRNE(
2479 1 ITYP ,KEY ,ITAB ,ITABM1,IXT ,
2480 3 NIXT ,IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2481 4 NVART ,VART ,NUMELT,VARTG ,NVARTG ,
2482 5 IVARTG,NSELT,NVALL,IMERGE,ITHVAR ,ITHERM_FE,
2483 6 FLAGABF,NVARABF,RFI,LSUBMODEL,MAP_TABLES%ITRUSSM, NUMELT)
2484 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2485 ENDDO
2486 ENDIF
2487 !-------------------------------------------
2488 ! /TH/SPRI
2489 !-------------------------------------------
2490 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/spring' ,NTHSPRING)
2491 IF (NTHSPRING > 0) THEN
2492 CALL HM_OPTION_START(TRIM(THFILE)//'/spring')
2493 HM_NTHGRP = HM_NTHGRP + NTHSPRING
2494 DO I = 1, NTHSPRING
2495 TITR = ''
2496 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
2497 IGS = IGS+1
2498 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2499 ITHGRP(1,IGS) = ID
2500 NVALL = 16
2501 ITYP = 6
2502 KEY = 'spring '
2503 CALL HM_READ_THGRNE(
2504 1 ITYP ,KEY ,ITAB ,ITABM1,IXR ,
2505 3 NIXR ,IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2506 4 NVARR ,VARR ,NUMELR,VARRG ,NVARRG ,
2507 5 IVARRG,NSELR,NVALL,IMERGE,ITHVAR ,ITHERM_FE,
2508 6 FLAGABF,NVARABF,RFI,LSUBMODEL,MAP_TABLES%ISPRINGM, NUMELR)
2509 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2510 ENDDO
2511 ENDIF
2512 !-------------------------------------------
2513 ! /TH/BRIC
2514 !-------------------------------------------
2515 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/bric' ,NTHBRIC)
2516 IF (NTHBRIC > 0) THEN
2517 CALL HM_OPTION_START(TRIM(THFILE)//'/bric')
2518 HM_NTHGRP = HM_NTHGRP + NTHBRIC
2519 DO I = 1, NTHBRIC
2520 TITR = ''
2521 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
2522 IGS = IGS+1
2523 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2524 ITHGRP(1,IGS) = ID
2525 NVALL = 40
2526 ITYP = 1
2527 KEY = 'brick '
2528 CALL HM_READ_THGRNE(
2529 1 ITYP ,KEY ,ITAB ,ITABM1,IXS ,
2530 3 NIXS ,IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2531 4 NVARS ,VARS ,NUMELS,VARSG ,NVARSG ,
2532 5 IVARSG,NSELS,NVALL,IMERGE,ITHVAR ,ITHERM_FE,
2533 6 FLAGABF,NVARABF,RFI,LSUBMODEL,MAP_TABLES%ISOLM, NUMELS)
2534 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2535 ENDDO
2536 ENDIF
2537 !-------------------------------------------
2538 ! /TH/NODE
2539 !-------------------------------------------
2540 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/node' ,NTHNODE)
2541 IF (NTHNODE > 0) THEN
2542 CALL HM_OPTION_START(TRIM(THFILE)//'/node')
2543 HM_NTHGRP = HM_NTHGRP + NTHNODE
2544 DO I = 1, NTHNODE
2545 TITR = ''
2546 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
2547 IGS = IGS+1
2548 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2549 ITHGRP(1,IGS) = ID
2550 NVALL = 9
2551 ITYP = 0
2552 KEY = 'node '
2553 CALL HM_READ_THGRNE(
2554 1 ITYP ,KEY ,ITAB ,ITABM1 ,IXS ,
2555 3 NIXS ,IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2556 4 NVARN ,VARN ,0 ,VARNG ,NVARNG ,
2557 5 IVARNG,NSNOD ,NVALL ,IMERGE ,ITHVAR ,ITHERM_FE,
2558 6 FLAGABF,NVARABF,RFI,LSUBMODEL,IBID,IBID)
2559C decode les numeros de skews.
2560 NNE =ITHGRP(4,IGS)
2561 IAD2=ITHGRP(5,IGS)+NNE
2562 DO K=1,NNE
2563 ISK=ITHBUF(IAD2)
2564 DO 100 J=1,(NUMSKW+1)+NSUBMOD
2565 IF(ISK == ISKWN(4,J)) THEN
2566 ITHBUF(IAD2)=J-1
2567 GO TO 110
2568 ENDIF
2569 100 CONTINUE
2570C decode les numeros de frames.
2571 DO 101 J=(NUMSKW+1)+MIN(IUN,NSPCOND)*NUMSPH+1+NSUBMOD,
2572 . (NUMSKW+1)+MIN(IUN,NSPCOND)*NUMSPH+NSUBMOD+
2573 . (NUMFRAM+1)
2574 IF(ISK == ISKWN(4,J)) THEN
2575 ITHBUF(IAD2)=J-1
2576 GO TO 110
2577 ENDIF
2578 101 CONTINUE
2579 CALL ANCMSG(MSGID=434,
2580 . MSGTYPE=MSGERROR,
2581 . ANMODE=ANINFO_BLIND_1,
2582 . I1=ID,
2583 . C1=TITR,
2584 . I2=ISK)
2585 110 CONTINUE
2586 IAD2=IAD2+1
2587 ENDDO
2588 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2589 ENDDO
2590 ENDIF
2591 !-------------------------------------------
2592 ! /TH/SHEL
2593 !-------------------------------------------
2594 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/shel' ,NTHSHEL)
2595 IF (NTHSHEL > 0) THEN
2596 CALL HM_OPTION_START(TRIM(THFILE)//'/shel')
2597 HM_NTHGRP = HM_NTHGRP + NTHSHEL
2598 DO I = 1, NTHSHEL
2599 TITR = ''
2600 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
2601 IGS = IGS+1
2602 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2603 ITHGRP(1,IGS) = ID
2604 NVALL = 22
2605 ITYP = 3
2606 KEY = 'shell '
2607 CALL HM_READ_THGRNE(
2608 1 ITYP ,KEY ,ITAB ,ITABM1,IXC ,
2609 3 NIXC ,IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2610 4 NVARC ,VARC ,NUMELC,VARCG ,NVARCG ,
2611 5 IVARCG,NSELC,NVALL,IMERGE,ITHVAR ,ITHERM_FE,
2612 6 FLAGABF,NVARABF,RFI,LSUBMODEL,MAP_TABLES%ISH4NM, NUMELC)
2613C decode les numeros de skews.
2614 NNE =ITHGRP(4,IGS)
2615 IAD2=ITHGRP(5,IGS)+3*NNE
2616 DO K=1,NNE
2617 ISK=ITHBUF(IAD2)
2618 DO 200 J=1,(NUMSKW+1)+NSUBMOD
2619 IF(ISK == ISKWN(4,J)) THEN
2620 ITHBUF(IAD2)=J-1
2621 GO TO 210
2622 ENDIF
2623 200 CONTINUE
2624C frame numbers user -> internal
2625 DO 201 J=(NUMSKW+1)+MIN(IUN,NSPCOND)*NUMSPH+1+NSUBMOD,
2626 . (NUMSKW+1)+MIN(IUN,NSPCOND)*NUMSPH+NSUBMOD+
2627 . (NUMFRAM+1)
2628 IF(ISK == ISKWN(4,J)) THEN
2629 ITHBUF(IAD2)=J-1
2630 GO TO 210
2631 ENDIF
2632 201 CONTINUE
2633 CALL ANCMSG(MSGID=434,
2634 . MSGTYPE=MSGERROR,
2635 . ANMODE=ANINFO_BLIND_1,
2636 . I1=ID,
2637 . C1=TITR,
2638 . I2=ISK)
2639 210 CONTINUE
2640 IAD2=IAD2+1
2641 ENDDO
2642 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2643 ENDDO
2644 ENDIF
2645 !-------------------------------------------
2646 ! /TH/SH3N
2647 !-------------------------------------------
2648 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/sh3n' ,NTHSH3N)
2649 IF (NTHSH3N > 0) THEN
2650 CALL HM_OPTION_START(TRIM(THFILE)//'/sh3n')
2651 HM_NTHGRP = HM_NTHGRP + NTHSH3N
2652 DO I = 1, NTHSH3N
2653 TITR = ''
2654 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
2655 IGS = IGS+1
2656 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2657 ITHGRP(1,IGS) = ID
2658 NVALL = 22
2659 ITYP = 7
2660 KEY = 'sh3n '
2661 CALL HM_READ_THGRNE(
2662 1 ITYP ,KEY ,ITAB ,ITABM1,IXTG ,
2663 3 NIXTG ,IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2664 4 NVARC ,VARC ,NUMELTG,VARCG ,NVARCG ,
2665 5 IVARCG,NSELTG,NVALL,IMERGE,ITHVAR ,ITHERM_FE,
2666 6 FLAGABF,NVARABF,RFI,LSUBMODEL,MAP_TABLES%ISH3NM, NUMELTG)
2667C skew numbers user -> internal
2668 NNE =ITHGRP(4,IGS)
2669 IAD2=ITHGRP(5,IGS)+3*NNE
2670 DO K=1,NNE
2671 ISK=ITHBUF(IAD2)
2672 DO 300 J=1,(NUMSKW+1)+NSUBMOD
2673 IF(ISK == ISKWN(4,J)) THEN
2674 ITHBUF(IAD2)=J-1
2675 GO TO 310
2676 ENDIF
2677 300 CONTINUE
2678C frame numbers user -> internal
2679 DO 301 J=(NUMSKW+1)+MIN(IUN,NSPCOND)*NUMSPH+1+NSUBMOD,
2680 . (NUMSKW+1)+MIN(IUN,NSPCOND)*NUMSPH+NSUBMOD+
2681 . (NUMFRAM+1)
2682 IF(ISK == ISKWN(4,J)) THEN
2683 ITHBUF(IAD2)=J-1
2684 GO TO 310
2685 ENDIF
2686 301 CONTINUE
2687 CALL ANCMSG(MSGID=434,
2688 . MSGTYPE=MSGERROR,
2689 . ANMODE=ANINFO_BLIND_1,
2690 . I1=ID,
2691 . C1=TITR,
2692 . I2=ISK)
2693 310 CONTINUE
2694 IAD2=IAD2+1
2695 ENDDO
2696 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2697 ENDDO
2698 ENDIF
2699 !-------------------------------------------
2700 ! /TH/RBODY
2701 !-------------------------------------------
2702 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/rbody' ,NTHRBODY)
2703 IF (NTHRBODY > 0) THEN
2704 CALL HM_OPTION_START(TRIM(THFILE)//'/rbody')
2705 HM_NTHGRP = HM_NTHGRP + NTHRBODY
2706 DO I = 1, NTHRBODY
2707 TITR = ''
2708 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
2709 IGS = IGS+1
2710 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2711 ITHGRP(1,IGS) = ID
2712 NVALL = 15
2713 ITYP = 103
2714 KEY = 'rbody '
2715 CALL HM_READ_THGRKI_RBODY(
2716 1 ITYP ,KEY ,0,
2717 3 IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2718 4 NVARRB ,VARRB ,NRBODY ,VARRBG ,NVARRBG,
2719 5 IVARRBG,NSRBY ,NVALL ,ITHVAR ,FLAGABF,NVARABF,
2720 6 NOM_OPT,IGS ,NPBY ,LSUBMODEL )
2721 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2722 ENDDO
2723 ENDIF
2724 !-------------------------------------------
2725 ! /TH/MONVOL
2726 !-------------------------------------------
2727 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/monv' ,NTHMONVOL)
2728 IF (NTHMONVOL > 0) THEN
2729 CALL HM_OPTION_START(TRIM(THFILE)//'/monv')
2730 HM_NTHGRP = HM_NTHGRP + NTHMONVOL
2731 DO I = 1, NTHMONVOL
2732 TITR = ''
2733 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
2734 IGS = IGS+1
2735 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2736 ITHGRP(1,IGS) = ID
2737 NVALL = 9
2738 ITYP = 107
2739 KEY = 'monvol '
2740 CALL HM_READ_THGRKI(
2741 1 ITYP ,KEY ,PTR_NOPT_MONV,
2742 3 IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2743 4 NVARMV ,VARMV ,NVOLU ,VARMVG,NVARMVG,
2744 5 IVARMVG,NSNE,NVALL,ITHVAR,FLAGABF,NVARABF,
2745 6 NOM_OPT,IGS ,LSUBMODEL,LITHBUFMX )
2746 NVALL = 9
2747 IGS=IGS+1
2748 TITR1(1:5)='vent '
2749 DO J=1,nchartitle-5
2750 TITR1(J+5:J+5) = TITR(J:J)
2751 ENDDO
2752 CALL FRETITL(TITR1,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2753 ITHGRP(1,IGS)=ID_LIMIT%TH
2754 ID_LIMIT%TH = ID_LIMIT%TH + 1
2755 CALL HM_THGRKI_VENT(
2756 1 ITYP ,KEY ,PTR_NOPT_MONV,
2757 3 IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2758 4 NVARMVENT ,VARMVENT ,NVOLU ,VARMVENTG, NVARMVENTG,
2759 5 IVARMVENTG,NSNE ,NVALL ,ITHVAR , FLAGABF ,NVARABF,
2760 6 NOM_OPT ,IGS ,T_MONVOL ,NVARMVENT, LSUBMODEL)
2761 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2762 ENDDO
2763 ENDIF
2764 !-------------------------------------------
2765 ! /TH/PART
2766 !-------------------------------------------
2767 !HM_NTHPART: total number of /TH/PART in 0000.rad
2768 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/part' ,HM_NTHPART)
2769 IF (HM_NTHPART > 0) THEN
2770 CALL HM_OPTION_START(TRIM(THFILE)//'/part')
2771 HM_NTHGRP = HM_NTHGRP + HM_NTHPART
2772 NUMTHPART = NPART + NTHPART
2773 DO I = 1, HM_NTHPART
2774 CALL HM_READ_THGRPA(
2775 . IPART ,LIPART1 ,ITHBUF ,IAD ,IFI ,
2776 . VARPA ,NVARPA ,VARPAG ,NVARPAG ,NUMTHPART ,
2777 . IVARPAG ,PATHID ,TAGP ,IPARTH ,NPARTH ,
2778 . NVPARTH ,LSUBMODEL)
2779 ENDDO
2780 ENDIF
2781 !-------------------------------------------
2782 ! /TH/SUBS
2783 !-------------------------------------------
2784 !HM_NTHSUBS: total number of /TH/SUBS in 0000.rad
2785 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/subs' ,HM_NTHSUBS)
2786 IF (HM_NTHSUBS > 0) THEN
2787 CALL HM_OPTION_START(TRIM(THFILE)//'/subs')
2788 HM_NTHGRP = HM_NTHGRP + HM_NTHSUBS
2789 ! HM_NTHGRP: Number of /TH/ Group in 0000.rad
2790 DO I = 1, HM_NTHSUBS
2791 ! TAGS: Total number of SUBS for the TIME history output
2792 CALL HM_READ_THGRPA_SUB(
2793 . IAD ,IFI ,ITHBUF ,NVARPA ,VARPA ,
2794 . VARPAG ,NVARPAG ,IVARPAG ,PATHID ,
2795 . SUTHID ,TAGS ,SUBSET ,ITHFLAG ,LSUBMODEL)
2796 ENDDO
2797 ENDIF
2798 !-------------------------------------------
2799 ! /TH/QUAD
2800 !-------------------------------------------
2801 TH_STRAIN = 0 ! Time history strain flag for QUADS -> TH output STRAIN TENSOR
2802 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/quad' ,HM_NTHQUAD )
2803 IF (HM_NTHQUAD > 0) THEN
2804 CALL HM_OPTION_START(TRIM(THFILE)//'/quad')
2805 HM_NTHGRP = HM_NTHGRP + HM_NTHQUAD
2806 DO I = 1, HM_NTHQUAD
2807 TITR = ''
2808 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID = ID, OPTION_TITR = TITR)
2809 IGS = IGS+1
2810 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2811 ITHGRP(1,IGS) = ID
2812 NVALL = 40
2813 KEY = 'quad '
2814 CALL HM_READ_THGRNE(
2815 1 2 ,KEY ,ITAB ,ITABM1,IXQ ,
2816 3 NIXQ ,IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2817 4 NVARS ,VARS ,NUMELQ,VARSG ,NVARSG ,
2818 5 IVARSG,NSELQ,NVALL,IMERGE,ITHVAR ,ITHERM_FE,
2819 6 FLAGABF,NVARABF,RFI,LSUBMODEL,MAP_TABLES%IQUADM, NUMELQ)
2820 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2821 ENDDO
2822 ENDIF
2823 !-------------------------------------------
2824 ! /TH/SPHCEL
2825 !-------------------------------------------
2826 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/sphcel' ,HM_NTHSPHCEL )
2827 IF (HM_NTHSPHCEL > 0) THEN
2828 CALL HM_OPTION_START(TRIM(THFILE)//'/sphcel')
2829 HM_NTHGRP = HM_NTHGRP + HM_NTHSPHCEL
2830 DO I = 1, HM_NTHSPHCEL
2831 TITR = ''
2832 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID = ID, OPTION_TITR = TITR)
2833 IGS = IGS+1
2834 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2835 ITHGRP(1,IGS) = ID
2836 NVALL = 41
2837 KEY = 'sphcel '
2838 CALL HM_READ_THGRNE(
2839 1 51 ,KEY ,ITAB ,ITABM1,KXSP ,
2840 2 NISP ,IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
2841 4 NVARSPH,VARSPH,NUMSPH,VARSPG,NVARSPG ,
2842 5 IVARSPG,NSELSP,NVALL,IMERGE,ITHVAR ,ITHERM_FE,
2843 6 FLAGABF,NVARABF,RFI,LSUBMODEL,MAP_TABLES%ISPHM,NUMSPH)
2844 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2845 ENDDO
2846 ENDIF
2847 !-------------------------------------------
2848 ! /TH/CYL_JO
2849 !-------------------------------------------
2850 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/cyl_jo' ,hm_nthcyljo)
2851 IF (hm_nthcyljo > 0) THEN
2852 CALL hm_option_start(trim(thfile)//'/CYL_JO')
2853 hm_nthgrp = hm_nthgrp + hm_nthcyljo
2854 DO i = 1, hm_nthcyljo
2855 titr = ''
2856 CALL hm_option_read_key(lsubmodel, option_id = id, option_titr = titr)
2857 igs = igs+1
2858 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2859 ithgrp(1,igs) = id
2860 nvall = 6
2861 key = 'CYL_JO '
2862 CALL hm_read_thgrki(
2863 1 105 ,key ,ptr_nopt_joint,
2864 2 iad ,ifi ,ithgrp(1,igs) ,ithbuf ,
2865 3 nvarjo ,varjo ,njoint,varjog,nvarjog,
2866 4 ivarjog ,nsne ,nvall ,ithvar ,flagabf ,nvarabf ,
2867 5 nom_opt ,igs ,lsubmodel, lithbufmx)
2868 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2869 ENDDO
2870 ENDIF
2871 !-------------------------------------------
2872 ! /TH/FXBODY
2873 !-------------------------------------------
2874 CALL hm_option_count(trim(thfile)//'/FXBODY' ,hm_nthfxbody)
2875 IF (hm_nthfxbody > 0) THEN
2876 CALL hm_option_start(trim(thfile)//'/FXBODY')
2877 hm_nthgrp = hm_nthgrp + hm_nthfxbody
2878 DO i = 1, hm_nthfxbody
2879 titr = ''
2880 CALL hm_option_read_key(lsubmodel, option_id = id, option_titr = titr)
2881 igs = igs+1
2882 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2883 ithgrp(1,igs) = id
2884 nvall = 0
2885 key = 'FXBODY '
2886 CALL hm_read_thgrki(
2887 1 111 ,key ,ptr_nopt_fxby,
2888 2 iad ,ifi ,ithgrp(1,igs) ,ithbuf ,
2889 3 nvarfx ,varfx ,nfxbody ,varfxg ,nvarfxg,
2890 4 ivarfxg ,nsne ,nvall ,ithvar ,flagabf ,nvarabf ,
2891 5 nom_opt ,igs ,lsubmodel, lithbufmx )
2892 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2893 ENDDO
2894 ENDIF
2895 !-------------------------------------------
2896 ! /TH/FRAME
2897 !-------------------------------------------
2898 CALL hm_option_count(trim(thfile)//'/FRAME' ,hm_nthframe)
2899 IF (hm_nthframe > 0) THEN
2900 CALL hm_option_start(trim(thfile)//'/FRAME')
2901 hm_nthgrp = hm_nthgrp + hm_nthframe
2902 DO i = 1, hm_nthframe
2903 titr = ''
2904 CALL hm_option_read_key(lsubmodel, option_id = id, option_titr = titr)
2905 igs = igs+1
2906 CALL fretitl(titr,ithgrp(nithgr-ltitr+1,igs),ltitr)
2907 ithgrp(1,igs) = id
2908 nvall = 24
2909 key = 'FRAME '
2910 CALL hm_read_thgrki(
2911 1 110 ,key ,ptr_nopt_skw+numskw+1,
2912 2 iad ,ifi ,ithgrp(1,igs) ,ithbuf ,
2913 3 nvarfr ,varfr ,numfram+1 ,varfrg, nvarfrg ,
2914 4 ivarfrg ,nsne ,nvall ,ithvar ,flagabf ,nvarabf ,
2915 5 nom_opt ,igs ,lsubmodel, lithbufmx)
2916 IF(igs/=0) nthwa=max(nthwa,ithgrp(4,igs)*ithgrp(6,igs))
2917 ENDDO
2918 ENDIF
2919 !-------------------------------------------
2920 ! /TH/GAUGE
2921 !-------------------------------------------
2922 !HM_NTHGAUGE: total number of /TH/GAUGE in 0000.rad
2923 CALL hm_option_count(trim(thfile)//'/GAUGE' ,hm_nthgauge)
2924 IF (hm_nthgauge > 0) THEN
2925 CALL hm_option_start(trim(thfile)//'/gauge')
2926 ! HM_NTHGRP: Number of /TH/ Group in 0000.rad
2927 HM_NTHGRP = HM_NTHGRP + HM_NTHGAUGE
2928 DO I = 1, HM_NTHGAUGE
2929 TITR = ''
2930 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID = ID, OPTION_TITR = TITR)
2931 IGS = IGS+1
2932 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2933 ITHGRP(1,IGS) = ID
2934 NVALL = 3
2935 KEY = 'gauge '
2936 CALL HM_READ_THGRKI(
2937 1 113 ,KEY ,PTR_NOPT_GAU,
2938 2 IAD ,IFI ,ITHGRP(1,IGS) ,ITHBUF ,
2939 3 NVARGAU ,VARGAU ,NBGAUGE ,VARGAUG,NVARGAUG ,
2940 4 IVARGAUG ,NSNE ,NVALL ,ITHVAR ,FLAGABF ,NVARABF ,
2941 5 NOM_OPT ,IGS ,LSUBMODEL, LITHBUFMX)
2942 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2943 ENDDO
2944 ENDIF
2945 !-------------------------------------------
2946 ! /TH/SPH_FLOW
2947 !-------------------------------------------
2948 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/sph_flow' ,HM_NTHSPHFLOW)
2949 IF (HM_NTHSPHFLOW > 0) THEN
2950 CALL HM_OPTION_START(TRIM(THFILE)//'/sph_flow')
2951 HM_NTHGRP = HM_NTHGRP + HM_NTHSPHFLOW
2952 DO I = 1, HM_NTHSPHFLOW
2953 TITR = ''
2954 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID = ID, OPTION_TITR = TITR)
2955 IGS = IGS+1
2956 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2957 ITHGRP(1,IGS) = ID
2958 NVALL = 1
2959 KEY = 'sph_flow '
2960 CALL HM_READ_THGRKI(
2961 1 115 ,KEY ,PTR_NOPT_SPHIO,
2962 2 IAD ,IFI ,ITHGRP(1,IGS) ,ITHBUF ,
2963 3 NVARFLOW,VARFLOW,NSPHIO,VARFLOWG,NVARFLOWG ,
2964 5 IVARFLOWG,NSNE ,NVALL ,ITHVAR ,FLAGABF ,NVARABF ,
2965 5 NOM_OPT ,IGS ,LSUBMODEL,LITHBUFMX)
2966C error message if type of inout is not correct
2967 DO J=1,ITHGRP(4,IGS)
2968 INOUT = ITHBUF(ITHGRP(5,IGS)+J-1)
2969 IF (INOUT>0) THEN
2970 IF (ISPHIO(1,INOUT) == 1) THEN
2971 CALL ANCMSG(MSGID=1111,
2972 . MSGTYPE=MSGERROR,
2973 . ANMODE=ANINFO_BLIND_1,
2974 . I1=ITHGRP(1,IGS),
2975 . C1=TITR,
2976 . I2=NOM_OPT(1,PTR_NOPT_SPHIO+INOUT))
2977 ENDIF
2978 ENDIF
2979 END DO
2980 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
2981 ENDDO
2982 ENDIF
2983 !-------------------------------------------
2984 ! /TH/NSTRAND
2985 !-------------------------------------------
2986 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/nstrand' ,HM_NTHNSTRAND )
2987 IF (HM_NTHNSTRAND > 0) THEN
2988 CALL HM_OPTION_START(TRIM(THFILE)//'/nstrand')
2989 HM_NTHGRP = HM_NTHGRP + HM_NTHNSTRAND
2990 DO I = 1, HM_NTHNSTRAND
2991 TITR = ''
2992 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
2993 IGS = IGS+1
2994 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
2995 ITHGRP(1,IGS) = ID
2996 NVALL = 4
2997 ITYP = 100
2998 KEY = 'nstrand '
2999 CALL HM_READ_THGRNS(
3000 1 100 ,KEY ,ITAB ,ITABM1,KXX ,
3001 3 IXX ,IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
3002 4 NVARNS ,VARNS ,VARNSG ,NVARNSG ,IVARNSG,
3003 5 NSELR ,IVNS2R,NVALL ,ID,TITR,
3004 6 ITHVAR,FLAGABF,NVARABF,LSUBMODEL)
3005 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
3006 ENDDO
3007 ENDIF
3008 !-------------------------------------------
3009 ! /TH/SURF
3010 !-------------------------------------------
3011 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/surf' ,HM_NTHSURF )
3012 IF (HM_NTHSURF> 0) THEN
3013 CALL HM_OPTION_START(TRIM(THFILE)//'/surf')
3014 HM_NTHGRP = HM_NTHGRP + HM_NTHSURF
3015 DO I = 1, HM_NTHSURF
3016 TITR = ''
3017 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
3018 IGS = IGS+1
3019 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
3020 ITHGRP(1,IGS) = ID
3021 NVALL = 1
3022 ITYP = 116
3023 KEY = 'surf '
3024 CALL HM_READ_THGRSURF(
3025 1 ITYP ,KEY ,
3026 3 IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
3027 4 NVARSURF,VARSURF,NSURF,VARSURFG,NVARSURFG ,
3028 5 IVARSURFG,NSNE,NVALL,ITHVAR,FLAGABF,NVARABF,
3029 6 IGRSURF,IGS,LSUBMODEL)
3030 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
3031 ENDDO
3032 ENDIF
3033 !-------------------------------------------
3034 ! /TH/TRIA
3035 !-------------------------------------------
3036 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/tria' ,HM_NTHTRIA )
3037 IF (HM_NTHTRIA > 0) THEN
3038 CALL HM_OPTION_START(TRIM(THFILE)//'/tria')
3039 HM_NTHGRP = HM_NTHGRP + HM_NTHTRIA
3040 DO I = 1, HM_NTHTRIA
3041 TITR = ''
3042 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
3043 IGS = IGS+1
3044 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
3045 ITHGRP(1,IGS) = ID
3046 NVALL = 40
3047 ITYP = 117
3048 KEY = 'tria '
3049 CALL HM_READ_THGRNE(
3050 1 ITYP ,KEY ,ITAB ,ITABM1,IXTG ,
3051 3 NIXTG ,IAD ,IFI ,ITHGRP(1,IGS),ITHBUF ,
3052 4 NVARS ,VARS ,NUMELTG,VARSG ,NVARSG ,
3053 5 IVARSG,NSELTG,NVALL,IMERGE,ITHVAR ,ITHERM_FE,
3054 6 FLAGABF,NVARABF,RFI,LSUBMODEL,MAP_TABLES%ITRIAM, NUMELTG)
3055 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
3056 ENDDO
3057 ENDIF
3058 !-------------------------------------------
3059 ! /TH/SLIPRING
3060 !-------------------------------------------
3061 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/slipring' ,HM_NTHSLIPRING)
3062 IF (HM_NTHSLIPRING > 0) THEN
3063 CALL HM_OPTION_START(TRIM(THFILE)//'/slipring')
3064 HM_NTHGRP = HM_NTHGRP + HM_NTHSLIPRING
3065 DO I = 1, HM_NTHSLIPRING
3066 TITR = ''
3067 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID = ID, OPTION_TITR = TITR)
3068 IGS = IGS+1
3069 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
3070 ITHGRP(1,IGS) = ID
3071 NVALL = 6
3072 ITYP = 118
3073 KEY = 'slipring '
3074 CALL HM_READ_THGRKI(
3075 1 ITYP ,KEY ,PTR_NOPT_SLIPRING,
3076 2 IAD ,IFI ,ITHGRP(1,IGS) ,ITHBUF ,
3077 3 NVARSLIP ,VARSLIP ,NSLIPRING ,VARSLIPG,NVARSLIPG ,
3078 4 IVARSLIPG,NSNE ,NVALL ,ITHVAR ,FLAGABF ,NVARABF ,
3079 5 NOM_OPT ,IGS ,LSUBMODEL, LITHBUFMX)
3080 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
3081 ENDDO
3082 ENDIF
3083 !-------------------------------------------
3084 ! /TH/RETRACTOR
3085 !-------------------------------------------
3086 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/retractor' ,HM_NTHRETRACTOR)
3087 IF (HM_NTHRETRACTOR > 0) THEN
3088 CALL HM_OPTION_START(TRIM(THFILE)//'/retractor')
3089 HM_NTHGRP = HM_NTHGRP + HM_NTHRETRACTOR
3090 DO I = 1, HM_NTHRETRACTOR
3091 TITR = ''
3092 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID = ID, OPTION_TITR = TITR)
3093 IGS = IGS+1
3094 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
3095 ITHGRP(1,IGS) = ID
3096 NVALL = 3
3097 ITYP = 119
3098 KEY = 'retractor '
3099 CALL HM_READ_THGRKI(
3100 1 ITYP ,KEY ,PTR_NOPT_RETRACTOR,
3101 2 IAD ,IFI ,ITHGRP(1,IGS) ,ITHBUF ,
3102 3 NVARRET ,VARRET ,NRETRACTOR ,VARRETG,NVARRETG ,
3103 4 IVARRETG ,NSNE ,NVALL ,ITHVAR ,FLAGABF ,NVARABF ,
3104 5 NOM_OPT ,IGS ,LSUBMODEL, LITHBUFMX)
3105 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
3106 ENDDO
3107 ENDIF
3108 !-------------------------------------------
3109 ! /TH/SENSOR
3110 !-------------------------------------------
3111 CALL HM_OPTION_COUNT(TRIM(THFILE)//'/sensor' ,HM_NTHSENS )
3112 IF (HM_NTHSENS > 0) THEN
3113 CALL HM_OPTION_START(TRIM(THFILE)//'/sensor')
3114 HM_NTHGRP = HM_NTHGRP + HM_NTHSENS
3115 DO I = 1, HM_NTHSENS
3116 TITR = ''
3117 CALL HM_OPTION_READ_KEY(LSUBMODEL, OPTION_ID=ID, OPTION_TITR=TITR)
3118 IGS = IGS+1
3119 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
3120 ITHGRP(1,IGS) = ID
3121 NVALL = 1
3122 ITYP = 120
3123 KEY = 'sensor '
3124 CALL HM_READ_THGRSENS(SENSORS,
3125 . ITYP ,KEY ,IGS ,LITHBUFMX ,ITHBUF ,
3126 . IAD ,IFI ,ITHGRP(1,IGS),ITHVAR,NVALL ,
3127 . NVARSENS ,NVARSENSG,VARSENS ,VARSENSG ,IVARSENSG ,
3128 . NSNE ,NVARABF ,LSUBMODEL )
3129c
3130 IF(IGS/=0) NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
3131 ENDDO
3132 ENDIF
3133 !-------------------------------------------
3134 ! /TH/CHECKSUM ( activated automaticaly if /CHECKSUM is used )
3135 !-------------------------------------------
3136 HM_NTHCHECKSUM=CHECKSUM%checksum_count
3137 IF (HM_NTHCHECKSUM > 0) THEN
3138 HM_NTHGRP = HM_NTHGRP + 1
3139 IGS = IGS+1
3140 CHECKSUM_TITLE = ''
3141 DECK_CHECKSUM = ''
3142
3143 TITR = 'checksum '
3144 CALL FRETITL(TITR,ITHGRP(NITHGR-LTITR+1,IGS),LTITR)
3145
3146
3147 ITHGRP(1,IGS) = 100000
3148 NVALL = 1
3149 ITYP = 121
3150 KEY = 'checksum '
3151 CALL HM_READ_THCHECKSUM(
3152 1 121 ,KEY ,IAD ,IFI ,NSNE ,
3153 2 ITHGRP(1:,IGS) ,ITHBUF ,NVARCHECKSUM ,VARCHECKSUM ,LITHBUFMX,
3154 3 CHECKSUM ,NSUBDOM ,NITHGR ,LTITR ,IPRI )
3155c
3156 NTHWA=MAX(NTHWA,ITHGRP(4,IGS)*ITHGRP(6,IGS))
3157 ENDIF
3158C
3159 DO N=1,NPART
3160 IF(IPART(8,N)>0) NSMAT=NSMAT+1
3161 ENDDO
3162C=======================================================================
3163C Tag des groupes avec des parts dans les TH PART et SUBSET non DEF
3164C=======================================================================
3165 DO NG=1,NGROUP
3166 NEL=IPARG(2,NG)
3167 NFT=IPARG(3,NG)
3168 ITY=IPARG(5,NG)
3169 IF (ITY == 1) THEN
3170 IPARTEL_P => IPARTS ! 3D solids
3171 ELSEIF (ITY == 2) THEN
3172 IPARTEL_P => IPARTQ ! Quads
3173 ELSEIF (ITY == 3) THEN
3174 IPARTEL_P => IPARTC ! Shells
3175 ELSEIF (ITY == 4) THEN
3176 IPARTEL_P => IPARTT ! Truss
3177 ELSEIF (ITY == 5) THEN
3178 IPARTEL_P => IPARTP ! Beams
3179 ELSEIF (ITY == 6) THEN
3180 IPARTEL_P => IPARTR ! Springs
3181.OR. ELSEIF (ITY == 7 ITY == 8) THEN
3182 IPARTEL_P => IPARTG ! SH3N or Tria
3183 ELSEIF (ITY == 51) THEN
3184 IPARTEL_P => IPARTSP ! Sph
3185 ELSEIF (ITY == 100) THEN
3186 IPARTEL_P => IPARTX ! Xelems
3187 ELSEIF (ITY == 101) THEN
3188 IPARTEL_P => IPARTIG3D ! Ig3d
3189 ENDIF
3190 DO J=NFT+1,NFT+NEL
3191 IF(IPARG(80,NG) == 1) EXIT
3192 NUMPART = IPARTEL_P(J)
3193 IF(IPARTH(NVPARTH,NUMPART)>0) THEN
3194 NVAR =IPARTH(NVPARTH,NUMPART)
3195 IAD =IPARTH(NVPARTH+1,NUMPART)
3196 DO N=IAD,IAD+NVAR-1
3197 K=ITHBUF(N)
3198.AND. IF(K>7K/=28) THEN
3199 IPARG(80,NG) = 1
3200 EXIT
3201 ENDIF
3202 ENDDO
3203 ENDIF
3204 ENDDO
3205 ENDDO
3206C=======================================================================
3207C print groupe de PART et SUBSET
3208C=======================================================================
3209 CALL THPRIN(
3210 1 'part ',IPART ,LIPART1 ,4 ,
3211 2 NVPARTH ,ITHBUF ,VARPA ,NPART + NTHPART )
3212 CALL THPRIN_SUB(
3213 1 'subset ',SUBSET ,ITHFLAG ,1 ,
3214 2 NVSUBTH ,ITHBUF ,VARPA ,NSUBS )
3215C
3216C LES GROUPES TH PART ET SUBSET NE SONT PLUS UTILISES
3217C
3218 NTHGRP2=IGS
3219C------------------------------------------------
3220C Recherce des parts ou subsets dans plusieurs TH
3221C------------------------------------------------
3222 DO I=1,TAGP
3223 DO J=1, TAGP
3224.AND. IF((I < J)(PATHID(I) == PATHID(J)))THEN
3225 CALL ANCMSG(MSGID=553,
3226 . MSGTYPE=MSGERROR,
3227 . ANMODE=ANINFO_BLIND_1,
3228 . C1='part',
3229 . I1=PATHID(I))
3230 ENDIF
3231 ENDDO
3232 ENDDO
3233 DO I=1,TAGS
3234 DO J=1, TAGS
3235.AND. IF((I < J)(SUTHID(I) == SUTHID(J)))THEN
3236 CALL ANCMSG(MSGID=553,
3237 . MSGTYPE=MSGERROR,
3238 . ANMODE=ANINFO_BLIND_1,
3239 . C1='subset',
3240 . I1=SUTHID(I))
3241 ENDIF
3242 ENDDO
3243 ENDDO
3244C-------------------------------------
3245C Recherche des ID doubles
3246C-------------------------------------
3247 CALL UDOUBLE(ITHGRP,NITHGR,NTHGRP2,MESS,0,BID)
3248C
3249C-------------------------------------
3250C Write th.nms1 file
3251C-------------------------------------
3252
3253 IF(DUMP_THNMS1_FILE == 1) THEN
3254 CALL WRITE_THNMS1(
3255 1 NVARN1 ,NVARN1A ,NVARN2 ,NVARNPINCH ,NVARS1 ,
3256 2 NVARS2 ,NVARS3 ,NVARS4 ,NVARS5 ,NVARS6 ,
3257 3 NVARS7 ,NVARS8 ,NVARS9 ,NVARS10 ,NVARSNLOC ,
3258 4 NVARP ,NVARR ,NVART ,NVARNS ,NVARSPH ,
3259 5 NVARIN ,NVARRW ,NVARRB ,NVARMV ,NVARSE ,
3260 5 NVARAC ,NVARJO ,NVARMVENT ,NVARPA ,NVARFX ,
3261 6 NVARGAU ,NVARFR ,NVARSLIP ,NVARRET ,NVARCLUS ,
3262 7 NVARFLOW ,NVARSURF ,NVARC ,NVARSENS ,
3263 8 VARN1_TITLE ,VARN1A_TITLE ,VARN2_TITLE ,
3264 9 VARNPINCH_TITLE,VARP_TITLE ,VARR_TITLE ,VART_TITLE ,
3265 A VARS1_TITLE ,VARS2_TITLE ,VARS3_TITLE ,VARS4_TITLE ,VARS5_TITLE ,
3266 B VARS6_TITLE ,VARS7_TITLE ,VARS8_TITLE ,VARS9_TITLE ,VARSNLOC_TITLE ,
3267 C VARC_TITLE ,VARS10_TITLE ,
3268 D VARNS_TITLE ,VARSPH_TITLE ,VARIN_TITLE ,
3269 E VARRW_TITLE ,VARRB_TITLE ,VARMV_TITLE ,VARSE_TITLE ,VARAC_TITLE ,
3270 F VARJO_TITLE ,VARMVENT_TITLE,VARPA_TITLE ,VARFX_TITLE ,VARGAU_TITLE ,
3271 G VARFR_TITLE ,VARSLIP_TITLE ,VARRET_TITLE ,VARCLUS_TITLE ,VARFLOW_TITLE ,
3272 H VARSURF_TITLE ,VARSENS_TITLE ,
3273 I VARN1 ,VARN1A ,VARN2 ,VARNPINCH ,
3274 J VARP ,VARR ,VART ,VARS1 ,VARS2 ,
3275 K VARS3 ,VARS4 ,VARS5 ,VARS6 ,VARS7 ,
3276 L VARS8 ,VARS9 ,VARS10 ,VARSNLOC ,
3277 M VARC ,
3278 N VARNS ,VARSPH ,VARIN ,
3279 O VARRW ,VARRB ,VARMV ,VARSE ,VARAC ,
3280 P VARJO ,VARMVENT ,VARPA ,VARFX ,VARGAU ,
3281 Q VARFR ,VARSLIP ,VARRET ,VARCLUS ,VARFLOW ,
3282 R VARSURF ,VARSENS )
3283 ENDIF
3284C-------------------------------------
3285
3286
3287
3288 ! Deallocate Variables
3289 DEALLOCATE(VARS6)
3290 DEALLOCATE(VARS7)
3291 DEALLOCATE(VARSG2)
3292 DEALLOCATE(VARSG3)
3293 DEALLOCATE(VARS)
3294 DEALLOCATE(VARC)
3295 DEALLOCATE(VARSG)
3296 DEALLOCATE(IVARSG)
3297
3298 DEALLOCATE(VARN1_TITLE,VARN1A_TITLE,VARN2_TITLE,VARNPINCH_TITLE,
3299 . VARP_TITLE,VARR_TITLE,VART_TITLE,
3300 . VARS1_TITLE,VARS2_TITLE,VARS3_TITLE,VARS4_TITLE,
3301 . VARS5_TITLE,VARS6_TITLE,VARS7_TITLE,VARS8_TITLE,VARS9_TITLE,
3302 . VARC_TITLE,VARS10_TITLE,
3303 . VARNS_TITLE,VARSPH_TITLE,
3304 . VARIN_TITLE,VARRW_TITLE,VARRB_TITLE,
3305 . VARMV_TITLE,VARSE_TITLE,VARAC_TITLE,
3306 . VARJO_TITLE,VARMVENT_TITLE,
3307 . VARPA_TITLE,VARFX_TITLE,
3308 . VARGAU_TITLE,VARFR_TITLE,
3309 . VARSLIP_TITLE,VARRET_TITLE,
3310 . VARCLUS_TITLE,VARFLOW_TITLE,
3311 . VARSURF_TITLE,VARSENS_TITLE)
3312
3313 DEALLOCATE(VARN)
3314
3315
3316 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:251
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_option_read_mod::hm_option_read_key
Definition hm_option_read_mod.F:410
max
#define max(a, b)
Definition macros.h:21
id
initmumps id
Definition multiplerhs_example.m:4
groupdef_mod
Definition groupdef_mod.F:662
hm_option_read_mod
Definition hm_option_read_mod.F:408
mapping_option_mod
Definition dichotomy_mod.F:36
message_mod
Definition message_mod.F:1249
monvol_struct_mod
Definition monvol_struct_mod.F:66
names_and_titles_mod
Definition names_and_titles_mod.F:997
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1003
names_and_titles_mod::ncharline
integer, parameter ncharline
Definition names_and_titles_mod.F:1001
submodel_mod
Definition submodel_mod.F:712
fretitl
subroutine fretitl(titr, iasc, l)
Definition freform.F:620