From 29d036c6d8214556c2a0669beb22bd2d11c6574d Mon Sep 17 00:00:00 2001
From: Antoine Hoffmann <antoine.hoffmann@epfl.ch>
Date: Thu, 9 Jun 2022 11:01:24 +0200
Subject: [PATCH] Error on the weights for simpson rule
---
src/grid_mod.F90 | 8 +++-----
1 file changed, 3 insertions(+), 5 deletions(-)
diff --git a/src/grid_mod.F90 b/src/grid_mod.F90
index 77ea54e5..f58cbb86 100644
--- a/src/grid_mod.F90
+++ b/src/grid_mod.F90
@@ -480,12 +480,10 @@ CONTAINS
! Weitghs for Simpson rule
ALLOCATE(zweights_SR(izs:ize))
DO iz = izs,ize
- IF((iz .EQ. 1) .OR. (iz .EQ. Nz)) THEN
- zweights_SR(iz) = 1._dp
- ELSEIF(MODULO(iz-1,2)) THEN
- zweights_SR(iz) = 4._dp
- ELSE
+ IF(MODULO(iz,2) .EQ. 1) THEN ! odd iz
zweights_SR(iz) = 2._dp
+ ELSE ! even iz
+ zweights_SR(iz) = 4._dp
ENDIF
ENDDO
END SUBROUTINE set_zgrid
--
GitLab