! -*- mode: f90 -*-

!===================
&BasicParam
!===================

! Flow type (1=Lock-exchange, 2=TGV, 3=Channel, 4=Periodic hill, 5=Cylinder, 6=dbg-schemes)
itype = 3

! Domain decomposition
p_row=0               ! Row partition
p_col=0               ! Column partition

! Mesh
nx=128                ! X-direction nodes
ny=21                 ! Y-direction nodes
nz=64                 ! Z-direction nodes
istret = 0            ! y mesh refinement (0:no, 1:center, 2:both sides, 3:bottom)
beta = 0.259065151    ! Refinement parameter (beta)

! Domain
xlx = 8.              ! Lx (Size of the box in x-direction)
yly = 2.              ! Ly (Size of the box in y-direction)
zlz = 4.              ! Lz (Size of the box in z-direction)

! Flow parameters
re  = 4200.           ! nu=1/re (Kinematic Viscosity)
cpg = F               ! if cpg=T, then re is friction Reynolds number

! Time stepping
dt = 0.005            ! Time step
ifirst = 1            ! First iteration
ilast = 100        ! Last iteration

! Enable modelling tools
ilesmod=0             ! if 0 then DNS
numscalar = 0         ! How many scalars? (Set to zero to disable scalars)
iibm=0                ! Flag for immersed boundary method

! Boundary and initial conditions
iin = 1               ! Inflow conditions (1: classic, 2: turbinit, 4: SEM)
u1  = 1.              ! u1 (max velocity) (for inflow condition)
u2  = 1.              ! u2 (min velocity) (for inflow condition)
init_noise = 0.125    ! Turbulence intensity (1=100%) !! Initial condition
inflow_noise = 0.0    ! Turbulence intensity (1=100%) !! Inflow condition
idir_stream  = 1      ! Index of the streamwise direction (1=X, 3=Z, 2 not available)

nclx1 = 0
nclxn = 0
ncly1 = 2
nclyn = 2
nclz1 = 0
nclzn = 0

/End

!====================
&NumOptions
!====================

! Spatial derivatives
ifirstder = 4         ! (1->2nd central, 2->4th central, 3->4th compact, 4-> 6th compact)
isecondder = 5        ! (1->2nd central, 2->4th central, 3->4th compact, 4-> 6th compact, 5->hyperviscous 6th)
ipinter = 3           ! interpolation scheme (1: classic, 2: optimized, 3: optimized agressive)

! Time scheme
iimplicit = 2
itimescheme = 3       ! Time integration scheme (1->Euler,2->AB2, 3->AB3, 4->AB4,5->RK3,6->RK4, 7-->CN2+AB3)

! Dissipation control
nu0nu = 4.0             ! Ratio between hyperviscosity/viscosity at nu
cnu = 0.44               ! Ratio between hypervisvosity at k_m=2/3pi and k_c= pi

/End

!=================
&InOutParam
!=================

! Basic I/O
irestart = 0          ! Read initial flow field ?
icheckpoint = 5000    ! Frequency for writing backup file
ioutput = 10000       ! Frequency for visualization
ilist = 1            ! Frequency for writing to screen        
nvisu = 1             ! Size for visualisation collection
nprobes = 4           ! Number of probes

/End

!=================
&Statistics
!=================

wrotation = 0.12      ! rotation speed to trigger turbulence
spinup_time = 5000    ! number of time steps with a rotation to trigger turbulence
nstat = 1             ! Size arrays for statistic collection
initstat = 40000      ! Time steps after which statistics are collected

/End

!=================
&ProbesParam
!=================

flag_all_digits = T     ! By default, only 6 digits, 16 when True
flag_extra_probes = T   ! By default, gradients are not monitored
xyzprobes(1,1) = 0.1    ! X position of probe 1
xyzprobes(2,1) = 0.1    ! Y position of probe 1
xyzprobes(3,1) = 0.1    ! Z position of probe 1
xyzprobes(1,2) = 0.1    ! X position of probe 2
xyzprobes(2,2) = 0.2    ! Y position of probe 2
xyzprobes(3,2) = 0.1    ! Z position of probe 2
xyzprobes(1,3) = 0.1    ! X position of probe 3
xyzprobes(2,3) = 0.3    ! Y position of probe 3
xyzprobes(3,3) = 0.1    ! Z position of probe 3
xyzprobes(1,4) = 0.1    ! X position of probe 4
xyzprobes(2,4) = 0.4    ! Y position of probe 4
xyzprobes(3,4) = 0.1    ! Z position of probe 4

/End

!########################
! OPTIONAL PARAMETERS
!#######################

!================
&ScalarParam
!================

Sc(1) = 1.0           ! Schmidt number

nclxS1 = 0
nclxSn = 0
nclyS1 = 2
nclySn = 2
nclzS1 = 0
nclzSn = 0

/End

!================
&LESModel
!================

iles = 0             ! LES Model (1: Phys Smag, 2: Phys WALE, 3: Phys dyn. Smag, 4: iSVV, 5: dyn SEV)
smagcst = 0.14        ! Smagorinsky constant
walecst = 0.5         ! WALES Model Coefficient
iwall = 0             ! Enable wall modeling

/End

!================
&WallModel
!================

SmagWallDamp = 0      ! Smagorinsky damping function for ABL if 1

/End

&CASE
/End