From 5f57067db3ddae5f0f76a2c37af552fc50892f3c Mon Sep 17 00:00:00 2001 From: gdecker1 Date: Tue, 11 Jul 2023 11:12:16 -0700 Subject: [PATCH] Added Bonan Fig 6.3 to docs * Added Bonan phi and psi plots to docs * Add Bonan to UniversalFunctions.md * Added suggested description to docs * Removed reference to assets folder --- docs/src/UniversalFunctions.md | 9 ++++++++- docs/src/plot_universal_functions.jl | 7 ++++++- 2 files changed, 14 insertions(+), 2 deletions(-) diff --git a/docs/src/UniversalFunctions.md b/docs/src/UniversalFunctions.md index 5db34e3d..354d3b2e 100644 --- a/docs/src/UniversalFunctions.md +++ b/docs/src/UniversalFunctions.md @@ -1,6 +1,6 @@ # Universal Functions -UniversalFunctions.jl provides universal functions for SurfaceFluxes.jl. Here, we reproduce some plots from literature, specifically from Gryanik et al. 2020, and Businger. +UniversalFunctions.jl provides universal functions for SurfaceFluxes.jl. The functions are defined in `src/UniversalFunctions.jl` and the plots are generated in `docs/src/plot_universal_functions.jl`. Here, we reproduce some plots from literature, specifically from Gryanik et al. 2020, Businger 1971, and Bonan 2019. Note that Bonan uses the forms of $\phi$ and $\psi$ from Dyer and Hicks 1970; Dyer 1974; Brutsaert 1982, pp. 68–71; Garratt 1992, pp. 52–54. ```@example include("plot_universal_functions.jl") @@ -24,3 +24,10 @@ include("plot_universal_functions.jl") ![](Businger_phi_h.svg) ![](Businger_phi_m.svg) + +## Figs 1,2 (Bonan) + +![](Bonan_phi_h.svg) +![](Bonan_phi_m.svg) +![](Bonan_psi_h.svg) +![](Bonan_psi_m.svg) diff --git a/docs/src/plot_universal_functions.jl b/docs/src/plot_universal_functions.jl index 85e535d8..5173898a 100644 --- a/docs/src/plot_universal_functions.jl +++ b/docs/src/plot_universal_functions.jl @@ -59,10 +59,10 @@ function save_ψ_figs(ufts, ζ; ylims = nothing, fig_prefix = "", xaxis = :ident end +# Gryanik Plots save_ϕ_figs(ufts, FT(0):FT(0.01):FT(15); ylims = (0, 30), fig_prefix = "Gryanik12") save_ψ_figs(ufts, FT(0):FT(0.01):FT(15); ylims = (-25, 0), fig_prefix = "Gryanik12") - save_ϕ_figs( ufts, 10 .^ (FT(-3):0.1:FT(2)); @@ -73,4 +73,9 @@ save_ϕ_figs( ) +# Businger Plots save_ϕ_figs(ufts, FT(-2.5):FT(0.01):FT(2); ylims = (-1, 8), fig_prefix = "Businger") + +# Bonan Plots +save_ϕ_figs(ufts, FT(-2):FT(0.01):FT(1); ylims = (0, 6), fig_prefix = "Bonan") +save_ψ_figs(ufts, FT(-2):FT(0.01):FT(1); ylims = (-5, 4), fig_prefix = "Bonan")