w_out.f90

自己编的地震图用mapgis拾取后的时深转换程序

real*8 :: x,y,b(1:3),z,v,dx
integer m,a(1:6),n,p
character :: head*8,filename*20,out*23,outdata*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


!---------------参数设置----------------
v=2500
dx=50
!------------------------------------------

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

  p=j/10
  if(j<10) then
  outdata="outdata"//char(j+47)
  else 
  
  outdata=char(p+48)//char(j+48-p*10)
  outdata="outdata"//outdata
  endif

 ! open(3,FILE=outdata,status='unknown')
  
  do k=1,n
    read(1,'(2f)') x,y
   !-------------------------处理-------------
    y=-y*v/2
    x=dx*x

   !---------------------------------------
	write(2,200) x,y
 	!write(3,200) x,y
  enddo
  !close(3)
  
  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(f16.4,',',f16.4)
300  FORMAT(I2,',',f13.6)
CLOSE(1)
close(2)
END