LES Intercomparison Study for Neutral Boundary Layers: Description
Last updated: May 2026
References
Andrén, A., Brown, A. R., Graf, J., Mason, P. J., Moeng, C.-H., Nieuwstadt, F. T. M., & Schumann, U. (1994). Large-eddy simulation of a neutrally stratified boundary layer: a comparison of four codes. Quart. J. Roy. Meteorol. Soc., 120, 1457–1484.
Anderson, W. C., Basu, S., & Letchford, C. W. (2007). Comparison of dynamic subgrid-scale models for simulations of neutrally buoyant shear-driven atmospheric boundary layer flows. Environmental Fluid Mechanics, 7, 195-215.
Case Description
The NBL_A94 case simulates a (truly) neutrally stratified Ekman layer with zero surface heat flux. The flow is driven by a constant geostrophic wind. The original intercomparison case is slightly modified to have realistic simulated boundary layer height.
Parameter |
Value |
|---|---|
Reference run |
|
Domain |
4000 m x 4000 m x 4000 m |
Reference grid |
\(128^3\) |
Run matrix |
\(64^3\), \(128^3\), \(256^3\), \(384^3\) x LASDD-SM, LASDD-WL, LAD-SM, LAD-WL x SP, DP |
Simulation time |
300000 s (83.3333 h) |
Reference time step |
1 s |
Geostrophic wind |
\(U_g = 5\) m s\(^{-1}\), \(V_g = 0\) m s\(^{-1}\) |
Coriolis parameter |
\(f = 1.0 \times 10^{-4}\) s\(^{-1}\) |
Roughness lengths |
\(z_{0m} = 0.01\) m, \(z_{0T} = 0.01\) m |
Thermal surface option |
|
Sensible heat flux |
0 K m s\(^{-1}\) |
Reference temperature |
\(T_0 = 300\) K |
Inversion strength |
0 / 1000 K m\(^{-1}\) |
SGS model |
LASDD-SM ( |
Filter-to-grid ratio |
|
Damping layer |
above \(z = 2500\) m |
Output statistics interval |
60 s |
3D output interval |
SimTime s |
Geostrophic forcing option |
|
Advection forcing option |
|
Moisture option |
|
Boundary-layer type |
Neutral Ekman layer |
Reference averaging window |
hours 75-83.3 |
The representative configuration below is copied from examples/NBL_A94/runs/128x128x128_LASDD_SM_SP/Config.py.
Reference Configuration
[1]:
# Copyright (C) 2025 Sukanta Basu
#
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program. If not, see <https://www.gnu.org/licenses/>.
"""
File: Config.py
===============
:Author: Sukanta Basu
:AI Assistance: Claude.AI (Anthropic) is used for documentation,
code restructuring, and performance optimization
:Date: 2026-05-20
:Description: Neutral BL benchmark after Andren et al. (1994).
Grid: 128x128x128, SGS: LASDD-SM,
Precision: single.
"""
# ============================================================
# Imports
# ============================================================
import numpy as np
# ============================================================
# User Input
# ============================================================
# ------------------------------------------------------------
# Platform options
# ------------------------------------------------------------
use_double_precision = False
# 0: use CPU, 1: use GPU
optGPU = 1
GPU_ID = 0
# ------------------------------------------------------------
# Domain configuration
# ------------------------------------------------------------
# Domain size (m)
l_x = 4000
l_y = 4000
l_z = 4000
# Number of grid points
nx = 128
ny = 128
nz = 128
# ------------------------------------------------------------
# Time integration configuration
# ------------------------------------------------------------
# Change this if it is a restart run
istep = 1
# Time stepping and simulation time
dt = 1.0 # unit: sec
SimTime = 300000 # unit: sec
# Galilean transformation (m/s)
Ugal = 3
# ------------------------------------------------------------
# Surface configuration
# ------------------------------------------------------------
# optSurfFlux: 0 = homogeneous, 1 = heterogeneous
optSurfFlux = 0
# optSurfBC: 0 = constant flux
# 1 = time-varying flux (from SurfaceBCFile)
# 2 = time-varying surface temperature (from SurfaceBCFile)
optSurfBC = 0
# Roughness lengths (m)
z0m = 0.01
z0T = 0.01
# Screen-level temperature reference height (m); 0 = use z0T
zTemperature = 0.0
# SensibleHeatFlux: not used when optSurfBC >= 1; set to 0 as placeholder
SensibleHeatFlux = 0.0 # K m/s
# Path to surface BC file (relative to run directory)
SurfaceBCFile = 'input/SurfaceBC.npz'
# ------------------------------------------------------------
# Forcing configuration
# ------------------------------------------------------------
# Geostrophic wind option:
# 0 = constant Ug2, Vg2 from Config
# 1 = time + height varying, loaded from GeoWindFile
optGeoWind = 0
# Constant geostrophic wind (m/s)
Ug2 = 5
Vg2 = 0
# Path to geostrophic wind file (relative to run directory)
GeoWindFile = 'input/GeoWind.npz'
# Coriolis parameter (1/s)
f_coriolis = 0.0001
# Potential temperature lapse rate above domain top (K/m)
inversion = 0/1000
# Buoyancy calculation: 0 = use reference T_0, 1 = use local THv
optBuoyancy = 1
# Reference temperature (K)
T_0 = 300
# ------------------------------------------------------------
# Subgrid-scale configuration
# ------------------------------------------------------------
# SGS model: 1 = LASDD-SM, 2 = LASDD-WL, 3 = LAD-SM, 4 = LAD-WL
optSgs = 1
# Dynamic SGS update frequency (every N steps)
dynamicSGS_call_time = 1
# Filter to grid ratio (FGR=1: implicit + dealiasing; FGR>=2: explicit)
FGR = 2
# Initial SGS coefficients (used before first dynamic update)
Cs2 = 0.1 ** 2 # SM models: initial Cs^2
Cwl = 0.1 ** 2 # WL models: initial C_WL
Cs2PrRatio = Cs2 / 1.0
CwlPrRatio = Cwl / 1.0
# ------------------------------------------------------------
# Damping layer configuration
# ------------------------------------------------------------
optDamping = 1 # 1: activate Rayleigh damping
z_damping = 2500 # unit: m
RelaxTime = 300 # unit: s
# ------------------------------------------------------------
# Statistics computation
# ------------------------------------------------------------
SampleInterval_sec = 10.0 # collect a sample every N s
OutputInterval_sec = 60.0 # output averaged stats every N s
Output3DInterval_sec = SimTime # output 3D fields only at the end of the run
# ------------------------------------------------------------
# Large-scale advection forcing
# ------------------------------------------------------------
# 0: none, 1: time/height-varying (from AdvectionFile)
optAdvection = 0
AdvectionFile = 'input/AdvForcing.npz'
# ------------------------------------------------------------
# Moisture configuration
# ------------------------------------------------------------
# 0: dry run, 1: prognostic specific humidity Q
optMoisture = 0
# Screen-level moisture reference height (m); 0 = use z0T
zMoisture = 0.0
# Surface moisture flux (kg/kg m/s); used when optMoistureSurfBC = 0
MoistureFlux = 0.0
# 0: constant flux, 1: time-varying flux, 2: time-varying surface Q
optMoistureSurfBC = 0
MoistureSurfaceBCFile = 'input/MoistureSurfaceBC.npz'
# Specific humidity lapse rate above domain top (kg/kg/m); 0 = zero gradient
q_inversion = 0.0
# Pressure solver: 0 = LU (original), 1 = Thomas (tridiagonal, faster)
optPressureSolver = 1