Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

- same as previous commit #4

Merged
merged 1 commit into from
Mar 26, 2014
Merged

- same as previous commit #4

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
7 changes: 4 additions & 3 deletions SOLVER/seismograms.f90
View file
Open in desktop
Expand Up @@ -444,7 +444,7 @@ subroutine prepare_from_recfile_seis
allocate(reclat(num_rec_glob),reclon(num_rec_glob),recelevation(num_rec_glob),recbury(num_rec_glob))
allocate(receiver_name(num_rec_glob))
open(unit=34,file='STATIONS',iostat=ierror,status='old',action='read',position='rewind')
open(unit=30,file=datapath(1:lfdata)//'/receiver_names.dat')
if (mynum==0) open(unit=30,file=datapath(1:lfdata)//'/receiver_names.dat')

do i=1,num_rec_glob
read(34,*)rec_name(i),rec_network(i),reclat(i),reclon(i),recelevation(i),recbury(i)
Expand All @@ -455,9 +455,10 @@ subroutine prepare_from_recfile_seis
endif
recfile_readth(i) = 90.d0 - reclat(i)
receiver_name(i) = trim(rec_name(i))//'_'//trim(rec_network(i))
write(30,*)trim(receiver_name(i)),recfile_readth(i),recfile_readph(i)
if (mynum==0) write(30,*)trim(receiver_name(i)),recfile_readth(i),recfile_readph(i)
enddo
close(34); close(30)
close(34)
if (mynum==0) close(30)


case default
Expand Down