Notation

This page summarizes the mathematical and code notation used in NumericalEarth.jl, following the conventions established in Breeze.jl.

How the notation works

Variable names are built by combining a base symbol with superscripts and, occasionally, a short plain-text tag.

Base symbols are single characters (often script letters) that identify the physical category of a quantity — for example, 𝒬 for heat flux, ℐ for radiative intensity, J for mass flux, and τ for kinematic momentum flux.

Superscripts refine the meaning in several ways:

  • Phase or species: ᵛ (vapor), ˡ (liquid), ⁱ (ice), ᶜ (condensate)
  • Component: ᵃᵗ (atmosphere), ᵒᶜ (ocean), ˢⁱ (sea ice), ˡᵃ (land)
  • Interface pair: ᵃᵒ (atm–ocean), ᵃⁱ (atm–ice), ⁱᵒ (ice–ocean)
  • Direction: ˣ / ʸ (spatial), ˢʷ / ˡʷ (shortwave / longwave)
  • Process: ⁱⁿ (interface), ᶠʳᶻ (frazil)

Component superscripts are used only in cross-component context — wherever a variable appears alongside variables from other components, as in interface computations or coupled-model discussions (Tˡᵃ next to Tᵃᵗ and Tⁱⁿ). Within a single component's own namespace the bare symbol is used: the land model's prognostic state is (; T, M) (as returned by prognostic_fields(land)), just as the ocean's is (; u, v, w, T, S).

Modifier arrows ꜜ (\^downarrow) and ꜛ (\^uparrow) denote downwelling and upwelling directions in radiative fluxes.

Subscripts encode radiative process (ₜ transmitted, ₐ absorbed) and the similarity-theory scale ★.

For example, 𝒬ᵛ is the latent (vapor) heat flux, ℐꜜˢʷ is the downwelling shortwave radiative intensity, and τˣ is the zonal kinematic momentum flux.

In Julia code, superscripts are entered with Unicode (e.g. \scrQ<tab> → 𝒬, then \^v<tab> → ᵛ). The modifier arrows ꜜ and ꜛ are entered with \^downarrow<tab> and \^uparrow<tab>.

Base symbols

MathCodeTab completionMeaning
$\mathcal{Q}$𝒬\scrQHeat flux (W m⁻²)
$\mathscr{I}$ℐ\scrIRadiative intensity (W m⁻²)
$J$JMass flux (kg m⁻² s⁻¹)
$\tau$τ\tauKinematic momentum flux (m² s⁻²)
$\mathcal{L}$ℒ\scrLLatent heat (J kg⁻¹)
$M$MLayer-integrated mass per area (kg m⁻²)

Note: $\tau^x$ (τˣ) is the kinematic momentum flux (stress divided by density). The mass-weighted stress is $\rho \tau^x$ (ρτˣ, in N m⁻²).

These base symbols are combined with superscript and subscript labels (documented below) to form specific variable names.

Superscript and subscript labels

Superscripts and subscripts are used systematically to label physical quantities. Superscripts generally denote the type or phase of a quantity, while subscripts denote the component or location.

Superscript labels

