AFiD/AuxiliaryRoutines.F90 at master · PhysicsofFluids/AFiD · 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
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!                                                         !
!    FILE: AuxiliaryRoutines.F90                          !
!    CONTAINS: subroutines Allocate*,Destroy*             !
!                                                         !
!    PURPOSE: Auxiliary routines used for memory allocs   !
!     and memory freeing                                  !
!                                                         !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

      module AuxiliaryRoutines
      contains

        subroutine AllocateReal1DArray(var,st1,en1)
        use decomp_2d
        implicit none
        integer, intent(in) :: st1,en1
        real,allocatable,dimension(:),intent(inout) :: var
        integer :: alloc_stat, errorcode

        if (.not. allocated(var)) allocate(var(st1:en1), stat=alloc_stat)

        if (alloc_stat /= 0) then
           call decomp_2d_abort(errorcode, &
                'Memory allocation failed when creating new arrays')
        end if

        var =0.0d0

        return
        end subroutine AllocateReal1DArray

!===========================================================================

        subroutine AllocateInt1DArray(var,st1,en1)
        use decomp_2d
        implicit none
        integer, intent(in) :: st1,en1
        integer,allocatable,dimension(:),intent(inout) :: var
        integer :: alloc_stat, errorcode

        if (.not. allocated(var)) allocate(var(st1:en1), stat=alloc_stat)

        if (alloc_stat /= 0) then
           errorcode = 8
           call decomp_2d_abort(errorcode, &
                'Memory allocation failed when creating new arrays')
        end if

        var =0

        return
        end subroutine AllocateInt1DArray

!===========================================================================

        subroutine AllocateReal2DArray(var,st1,en1,st2,en2)
        use decomp_2d
        implicit none
        integer, intent(in) :: st1,en1,st2,en2
        real,allocatable,dimension(:,:),intent(inout) :: var
        integer :: alloc_stat, errorcode

        if (.not. allocated(var)) allocate(var(st1:en1,st2:en2), stat=alloc_stat)

        if (alloc_stat /= 0) then
           errorcode = 8
           call decomp_2d_abort(errorcode, &
                'Memory allocation failed when creating new arrays')
        end if

        var =0.0

        return
        end subroutine AllocateReal2DArray

!===========================================================================

        subroutine AllocateReal3DArray(var,st1,en1,st2,en2,st3,en3)
        use decomp_2d
        implicit none
        integer, intent(in) :: st1,en1,st2,en2,st3,en3
        real,allocatable,dimension(:,:,:),intent(inout) :: var
        integer :: alloc_stat, errorcode

        if (.not. allocated(var)) allocate(var(st1:en1,st2:en2,st3:en3), stat=alloc_stat)

        if (alloc_stat /= 0) then
           errorcode = 8
           call decomp_2d_abort(errorcode, &
                'Memory allocation failed when creating new arrays')
        end if

        var =0.0

        return
        end subroutine AllocateReal3DArray

!===========================================================================

        subroutine DestroyReal1DArray(var)
        use decomp_2d
        implicit none
        real,allocatable,dimension(:),intent(inout) :: var

        if (allocated(var)) deallocate(var)

        return
        end subroutine DestroyReal1DArray

!===========================================================================

        subroutine DestroyInt1DArray(var)
        use decomp_2d
        implicit none
        integer,allocatable,dimension(:),intent(inout) :: var

        if (allocated(var)) deallocate(var)

        return
        end subroutine DestroyInt1DArray

!===========================================================================

        subroutine DestroyReal2DArray(var)
        use decomp_2d
        implicit none
        real,allocatable,dimension(:,:),intent(inout) :: var

        if (allocated(var)) deallocate(var)

        return
        end subroutine DestroyReal2DArray

!===========================================================================

        subroutine DestroyReal3DArray(var)
        use decomp_2d
        implicit none
        real,allocatable,dimension(:,:,:),intent(inout) :: var

        if (allocated(var)) deallocate(var)

        return
        end subroutine DestroyReal3DArray

      end module AuxiliaryRoutines