43
44
45
46 USE elbufdef_mod
47 use element_mod , only : nixt
48
49
50
51#include "implicit_f.inc"
52
53
54
55#include "mvsiz_p.inc"
56#include "com04_c.inc"
57
58
59
60 INTEGER NCT(NIXT,*),OFFSET,JFT,JLT,NEL ,IKGEO,
61 . ETAG(*),IDDL(*) ,NDOF(*) ,IADK(*) ,JDIK(*)
63 . pm(*), geo(*), x(*),
64 . ke11(9,mvsiz),ke12(9,mvsiz),ke22(9,mvsiz),off(mvsiz),
65 . k_diag(*) ,k_lt(*)
66
67 TYPE (ELBUF_STRUCT_), TARGET:: ELBUF_TAB
68
69
70
71
72 INTEGER MAT(MVSIZ),NGL(MVSIZ)
74 . al(mvsiz),kx(mvsiz),ex(mvsiz),ey(mvsiz),ez(mvsiz)
75
76 TYPE(G_BUFEL_) ,POINTER :: GBUF
77
78
79
80 gbuf => elbuf_tab%GBUF
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99 CALL tcoork3(jft ,jlt ,x ,nct ,mat ,
100 1 ngl ,al ,ex ,ey ,ez ,
101 2 gbuf%OFF,off )
102 CALL tmat3(jft ,jlt ,pm ,mat ,gbuf%AREA,
103 . al ,kx )
104
105
106
107 CALL r4sumg3(jft ,jlt ,al ,kx ,ex ,
108 1 ey ,ez ,ke11 ,ke12 ,ke22 )
109 IF (ikgeo>0)
CALL tkeg3(jft ,jlt ,
110 1 gbuf%FOR ,al ,ke11 ,ke12 ,ke22 )
112 1 nct ,nel ,iddl ,ndof ,k_diag ,
113 2 k_lt ,iadk ,jdik ,ke11 ,ke12 ,
114 3 ke22 ,off )
115 IF (neig>0)
CALL peoff(
116 1 nixt, jft, jlt, nct , etag,
117 2 off )
118
119 RETURN
subroutine assem_p(nixpl, nd, ixp, nel, iddl, ndof, k_diag, k_lt, iadk, jdik, kc11, kc12, kc22, off)
subroutine peoff(nixpl, jft, jlt, ixp, etag, off)
subroutine r4sumg3(jft, jlt, al, kx, r11, r21, r31, ke11, ke12, ke22)
subroutine tcoork3(jft, jlt, x, ncc, mat, ngl, al, r11, r21, r31, offg, off)
subroutine tkeg3(jft, jlt, fx, al, ke11, ke12, ke22)
subroutine tmat3(jft, jlt, pm, mat, area, al, kx)