LabelCodeMeaningExample
$v$ᵛwater vapor$\mathcal{Q}^v$ (latent heat flux)
$T$ᵀtemperature / sensible$\mathcal{Q}^T$ (sensible heat flux)
$\mathrm{rn}$ʳⁿrain$J^{\mathrm{rn}}$ (rainfall)
$\mathrm{sn}$ˢⁿsnow$J^{\mathrm{sn}}$ (snowfall)
$S$ˢsalinity$J^S$ (salinity flux)
$w$ʷfreshwater$J^w$ (freshwater volume flux per unit area)
$h$ʰheight (roughness element)$\sigma^h$ (height standard deviation)
$i$ⁱice$\mathcal{L}^i$ (latent heat of sublimation)
$\ell$ˡliquid$\mathcal{L}^\ell$ (latent heat of vaporization)
$p$ᵖconstant pressure$c^{pm}$ (moist isobaric heat capacity)
$m$ᵐmixture (moist air)$c^{pm}$ (moist isobaric heat capacity)
$d$ᵈdry (air)$c^{pd}$ (dry air heat capacity)
$D$ᴰdrag$C^D$ (drag coefficient)
$\mathrm{in}$ⁱⁿinterface$T^{\mathrm{in}}$ (interface temperature)
$\mathrm{frz}$ᶠʳᶻfrazil$\mathcal{Q}^{\mathrm{frz}}$ (frazil heat flux)
$x$ˣzonal / x-direction$\tau^x$ (zonal kinematic stress)
$y$ʸmeridional / y-direction$\tau^y$ (meridional kinematic stress)
$\mathrm{at}$ᵃᵗatmosphere$\rho^{\mathrm{at}}$ (air density)
$\mathrm{ae}$ᵃᵉair entry$\alpha^{\mathrm{ae}}$ (inverse air-entry head)
$\mathrm{oc}$ᵒᶜocean$\rho^{\mathrm{oc}}$ (ocean reference density)
$\mathrm{si}$ˢⁱsea ice$h^{\mathrm{si}}$ (sea ice thickness)
$\mathrm{la}$ˡᵃland$M^{\mathrm{la}}$ (land water mass per area)
$\mathrm{ao}$ᵃᵒatmosphere–ocean interface$\mathcal{Q}^{\mathrm{ao}}$ (atm–ocean heat flux)
$\mathrm{ai}$ᵃⁱatmosphere–ice interface$\mathcal{Q}^{\mathrm{ai}}$ (atm–ice heat flux)
$\mathrm{io}$ⁱᵒice–ocean interface$\mathcal{Q}^{\mathrm{io}}$ (ice–ocean heat flux)
$\mathrm{sw}$ˢʷshortwave$\mathscr{I}$ ꜜ ${}^{\mathrm{sw}}$ (downwelling shortwave)
$\mathrm{lw}$ˡʷlongwave$\mathscr{I}$ ꜜ ${}^{\mathrm{lw}}$ (downwelling longwave)

Modifier arrows

SymbolCodeTab completionMeaning
ꜜꜜ\^downarrowdownwelling
ꜛꜛ\^uparrowupwelling

Subscript labels

LabelCodeMeaningExample
$t$ₜtransmitted$\mathscr{I}_{t}^{\mathrm{sw}}$ (transmitted shortwave)
$a$ₐabsorbed$\mathscr{I}_{a}^{\mathrm{lw}}$ (absorbed longwave)
$\star$★similarity theory scale$u_\star$ (friction velocity)
$b$bboundary value, or background$T_b$ (interface temperature), $u_b$ (background bottom velocity)

Atmosphere state variables

MathCodePropertyDescription
$T$TtemperatureAir temperature (K)
$p$ppressureAir pressure (Pa)
$q$qspecific humidityMass mixing ratio of water vapor (kg kg⁻¹)
$u$uzonal velocityEastward wind component (m s⁻¹)
$v$vmeridional velocityNorthward wind component (m s⁻¹)
$\mathscr{I}_\downarrow^{\mathrm{sw}}$ℐꜜˢʷdownwelling shortwaveDownwelling shortwave radiation (W m⁻²)
$\mathscr{I}_\downarrow^{\mathrm{lw}}$ℐꜜˡʷdownwelling longwaveDownwelling longwave radiation (W m⁻²)
$h_{b\ell}$h_bℓboundary layer heightAtmospheric boundary layer height (m)
$pᵛ⁺$$pᵛ⁺$saturation vapor pressureVapor pressure at saturation (Pa)
$qᵛ⁺$qᵛ⁺saturation specific humiditySpecific humidity at saturation, $q^{v+}(T)$ (kg kg⁻¹)
$qˢ$qsurface specific humiditySpecific humidity at the interface; set by the humidity model (β·qᵛ⁺ for FractionalHumidity, a vapor-flux balance for SkinHumidity) (kg kg⁻¹)

Land state variables and parameters

Bare symbols below are the land model's internal names; in cross-component context they take the ˡᵃ superscript (Tˡᵃ, Mˡᵃ, Mˡᵃ⁺) per the component-superscript rule above.

