Listing of file='PLOT3D.FOR;01' on disk='vmedia/chip_61-sector.ccvf'
PROGRAM PLOT3D
DIMENSION P(3,18),G(3,18),T(3,3),I1(3),I2(3),
1 DEG(3),TI(3,3)
LOGICAL CLEAR,COLOR,P2,P242,P255
DATA P/0.0,0.0,0.0,30.0,0.0,0.0,30.0,30.0,0.0,
1 30.0,30.0,30.0,30.0,30.0,0.0,0.0,30.0,0.0,
2 0.0,0.0,0.0,0.0,0.0,30.0,0.0,30.0,30.0,
3 30.0,30.0,30.0,0.0,30.0,30.0,0.0,30.0,0.0,
4 0.0,0.0,0.0,0.0,0.0,30.0,30.0,0.0,30.0,
5 30.0,30.0,30.0,30.0,0.0,30.0,30.0,0.0,0.0/
6 I1,I2/1,1,2,2,3,3/
7 TI/1.0,0.0,0.0,0.0,1.0,0.0,0.0,0.0,1.0/
8 IB,NOPTS,NS,ANGLE,DX,DY/16,18,1,10.0,60.0,60.0/
9 P2,P242,P255,CLEAR/2,242,255,12/
CALL PLTOBJ (NS,NOPTS,P,3,18,CLEAR,DX,DY,P2,P242,P255)
1 DO 30 IR=1,3
DO 30 JR=1,6
COLOR=IB+JR
WRITE(3)COLOR
CALL INIMAT (T,3,3,TI,3,3)
CALL ROTMAT (T,3,3,I1,3,I2,3,IR,ANGLE)
CALL MLTMAT (NOPTS,P,3,18,T,3,3)
CALL PLTOBJ (NS,NOPTS,P,3,18,CLEAR,DX,DY,P2,P242,P255)
30 CONTINUE
GO TO 1
END
SUBROUTINE PLTOBJ(NDS,NODPTS,DP,M1,M2,DCLEAR,DDX,DDY,
1 PD2,PD242,PD255)
DIMENSION DP(M1,M2)
LOGICAL LX,LY,PD2,PD242,PD255,DCLEAR
WRITE(3)DCLEAR
LX=DP(1,NDS)+DDX
LY=DP(2,NDS)+DDY
WRITE(3)PD2,LX,LY,PD242
NSP=NDS+1
DO 10 I=NSP,NODPTS
LX=DP(1,I)+DDX
LY=DP(2,I)+DDY
WRITE(3)LX,LY
10 CONTINUE
WRITE(3)PD255
RETURN
END
SUBROUTINE ROTMAT (DT,N3,N4,ID1,N5,ID2,N6,II,DANGLE)
DIMENSION DT(N3,N4),ID1(N5),ID2(N6)
RADIAN=(3.14159*DANGLE)/180.0
K=ID1(II)
J=ID2(II)
DT(K,K)=COS(RADIAN)
DT(K,J)=-SIN(RADIAN)
DT(J,K)=SIN(RADIAN)
DT(J,J)=COS(RADIAN)
RETURN
END
SUBROUTINE INIMAT (DT,M1,M2,DTI,M3,M4)
DIMENSION DT(M1,M2),DTI(M3,M4)
DO 10 I=1,M1
DO 10 J=1,M4
10 DT(I,J)=DTI(I,J)
RETURN
END
SUBROUTINE MLTMAT (NODPTS,DP,M1,M2,DT,M3,M4)
DIMENSION DP(M1,M2),DT(M3,M4)
DO 10 I=1,NODPTS
PX=DP(1,I)
PY=DP(2,I)
PZ=DP(3,I)
DP(1,I)=PX*DT(1,1)+PY*DT(1,2)+PZ*DT(1,3)
DP(2,I)=PX*DT(2,1)+PY*DT(2,2)+PZ*DT(2,3)
DP(3,I)=PX*DT(3,1)+PY*DT(3,2)+PZ*DT(3,3)
10 CONTINUE
RETURN
END
CONTINU