33
34
35
36#include "implicit_f.inc"
37
38
39
40#include "com01_c.inc"
41#include "param_c.inc"
42
43
44
45 INTEGER NODFT,NODLT,ICODT(*),ICODR(*),ISKEW(*)
47 . a(3,*), ar(3,*), ms(*), v(3,*), vr(3,*),skew(lskew,*)
48
49
50
51 INTEGER I, N, K, ISK, LCOD, NINDX,
52 . INDX(1024),L,J
53
54 DO 420 i = nodft,nodlt,1024
55 nindx = 0
56 DO 400 n=i,
min(nodlt,i+1023)
57 IF (icodt(n)/=0)THEN
58 nindx = nindx + 1
59 indx(nindx) = n
60 ENDIF
61 400 CONTINUE
62 CALL bcs1v(nindx,indx,iskew,icodt,
63 . a ,skew,v )
64 420 CONTINUE
65
66 IF(iroddl/=0) THEN
67 DO 440 i = nodft,nodlt,1024
68 nindx = 0
69 DO 430 n=i,
min(nodlt,i+1023)
70 IF (icodr(n)/=0)THEN
71 nindx = nindx + 1
72 indx(nindx) = n
73 ENDIF
74 430 CONTINUE
75 CALL bcs1v(nindx,indx,iskew,icodr,
76 . ar ,skew,vr )
77 440 CONTINUE
78 ENDIF
79
80 RETURN
subroutine bcs1v(nindx, indx, iskew, icodt, a, skew, v)