MathCodePropertyDescription
$T$temperatureground temperaturePrognostic land-column temperature (K)
$M$water_storageland waterPrognostic land water mass per area (kg m⁻²)
$M^{+}$maximum_water_storagemaximum land waterBucket capacity; soil-science "field capacity" (kg m⁻²)
$𝒮$saturationsurface saturationContinuous land surface saturation $\mathrm{clamp}(M/M⁺, 0, 1)$; the interface humidity models derive their availability $β$ from it (–)
$𝒮ᶜ$critical_saturationcritical saturationSaturation above which the surface evaporates at full efficiency, for CriticalSaturation (–)
$𝒮ᶜ$dry_layer_onset_saturationdry-layer onset saturationSaturation below which a dry surface layer forms, for StorageBasedDryLayerDepth; shares the symbol $𝒮ᶜ$ with critical_saturation above (–)
$T^{\mathrm{deep}}$deep_temperaturedeep climatological temperaturePrescribed deep/climatological target temperature for force-restore (K)
$τ^{\mathrm{deep}}$deep_time_scaledeep-restore time scaleTime scale of surface relaxation toward $T^{\mathrm{deep}}$ (s)
$c^{\mathrm{dry}}$dry_heat_capacitydry areal heat capacityAreal heat capacity of the water-free slab; a Number or an AbstractField (J m⁻² K⁻¹)
$c^l$liquid_heat_capacityliquid heat capacitySpecific heat capacity of slab liquid water (J kg⁻¹ K⁻¹)
$c^{\mathrm{la}}$–land areal heat capacityWater-dependent areal heat capacity $c^{\mathrm{la}} = c^{\mathrm{dry}} + c^l M$ (J m⁻² K⁻¹)
$d$surface_thicknesssurface thicknessThickness of the dry surface layer through which soil vapor diffuses, for SkinHumidity (m)
$κ^q$vapor_diffusivitysoil vapor diffusivityVapor mass diffusivity in the surface soil layer, for SkinHumidity (kg m⁻¹ s⁻¹)
$\chi^{\mathrm{sand}}$sandsoil sand fractionMass fraction of sand grains in the mineral (non-organic) solid matrix (kg kg⁻¹)
$\chi^{\mathrm{silt}}$siltsoil silt fractionMass fraction of silt grains in the mineral (non-organic) solid matrix (kg kg⁻¹)
$\chi^{\mathrm{clay}}$claysoil clay fractionMass fraction of clay grains in the mineral (non-organic) solid matrix (kg kg⁻¹)
$\chi^{\mathrm{soc}}$SOCsoil organic carbon concentrationMass fraction of organic carbon in the solid matrix (kg kg⁻¹)
$\rho^{\mathrm{soil}}$ρ_soilsoil bulk dry densityBulk dry density of the soil within each vertical layer (kg m⁻³)
$\rho^{\mathrm{soc}}$ρ_socsoil organic carbon densityBulk density of organic material within each vertical layer (kg m⁻³)

Variably-saturated slab land

Symbols introduced by VariablySaturatedHydrology, WaterCoupledEnergy, and DryLayerHumidity.

