LES Intercomparison Study for Convective Boundary Layers: Description
Last updated: May 2026
Reference
Nieuwstadt, F. T., Mason, P. J., Moeng, C. H., & Schumann, U. (1993). Large-eddy simulation of the convective boundary layer: A comparison of four computer codes. In Turbulent Shear Flows 8: Selected Papers from the Eighth International Symposium on Turbulent Shear Flows, Munich, Germany, September 9–11, 1991 (pp. 343-367). Berlin, Heidelberg: Springer Berlin Heidelberg.
Case Description
This intercomparison case (henceforth CBL_N91) simulates a (free) convective boundary layer driven by a constant surface sensible heat flux. The flow develops a deepening, well-mixed convective layer with no geostrophic wind.
Parameter |
Value |
|---|---|
Reference run |
|
Domain |
6400 m x 6400 m x 5000 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 |
12600 s (3.5 h) |
Reference time step |
0.25 s |
Geostrophic wind |
\(U_g = 0\) m s\(^{-1}\), \(V_g = 0\) m s\(^{-1}\) |
Coriolis parameter |
\(f = 1.0 \times 10^{-4}\) s\(^{-1}\) |
Roughness lengths |
\(z_{0m} = 0.16\) m, \(z_{0T} = 0.16\) m |
Thermal surface option |
|
Sensible heat flux |
0.06 K m s\(^{-1}\) |
Reference temperature |
\(T_0 = 300\) K |
Inversion strength |
3 / 1000 K m\(^{-1}\) |
SGS model |
selected by the user |
Filter-to-grid ratio |
|
Damping layer |
above \(z = 3300\) m |
Output statistics interval |
60 s |
3D output interval |
SimTime s |
Geostrophic forcing option |
|
Advection forcing option |
|
Moisture option |
|
Surface forcing |
Constant sensible heat flux |
Reference averaging window |
hours 3-3.5 |
The representative configuration below is copied from examples/CBL_N91/runs/128x128x128_LASDD_SM_DP/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: Convective BL benchmark after Nieuwstadt et al. (1991).
Grid: 128x128x128, SGS: LASDD-SM,
Precision: double.
"""
# ============================================================
# Imports
# ============================================================
import numpy as np
# ============================================================
# User Input
# ============================================================
# ------------------------------------------------------------
# Platform options
# ------------------------------------------------------------
use_double_precision = True
# 0: use CPU, 1: use GPU
optGPU = 1
GPU_ID = 0
# ------------------------------------------------------------
# Domain configuration
# ------------------------------------------------------------
# Domain size (m)
l_x = 6400
l_y = 6400
l_z = 5000
# 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 = 0.25 # unit: sec
SimTime = 3.5*3600 # unit: sec
# Galilean transformation (m/s)
Ugal = 0
# ------------------------------------------------------------
# Surface configuration
# ------------------------------------------------------------
# optSurfFlux: 0 = homogeneous, 1 = heterogeneous
optSurfFlux = 1
# optSurfBC: 0 = constant flux
# 1 = time-varying flux (from SurfaceBCFile)
# 2 = time-varying surface temperature (from SurfaceBCFile)
optSurfBC = 0
# Roughness lengths (m)
z0m = 0.16
z0T = 0.16
# 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.06 # 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 = 0
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 = 3/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: 0 = Static SM, 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 = 3300 # 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