Skip to content
GitLab
Explore
Sign in
Primary navigation
Search or go to…
Project
Gyacomo
Manage
Activity
Members
Labels
Plan
Issues
Issue boards
Milestones
Wiki
Code
Merge requests
Repository
Branches
Commits
Tags
Repository graph
Compare revisions
Snippets
Build
Pipelines
Jobs
Pipeline schedules
Artifacts
Deploy
Releases
Package registry
Model registry
Operate
Environments
Terraform modules
Monitor
Incidents
Analyze
Value stream analytics
Contributor analytics
CI/CD analytics
Repository analytics
Model experiments
Help
Help
Support
GitLab documentation
Compare GitLab plans
Community forum
Contribute to GitLab
Provide feedback
Keyboard shortcuts
?
Snippets
Groups
Projects
Show more breadcrumbs
Antoine Cyril David Hoffmann
Gyacomo
Commits
71e00e53
Commit
71e00e53
authored
4 years ago
by
Antoine Cyril David Hoffmann
Browse files
Options
Downloads
Patches
Plain Diff
cleaned unused variables
parent
8f1d3b41
No related branches found
No related tags found
No related merge requests found
Changes
1
Hide whitespace changes
Inline
Side-by-side
Showing
1 changed file
src/fourier_grid_mod.F90
+13
-17
13 additions, 17 deletions
src/fourier_grid_mod.F90
with
13 additions
and
17 deletions
src/fourier_grid_mod.F90
+
13
−
17
View file @
71e00e53
...
@@ -7,7 +7,7 @@ MODULE fourier_grid
...
@@ -7,7 +7,7 @@ MODULE fourier_grid
! GRID Namelist
! GRID Namelist
INTEGER
,
PUBLIC
,
PROTECTED
::
pmaxe
=
2
! The maximal electron Hermite-moment computed
INTEGER
,
PUBLIC
,
PROTECTED
::
pmaxe
=
2
! The maximal electron Hermite-moment computed
INTEGER
,
PUBLIC
,
PROTECTED
::
jmaxe
=
2
! The maximal electron Laguerre-moment computed
INTEGER
,
PUBLIC
,
PROTECTED
::
jmaxe
=
2
! The maximal electron Laguerre-moment computed
INTEGER
,
PUBLIC
,
PROTECTED
::
pmaxi
=
2
! The maximal ion Hermite-moment computed
INTEGER
,
PUBLIC
,
PROTECTED
::
pmaxi
=
2
! The maximal ion Hermite-moment computed
INTEGER
,
PUBLIC
,
PROTECTED
::
jmaxi
=
2
! The maximal ion Laguerre-moment computed
INTEGER
,
PUBLIC
,
PROTECTED
::
jmaxi
=
2
! The maximal ion Laguerre-moment computed
INTEGER
,
PUBLIC
,
PROTECTED
::
nkr
=
10
! Number of total internal grid points in kr
INTEGER
,
PUBLIC
,
PROTECTED
::
nkr
=
10
! Number of total internal grid points in kr
...
@@ -20,7 +20,6 @@ MODULE fourier_grid
...
@@ -20,7 +20,6 @@ MODULE fourier_grid
! Indices of s -> start , e-> END
! Indices of s -> start , e-> END
INTEGER
,
PUBLIC
,
PROTECTED
::
ips_e
,
ipe_e
,
ijs_e
,
ije_e
INTEGER
,
PUBLIC
,
PROTECTED
::
ips_e
,
ipe_e
,
ijs_e
,
ije_e
INTEGER
,
PUBLIC
,
PROTECTED
::
ips_i
,
ipe_i
,
ijs_i
,
ije_i
INTEGER
,
PUBLIC
,
PROTECTED
::
ips_i
,
ipe_i
,
ijs_i
,
ije_i
INTEGER
,
PUBLIC
,
PROTECTED
::
ns_e
,
ne_e
,
ns_i
,
ne_i
!start/end indices for coll. mat.
INTEGER
,
PUBLIC
,
PROTECTED
::
ikrs
,
ikre
,
ikzs
,
ikze
INTEGER
,
PUBLIC
,
PROTECTED
::
ikrs
,
ikre
,
ikzs
,
ikze
! Toroidal direction
! Toroidal direction
...
@@ -30,7 +29,7 @@ MODULE fourier_grid
...
@@ -30,7 +29,7 @@ MODULE fourier_grid
! Grids containing position in fourier space
! Grids containing position in fourier space
REAL
(
dp
),
DIMENSION
(:),
ALLOCATABLE
,
PUBLIC
,
PROTECTED
::
kzarray
REAL
(
dp
),
DIMENSION
(:),
ALLOCATABLE
,
PUBLIC
,
PROTECTED
::
kzarray
REAL
(
dp
),
DIMENSION
(:),
ALLOCATABLE
,
PUBLIC
,
PROTECTED
::
krarray
REAL
(
dp
),
DIMENSION
(:),
ALLOCATABLE
,
PUBLIC
,
PROTECTED
::
krarray
! Grid containing the polynomials degrees
! Grid containing the polynomials degrees
INTEGER
,
DIMENSION
(:),
ALLOCATABLE
,
PUBLIC
,
PROTECTED
::
parray_e
INTEGER
,
DIMENSION
(:),
ALLOCATABLE
,
PUBLIC
,
PROTECTED
::
parray_e
INTEGER
,
DIMENSION
(:),
ALLOCATABLE
,
PUBLIC
,
PROTECTED
::
parray_i
INTEGER
,
DIMENSION
(:),
ALLOCATABLE
,
PUBLIC
,
PROTECTED
::
parray_i
...
@@ -47,10 +46,10 @@ CONTAINS
...
@@ -47,10 +46,10 @@ CONTAINS
SUBROUTINE
set_krgrid
SUBROUTINE
set_krgrid
USE
prec_const
USE
prec_const
IMPLICIT
NONE
IMPLICIT
NONE
INTEGER
::
ikr
INTEGER
::
ikr
! Start and END indices of grid
! Start and END indices of grid
ikrs
=
1
ikrs
=
1
ikre
=
nkr
ikre
=
nkr
! Grid spacings, precompute some inverses
! Grid spacings, precompute some inverses
IF
(
nkr
.GT.
1
)
THEN
! To avoid case with 0 intervals
IF
(
nkr
.GT.
1
)
THEN
! To avoid case with 0 intervals
deltakr
=
(
krmax
-
krmin
)
/
REAL
(
nkr
-1
,
dp
)
deltakr
=
(
krmax
-
krmin
)
/
REAL
(
nkr
-1
,
dp
)
...
@@ -70,7 +69,7 @@ CONTAINS
...
@@ -70,7 +69,7 @@ CONTAINS
INTEGER
::
ikz
INTEGER
::
ikz
! Start and END indices of grid
! Start and END indices of grid
ikzs
=
1
ikzs
=
1
ikze
=
nkz
ikze
=
nkz
! Grid spacings, precompute some inverses
! Grid spacings, precompute some inverses
IF
(
nkz
.GT.
1
)
THEN
IF
(
nkz
.GT.
1
)
THEN
deltakz
=
(
kzmax
-
kzmin
)
/
REAL
(
nkz
-1
,
dp
)
deltakz
=
(
kzmax
-
kzmin
)
/
REAL
(
nkz
-1
,
dp
)
...
@@ -89,7 +88,7 @@ CONTAINS
...
@@ -89,7 +88,7 @@ CONTAINS
IMPLICIT
NONE
IMPLICIT
NONE
INTEGER
::
ip
,
ij
INTEGER
::
ip
,
ij
ips_e
=
1
;
ipe_e
=
pmaxe
+
1
ips_e
=
1
;
ipe_e
=
pmaxe
+
1
ips_i
=
1
;
ipe_i
=
pmaxi
+
1
ips_i
=
1
;
ipe_i
=
pmaxi
+
1
ALLOCATE
(
parray_e
(
ips_e
:
ipe_e
))
ALLOCATE
(
parray_e
(
ips_e
:
ipe_e
))
ALLOCATE
(
parray_i
(
ips_i
:
ipe_i
))
ALLOCATE
(
parray_i
(
ips_i
:
ipe_i
))
DO
ip
=
ips_e
,
ipe_e
;
parray_e
(
ip
)
=
ip
-1
;
END
DO
DO
ip
=
ips_e
,
ipe_e
;
parray_e
(
ip
)
=
ip
-1
;
END
DO
...
@@ -102,9 +101,6 @@ CONTAINS
...
@@ -102,9 +101,6 @@ CONTAINS
DO
ij
=
ijs_e
,
ije_e
;
jarray_e
(
ij
)
=
ij
-1
;
END
DO
DO
ij
=
ijs_e
,
ije_e
;
jarray_e
(
ij
)
=
ij
-1
;
END
DO
DO
ij
=
ijs_i
,
ije_i
;
jarray_i
(
ij
)
=
ij
-1
;
END
DO
DO
ij
=
ijs_i
,
ije_i
;
jarray_i
(
ij
)
=
ij
-1
;
END
DO
ns_e
=
bare
(
ips_e
,
ijs_e
);
ne_e
=
bare
(
ipe_e
,
ije_e
)
ns_i
=
bari
(
ips_i
,
ijs_i
);
ne_i
=
bari
(
ipe_i
,
ije_i
)
END
SUBROUTINE
set_pj
END
SUBROUTINE
set_pj
SUBROUTINE
fourier_grid_readinputs
SUBROUTINE
fourier_grid_readinputs
...
@@ -146,17 +142,17 @@ CONTAINS
...
@@ -146,17 +142,17 @@ CONTAINS
CALL
attach
(
fidres
,
TRIM
(
str
),
"kzmax"
,
kzmax
)
CALL
attach
(
fidres
,
TRIM
(
str
),
"kzmax"
,
kzmax
)
END
SUBROUTINE
fourier_grid_outputinputs
END
SUBROUTINE
fourier_grid_outputinputs
FUNCTION
bare
(
ip
,
ij
)
FUNCTION
bare
(
p_
,
j_
)
IMPLICIT
NONE
IMPLICIT
NONE
INTEGER
::
bare
,
ip
,
ij
INTEGER
::
bare
,
p_
,
j_
bare
=
(
jmaxe
+1
)
*
i
p
+
i
j
+
1
bare
=
(
jmaxe
+1
)
*
p
_
+
j
_
+
1
END
FUNCTION
END
FUNCTION
FUNCTION
bari
(
ip
,
ij
)
FUNCTION
bari
(
p_
,
j_
)
IMPLICIT
NONE
IMPLICIT
NONE
INTEGER
::
bari
,
ip
,
ij
INTEGER
::
bari
,
p_
,
j_
bari
=
bare
(
pmaxe
,
jmaxe
)
+
(
jmaxi
+1
)
*
i
p
+
i
j
+
1
bari
=
(
jmaxi
+1
)
*
p
_
+
j
_
+
1
END
FUNCTION
END
FUNCTION
...
...
This diff is collapsed.
Click to expand it.
Preview
0%
Loading
Try again
or
attach a new file
.
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Save comment
Cancel
Please
register
or
sign in
to comment