MathCodePropertyDescription
$h^{\mathrm{la}}$slab_depthdepth of prognostic landVertical thickness of the integrated land slab, from $z_b$ to $z_s$ (m)
$\nu$porositysoil porosityTotal pore fraction (–)
$\theta^l$–pore liquid fractionPhysical liquid-filled pore fraction; surface physics consumes this (–)
$\vartheta^l$–augmented liquid fractionConservative storage variable $= \theta^l + \max(\Pi, 0)/h^{\mathrm{ss}}$; allows $M > M⁺$ saturated overflow (–)
$\theta^r$residual_liquid_fractionresidual liquid fractionMinimum liquid-filled pore fraction (–)
$𝒮$saturationeffective saturationEffective (relative) saturation $𝒮 = \mathrm{clamp}\!\left((\theta^l - \theta^r)/(\nu - \theta^r),\, 0,\, 1\right)$; the humidity availability and the front depth $\delta^v$ derive from it (–)
$h^{\mathrm{ss}}$storage_heightstorage heightSaturated storage height — the head built per unit fractional over-saturation; reciprocal of the specific storage ($1/S_s$) (m)
$\Pi$–soil pressure headMatric/pressure head; $\Pi \le 0$ unsaturated, $\Pi > 0$ saturated overflow (m)
$\Pi^d$deep_pressure_headdeep pressure headPressure head of the deep reservoir below the slab, passed to the deep-flux closure (m)
$h$–hydraulic head$h = z + \Pi$ (m)
$K$–hydraulic conductivityDarcy conductivity (m s⁻¹)
$K_0$matching_point_conductivitymatching-point conductivityConductivity the Mualem–van Genuchten curve reaches at $𝒮 = 1$; a matrix value when it comes from a pedotransfer function, excluding macropore flow (m s⁻¹)
$K^{+}$saturated_conductivitysaturated hydraulic conductivityMacropore-inclusive Darcy conductivity of saturated soil, for an infiltration cap (m s⁻¹)
$\alpha^{\mathrm{ae}}$inverse_air_entry_headinverse air-entry headReciprocal of the air-entry pressure head, the suction at which the largest pores begin to empty (m⁻¹)
$\mathscr{n}$pore_size_uniformitypore-size uniformityvan Genuchten retention exponent; large $\mathscr{n}$ drains over a narrow band of suctions, $\mathscr{n} \to 1$ over orders of magnitude. Sets $\mathscr{m} = 1 - 1/\mathscr{n}$ (–)
$\eta^K$ηᴷ, pore_connectivity_exponentpore-connectivity exponentMualem exponent on saturation in $K(\mathcal S)$ (–)
$T_0$reference_temperatureviscosity reference temperatureTemperature at which $K_0$ holds and the viscosity factor $\Theta(T)$ on $K$ is unity (K)
$J^{Es}$surface_energy_fluxsurface energy fluxSigned surface energy flux, positive upward (out of the slab) (W m⁻²)
$J^{lb}$deep_liquid_fluxdeep-boundary liquid fluxLiquid mass flux across the slab bottom, positive upward (into the slab, capillary rise / groundwater return); drainage is $J^{lb} < 0$ (kg m⁻² s⁻¹)
$J^{ls}$surface_liquid_fluxsurface liquid fluxLiquid mass flux at the surface $J^{ls} = -P^l + R^{\mathrm{sfc}}$, positive upward (out of the slab); infiltration is $J^{ls} < 0$ (kg m⁻² s⁻¹)
$R^{\mathrm{sfc}}$surface_runoffsurface runoffLiquid input rejected at the surface, $\ge 0$ (kg m⁻² s⁻¹)
$R^{\mathrm{lat}}$subsurface_runoffsubsurface runoffLateral storage export, $\ge 0$ (kg m⁻² s⁻¹)
$\kappa^T$thermal_conductivitythermal conductivityEffective ground thermal conductivity (W m⁻¹ K⁻¹)
$\Lambda^{\mathrm{deep}}$deep_conductancedeep energy conductanceForce-restore deep energy conductance (W m⁻² K⁻¹); see also $τ^{\mathrm{deep}}$
$T_r$reference_temperaturereference temperatureReference temperature for internal energy $e^l(T) = c^l (T - T_r)$ (K)
$T^{\mathrm{in}}$–interface temperatureAtmosphere-facing skin temperature, $T^{\mathrm{in}}$ (K)
$q^{\mathrm{in}}$–interface specific humidityAtmosphere-facing skin humidity, $q^{\mathrm{in}}$ (kg kg⁻¹)
$T^e$–dry-layer temperatureDiagnostic temperature at the dry layer (K)
$q^e$–dry-layer specific humidityVapor source humidity at the dry layer (kg kg⁻¹)
$\delta^v$dry_layer_depthdry-layer depthDry-layer thickness through which vapor diffuses, diagnostic of $𝒮$ (m)
$\chi$–blend coefficient$\chi = \mathrm{clamp}(\delta^v/\ell^T, 0, 1)$; weights $T^e$ between $T^{\mathrm{in}}$ and $T^{\mathrm{la}}$ (–)
$\eta$dry_layer_exponentfront-depth exponentExponent in $\delta^v = \delta^v_{max}[1 - \min(𝒮/𝒮^c, 1)]^\eta$ (–)
$\ell^T$thermal_exchange_depth, exchange_depththermal exchange depthDepth over which $\Lambda^{\mathrm{in}} = \kappa^T/\ell^T$ couples $T^{\mathrm{la}}$ to $T^{\mathrm{in}}$ (m)
$D^v$molecular_diffusivityvapor diffusivity in airMolecular vapor diffusivity in air (m² s⁻¹)
$w^d$–dry-layer piston velocity$w^d = D^v_{eff}/\max(\delta^v, \delta^v_{min})$ (m s⁻¹)

