C MAKELDT.F - TO MAKE A STANDARD LDT FROM A GEOSOFT XYZ FILE CHARACTER*138 BLINE REAL*8 ANUM INTEGER*4 LINENUM C OPEN(UNIT=3,FILE='COMMag.xyz',ACTION='READ') OPEN(UNIT=4,FILE='magall.ldt') C LINENUM=0 C START WITH THIS 10 READ(3,100,END=30) BLINE 100 FORMAT(A138) IF (BLINE(1:1) .EQ. '/') GOTO 10 C SKIP LINES STARTING WITH / IF (BLINE(1:1) .EQ. 'L') THEN READ(BLINE,101) ANUM 101 FORMAT(5X,F7.1) LINENUM=ANUM*10.0 GOTO 10 ENDIF IF (BLINE(1:1) .EQ. 'T') THEN READ(BLINE,101) ANUM LINENUM=ANUM*10.0 LINENUM=LINENUM+90000 GOTO 10 ENDIF C WORK OUT IF LINE OR TIE C IF (BLINE(13:13) .EQ. ' ') GOTO 10 IF (BLINE(13:13) .EQ. ' ') PAUSE ' COORDS' C NO COORDS C ADD SPECIAL CASE FOR DUNDAS MAGS WHERE LINE 13.2 (132) C HAS NULLS FOR EASTING SO NEED TO DELETE THESE AND BOOST C THE LINE NUMBER C IF (LINENUM .EQ. 132 .AND. BLINE(3:3) .NE. '3') THEN C LINENUM=LINENUM+1 C GOTO 10 C ENDIF C IF (BLINE(3:3) .NE. '3') GOTO 10 C INCREMENT LINE NUMBER AND ALL OTHER ERRORS ARE IN THIS LINE WRITE(4,200) LINENUM,BLINE(1:54),BLINE(56:57),BLINE(59:138) 200 FORMAT(I10,A54,A2,A70) GOTO 10 C 30 CONTINUE CLOSE(UNIT=4,STATUS='KEEP') END