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

Removed debug statements

parent 4f1b99d5
No related branches found
No related tags found
No related merge requests found
...@@ -76,7 +76,6 @@ call mpi_comm_size (mpi_comm_world,nproc,ierr) ...@@ -76,7 +76,6 @@ call mpi_comm_size (mpi_comm_world,nproc,ierr)
call mpi_comm_rank (mpi_comm_world,iproc,ierr) call mpi_comm_rank (mpi_comm_world,iproc,ierr)
if (iproc.eq.0) then if (iproc.eq.0) then
print *,'Entered write_global_output'
call int_to_char (cistep,4,istep) call int_to_char (cistep,4,istep)
select case (outputtype) select case (outputtype)
case ('debug') case ('debug')
...@@ -88,7 +87,6 @@ if (iproc.eq.0) then ...@@ -88,7 +87,6 @@ if (iproc.eq.0) then
case default case default
call stop_run ('pb with argument in write_global_output$') call stop_run ('pb with argument in write_global_output$')
end select end select
print *,'Case chosen'
! write array lengths ! write array lengths
write (9) osolve%octree(3), & write (9) osolve%octree(3), &
osolve%nnode, & osolve%nnode, &
...@@ -97,7 +95,6 @@ if (iproc.eq.0) then ...@@ -97,7 +95,6 @@ if (iproc.eq.0) then
osolve%nlsf, & osolve%nlsf, &
cl%np, & cl%np, &
current_time current_time
print *,'Array lengths written'
! write info on octree solve nodes (x,y,z,u,v,w,lsf,temp) ! write info on octree solve nodes (x,y,z,u,v,w,lsf,temp)
write (9) (osolve%x(i), & write (9) (osolve%x(i), &
osolve%y(i), & osolve%y(i), &
...@@ -115,7 +112,6 @@ if (iproc.eq.0) then ...@@ -115,7 +112,6 @@ if (iproc.eq.0) then
osolve%kfix((i-1)*3+3), & osolve%kfix((i-1)*3+3), &
osolve%kfixt(i), & osolve%kfixt(i), &
i=1,osolve%nnode) i=1,osolve%nnode)
print *,'osolve nodal values written'
! write icon array ! write icon array
write (9) ((osolve%icon(k,i),k=1,8), & write (9) ((osolve%icon(k,i),k=1,8), &
! Line below uncommented by dwhipp - 12/09 ! Line below uncommented by dwhipp - 12/09
...@@ -128,13 +124,10 @@ if (iproc.eq.0) then ...@@ -128,13 +124,10 @@ if (iproc.eq.0) then
osolve%is_plastic(i), & osolve%is_plastic(i), &
ov%whole_leaf_in_fluid(i), & ov%whole_leaf_in_fluid(i), &
i=1,osolve%nleaves) i=1,osolve%nleaves)
print *,'osolve elemental values written'
! write octree information ! write octree information
write (9) (osolve%octree(i),i=1,osolve%octree(3)) write (9) (osolve%octree(i),i=1,osolve%octree(3))
print *,'octree info written'
! write bad face information ! write bad face information
write (9) ((osolve%iface(k,i),k=1,9),i=1,osolve%nface) write (9) ((osolve%iface(k,i),k=1,9),i=1,osolve%nface)
print *,'bad faces written'
! write void information ! write void information
write (9) (vo%node(i), & write (9) (vo%node(i), &
vo%leaf(i), & vo%leaf(i), &
...@@ -142,17 +135,13 @@ if (iproc.eq.0) then ...@@ -142,17 +135,13 @@ if (iproc.eq.0) then
vo%rtf(i), & vo%rtf(i), &
vo%influid(i), & vo%influid(i), &
i=1,osolve%nnode) i=1,osolve%nnode)
print *,'void written'
! write bad faces information ! write bad faces information
write (9) (vo%face(i),i=1,osolve%nface) write (9) (vo%face(i),i=1,osolve%nface)
print *,'bad faces written'
! write surface information (r,s,x,y,z,xn,yn,zn) ! write surface information (r,s,x,y,z,xn,yn,zn)
do is=1,osolve%nlsf do is=1,osolve%nlsf
print *,'In surface loop'
write (9) surface(is)%nsurface, & write (9) surface(is)%nsurface, &
surface(is)%activation_time, & surface(is)%activation_time, &
surface(is)%nt surface(is)%nt
print *,'Basic info written'
select case (outputtype) select case (outputtype)
case ('debug') case ('debug')
write (9) (surface(is)%r(i), & write (9) (surface(is)%r(i), &
...@@ -178,9 +167,7 @@ if (iproc.eq.0) then ...@@ -178,9 +167,7 @@ if (iproc.eq.0) then
surface(is)%w(i), & surface(is)%w(i), &
i=1,surface(is)%nsurface) i=1,surface(is)%nsurface)
end select end select
print *,'Surface info written'
write (9) (surface(is)%icon(1:3,i),i=1,surface(is)%nt) write (9) (surface(is)%icon(1:3,i),i=1,surface(is)%nt)
print *,'Surface connectivity written'
enddo enddo
! write cloud information ! write cloud information
write (9) (cl%x(i), & write (9) (cl%x(i), &
...@@ -195,21 +182,16 @@ if (iproc.eq.0) then ...@@ -195,21 +182,16 @@ if (iproc.eq.0) then
cl%press(i), & cl%press(i), &
cl%tag(i), & cl%tag(i), &
i=1,cl%np) i=1,cl%np)
print *,'Cloud written'
if (params%isobc) then if (params%isobc) then
! write isostasy basal displacement array - dwhipp 11/09 ! write isostasy basal displacement array - dwhipp 11/09
write (9) 2**params%levelmax_oct write (9) 2**params%levelmax_oct
print *,'isobc size written'
write (9) ((zi%zisodisp(i,j)+surface(surface(is)%nsurface)%sp01,& write (9) ((zi%zisodisp(i,j)+surface(surface(is)%nsurface)%sp01,&
j=1,2**params%levelmax_oct+1),i=1,2**params%levelmax_oct+1) j=1,2**params%levelmax_oct+1),i=1,2**params%levelmax_oct+1)
print *,'zisodisp written'
end if end if
close (9) close (9)
print *,'And we out...'
! do i=1,osolve%nnode ! do i=1,osolve%nnode
! if (abs(osolve%x(i)-0.5d0)<1.d-6 .and. osolve%z(i)<0.08d0) then ! if (abs(osolve%x(i)-0.5d0)<1.d-6 .and. osolve%z(i)<0.08d0) then
! write (1234,'(5f30.15)') osolve%y(i),osolve%z(i),osolve%v(i),osolve%w(i),ov%temporary_nodal_pressure(i) ! write (1234,'(5f30.15)') osolve%y(i),osolve%z(i),osolve%v(i),osolve%w(i),ov%temporary_nodal_pressure(i)
......
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