1 SUBROUTINE EQSPC1(X,S,NP,XE,SE,NE,DSEDX) 2C 3C ----COMPUTE S AND DS/DX ON EQUAL SPACED X 4C 5 DIMENSION ROUTID(2) 6 DIMENSION CS(6),X(1),S(1),XE(1),SE(1),DSEDX(1) 7 DIMENSION VAR(4),LGH(4) 8 DATA ROUTID/4HEQSP,4HC1 / 9C 10 IF(NP.EQ.1)GO TO 1020 11 IS=0 12 FNE=NE-1 13 XIN=(X(NP)-X(1))/FNE 14 XE(1)=X(1) 15 SE(1)=S(1) 16 XE(NE)=X(NP) 17 SE(NE)=S(NP) 18 NN=NE-1 19 DO 1000 I=2,NN 20 XE(I)=XE(I-1)+XIN 21 VAR(1)=XE(I) 22 LGH(1)=NP 23 CALL INTERX(1,X,VAR,LGH,S,SE(I),NP,NP,0,0,0,0,0,0,0,0, 24 1 4HSE ,1,ROUTID) 25 1000 CONTINUE 26 DO 1010 I=1,NE 27 CALL TBFUNX(XE(I),DUMM,DSEDX(I),NE,XE,SE,CS,IS,MI,NG,0,0, 28 1 4HDSED,1,ROUTID) 29 1010 CONTINUE 30 RETURN 31 1020 DO 1030 I=1,NE 32 XE(I)=X(1) 33 SE(I)=S(1) 34 1030 DSEDX(I)=0.0 35 RETURN 36 END 37