Skip to content
Snippets Groups Projects
Commit fa985e6d authored by holger's avatar holger
Browse files

deleted: src/writedat_fluid.f90

parent 52d2b5dc
No related branches found
No related tags found
No related merge requests found
!####################################################################
!> @author Holger Grosshans
subroutine writedat_ufluid_xz
use var
character(70) :: filename,rowfmt
integer :: i,j,l
100 format(es11.4e2)
write(rowfmt,'(a,i4,a)') '(',(imax-imin+2),'(1x,es11.4e2))'
j=int(jj/2)
write(filename,'(a,i3.3,a,i6.6,a)') &
'results/fluid_u_xz_p',myid,'_',nt,'.dat'
open(10,file=filename)
write(10,'(a,es14.6e2,a,x,a)') '# u: t=',t,' s --',version
write(10,'(a11)',advance='no') '# z/x'
write(10,fmt=rowfmt) (xf(i),i=imin,imax+1)
do l=lmin-1,lmax+1
write(10,fmt=100,advance='no') zc(l)
write(10,fmt=rowfmt) (u(i,j,l),i=imin,imax+1)
enddo
close(10)
return
end
!####################################################################
!> @author Holger Grosshans
subroutine writedat_vfluid_xz
use var
character(70) :: filename,rowfmt
integer :: i,j,l
100 format(es11.4e2)
write(rowfmt,'(a,i4,a)') '(',(imax-imin+2),'(1x,es11.4e2))'
j=int(jj/2)
write(filename,'(a,i3.3,a,i6.6,a)') &
'results/fluid_v_xz_p',myid,'_',nt,'.dat'
open(10,file=filename)
write(10,'(a,es14.6e2,a,x,a)') '# v: t=',t,' s --',version
write(10,'(a11)',advance='no') '# z/x'
write(10,fmt=rowfmt) (xc(i),i=imin,imax+1)
do l=lmin-1,lmax+1
write(10,fmt=100,advance='no') zc(l)
write(10,fmt=rowfmt) (v(i,j,l),i=imin,imax+1)
enddo
close(10)
return
end
!####################################################################
!> @author Holger Grosshans
subroutine writedat_wfluid_xz
use var
character(70) :: filename,rowfmt
integer :: i,j,l
100 format(es11.4e2)
write(rowfmt,'(a,i4,a)') '(',(imax-imin+2),'(1x,es11.4e2))'
j=int(jj/2)
write(filename,'(a,i3.3,a,i6.6,a)') &
'results/fluid_w_xz_p',myid,'_',nt,'.dat'
open(10,file=filename)
write(10,'(a,es14.6e2,a,x,a)') '# w: t=',t,' s --',version
write(10,'(a11)',advance='no') '# z/x'
write(10,fmt=rowfmt) (xc(i),i=imin,imax+1)
do l=lmin-1,lmax+1
write(10,fmt=100,advance='no') zf(l)
write(10,fmt=rowfmt) (w(i,j,l),i=imin,imax+1)
enddo
close(10)
return
end
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment