Skip to content

Commit

Permalink
first commit
Browse files Browse the repository at this point in the history
  • Loading branch information
@tjz2026
tjz2026 committed Jan 9, 2016
1 parent 4d1a459 commit 1543dce
Show file tree
Hide file tree
Showing 25 changed files with 9,954 additions and 0 deletions.
77 changes: 77 additions & 0 deletions 2Dwall_tecplot.f90
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,77 @@
!!!! this program is used for gathering the density profile from each processor
!and write the profile in the format of tecplot .

program main
implicit none
integer :: ii,K_i,K_j,K_k,K,aaa,bbb,jj
real*8 :: dx,dy,dz
real*8,allocatable :: RHO_A(:),RHO_B(:),RHO_C(:)
real*8,allocatable :: R_X(:),R_Y(:),R_Z(:)
integer,parameter :: SIDEx=32,SIDEy=32,local_nx=8,filenumber=4
character(len=30)::aa,bb


allocate(RHO_A(0:local_nx*SIDEy-1))
allocate(RHO_B(0:local_nx*SIDEy-1))
allocate(RHO_C(0:local_nx*SIDEy-1))
allocate(R_X(0:local_nx-1))
allocate(R_Y(0:SIDEy-1))





open(unit=45,file='dfpde_wall_sum.dat',status='new')
WRITE (45,2501)
WRITE (45,2502) SIDEx,SIDEy
close(45)

! open(unit=45,file='Pu_total.dat',status='new')
! WRITE (45,2501)
! WRITE (45,2502) SIDEz,SIDEy,SIDEx
! close(45)
do ii=0,filenumber-1
aaa=ii
write(aa,*) aaa
bbb=1
write(bb,*) bbb

