Skip to content
Snippets Groups Projects
Commit 29d036c6 authored by Antoine Cyril David Hoffmann's avatar Antoine Cyril David Hoffmann :seedling:
Browse files

Error on the weights for simpson rule

parent 0d7d13ac
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 3 additions and 5 deletions
src/grid_mod.F90
+ 3
− 5
View file @ 29d036c6
...@@ -480,12 +480,10 @@ CONTAINS ...@@ -480,12 +480,10 @@ CONTAINS
! Weitghs for Simpson rule ! Weitghs for Simpson rule
ALLOCATE(zweights_SR(izs:ize)) ALLOCATE(zweights_SR(izs:ize))
DO iz = izs,ize DO iz = izs,ize
IF((iz .EQ. 1) .OR. (iz .EQ. Nz)) THEN IF(MODULO(iz,2) .EQ. 1) THEN ! odd iz
zweights_SR(iz) = 1._dp
ELSEIF(MODULO(iz-1,2)) THEN
zweights_SR(iz) = 4._dp
ELSE
zweights_SR(iz) = 2._dp zweights_SR(iz) = 2._dp
ELSE ! even iz
zweights_SR(iz) = 4._dp
ENDIF ENDIF
ENDDO ENDDO
END SUBROUTINE set_zgrid END SUBROUTINE set_zgrid
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment