From c6dca2fbb85f1b864928065cd784392fe84988ee Mon Sep 17 00:00:00 2001 From: Paul Bartholomew Date: Sat, 5 Oct 2024 16:39:01 +0100 Subject: [PATCH] The volume averaging denominator overflows in TGV postprocessing This is resolved by using 64-bit integers, placing the mesh limit at 2M cubed domains which should be sufficient for a while --- src/Case-TGV.f90 | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/src/Case-TGV.f90 b/src/Case-TGV.f90 index f55014377..910f3df01 100644 --- a/src/Case-TGV.f90 +++ b/src/Case-TGV.f90 @@ -188,6 +188,8 @@ end subroutine init_post !############################################################################ subroutine postprocess_tgv(ux1,uy1,uz1,phi1,ep1) + use iso_fortran_env + USE decomp_2d_io USE variables, only: nx,ny,nz USE MPI @@ -206,7 +208,8 @@ subroutine postprocess_tgv(ux1,uy1,uz1,phi1,ep1) real(mytype), allocatable, dimension(:,:,:) :: bx2,by2,bz2,bx3,by3,bz3 real(mytype) :: eek, enst, eps, eps2 real(mytype) :: eem,omegam,jmax,disb,omegab - integer :: nxc, nyc, nzc, xsize1, xsize2, xsize3 + integer :: xsize1, xsize2, xsize3 + integer(int64) :: nxc, nyc, nzc integer :: i,j,k,is,code,nvect1 character(len=30) :: filename