w.f90

地震数据时间深度转换

real*8 :: x(1:100),y(1:100),b(1:3),z
integer m,a(1:6),n
character :: head*8,filename*20,out*23
read(*,*)filename
out="out"//filename
open(1,FILE=filename,status='old')
open(2,FILE=out,status='unknown')

read(1,'(8A/)') head
read(1,'(I2)') m
!WRITE(2,'(8A/)') head
!write(2,'(I2)') m


do j=1,m-1

  read(1,100)(a(i),i=1,3),(b(i),i=1,3),(a(i),i=4,6)
  !write(2,400)a(1),a(2),a(3),b(1),b(2),b(3),a(4),a(5),a(6)
  read(1,'(I)') n
  !write(2,'(I2)') n

 
  
  do k=1,n
    read(1,'(2f)') x(k),y(k)
   !---------------------------------------
    y(k)=-y(k)



   !---------------------------------------
    !write(2,200) x,y
  enddo
  write(2,'(f13.6$)') (x(k),k=1,n)
  write(2,'(/)')
  write(2,'(f13.6$)') (y(k),k=1,n)
  
  read(1,'(I,f)')l,z
  !write(2,300)l,z

enddo

100  FORMAT(3I,3f,3I)
400  FORMAT(I1,',',I1,',',I3,',',f8.6,',',f9.6,',',f9.6,',',I1,',',I1,',',I1,',')
200  FORMAT(f13.6,',',f13.6)
!300  FORMAT(I2,',',f13.6)
CLOSE(1)
close(2)
END