From 44bd301d6b65319d61290ae80cb554702933accd Mon Sep 17 00:00:00 2001
From: Antoine Cyril David Hoffmann <ahoffman@spcpc606.epfl.ch>
Date: Wed, 28 Oct 2020 10:38:58 +0100
Subject: [PATCH] Condition for Nkr=1 was wrong
---
src/grid_mod.F90 | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/grid_mod.F90 b/src/grid_mod.F90
index afd34c69..125fa4a9 100644
--- a/src/grid_mod.F90
+++ b/src/grid_mod.F90
@@ -90,7 +90,7 @@ CONTAINS
IF ( Nr .GT. 1 ) THEN ! To avoid case with 0 intervals
deltar = Lr / REAL(Nr-1,dp)
ELSE
- deltar = 0;
+ deltar = 1;
ENDIF
! Discretized r positions ordered as dx*(0 1 2 -3 -2 -1)
ALLOCATE(rarray(irs:ire))
@@ -116,7 +116,7 @@ CONTAINS
IF ( Nz .GT. 1 ) THEN ! To avoid case with 0 intervals
deltaz = Lz / REAL(Nz-1,dp)
ELSE
- deltaz = 0;
+ deltaz = 1;
ENDIF
! Discretized r positions ordered as dx*(0 1 2 -3 -2 -1)
ALLOCATE(zarray(irs:ire))
--
GitLab