diff --git a/code/m_parameters.f90 b/code/m_parameters.f90 index 7ec296b..c300213 100644 --- a/code/m_parameters.f90 +++ b/code/m_parameters.f90 @@ -4,6 +4,8 @@ module m_parameters integer :: nxp,nyp,nzp + real :: hxp,hyp,hzp + real :: l_0 integer :: startnum,endnum,skipnum,shiftnum @@ -22,6 +24,8 @@ module m_parameters integer, allocatable :: omit_t(:,:) + real, parameter :: pi=3.14159265358979323846 + real, parameter :: me=1.00e-20 contains @@ -109,6 +113,11 @@ subroutine read_intro prp=scp/lep +! hyp=l_0*pi/REAL(nyp) + hyp=l_0*pi/REAL(nyp-1) ! kwon + hxp=hyp + hzp=hyp + close (100) close (101) diff --git a/code/post.f90 b/code/post.f90 index 526dcf6..df0e6b7 100644 --- a/code/post.f90 +++ b/code/post.f90 @@ -9,9 +9,6 @@ INTEGER :: countnum INTEGER :: num_, dummyu_ ! hybrid REAL :: tnow - REAL :: hxp,hyp,hzp - REAL, PARAMETER :: pi=3.14159265358979323846 - REAL, PARAMETER :: me=1.00e-20 REAL, DIMENSION(:), ALLOCATABLE :: favg_ndata !REAL, DIMENSION(:,:), ALLOCATABLE :: SMF,SMC,FMS,KMS REAL, DIMENSION(:,:), ALLOCATABLE :: SMF,SMC @@ -1580,10 +1577,6 @@ !ALLOCATE(IRE9(14,nxp),STAT=ierr) ; IRE9=0. -! hyp=l_0*pi/REAL(nyp) - hyp=l_0*pi/REAL(nyp-1) ! kwon - hxp=hyp; hzp=hyp - WRITE(*,'(a6,i3,a8,i3,a8,i3)') ' NX = ',nxp,' / NY = ',nyp,' / NZ = ',nzp WRITE(*,*) ' Preparing memory space for COMPACT SCHEME' CALL ludcmp(nxp,nyp,nzp,1,0,0) ! 1,1,0