Ocean state variables

MathCodePropertyDescription
$T$TtemperatureOcean potential temperature (ᵒC or K)
$S$SsalinityPractical salinity (g kg⁻¹)
$u$uzonal velocityEastward ocean velocity (m s⁻¹)
$v$vmeridional velocityNorthward ocean velocity (m s⁻¹)
$u_b$ubbackground bottom velocityUnresolved velocity, mostly tidal, added in quadrature to the resolved speed in the bottom drag (m s⁻¹)
$\rho^{\mathrm{oc}}$ρᵒᶜreference densityOcean reference density (kg m⁻³)
$c^{\mathrm{oc}}$cᵒᶜheat capacityOcean heat capacity (J kg⁻¹ K⁻¹)

Sea ice state variables

MathCodePropertyDescription
$h^{\mathrm{si}}$hˢⁱice thicknessSea ice thickness (m)
$\aleph$ℵice concentrationAreal fraction of ice cover (–)
$S^{\mathrm{si}}$Sˢⁱice salinitySea ice bulk salinity (g kg⁻¹)

Radiation properties

MathCodePropertyDescription
$\sigma$σStefan–Boltzmann constant(W m⁻² K⁻⁴)
$\alpha$αalbedoSurface reflectivity (–)
$\epsilon$ϵemissivitySurface emissivity (–)

Similarity theory / surface layer

MathCodePropertyDescription
$u_\star$u★friction velocitySurface friction velocity (m s⁻¹)
$\theta_\star$θ★temperature scaleFlux characteristic temperature (K)
$q_\star$q★humidity scaleFlux characteristic specific humidity (kg kg⁻¹)
$b_\star$b★buoyancy scaleFlux characteristic buoyancy (m s⁻²)
$L_\star$L★Obukhov lengthMonin–Obukhov length scale (m)
$C^D$Cᴰdrag coefficientBulk transfer coefficient for momentum (–)
$\psi$ψstability functionIntegrated stability correction (–)
$\Psi$Ψinterface stateAggregate interface state (an AbstractInterfaceState) carried through the similarity-theory fixed-point solver compute_interface_state
$\zeta$ζstability parameter$z / L_\star$ (–)
$\Delta h^d$Δhᵈdisplaced profile heightHeight $\Delta h - d$ above the zero-plane displacement at which the similarity profiles are evaluated (m)
$\ell$ℓroughness lengthAerodynamic roughness length (m)
$\ell^\mathrm{m}$ℓᵐmomentum roughness lengthAerodynamic momentum roughness length (m)
$\ell^\mathrm{s}$ℓˢscalar roughness lengthAerodynamic scalar roughness length (m)
$\varkappa$ϰvon Kármán constant$\approx 0.4$ (–)

Note the case distinction: lowercase $\psi$ (ψ) is the stability function, while capital $\Psi$ (Ψ) is the aggregate interface-state object.

Surface morphometry

Roughness elements — buildings or plant canopy — that set the aerodynamic parameters $\ell^\mathrm{m}$ and $d$ of a land surface.

