Skip to content

Commit 2eb4c49

Browse files
committed
more doc updates
1 parent 080e32d commit 2eb4c49

File tree

1 file changed

+14
-7
lines changed

1 file changed

+14
-7
lines changed

docs/src/perturbed_dycoms_rf02.md

Lines changed: 14 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -115,14 +115,20 @@ Finally, we need to update type_getters.jl with our new keyword arguments as wel
115115
With this step complete, we are now ready to pass a new keyword argument that we defined ourselves to modify or change the initial conditions. The scripts currently contained in /examples/perturbed_dycoms_rf02 contain code for running the ensemble as well as post-processing/plotting capabilities.
116116

117117
# Creating & Running Ensembles
118+
All example scripts referenced are stored in examples/perturbed_dycoms_rf02:
119+
120+
```
121+
script_folder = joinpath(pkgdir(CA), "examples", "perturbed_dycoms_rf02")
122+
```
123+
118124
The script `pipeline.jl` handles this entire step without parallelized runs, but the steps can also be performed individually.
119125

120126
The recommended workflow is the create multiple copies of the .yml files and label them according to the initial conditions used. For example, the default file might look like `prognostic_edmfx_dycoms_rf02_column_qtot0_9.45_theta0_288.3_thetai_295.0_zi_795.0_prescribedN_1.0e8.yml`.
121127

122128
We can use `make_yaml.jl` to generate multiple different configurations with ease. To do so, provide the function `make_yamls()` with a default file for which to copy, and an output path. This script uses nested for loops to generate every possible combination of parameters, and they can be altered directly in the file. This might look like:
123129

124130
```
125-
include(joinpath(pkgdir(CA),"LWP_N_scripts", "make_yaml.jl"))
131+
include(joinpath(script_folder, "make_yaml.jl"))
126132
127133
default_data_path = default_prog_2M
128134
output_dir = prognostic_2M_config
@@ -132,20 +138,21 @@ make_yamls(default_data_path, output_dir, is_prog=true)
132138
To use these configuration files, we can either parallelize or run them serially. If parallelization is desired, run:
133139

134140
```
135-
julia -p <num_processors> parallel_driver.jl
141+
julia -p <num_processors> ./examples/perturbed_dycoms_rf02/parallel_driver.jl
136142
```
137143

138144
Be sure that the paths, which are handled by `all_paths.jl`, are referring to the right set of simulations. This example focuses in on prognostic EDMF+2M simulations, but the scripts have implementation for diagnostic EDMF as well as 1-moment microphysics as well.
139145

140146
# Visualizing Ensemble Output
141-
After the simulations are done being ran, we can use `process_plot_outputs.jl` to visualize changes in the liquid water path (LWP) and cloud droplet number concentration (N), in order to recreate the figure above. An example code block to recreate the plot is shown below, but additionally plotting capabilities are demonstrated in `LWP_N_Experiments.ipynb`.
147+
After the simulations are done being ran, we can use `process_plot_outputs.jl` to visualize changes in the liquid water path (LWP) and cloud droplet number concentration (N), in order to recreate the figure above.
142148

143149
```
144-
include(joinpath(pkgdir(CA), "LWP_N_scripts", "process_plot_outputs.jl"))
145-
output_dir = output_2M
150+
include(joinpath(script_folder, "process_plot_outputs.jl"))
151+
152+
output_dir = "./output" #output_2M
146153
all_outputs = make_edmf_vec(output_dir)
147154
filtered = filter_runs(all_outputs)
148-
sampled = StatsBase.sample(filtered, 40; replace = false)
155+
# sampled = StatsBase.sample(filtered, 40; replace=false)
149156
150-
fig = plot_edmf(sampled, is_1M = false, is_time = false, save = false)
157+
fig = plot_edmf(filtered, is_1M=false, is_time=false)
151158
```

0 commit comments

Comments
 (0)