33
34
35
36#include "implicit_f.inc"
37
38
39
40#include "com04_c.inc"
41#include "param_c.inc"
42#include "task_c.inc"
43
44
45
47 . x(3,*),rwbuf(nrwlp,*)
48 INTEGER NPRW(*)
49
50
51
52 INTEGER J,N,K,,N3,N4,ITYP,I3000,NXX,NYY,NZZ
54 i3000 = 3000
55 IF (ispmd/=0) GOTO 100
56 DO n=1,nsect+nrwall
60 ENDDO
61
62 k=1
63 DO n=1,nrwall
64 n2=n +nrwall
65 n3=n2+nrwall
66 n4=n3+nrwall
67 ityp= nprw(n4)
68 IF(iabs(ityp)==1.OR.ityp==4)THEN
69 nxx = three1000*rwbuf(1,n)
70 nyy = three1000*rwbuf(2,n)
71 nzz = three1000*rwbuf(3,n)
72 DO j=1,4
76 ENDDO
77 ELSEIF(ityp==2)THEN
78 CALL dnwalc(x,rwbuf(1,n),nprw(n3))
79 ELSEIF(ityp==3)THEN
80 CALL dnwals(x,rwbuf(1,n),nprw(n3))
81 ENDIF
82 k=k+nprw(n)
83 IF(nprw(n4)==-1)k=k+nint(rwbuf(8,n))
84 ENDDO
85
86 100 CONTINUE
87 RETURN
subroutine dnwalc(x, rwl, msr)
subroutine dnwals(x, rwl, msr)
void write_s_c(int *w, int *len)