MathCodePropertyDescription
$h$helement heightMean height of the roughness elements (m)
$h^\mathrm{max}$hᵐᵃˣmaximum element heightTallest roughness element in the cell (m)
$\sigma^h$σʰheight standard deviationSpread of the roughness-element heights (m)
$\lambda^p h$gross_building_heightgross building heightMean element height over all cells, built or not, $\lambda^p h$ — the digital-surface lift (m)
$d$dzero-plane displacementDisplacement height of the roughness sublayer (m)
$d^h$dʰdisplacement ratioPacking displacement $d/h$ of the obstacle-array fit, before the height-spread correction (–)
$\lambda^p$λᵖplan-area indexGround-plan area covered by the elements, per unit ground area (–)
$\lambda^f$λᶠfrontal-area indexWindward face area of the elements, per unit ground area (–)
$\lambda^{p\,\mathrm{min}}$λᵖᵐⁱⁿplan-area index floorBelow this $\lambda^p$ a cell reduces to bare soil (–)
$\ell^\mathrm{soil}$ℓˢᵒⁱˡbare-soil roughness lengthMomentum roughness length where the plan-area index vanishes (m)
$A$Aarray constantPacking constant of the Macdonald displacement fit (–)
$\beta$βdrag correction factorLumped drag correction of the Macdonald drag partition (–)
$X$Xdisplacement parameter$(\sigma^h + h)/h^\mathrm{max} \in [0, 1]$, argument of the Kanda displacement fit (–)
$Y$Yheight-spread parameter$\lambda^p \sigma^h / h$, argument of the Kanda roughness rescaling (–)

The superscripts on $\lambda^p$ and $\lambda^f$ distinguish them from the longitude $\lambda$.

Two symbols here are shared with other land parameters and mean something different. $\beta$ is the drag correction of the Macdonald drag partition, not the moisture availability $\beta$ of the hydrology and interface-humidity models; it keeps the symbol of its source parametrization. Likewise $d$ is the zero-plane displacement, not the surface_thickness $d$ of SkinHumidity.

Radiative fluxes

MathCodePropertyDescription
$\mathscr{I}_\downarrow^{\mathrm{sw}}$ℐꜜˢʷdownwelling shortwaveDownwelling shortwave radiation (W m⁻²)
$\mathscr{I}_\downarrow^{\mathrm{lw}}$ℐꜜˡʷdownwelling longwaveDownwelling longwave radiation (W m⁻²)
$\mathscr{I}_\uparrow^{\mathrm{lw}}$ℐꜛˡʷupwelling longwaveEmitted longwave radiation (W m⁻²)

| $\mathscr{I}_{t}^{\mathrm{sw}}$ | ℐₜˢʷ | transmitted shortwave | Shortwave passing through the surface, $(1-\alpha) \mathscr{I}_\downarrow^{\mathrm{sw}}$ (W m⁻²) | | $\mathscr{I}_{a}^{\mathrm{lw}}$ | ℐₐˡʷ | absorbed longwave | Longwave absorbed at the surface, $\epsilon \mathscr{I}_\downarrow^{\mathrm{lw}}$ (W m⁻²) |

Radiative fluxes use $\mathscr{I}$ (ℐ, for "intensity") with a modifier arrow (ꜜ/ꜛ for downwelling/upwelling) and superscript band (ˢʷ/ˡʷ). Derived radiative quantities use a subscript process label (ₜ, ₐ) with a superscript band.

Turbulent interface fluxes

MathCodePropertyDescription
$\mathcal{Q}^v$𝒬ᵛlatent heat fluxTurbulent latent heat flux (W m⁻²)
$\mathcal{Q}^T$𝒬ᵀsensible heat fluxTurbulent sensible heat flux (W m⁻²)
$J^v$Jᵛwater vapor fluxTurbulent mass flux of water vapor (kg m⁻² s⁻¹)
$\tau^x$τˣzonal kinematic stressKinematic zonal momentum flux (m² s⁻²)
$\tau^y$τʸmeridional kinematic stressKinematic meridional momentum flux (m² s⁻²)
$\rho \tau^x$ρτˣzonal wind stressMass-weighted zonal stress (N m⁻²)
$\rho \tau^y$ρτʸmeridional wind stressMass-weighted meridional stress (N m⁻²)

Net ocean fluxes

MathCodePropertyDescription
$J^T$Jᵀtemperature fluxNet ocean temperature flux (K m s⁻¹)
$J^S$Jˢsalinity fluxNet ocean salinity flux (g kg⁻¹ m s⁻¹)
$\mathcal{Q}^{\mathrm{frz}}$𝒬ᶠʳᶻfrazil heat fluxHeat released by frazil ice formation (W m⁻²)

Net surface freshwater fluxes

