The ABL case user-file

Contains: user specified routines for the ABL case. More...

Functions/Subroutines
subroutine	uservp (ix, iy, iz, ieg)
	Set variable properties, does not call any subroutines / functions. More...
subroutine	userf (ix, iy, iz, eg)
subroutine	userq (ix, iy, iz, eg)
subroutine	userbc (ix, iy, iz, iside, eg)
subroutine	useric (ix, iy, iz, eg)
	Set up initial conditions. More...
subroutine	userchk ()
	Compute the turbulent stress tensors and write statistics. More...
subroutine	userqtl
subroutine	usrdat ()
subroutine	usrdat2 ()
subroutine	usrdat3 ()
subroutine	eddy_visc (e)
	Compute eddy viscosity using dynamic smagorinsky model. More...
subroutine	set_grid_spacing (dg2)
	Compute D^2, the grid spacing used in the DS sgs model. More...
subroutine	set_obj
	Define objects for surface integrals. More...
subroutine	set_forcing (f_new, u, idir)
	Compute inflow / outflow conditions a.k.a. driving force. Ubar=1. More...
subroutine	coriolis (ix, iy, iz, ieg)
	Compute Coriolis acceleration \[ F_x, F_z = - f\hat{e}_y \times \vec{u} \] See https://lists.mcs.anl.gov/pipermail/nek5000-users/2014-May/002798.html. More...

Detailed Description

Contains: user specified routines for the ABL case.

• uservp : variable properties
• userf : local acceleration term for fluid
• userq : local source term for scalars
• userbc : boundary conditions
• useric : initial conditions
• userchk : general purpose routine for checking errors etc.
• userqtl : thermal divergence for lowMach number flows
• usrdat : modify element vertices
• usrdat2 : modify mesh coordinates
• usrdat3 : general purpose routine for initialization

Function/Subroutine Documentation

coriolis()

subroutine coriolis	(	integer	ix,
		integer	iy,
		integer	iz,
		integer	ieg
	)

Compute Coriolis acceleration

\[ F_x, F_z = - f\hat{e}_y \times \vec{u} \]

See https://lists.mcs.anl.gov/pipermail/nek5000-users/2014-May/002798.html.

Parameters

[in]	ix,iy,iz	GLL point index
[in]	ieg	global element number

eddy_visc()

subroutine eddy_visc

(

integer

e

)

Compute eddy viscosity using dynamic smagorinsky model.

Here is the call graph for this function:

set_forcing()

subroutine set_forcing	(		f_new,
			u,
			idir
	)

Compute inflow / outflow conditions a.k.a. driving force. Ubar=1.

Here is the call graph for this function:

set_grid_spacing()

subroutine set_grid_spacing

(

real, dimension(nx1,ny1,nz1,nelv)

dg2

)

Compute D^2, the grid spacing used in the DS sgs model.

Here is the call graph for this function:

set_obj()

subroutine set_obj

(

)

Define objects for surface integrals.

Here is the call graph for this function:

userbc()

subroutine userbc	(	integer	ix,
		integer	iy,
		integer	iz,
		integer	iside,
		integer	eg
	)

userchk()

subroutine userchk

(

)

Compute the turbulent stress tensors and write statistics.

Here is the call graph for this function:

userf()

subroutine userf	(	integer	ix,
		integer	iy,
		integer	iz,
		integer	eg
	)

useric()

subroutine useric	(	integer	ix,
		integer	iy,
		integer	iz,
		integer	eg
	)

Set up initial conditions.

Use standard subroutine RANDOM_NUMBER. See https://gcc.gnu.org/onlinedocs/gfortran/RANDOM_005fNUMBER.html

userq()

subroutine userq	(	integer	ix,
		integer	iy,
		integer	iz,
		integer	eg
	)

userqtl()

subroutine userqtl

(

)

Here is the call graph for this function:

uservp()

subroutine uservp	(	integer	ix,
		integer	iy,
		integer	iz,
		integer	ieg
	)

Set variable properties, does not call any subroutines / functions.

Parameters

ix	x-index
iy	y-index
iz	z-index
ieg	element index

usrdat()

subroutine usrdat

(

)

Here is the call graph for this function:

usrdat2()

subroutine usrdat2

(

)

usrdat3()

subroutine usrdat3

(

)