!open(unit=33,file='Wall_sum' // trim(adjustl(aa)) // '.dat',status='old')
open(unit=33,file='Wall_sum' // trim(adjustl(aa)) // '.dat',status='old')
k=0
do k_i=0,Local_nx-1
do k_j=0,SIDEy-1
read(33,*) R_X(K_i),R_Y(K_j),RHO_A(k)
k=k+1
enddo
enddo
close(33)

open(unit=45,file='dfpde_wall_sum.dat',status='old',position='append')

k=0
do k_i=0,Local_nx-1
do k_j=0,SIDEy-1
write(45,2504) R_X(K_i),R_Y(K_j),RHO_A(k)
k=k+1
enddo
enddo
close(45)
enddo


2501 FORMAT ('VARIABLES =',1X,'"X"',1X,'"Y"',1X,'"A"')
2502 FORMAT ('ZONE',1x,'I=',I3,1X,'J=',I3,1X,'F=POINT')
2503 FORMAT (I4,2X,I4,2X,I4,2X,F9.5,2X,F9.5)
2504 FORMAT (F9.5,2X,F9.5,2X,F9.5,2X)

deallocate(RHO_A)
deallocate(RHO_B)
deallocate(RHO_C)
deallocate(R_X)
deallocate(R_Y)

end



125 changes: 125 additions & 0 deletions 3D_tecplot.f90
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,125 @@
!!!! this program is used for gathering the density profile from each processor
!and write the profile in the format of tecplot .

program main
implicit none
integer :: ii,K_i,K_j,K_k,K,aaa
real*8 :: dx,dy,dz
real*8,allocatable :: RHO_A(:),RHO_B(:),RHO_C(:)
real*8,allocatable :: R_X(:),R_Y(:),R_Z(:)
integer,parameter :: SIDEx=32,SIDEy=32,SIDEz=32,local_nx=4,filenumber=8
character(len=30)::aa


allocate(RHO_A(0:local_nx*SIDEy*SIDEz-1))
allocate(RHO_B(0:local_nx*SIDEy*SIDEz-1))
allocate(RHO_C(0:local_nx*SIDEy*SIDEz-1))
allocate(R_X(0:local_nx-1))
allocate(R_Y(0:SIDEy-1))
allocate(R_Z(0:SIDEz-1))





open(unit=45,file='dfpde_total.dat',status='new')
WRITE (45,2501)
WRITE (45,2502) SIDEz,SIDEy,SIDEx
close(45)

! open(unit=45,file='Pu_total.dat',status='new')
! WRITE (45,2501)
! WRITE (45,2502) SIDEz,SIDEy,SIDEx
! close(45)

do ii=0,filenumber-1
aaa=ii
write(aa,*) aaa
open(unit=33,file='RHO_total' //trim(adjustl(aa)) // '.dat',status='old')
k=0
do k_i=0,Local_nx-1
do k_j=0,SIDEy-1
do K_k=0,SIDEz-1
read(33,*) R_X(K_i),R_Y(K_j),R_Z(K_k),RHO_A(k),RHO_B(k),RHO_C(k)
k=k+1
enddo
enddo
enddo
close(33)

open(unit=45,file='dfpde_total.dat',status='old',position='append')

k=0
do k_i=0,Local_nx-1
do k_j=0,SIDEy-1
do K_k=0,SIDEz-1
write(45,2504) R_X(K_i),R_Y(K_j),R_Z(K_k),RHO_A(k),RHO_B(k),1.0-RHO_C(k)
k=k+1
enddo
enddo
enddo
close(45)
enddo

! do ii=0,filenumber-1
! aaa=ii
! write(aa,*) aaa
! open(unit=33,file='Pu' //trim(adjustl(aa)) // '.dat',status='old')
! k=0
! do k_i=0,Local_nx-1
! do k_j=0,SIDEy-1
! do K_k=0,SIDEz-1
! read(33,*) R_X(K_i),R_Y(K_j),R_Z(K_k),RHO_A(k),RHO_B(k),RHO_C(k)
! k=k+1
! enddo
! enddo
! enddo
! close(33)

! open(unit=45,file='Pu_total.dat',status='old',position='append')
!
! k=0
! do k_i=0,Local_nx-1
! do k_j=0,SIDEy-1
! do K_k=0,SIDEz-1
! write(45,2504) R_X(K_i),R_Y(K_j),R_Z(K_k),RHO_A(k),RHO_B(k),RHO_C(k)
! k=k+1
! enddo
! enddo
! enddo
! close(45)







enddo











2501 FORMAT ('VARIABLES =',1X,'"X"',1X,'"Y"',1X,'"Z"',1X,'"A"',1X,'"B"',1X,' "C"')
2502 FORMAT ('ZONE',1x,'I=',I3,1X,'J=',I3,1X,'K=',I3,1X,'F=POINT')
2503 FORMAT (I4,2X,I4,2X,I4,2X,F9.5,2X,F9.5)
2504 FORMAT (F9.5,2X,F9.5,2X,F9.5,2X,F9.5,2X,F9.5,2X,F9.5)

deallocate(RHO_A)
deallocate(RHO_B)
deallocate(RHO_C)
deallocate(R_X)
deallocate(R_Y)
deallocate(R_Z)

end



128 changes: 128 additions & 0 deletions 3Dwall_tecplot.f90
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,128 @@
!!!! this program is used for gathering the density profile from each processor
!and write the profile in the format of tecplot .

program main
implicit none
integer :: ii,K_i,K_j,K_k,K,aaa,bbb,jj
real*8 :: dx,dy,dz
real*8,allocatable :: RHO_A(:),RHO_B(:),RHO_C(:)
real*8,allocatable :: R_X(:),R_Y(:),R_Z(:)
integer,parameter :: SIDEx=64,SIDEy=64,SIDEz=8,local_nx=16,filenumber=4
character(len=30)::aa,bb


allocate(RHO_A(0:local_nx*SIDEy*SIDEz-1))
allocate(RHO_B(0:local_nx*SIDEy*SIDEz-1))
allocate(RHO_C(0:local_nx*SIDEy*SIDEz-1))
allocate(R_X(0:local_nx-1))
allocate(R_Y(0:SIDEy-1))
allocate(R_Z(0:SIDEz-1))





open(unit=45,file='dfpde_wall_sum.dat',status='new')
WRITE (45,2501)
WRITE (45,2502) SIDEx,SIDEy,SIDEz
close(45)

! open(unit=45,file='Pu_total.dat',status='new')
! WRITE (45,2501)
! WRITE (45,2502) SIDEz,SIDEy,SIDEx
! close(45)
do ii=0,filenumber-1
aaa=ii
write(aa,*) aaa
bbb=1
write(bb,*) bbb

!open(unit=33,file='Wall_sum' // trim(adjustl(aa)) // '.dat',status='old')
open(unit=33,file='Wall' //trim(adjustl(bb))// 'on' // trim(adjustl(aa)) // '.dat',status='old')
k=0
do k_i=0,Local_nx-1
do k_j=0,SIDEy-1
do K_k=0,SIDEz-1
read(33,*) R_X(K_i),R_Y(K_j),R_Z(K_k),RHO_A(k)
k=k+1
enddo
enddo
enddo
close(33)

open(unit=45,file='dfpde_wall_sum.dat',status='old',position='append')

k=0
do k_i=0,Local_nx-1
do k_j=0,SIDEy-1
do K_k=0,SIDEz-1
write(45,2504) R_X(K_i),R_Y(K_j),R_Z(K_k),RHO_A(k)
k=k+1
enddo
enddo
enddo
close(45)
enddo

! do ii=0,filenumber-1
! aaa=ii
! write(aa,*) aaa
! open(unit=33,file='Pu' //trim(adjustl(aa)) // '.dat',status='old')
! k=0
! do k_i=0,Local_nx-1
! do k_j=0,SIDEy-1
! do K_k=0,SIDEz-1
! read(33,*) R_X(K_i),R_Y(K_j),R_Z(K_k),RHO_A(k),RHO_B(k),RHO_C(k)
! k=k+1
! enddo
! enddo
! enddo
! close(33)

! open(unit=45,file='Pu_total.dat',status='old',position='append')
!
! k=0
! do k_i=0,Local_nx-1
! do k_j=0,SIDEy-1
! do K_k=0,SIDEz-1
! write(45,2504) R_X(K_i),R_Y(K_j),R_Z(K_k),RHO_A(k),RHO_B(k),RHO_C(k)
! k=k+1
! enddo
! enddo
! enddo
! close(45)







! enddo











2501 FORMAT ('VARIABLES =',1X,'"X"',1X,'"Y"',1X,'"Z"',1X,'"A"')
2502 FORMAT ('ZONE',1x,'I=',I3,1X,'J=',I3,1X,'K=',I3,1X,'F=POINT')
2503 FORMAT (I4,2X,I4,2X,I4,2X,F9.5,2X,F9.5)
2504 FORMAT (F9.5,2X,F9.5,2X,F9.5,2X,F9.5,2X,F9.5,2X,F9.5)

deallocate(RHO_A)
deallocate(RHO_B)
deallocate(RHO_C)
deallocate(R_X)
deallocate(R_Y)
deallocate(R_Z)

end



Loading

0 comments on commit 1543dce

Please sign in to comment.