-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathsaveBox.f90
53 lines (45 loc) · 1.96 KB
/
saveBox.f90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
!############################################################################
!
subroutine saveBox (ux1,uy1,uz1,pp3,nx1,nx2,ny1,ny2,nz1,nz2)
!
!############################################################################
USE param
USE variables
USE decomp_2d
USE decomp_2d_io
implicit none
TYPE(DECOMP_INFO) :: phG
real(mytype),dimension(xsize(1),xsize(2),xsize(3)) :: ux1,uy1,uz1,pp3
integer :: code,icomplet,ip,jp,kp,nx1,nx2,ny1,ny2,nz1,nz2
character(len=20) nfichier,nfichier1
character(len=20) :: filename
if (itime==ifirst ) then
call system('rm -rf subSave ')
call system('mkdir subSave')
call system('mkdir ./subSave/ux')
call system('mkdir ./subSave/uy')
call system('mkdir ./subSave/uz')
call system('mkdir ./subSave/pp')
endif
!! ###########################################################################
!! ------------------- UX ---------------------------------------------------
1 format('./subSave/ux/ux',I5.5)
write(filename, 1) (itime)
call decomp_2d_write_subdomain(1,ux1,nx1,nx2,ny1,ny2,nz1,nz2,filename)
!! ###########################################################################
!! ------------------- UY ---------------------------------------------------
2 format('./subSave/uy/uy',I5.5)
write(filename, 2) (itime)
call decomp_2d_write_subdomain(1,uy1,nx1,nx2,ny1,ny2,nz1,nz2,filename)
!! ###########################################################################
!! ------------------- UZ ---------------------------------------------------
3 format('./subSave/uz/uz',I5.5)
write(filename, 3) (itime)
call decomp_2d_write_subdomain(1,uz1,nx1,nx2,ny1,ny2,nz1,nz2,filename)
!! ###########################################################################
!! ------------------- Prerssure ---------------------------------------------------
4 format('./subSave/pp/pp',I5.5)
write(filename, 4) (itime)
call decomp_2d_write_subdomain(1,pp3,nx1,nx2,ny1,ny2,nz1,nz2,filename)
!############################################################################
end subroutine saveBox