MathCodePropertyDescription
$J^{\mathrm{rn}}$Jʳⁿrain freshwater fluxRain mass flux at the surface (kg m⁻² s⁻¹)
$J^{\mathrm{sn}}$Jˢⁿsnow freshwater fluxSnow mass flux at the surface (kg m⁻² s⁻¹)
$J^{\mathrm{ib}}$Jⁱᵇiceberg freshwater fluxIceberg calving mass flux into the ocean (kg m⁻² s⁻¹)

Thermodynamic properties

MathCodePropertyDescription
$\mathcal{L}^\ell$ℒˡlatent heat of vaporizationLiquid-phase latent heat (J kg⁻¹)
$\mathcal{L}^i$ℒⁱlatent heat of sublimationIce-phase latent heat (J kg⁻¹)
$\mathcal{L}^f$ℒᶠlatent heat of fusionHeat the ocean supplies to melt snowfall and icebergs (J kg⁻¹)
$c^{pm}$cᵖᵐmoist air heat capacityMoist isobaric specific heat (J kg⁻¹ K⁻¹)
$c^{pd}$cᵖᵈdry air heat capacityDry-air isobaric specific heat (J kg⁻¹ K⁻¹)
$\rho^{\mathrm{at}}$ρᵃᵗair densityAtmospheric air density (kg m⁻³)
$\varepsilon^{\mathrm{dv}}$$εᵈᵛ$vapor / dry-air gas-constant ratio$εᵈᵛ = R_v / R_d$ (so $(εᵈᵛ)^{-1} = R_d / R_v ≈ 0.622$ is the conventional ε in $q = ε e / p$) (–)

CF standard name mapping

The following table maps code variable names to their CF standard names where applicable.

CodeCF standard name
T (atm)air_temperature
T (ocn)sea_water_potential_temperature
Ssea_water_practical_salinity
u (atm)eastward_wind
v (atm)northward_wind
qspecific_humidity
pair_pressure
ℐꜜˢʷsurface_downwelling_shortwave_flux_in_air
ℐꜜˡʷsurface_downwelling_longwave_flux_in_air
𝒬ᵛsurface_upward_latent_heat_flux
𝒬ᵀsurface_upward_sensible_heat_flux
Jᵛwater_evapotranspiration_flux
ρτˣsurface_downward_eastward_stress
ρτʸsurface_downward_northward_stress
hˢⁱsea_ice_thickness
ℵsea_ice_area_fraction

Typing Unicode symbols in Julia

Most symbols can be entered in the Julia REPL and in editors with Julia support by typing a LaTeX-like abbreviation followed by <tab>. The table below collects the less obvious completions used in this notation.

SymbolTab completionDescription
𝒬\scrQScript Q (heat flux)
ℐ\scrIScript I (radiative intensity)
ℒ\scrLScript L (latent heat)
𝓃\scrnScript n (pore-size uniformity)
𝓂\scrmScript m (van Genuchten 1 - 1/𝓃)
τ\tauTau (kinematic stress)
ρ\rhoRho (density)
σ\sigmaSigma (Stefan–Boltzmann constant)
α\alphaAlpha (albedo)
ϵ\epsilonEpsilon (emissivity)
ℵ\alephAleph (ice concentration)
ϰ\varkappaVarkappa (von Kármán constant)
Ψ\PsiCapital Psi (interface state)
★\bigstarStar (similarity-theory scale)
ꜜ\^downarrowModifier down arrow (downwelling)
ꜛ\^uparrowModifier up arrow (upwelling)
ᵛ\^vSuperscript v
ᵀ\^TSuperscript T
ˢ\^sSuperscript s
ʷ\^wSuperscript w
ʰ\^hSuperscript h
ⁱ\^iSuperscript i
ˡ\^lSuperscript l
ᵖ\^pSuperscript p
ᵐ\^mSuperscript m
ᵈ\^dSuperscript d
ᴰ\^DSuperscript D
ˣ\^xSuperscript x
ʸ\^ySuperscript y
ᵃ\^aSuperscript a
ᵗ\^tSuperscript t
ᵒ\^oSuperscript o
ᶜ\^cSuperscript c
ⁿ\^nSuperscript n
ᶠ\^fSuperscript f
ʳ\^rSuperscript r
ᶻ\^zSuperscript z
ₜ\_tSubscript t (transmitted)
ₐ\_aSubscript a (absorbed)
ₚ\_pSubscript p (penetrating)