FixedPivot2D/src/MAIN.f90 at master · HugoMVale/FixedPivot2D · GitHub

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
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
PROGRAM FIXEDPIVOT2D

!----------------------------------------------------------------------------------------------
!
!DESCRIPTION:
!
!This program is a numerical implementation of the extended fixed pivot method to solve the PBE
!for 2-component aggregation, as described in https://doi.org/10.1021/ie050179s.
!
!DISCLAIMER:
!
!Some parts of the code were modified in 2019 to improve its readability and use. During this
!process some bugs may also have been introduced.
!
!KEY STEPS FOR USE:
!
!1. Define the desired coagulation kernel
!   Function COAGKERNEL in file INPUTS.f90
!
!2. Define initial PSD
!   Function PSDINI in file INPUTS.f90
!
!3. Define end integration time
!   Variable: time_final just below
!
!4. Define grid range
!   Variables: xmin, xmax, ymin, ymax just below
!
!5. Define grid type
!   Variables: meshx, meshy just below
!   See subroutine GRID in file GRID.f90 for more information.
!
!6. Define number of grid cells
!   Variables: Mx, My in module GLOBAL in file GLOBAL.f90
!
!7. Update folder path
!   Subroutine OUTPUT in file OUTPUT.f90
!
!8. Build and run
!   DO not forget to include DLSODE or ODEPACK
!   https://computing.llnl.gov/casc/odepack/
!
!----------------------------------------------------------------------------------------------

USE GLOBAL

IMPLICIT NONE

!\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\

INTEGER :: meshx, meshy

!\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\

!	*** END-TIME SPECIFICATION ***

time_final = 100.0d0

!	*** GRID SPECIFICATION ***

xmin = 1.0d-3
xmax = 1.0d4

ymin = 1.0d-3
ymax = 1.0d4

! For options, see GRID in GRID.f90
meshx=2
meshy=2


!	*** START OF CALCULATION SEQUENCE ***

!Start-of-simulation message
CALL OUTPUT (1)

!Build grid (different options available)
CALL GRID(meshx,meshy)

!Compute possible particle combinations and weights
CALL COMBINATIONS

!Set initial PSD
CALL AVERAGECELL

!Evaluate the coagulation matrix (for systems where coagulation kernel is time independent)
CALL COAGMATRIX

!Integrate discretized PBE
CALL INTEGRATION

!End-of-simulation message
CALL OUTPUT(4)


END PROGRAM FIXEDPIVOT2D
!##############################################################################################