ok.so heres the full subroutine progamme:
SUBROUTINE PASSIN(X,Y,Z,E,PR,WD,NDF,NSN,NMP,NEQ)
IMPLICIT REAL*8(A-H,O-Z)
COMMON/PRECI/ITWO
COMMON/TAPES/ISTRES,NDARAY,IPR
DIMENSION X(NSN),Y(NSN),Z(NSN),E(NMP),PR(NMP),WD(NMP)
DIMENSION NDF(6,NSN),JF(6)
DIMENSION NCYLT(1000)
FLOAT(I)=DBLE(I)
MN=0
100 READ(5,33)NN,(JF(I),I=1,6),X(NN),Y(NN),Z(NN),NI,NCYL
NCYLT(NN)=NCYL
N=MN+NI
MN=MN+1
110 DO 120 I=1,6
NDF(I,NN)=JF(I)
120 CONTINUE
IF (NI.EQ.0) GO TO 130
IF (NN-MN)130,125,140
125 CONTINUE
IF(NSN-NN)170,170,100
130 MN=NN
GO TO 125
C
C AUTOMATIC GENERATION OF NODAL DATA
C
140 NX=(NN-N+NI)/NI
XD=(X(NN)-X(N-NI))/FLOAT(NX)
YD=(Y(NN)-Y(N-NI))/FLOAT(NX)
ZD=(Z(NN)-Z(N-NI))/FLOAT(NX)
MN=NN
150 X(N)=X(N-NI)+XD
Y(N)=Y(N-NI)+YD
Z(N)=Z(N-NI)+ZD
NCYLT(N)=NCYLT(NN)
DO 160 I=1,6
NDF(I,N)=JF(I)
160 CONTINUE
N=N+NI
IF(N.LT.NN) GO TO 150
IF(NSN-NN)170,170,100
170 CONTINUE
DO 180 I=1,NSN
IF(NCYLT(I).EQ.0) GO TO 180
THETA=Y(I)*3.14159/180.0
Y(I)=X(I)*DSIN(THETA)
Y(I)=X(I)*DCOS(THETA)
180 CONTINUE
IF(IPR.EQ.0)WRITE(6,44)(I,(NDF(J,I),J=1,6),X(I),Y(I),Z(I),I=1,NSN)
C CONVERT '0' & '1' OF 'NDF' ARRAY TO EQUATION NUMBERS & '0'S.
NEQ=0
DO 30 N=1,NSN
DO 30 I=1,6
IF(NDF(I,N))10,20,10
20 NEQ=NEQ+1
NDF(I,N)=NEQ
GO TO 30
10 NDF(I,N)=0
30 CONTINUE
IF(IPR.EQ.0)WRITE(6,77)(I,(NDF(J,I),J=1,6),I=1,NSN)
C READ & PRINT MATERIAL PROPERTIES
DO 1 J=1,NMP
1 READ(5,55)I,E(I),PR(I),WD(I)
WRITE(6,66)(I,E(I),PR(I),WD(I),I=1,NMP)
[color=red:6a7731d8f8]& thanx for your valuable advice,[/color:6a7731d8f8]
& if u also find some problem with this subroutine plz let me know..
i will also now try to eliminate the error shown to me.