From cec45a8d7a95a7c250372081398f70a7171a5c58 Mon Sep 17 00:00:00 2001
From: Antoine Cyril David Hoffmann <ahoffman@spcpc606.epfl.ch>
Date: Tue, 7 Jul 2020 10:11:11 +0200
Subject: [PATCH] Fixed the delta k computation
---
src/fourier_grid_mod.F90 | 12 ++++++++++--
1 file changed, 10 insertions(+), 2 deletions(-)
diff --git a/src/fourier_grid_mod.F90 b/src/fourier_grid_mod.F90
index b2765095..fae1d38a 100644
--- a/src/fourier_grid_mod.F90
+++ b/src/fourier_grid_mod.F90
@@ -52,7 +52,11 @@ CONTAINS
ikrs = 1
ikre = nkr
! Grid spacings, precompute some inverses
- deltakr = (krmax - krmin) / REAL(nkr,dp)
+ IF ( nkr .GT. 1 ) THEN ! To avoid case with 0 intervals
+ deltakr = (krmax - krmin) / REAL(nkr-1,dp)
+ ELSE
+ deltakr = 0;
+ ENDIF
! Discretized kr positions
ALLOCATE(krarray(ikrs:ikre))
DO ikr = ikrs,ikre
@@ -68,7 +72,11 @@ CONTAINS
ikzs = 1
ikze = nkz
! Grid spacings, precompute some inverses
- deltakz = (kzmax - kzmin) / REAL(nkz,dp)
+ IF ( nkz .GT. 1 ) THEN
+ deltakz = (kzmax - kzmin) / REAL(nkz-1,dp)
+ ELSE
+ deltakz = 0;
+ ENDIF
! Discretized kz positions
ALLOCATE(kzarray(ikzs:ikze))
DO ikz = ikzs,ikze
--
GitLab