MODULE mesh
SAVE
! ---- the declarationf for
! ---- the upper/right (r) and lower/left (l) boundaries
! ---- of our (up to) 3D patch including the ghost-cells
! ---- The grid spacing is (dx,dy,dz)
! ---- The cells are centered at (xc,yc,zc)
DOUBLE PRECISION:: xl, xr
DOUBLE PRECISION:: yl, yr
DOUBLE PRECISION:: zl, zr
DOUBLE PRECISION:: dx, dy, dz
! ---- the declarationf for
! ---- the upper/right (r) and lower/left (l) boundaries
! ---- of the entire (up to) 3D domain excluding the ghost-cells
DOUBLE PRECISION:: xlg, xrg
DOUBLE PRECISION:: ylg, yrg
DOUBLE PRECISION:: zlg, zrg
END MODULE mesh
MODULE array_bounds
SAVE
! ---- The array indicies of the upper(hi) and lower (lo)
! ---- bounds of our domain, including ghost-cells
! ---- The array address of the physical domain
! ---- is (1:nx) x (1:ny) x (1:nz)
! ---- Most of these are set in Generic_SetArrayBounds.f90
INTEGER:: nx,ny,nz
INTEGER:: mx,my,mz
! ---- Bounds including ghost-cells for arrays
INTEGER:: ixlo
INTEGER:: ixhi
INTEGER:: iylo
INTEGER:: iyhi
INTEGER:: izlo
INTEGER:: izhi
! ---- Generic array bounds used when doing x,y or z sweeps.
INTEGER:: inlo
INTEGER:: inhi
! ---- Generic array bounds used when doing x,y or sweeps.
! ---- across derivatives
INTEGER:: idnlo
INTEGER:: idnhi
END MODULE array_bounds
MODULE method_parms
SAVE
! ---- The parameters that relate
! ---- the number of ghost cells to the size of a
! ---- given ENO stencil.
! ---- The parameter nvar is the number of
! ---- elements in the state vector -
! ---- This number should be >= 6 ( for the 3D problem)
! ---- 5 for the phsycial dependant variables, 1 for a
! ---- passive scalar used to determine gamma
! ---- more for other passive scalars
INTEGER :: dim
INTEGER :: ncomps
INTEGER :: nvars
INTEGER :: nscal
INTEGER :: enoOrder
INTEGER :: nghost
! ---- Flag set for using standard scheme
! ---- or optimized ( set to 1 for optimized)
INTEGER :: order
INTEGER :: optimized
INTEGER :: use_carbfix
INTEGER :: use_dcflag
INTEGER :: useExOutput
! ---- Determins the width of the stencil to use.
! ---- enoOder = (stencil - 1) /2 is the order of
! ---- accuracy of a canidate stencil inside WENO
INTEGER :: stencil
! ---- upper and lower bounds for center difference interp
INTEGER :: upb,lob
! ---- Selects the correct different solver
! ---- method = 0 is WENO
! ---- method = 1 is WENO-TCD
! ---- method = 2 is 'New Hybrid'
INTEGER :: method = 0
! ---- if not using Time Refinement set to 1.
INTEGER :: noTimeRefine = 0
! ---- Select the resolved scale viscous terms
! ---- useViscous = 0 no viscous terms
! ---- useViscous = 1 use viscous terms
! ---- ? this requires more ghost cells
INTEGER :: useViscous = 0
! ---- Determines if LES is being used in 3D
! ---- useLES = 0 no LES model
! ---- useLES = 1 use the LES model
! ---- this requires more ghost cells
INTEGER :: useLES = 0
! ---- Determine if the source term is to be used
! ---- useSource = 0 no source term
! ---- useSource = 1 source term - defined as weno_source
! in local problem directory
INTEGER :: useSource = 0
! ---- Are periodic boundary conditions used
! ---- in the x,y,z direction
INTEGER :: xper = 0
INTEGER :: yper = 0
INTEGER :: zper = 0
! these flags control the characteristic boundary condition type
INTEGER :: cbc_direction(6)
DATA cbc_direction /0, 0, 0, 0, 0, 0/
INTEGER :: cbc_ixlow, cbc_ixup
INTEGER :: cbc_iylow, cbc_iyup
INTEGER :: cbc_izlow, cbc_izup
! these flags control the boundary stencil
INTEGER :: bc_ixlow, bc_ixup
INTEGER :: bc_iylow, bc_iyup
INTEGER :: bc_izlow, bc_izup
! filtering variables
DOUBLE PRECISION :: alpha_eta1, alpha_eta2
END MODULE method_parms
subroutine cles_getiparam(name, value, ierr)
use method_parms
implicit none
include 'cles.i'
character*(*) name
integer value, ierr
character(LEN=12) cname
cname = name
ierr = CLES_RETURN_OK
select case (cname)
case ('nvars')
value = nvars
case ('ncomps')
value = ncomps
case ('dim')
value = dim
case ('nscal')
value = nscal
case ('stencil')
value = stencil
case ('method')
value = method
case ('notimerefine')
value = noTimeRefine
case ('order')
value = order
case ('nghost')
value = nghost
case ('optimized')
value = optimized
case ('use_carbfix')
value = use_carbfix
case ('enoorder')
value = enoOrder
case ('useles')
value = useLES
case ('useviscous')
value = useViscous
case ('usesource')
value = useSource
case default
ierr = CLES_RETURN_ERROR
end select
end subroutine cles_getiparam
MODULE Interp_coeffs
USE method_parms
SAVE
! ---- The declarations for the coefficients used
! ---- in calculating flux interpolants by WENO
! ---- make these too big, but static
DOUBLE PRECISION :: aw(0:10,0:10-1)
DOUBLE PRECISION :: cw(0:10)
DOUBLE PRECISION :: dw(0:10,10-1,0:10-1)
DOUBLE PRECISION :: ac, cc, dc, b2, b3
! ---- used in the Center Difference (Finite difference)
INTEGER :: tcd_bndry_width = 2
INTEGER :: tcd_bndry_points = 4
DOUBLE PRECISION :: tcd_center(3)
DOUBLE PRECISION :: tcd_flux(3)
! dimensions of tcd_bndry are max number of points (4) x max number of stencils (2)
DOUBLE PRECISION :: tcd_bndry(4,2)
END MODULE Interp_coeffs
module cles_interfaces
interface
subroutine cles_xLocation(i, x)
integer i
double precision x
end subroutine cles_xLocation
end interface
interface
subroutine cles_yLocation(i, x)
integer i
double precision x
end subroutine cles_yLocation
end interface
interface
subroutine cles_zLocation(i, x)
integer i
double precision x
end subroutine cles_zLocation
end interface
interface
subroutine cles_eqstate(q,ncomps,qt,nvars,n,useLES)
integer n, ncomps, nvars, useLES
DOUBLE PRECISION q(ncomps,n), qt(nvars,n)
end subroutine cles_eqstate
end interface
interface
subroutine cles_roe(ux0, ux1, ncomps, vx0, vx1, nvars, &
gamma, Rr, mu, nscal)
integer ncomps, nvars, nscal
double precision ux0(ncomps), ux1(ncomps)
double precision vx0(nvars), vx1(nvars)
double precision gamma, Rr, mu(*)
end subroutine cles_roe
end interface
interface
subroutine SetUpLES(version_)
INTEGER, INTENT(IN) :: version_
end subroutine SetUpLES
end interface
interface
subroutine InitializeLES(ux,vx)
use mesh
use array_bounds
use method_parms
DOUBLE PRECISION, INTENT(INOUT) :: ux(ncomps,ixlo:ixhi,iylo:iyhi,izlo:izhi)
DOUBLE PRECISION, INTENT(IN) :: vx(nvars,ixlo:ixhi,iylo:iyhi,izlo:izhi)
end subroutine InitializeLES
end interface
interface
subroutine AddSgsFlux(fxi,ux,vx,direction)
use mesh
use array_bounds
use method_parms
DOUBLE PRECISION, INTENT(IN) :: ux(ncomps,ixlo:ixhi,iylo:iyhi,izlo:izhi)
DOUBLE PRECISION, INTENT(IN) :: vx(nvars,ixlo:ixhi,iylo:iyhi,izlo:izhi)
DOUBLE PRECISION, INTENT(INOUT)::fxi(ncomps,ixlo:ixhi,iylo:iyhi,izlo:izhi,3)
INTEGER, INTENT(IN) :: direction
end subroutine AddSgsFlux
end interface
interface
integer function cleslog_active(mod, loc)
integer mod
integer loc
end function cleslog_active
end interface
interface
subroutine cleslog_log(slog)
character*(*) slog
end subroutine cleslog_log
end interface
interface
subroutine cleslog_log_flush()
end subroutine cleslog_log_flush
end interface
end module cles_interfaces