Diurnal Cycle over Wangara, Australia: Description
Last updated: May 2026
References
Clarke, R. H., Dyer, A. J., Brook, R. R., Reid, D. G., and Troup, A. J. (1971). The Wangara experiment: Boundary layer data. Division of Meteorological Physics Tech. Paper 19, 362 pp.
Yamada, T., and Mellor, G. (1975) A simulation of the Wangara atmospheric boundary layer data. Journal of the Atmos. Sciences, 32, 2309–2329.
Basu, S., Vinuesa, J.-F., and Swift, A. (2008). Dynamic LES modeling of a diurnal cycle. Journal of Applied Meteorology and Climatology, 47, 1156–1174.
Case Description
The Wangara diurnal-cycle case simulates a full 24-hour period over the Wangara field site in Australia, spanning daytime convective growth and nocturnal stable-layer collapse. The surface temperature and geostrophic forcing are time dependent.
Parameter |
Value |
|---|---|
Reference run |
|
Domain |
5000 m x 5000 m x 2000 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 |
86400 s (24 h) |
Reference time step |
0.5 s |
Coriolis parameter |
\(f = -8.26 \times 10^{-5}\) s\(^{-1}\) |
Roughness lengths |
\(z_{0m} = 0.01\) m, \(z_{0T} = 0.01\) m |
Thermal surface option |
|
Reference temperature |
\(T_0 = 278.5\) K |
Inversion strength |
0.001 K m\(^{-1}\) |
SGS model |
selected by the user |
Filter-to-grid ratio |
|
Damping layer |
above \(z = 1500\) m |
Output statistics interval |
60 s |
3D output interval |
SimTime s |
Geostrophic forcing option |
|
Advection forcing option |
|
Moisture option |
|
Surface boundary condition |
Observed screen-level temperature |
Reference output period |
24 h |
The representative configuration below is copied from examples/DC_Wangara/runs/128x128x128_LASDD_SM_DP/Config.py.
Reference Configuration
[ ]:
# 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: Wangara diurnal cycle after Basu et al. (2008).
Domain 5000x5000x2000 m; 24-hour simulation.
Cabauw-style geostrophic wind (time + height varying).
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 = 5000
l_y = 5000
l_z = 2000
# 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.5 # unit: sec
SimTime = 86400 # unit: sec
# Galilean transformation (m/s)
Ugal = 0
# ------------------------------------------------------------
# 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 = 2
# Roughness lengths (m)
z0m = 0.01
z0T = 0.01
# Screen-level temperature reference height (m); 0 = use z0T
zTemperature = 1.2
# 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 = 1
# Constant geostrophic wind (m/s)
Ug2 = -5.34
Vg2 = -0.43
# Path to geostrophic wind file (relative to run directory)
GeoWindFile = 'input/GeoWind.npz'
# Coriolis parameter (1/s)
f_coriolis = -8.26e-05
# Potential temperature lapse rate above domain top (K/m)
inversion = 0.001
# Buoyancy calculation: 0 = use reference T_0, 1 = use local THv
optBuoyancy = 1
# Reference temperature (K)
T_0 = 278.5
# ------------------------------------------------------------
# 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 = 1500 # unit: m
RelaxTime = 600 # 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