-
Notifications
You must be signed in to change notification settings - Fork 5
/
Copy pathsizes.f90
77 lines (58 loc) · 2.86 KB
/
sizes.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
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
!>
!! \brief This module contains basic size parameter definitions
!!
!! Module for C2-Ray
!!
!! \b Author: Garrelt Mellema
!!
!! \b Date: 2006-08-20
module sizes
! Module for C2-Ray
! Author: Garrelt Mellema
! Date: 2006-08-20
! This module is also accepted by the F compiler
! This module contains basic parameter definitions
! use precision
implicit none
private
!> Number of spatial dimensions
integer,parameter,public :: Ndim=3
!> Size of the mesh for spatial coordinate.
!for test problem
!integer,dimension(Ndim),parameter,public :: mesh=(/ 100, 100, 100 /)
!for 6912^3 particles, 13824^3 grid, 500/h Mpc box
integer,dimension(Ndim),parameter,public :: mesh=(/ 300, 300, 300 /)
!integer,dimension(Ndim),parameter,public :: mesh=(/ 600, 600, 600 /)
!integer,dimension(Ndim),parameter,public :: mesh=(/ 1200, 1200, 1200 /)
!for 4000^3 particles, 8000^3 grid, 244/h Mpc box
!integer,dimension(Ndim),parameter,public :: mesh=(/ 250, 250, 250 /)
!integer,dimension(Ndim),parameter,public :: mesh=(/ 500, 500, 500 /)
!for 1728^3 particles, 3456^3 grid, 47/h Mpc box
!integer,dimension(Ndim),parameter,public :: mesh=(/ 306, 306, 306 /)
!integer,dimension(Ndim),parameter,public :: mesh=(/ 612, 612, 612 /)
!for 5488^3 particles, 10976^3 grid, 425/h Mpc box
!integer,dimension(Ndim),parameter,public :: mesh=(/ 252, 252, 252 /)
!integer,dimension(Ndim),parameter,public :: mesh=(/ 504, 504, 504 /)
!integer,dimension(Ndim),parameter,public :: mesh=(/ 784, 784, 784 /)
!for 1024^3 particles, 2048^3 grid, 37/h Mpc box
!integer,dimension(Ndim),parameter,public :: mesh=(/ 256, 256, 256 /)
!for 1728^3 particles, 3456^3 grid, 64/h Mpc box
!integer,dimension(Ndim),parameter,public :: mesh=(/ 216, 216, 216 /)
!integer,dimension(Ndim),parameter,public :: mesh=(/ 432, 432, 432 /)
!integer,dimension(Ndim),parameter,public :: mesh=(/ 864, 864, 864 /)
!for 3072^3 particles, 6144^3 grid, 114/h Mpc box
!integer,dimension(Ndim),parameter,public :: mesh=(/ 256, 256, 256 /)
!integer,dimension(Ndim),parameter,public :: mesh=(/ 384, 384, 384 /)
!integer,dimension(Ndim),parameter,public :: mesh=(/ 512, 512, 512 /)
!for LG simulation: 1024^3 particles, 1024^3 grid, 64/h Mpc box
!integer,dimension(Ndim),parameter,public :: mesh=(/ 512, 512, 512 /)
!integer,dimension(Ndim),parameter,public :: mesh=(/ 256, 256, 256 /)
!integer,dimension(Ndim),parameter,public :: mesh=(/ 128, 128, 128 /)
!for old-style PMFAST simulation, 1624^3 particles, 3248^3 grid, 100/h or 35/h Mpc box
!integer,dimension(Ndim),parameter,public :: mesh=(/ 203, 203, 203 /)
!integer,dimension(Ndim),parameter,public :: mesh=(/406, 406, 406/)
!integer,dimension(Ndim),parameter,public :: mesh=(/812, 812, 812/)
integer,parameter,public :: meshx=mesh(1)
integer,parameter,public :: meshy=mesh(2)
integer,parameter,public :: meshz=mesh(3)
end module sizes