C ******************************************************************* C * * C * * C * * C * PROGRAM EVALUE * C * -------------- * C * * C * * C * THIS PROGRAM COMPUTES SEMI-MAJOR AXIS,SEMI-MINOR AXIS AND * C * ORIENTATION OF SEMI-MAJOR AXIS OF THE STRAIN ELLIPSE. * C * * C * * C * JOB CONTROL LANGUAGE * C * -------------------- * C * * C * * C * 1.0 INPUT DATA SETS * C * * C * STATION NUMBERS AND CO-ORDINATES * C * * C * //GO.FT01F001 DD DSN=A.M6068.ADJCRD31.DATA, * C * // UNIT=P3350,VOL=SER=USER11,DISP=SHR * C * * C * ELEMENTS OF SYMMETRIC STRAIN TENSOR MATRIX * C * //GO.FT15F001 DD DSN=A.M6068.TMATRIX2.DATA, * C * // UNIT=P3350,VOL=SER=USER11,DISP=SHR * C * * C * * C * 2.0 OUTPUT DATA SETS * C * * C * AXES AND ORIENTATION OF STRAIN ELLIPSES * C * * C * //GO.FT31F001 DD DSN=A.M6068.TEVALUE2.DATA, * C * // UNIT=DASD,VOL=SER=USER11,DISP=(NEW,CATLG,DELETE), * C * // DCB=(RECFM=FB,LRECL=33,BLKSIZE=3300),SPACE=(TRK,(10,10),RLSE)* C * //GO.SYSIN DD * * C * // * C * * C * * C * INPUT: * C * * C * CARD 1 FORMAT(I4) * C * NS=NUMBER OF STATIONS IN THE NETWORK * C * * C * ALL OTHER INPUT IS DONE FROM FILES * C * * C * N1,X,Y STATION NAMES AND CO-ORDINATES OF POINTS FROM EITHER * C * OF THE TWO ADJUSTMENTS. * C * FORMAT(I8,2F15.4) * C * * C * A= SYMMETRIC STRAIN TENSOR MATRIX IN THE FORM 1ST DIAGONAL * C * ELEMENT,OFF-DIAGONAL ELEMENT AND SECOND DIAGONAL ELEMENT * C * (DIRECT OUTPUT OF PROGRAM 'STRAIN') * C * FORMAT(1X,3F18.10) * C * * C * * C * OUTPUT: * C * D(1) AND D(2) ARE THE TWO AXES OF OF THE STRAIN ELLIPSES, * C * THAT IS, THE EIGEN VALUES OF STRAIN TENSOR MATRIX. * C * PHI=ORIENTATION ANGLE OF THE AXES OF STRAIN ELLIPSE * C * FORMAT(1X,2F16.10,F10.6) * C * (USED IN PROGRAM 'NETPLOT') * C * * C * WK=PERFORMANCE INDEX OF THE SUBROUTINE EIGRS. * C * IER=GIVES INDICATION OF ANY ERROR. * C * * C * * C * REFERENCES: * C * * C * IMS MANUAL * C * * C * * C * WRITTEN BY KHAGENDRA THAPA NOVEMBER,1979. * C * * C * * C ******************************************************************* C C IMPLICIT REAL*8(A-H,O-Z) REAL*4 PHI DIMENSION A(3),D(2),Z(2,2),WK(5),P(2) N=2 IZ=2 IJOB=2 C C READ NUMBER OF STATION NUMBERS C READ(5,126)NS 126 FORMAT(I4) DO 77 I=1,NS C C READ STATION NUMBER AND CO-ORDINATES OBTAINED FROM EITHER OF C THE TWO ADJUSTMENTS C READ(1,100)N1,X,Y C C READ ELEMENTS OD STRAIN TENSOR MATRIX IN UNIT 15(THIS IS C AN OUTPUT OF PROGRAM STRAIN) C READ(15,455)A 455 FORMAT(1X,3F18.10) C C CALL SUBROUTINE EIGRS TO COMPUTE THE EIGEN VALUES AND C EIGENVECTORS OF THE STRAIN TENSOR C CALL EIGRS(A,N,IJOB,D,Z,IZ,WK,IER) 100 FORMAT(I8,2F15.4) P(1)=DABS(D(1)) P(2)=DABS(D(2)) WRITE(6,25)N1 25 FORMAT(7X,'THE EIGENVALUES AND ORIENTATION ANGLE PHI'/, @7X,'FOR STATION NUMBER',I8/) PI=DARCOS(-1.D0) IF(P(1).EQ.0.0.OR.P(2).EQ.0.0)GOTO 99 C C DECIDE WHICH OF THE EIGEN VALUES IS GREATER C IF(P(1).GT.P(2)) GOTO 29 C C COMPUTE ORIENTATION ANGLE FOR THE SEMI-MAJOR AXIS OF C THE STRAIN ELLIPSE C ANGL=-(A(1)-D(2))/A(2) PHI=DATAN(ANGL) GOTO 69 29 ANGL=-(A(1)-D(1))/A(2) PHI=DATAN(ANGL) GOTO 69 99 PHI=0.0 69 CONTINUE D(1)=D(1)*1.0D5 D(2)=D(2)*1.0D5 IF(P(2).GT.P(1))GOTO 177 C C WRITE THE VALUES OF AXES OF STRAIN ELLIPSE AND ORIENTATION ANGLE C ON FILE IN UNIT 31 AS IN INPUT FOR PROGRAM NETPLOT C WRITE(31,140)(D(J),J=1,2),PHI WRITE(6,141)(D(J),J=1,2),PHI GOTO 48 177 WRITE(31,140)D(2),D(1),PHI WRITE(6,141)D(2),D(1),PHI 48 CONTINUE 140 FORMAT(1X,2F16.10,F10.6) WRITE(6,123) 141 FORMAT(7X,2F14.8,3X,F10.6/) 123 FORMAT(7X,'PRINT THE NUMBERS WHICH INDICATE THE RELIABILITY'/, @ 7X,'OF THE COMPUTATION OF EIGENVALUES AND EIGENVECTORS'/) WRITE(6,136)IER,WK(1) 136 FORMAT(7X,I8,F10.6/) WRITE(6,121)N1 121 FORMAT(7X,'THE EIGENVECTORS FOR STATION #',I8/) WRITE(6,145)((Z(J,K),K=1,2),J=1,2) 145 FORMAT(7X,2F12.8,7X,F12.8/) 77 CONTINUE STOP END