ccx_2.18/src/elements.f

32       character*8 lakon(*),label  local
52       label='        '
136         read(textpart(i)(6:13),'(a8)') label
140         if((label(1:2).eq.'DC').and.(label(1:7).ne.'DCOUP3D')) then
141           label(1:7)=label(2:8)
142           label(8:8)=' '
146         elseif(label(1:3).eq.'B21') then
147           label='B31     '
148         elseif(label(1:4).eq.'T2D2') then
149           label='T3D2    '
155         if((label.eq.'C3D20   ').or.
156      &       (label.eq.'CPE8    ').or.
157      &       (label.eq.'CPS8    ').or.
158      &       (label.eq.'CAX8    ').or.
159      &       (label.eq.'M3D8    ').or.
160      &       (label.eq.'S8      ').or.
161      &       (label.eq.'T3D3    ').or.
162      &       (label.eq.'B32     ').or.
167      &       (label.eq.'C3D20R  ').or.
168      &       (label.eq.'CPE8R   ').or.
169      &       (label.eq.'CPS8R   ').or.
170      &       (label.eq.'CAX8R   ').or.
171      &       (label.eq.'M3D8R   ').or.
172      &       (label.eq.'S8R     ').or.
173      &       (label.eq.'B32R    ').or.
178      &       (label.eq.'C3D8    ').or.
179      &       (label.eq.'CPE4    ').or.
180      &       (label.eq.'CPS4    ').or.
181      &       (label.eq.'CAX4    ').or.
182      &       (label.eq.'M3D4    ').or.
183      &       (label.eq.'S4      ').or.
184      &       (label.eq.'T3D2    ').or.
185      &       (label.eq.'B31     ').or.
190      &       (label.eq.'C3D8R   ').or.
191      &       (label.eq.'CPE4R   ').or.
192      &       (label.eq.'CPS4R   ').or.
193      &       (label.eq.'CAX4R   ').or.
194      &       (label.eq.'M3D4R   ').or.
195      &       (label.eq.'S4R     ').or.
196      &       (label.eq.'B31R    ').or.
201      &       (label.eq.'C3D8I   ').or.
206      &       (label.eq.'C3D10   ').or.
207      &       (label.eq.'C3D10T  ').or.
211      &       (label.eq.'C3D4    ').or.
216      &       (label.eq.'C3D15   ').or.
217      &       (label.eq.'CPE6    ').or.
218      &       (label.eq.'CPS6    ').or.
219      &       (label.eq.'CAX6    ').or.
220      &       (label.eq.'M3D6    ').or.
221      &       (label.eq.'S6      ').or.
226      &       (label.eq.'C3D6    ').or.
227      &       (label.eq.'CPE3    ').or.
228      &       (label.eq.'CPS3    ').or.
229      &       (label.eq.'CAX3    ').or.
230      &       (label.eq.'M3D3    ').or.
231      &       (label.eq.'S3      ').or.
235      &       (label.eq.'GAPUNI  ').or.
239      &       (label.eq.'SPRINGA ').or.
240      &       (label.eq.'SPRING1 ').or.
241      &       (label.eq.'SPRING2 ').or.
245      &       (label.eq.'DASHPOTA'))
251       elseif((label.eq.'F3D8    ').or.
252      &       (label.eq.'F3D8R   ').or.
253      &       (label.eq.'F3D4    ').or.
254      &       (label.eq.'F3D6    ')) then
258       elseif(label(1:7).eq.'DCOUP3D') then
262       elseif(label(1:1).eq.'D') then
267       elseif(label(1:4).eq.'MASS') then
271       elseif(label(1:1).eq.'U') then
277         write(*,*) label,' is an unknown element type'
282       if(label(1:3).eq.'CAX') then
284       elseif((label(1:2).eq.'S3').or.
285      &       (label(1:2).eq.'S4').or.
286      &       (label(1:2).eq.'S6').or.
287      &       (label(1:2).eq.'S8').or.
288      &       (label(1:1).eq.'B')) then
309       if(label.eq.'        ') then
322       if(label(1:5).eq.'C3D8I') then
325       elseif(label(4:5).eq.'20') then
329       elseif((label(1:4).eq.'CPE8').or.(label(1:4).eq.'CPS8').or.
330      &       (label(1:4).eq.'CAX8').or.(label(1:2).eq.'S8').or.
331      &       (label(1:4).eq.'M3D8')) then
334       elseif((label(1:4).eq.'CPE6').or.(label(1:4).eq.'CPS6').or.
335      &       (label(1:4).eq.'CAX6').or.(label(1:2).eq.'S6').or.
336      &       (label(1:4).eq.'M3D6')) then
339       elseif((label(1:3).eq.'B32').or.
340      &       (label(1:4).eq.'T3D3')) then
343       elseif((label(1:4).eq.'B31 ').or.
344      &       (label(1:4).eq.'T3D2')) then
348       elseif(label(1:4).eq.'B31R') then
351       elseif(label(4:4).eq.'8') then
354       elseif(label(1:3).eq.'S4 ') then
358       elseif((label(1:4).eq.'CPE4').or.(label(1:4).eq.'CPS4').or.
359      &       (label(1:4).eq.'CAX4').or.(label(1:2).eq.'S4').or.
360      &       (label(1:4).eq.'M3D4')) then
363       elseif(label(4:5).eq.'10') then
366       elseif(label(4:4).eq.'4') then
369       elseif(label(4:5).eq.'15') then
372       elseif(label(4:4).eq.'6') then
375       elseif((label(1:4).eq.'CPE3').or.(label(1:4).eq.'CPS3').or.
376      &       (label(1:4).eq.'CAX3').or.(label(1:2).eq.'S3').or.
377      &       (label(1:4).eq.'M3D3')) then
380       elseif(label(1:8).eq.'DASHPOTA') then
381         label='EDSHPTA1'
384       elseif(label(1:7).eq.'DCOUP3D') then
387       elseif(label(1:1).eq.'D') then
390       elseif(label(1:1).eq.'G') then
391         label='ESPGAPA1'
394       elseif(label(1:7).eq.'SPRINGA') then
395         label='ESPRNGA1'
398       elseif(label(1:7).eq.'SPRING1') then
399         label='ESPRNG10'
402       elseif(label(1:7).eq.'SPRING2') then
403         label='ESPRNG21'
406       elseif(label(1:4).eq.'MASS') then
409       elseif(label(1:1).eq.'U') then
410         number=ichar(label(2:2))*256**3+
411      &       ichar(label(3:3))*256**2+
412      &       ichar(label(4:4))*256+
413      &       ichar(label(5:5))
419         write(label(6:6),'(a1)') char(intpoints)
420         write(label(7:7),'(a1)') char(ndof)
421         write(label(8:8),'(a1)') char(nope)
465         lakon(i)=label
524         if((label(1:2).eq.'CP').or.(label(1:2).eq.'CA')) then