150
151
152
154
155 IMPLICIT NONE
156
157
158
159#include "my_real.inc"
160
161
162
163 TYPE(STACK_PLY),INTENT(INOUT) :: STACK
164 INTEGER,INTENT(IN) :: IPART_STACK
165 INTEGER,INTENT(IN) :: NPT_STACK
166 INTEGER,INTENT(IN) :: NS_STACK
167
168
169
170 INTEGER IERR
171
172 IF(ipart_stack >0) THEN
173
174
175 ALLOCATE(stack%IGEO(4*npt_stack+2,ns_stack),stat=ierr)
176 IF(ierr/=0) THEN
177 CALL ancmsg(msgid=19,anmode=aninfo,c1=
'(/PROP/STACK/...)')
179 ENDIF
180 stack%IGEO = 0
181
182 ALLOCATE(stack%GEO(6*npt_stack+1,ns_stack),stat=ierr)
183 IF(ierr/=0) THEN
184 CALL ancmsg(msgid=19,anmode=aninfo,c1=
'(/PROP/STACK/...)')
186 ENDIF
187 stack%GEO = 0.0
188
189 ALLOCATE(stack%PM(20,ns_stack),stat=ierr)
190 IF(ierr/=0) THEN
191 CALL ancmsg(msgid=19,anmode=aninfo,c1=
'(/PROP/STACK/...)')
193 ENDIF
194 stack%PM = 0.0
195 ELSE
196 ALLOCATE(stack%IGEO(0,0),stack%GEO(0,0),stack%PM(0,0))
197 ENDIF
198
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)