Skip to content
Snippets Groups Projects
Commit 94a7a151 authored by Dave Whipp's avatar Dave Whipp
Browse files

Added read/echo of lsoutfile

parent 2038eb10
No related branches found
No related tags found
No related merge requests found
...@@ -816,6 +816,10 @@ end if ...@@ -816,6 +816,10 @@ end if
call mpi_bcast(params%nest,1,mpi_logical,0,mpi_comm_world,ierr) call mpi_bcast(params%nest,1,mpi_logical,0,mpi_comm_world,ierr)
if (params%nest) then if (params%nest) then
nest%lsoutfile='OUT/time_0001.bin'
if (iproc.eq.0) call scanfile (params%infile,'lsoutfile',nest%lsoutfile,ires)
call mpi_bcast(nest%lsoutfile,128,mpi_character,0,mpi_comm_world,ierr)
nest%sselemx=1.d0 nest%sselemx=1.d0
if (iproc==0) call scanfile (params%infile,'sselemx',nest%sselemx,ires) if (iproc==0) call scanfile (params%infile,'sselemx',nest%sselemx,ires)
call mpi_bcast(nest%sselemx,1,mpi_double_precision,0,mpi_comm_world,ierr) call mpi_bcast(nest%sselemx,1,mpi_double_precision,0,mpi_comm_world,ierr)
...@@ -1007,6 +1011,7 @@ if (params%debug.gt.0 .and. iproc.eq.0) then ...@@ -1007,6 +1011,7 @@ if (params%debug.gt.0 .and. iproc.eq.0) then
write(*,'(a,l1)') shift//'nest ',params%nest write(*,'(a,l1)') shift//'nest ',params%nest
if (params%nest) then if (params%nest) then
write(*,'(a)') shift//'--- Nest parameters ---' write(*,'(a)') shift//'--- Nest parameters ---'
write(*,'(a,a)') shift//'lsoutfile ',trim(nest%lsoutfile)
write(*,'(a,e11.4)') shift//'sselemx ',nest%sselemx write(*,'(a,e11.4)') shift//'sselemx ',nest%sselemx
write(*,'(a,e11.4)') shift//'sselemy ',nest%sselemy write(*,'(a,e11.4)') shift//'sselemy ',nest%sselemy
write(*,'(a,e11.4)') shift//'sselemz ',nest%sselemz write(*,'(a,e11.4)') shift//'sselemz ',nest%sselemz
...@@ -1176,6 +1181,7 @@ if (params%debug.gt.1) then ...@@ -1176,6 +1181,7 @@ if (params%debug.gt.1) then
write(threadinfo%Logunit,'(a32,l1)') 'nest ',params%nest write(threadinfo%Logunit,'(a32,l1)') 'nest ',params%nest
if (params%nest) then if (params%nest) then
write(threadinfo%Logunit,'(a)') '--- Nest parameters ---' write(threadinfo%Logunit,'(a)') '--- Nest parameters ---'
write(threadinfo%Logunit,'(a32,a)') 'lsoutfile ',trim(nest%lsoutfile)
write(threadinfo%Logunit,'(a32,e11.4)') 'sselemx ',nest%sselemx write(threadinfo%Logunit,'(a32,e11.4)') 'sselemx ',nest%sselemx
write(threadinfo%Logunit,'(a32,e11.4)') 'sselemy ',nest%sselemy write(threadinfo%Logunit,'(a32,e11.4)') 'sselemy ',nest%sselemy
write(threadinfo%Logunit,'(a32,e11.4)') 'sselemz ',nest%sselemz write(threadinfo%Logunit,'(a32,e11.4)') 'sselemz ',nest%sselemz
...@@ -1189,6 +1195,8 @@ endif ...@@ -1189,6 +1195,8 @@ endif
if(iproc.eq.0) call flush(8) if(iproc.eq.0) call flush(8)
if(params%debug.gt.1) call flush(threadinfo%logunit) if(params%debug.gt.1) call flush(threadinfo%logunit)
!if (iproc.eq.0) call system ('rm fort.8')
end subroutine read_input_file end subroutine read_input_file
!------------------------------------------------------------------------------| !------------------------------------------------------------------------------|
!------------------------------------------------------------------------------| !------------------------------------------------------------------------------|
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment