Skip to content
Snippets Groups Projects
Commit 423157f9 authored by Stefano Serafin's avatar Stefano Serafin
Browse files

improved documentation

parent d3af461f
No related branches found
No related tags found
No related merge requests found
%% Cell type:markdown id: tags: %% Cell type:markdown id: tags:
# References # References
For the data assimilation algorithms: For the data assimilation algorithms:
* Anderson, J. L., 2001: An ensemble adjustment Kalman filter for data assimilation. Mon. Weather Rev., 129, 2884–2903. [https://doi.org/10.1175/1520-0493(2001)129%3C2884:AEAKFF%3E2.0.CO;2] * Anderson, J. L., 2001: [An ensemble adjustment Kalman filter for data assimilation](https://doi.org/10.1175/1520-0493(2001)129%3C2884:AEAKFF%3E2.0.CO;2). Mon. Weather Rev., 129, 2884–2903.
* Anderson, J. L., 2003: A local least squares framework for ensemble filtering. Mon. Weather Rev., 131, 634–642. [https://doi.org/10.1175/1520-0493(2003)131%3C0634:ALLSFF%3E2.0.CO;2] * Anderson, J. L., 2003: [A local least squares framework for ensemble filtering](https://doi.org/10.1175/1520-0493(2003)131%3C0634:ALLSFF%3E2.0.CO;2). Mon. Weather Rev., 131, 634–642.
* Hunt, B. R., E. J. Kostelich, and I. Szunyogh, 2007: Efficient data assimilation for spatiotemporal chaos: A local ensemble transform Kalman filter. Physica D, 230, 112–126. * Hunt, B. R., E. J. Kostelich, and I. Szunyogh, 2007: [Efficient data assimilation for spatiotemporal chaos: A local ensemble transform Kalman filter](http://dx.doi.org/10.1016/j.physd.2006.11.008). Physica D, 230, 112–126.
[http://dx.doi.org/10.1016/j.physd.2006.11.008]
For the single-column model of the convective boundary layer (CBL): For the single-column model of the convective boundary layer (CBL):
* Troen, I. B., and L. Mahrt, 1986: A simple model of the atmospheric boundary layer; sensitivity to surface evaporation. Bound.-Layer Meteorol., 37, 129–148. [https://doi.org/10.1007/BF00122760] * Troen, I. B., and L. Mahrt, 1986: [A simple model of the atmospheric boundary layer; sensitivity to surface evaporation](https://doi.org/10.1007/BF00122760). Bound.-Layer Meteorol., 37, 129–148.
%% Cell type:markdown id: tags: %% Cell type:markdown id: tags:
# How to run a simple OSSE experiment # How to run a simple OSSE experiment
The default configuration files (in `json` format) are set to compute a single analysis that assimilates a single observation. For configuration details, see below. Simply run: The default configuration files (in `json` format) are set to compute a single analysis that assimilates a single observation. For configuration details, see below. Simply run:
%% Cell type:code id: tags: %% Cell type:code id: tags:
``` python ``` python
from ENDA import experiment from ENDA import experiment
import json import json
from PE_CBL_graphics import *
with open('./default_da.json', 'r') as fp: with open('./default_da.json', 'r') as fp:
settings = json.load(fp) settings = json.load(fp)
exp = experiment(settings) exp = experiment(settings)
``` ```
%% Output %% Output
Spinup : # ...done Spinup : # ...done
Nature run : # ...done Nature run : # ...done
Cycle 0000 : ########## ...done Cycle 0000 : ########## ...done
%% Cell type:markdown id: tags:
Now plot the results:
%% Cell type:code id: tags: %% Cell type:code id: tags:
``` python ``` python
from matplotlib import pyplot as p from matplotlib import pyplot as p
from PE_CBL_graphics import *
fig = p.figure(151) fig = p.figure(151)
fig.set_size_inches(6,4) fig.set_size_inches(6,4)
ax1 = fig.add_subplot(1,1,1) ax1 = fig.add_subplot(1,1,1)
ax1 = plot_CBL_assimilation(exp,_,ax = ax1) ax1 = plot_CBL_assimilation(exp,_,ax = ax1)
ax1.set_xlabel(r'$\theta$ (K)') ax1.set_xlabel(r'$\theta$ (K)')
ax1.set_ylabel(r'$z$ (m)') ax1.set_ylabel(r'$z$ (m)')
ax1.set_title(r'Ensemble state at $t=3600$ s') ax1.set_title(r'Ensemble state at $t=3600$ s')
ax1.legend(frameon=False) l=ax1.legend(frameon=False)
``` ```
%% Output %% Output
<matplotlib.legend.Legend at 0x1542d5630430>
%% Cell type:markdown id: tags: %% Cell type:markdown id: tags:
# Configuration of the CBL model # Configuration of the CBL model
See the file `default_cbl.json` See the file `default_cbl.json`
``` ```
{ {
"g": 9.806, "g": 9.806,
"f": 0.0001, "f": 0.0001,
"kvonk": 0.4, "kvonk": 0.4,
"z0": 0.1, "z0": 0.1,
"ug": 0, "ug": 0,
"vg": 10, "vg": 10,
"theta_0": 290, "theta_0": 290,
"gamma": 0.003, "gamma": 0.003,
"Hmax": 0.1, "Hmax": 0.1,
"H0_perturbation_ampl_init": 0.0, "H0_perturbation_ampl_init": 0.0,
"H0_perturbation_ampl_time": 0.0, "H0_perturbation_ampl_time": 0.0,
"exct": 0.3, "exct": 0.3,
"pfac": 1.5, "pfac": 1.5,
"Kmin": 0.1, "Kmin": 0.1,
"Kmax": 200, "Kmax": 200,
"is_bwind": false, "is_bwind": false,
"is_cgrad": true, "is_cgrad": true,
"is_cycle": false, "is_cycle": false,
"dz": 50, "dz": 50,
"ztop": 4000, "ztop": 4000,
"rnseed": 181612, "rnseed": 181612,
"perturb_ensemble_state": true, "perturb_ensemble_state": true,
"perturbations_type": "smooth", "perturbations_type": "smooth",
"perturbations_theta_amplitude": 0.1, "perturbations_theta_amplitude": 0.1,
"perturbations_uv_amplitude": 0.1, "perturbations_uv_amplitude": 0.1,
"perturbations_smooth_number": 11, "perturbations_smooth_number": 11,
"perturbations_symmetric": true, "perturbations_symmetric": true,
"error_growth_perturbations_amplitude": 0.1, "error_growth_perturbations_amplitude": 0.1,
"perturb_ensemble_parameters": true, "perturb_ensemble_parameters": true,
"parameter_number": 1, "parameter_number": 1,
"parameter_transform": "pit", "parameter_transform": "pit",
"parameter_ensemble_min": 0.5, "parameter_ensemble_min": 0.5,
"parameter_ensemble_max": 4.5, "parameter_ensemble_max": 4.5,
"parameter_true": 1.5, "parameter_true": 1.5,
"do_parameter_estimation": true, "do_parameter_estimation": true,
"parameter_inflation_rtps_alpha": 0.4, "parameter_inflation_rtps_alpha": 0.4,
"return_covariances_increments_and_innovations": true "return_covariances_increments_and_innovations": true
} }
``` ```
%% Cell type:markdown id: tags: %% Cell type:markdown id: tags:
# Configuration of the data assimilation cycle # Configuration of the data assimilation cycle
See the file `default_da.json` See the file `default_da.json`
``` ```
{ {
"cbl_settings_file": "./default_cbl.json", "cbl_settings_file": "./default_cbl.json",
"type": "OSSE", "type": "OSSE",
"tspinup": 10800, "tspinup": 10800,
"trun": 3600, "trun": 3600,
"tspinup_assim": 0, "tspinup_assim": 0,
"assimilation_interval": 3600, "assimilation_interval": 3600,
"randomize_obs": true, "randomize_obs": true,
"nobs": 1, "nobs": 1,
"obs_coordinates_min": 433, "obs_coordinates_min": 433,
"obs_coordinates_max": 433, "obs_coordinates_max": 433,
"obs_kinds": "theta", "obs_kinds": "theta",
"obs_error_sdev_generate": 0.2, "obs_error_sdev_generate": 0.2,
"obs_error_sdev_assimilate": 0.1, "obs_error_sdev_assimilate": 0.1,
"localization_cutoff": 400, "localization_cutoff": 400,
"nens": 50, "nens": 50,
"FILTER": "EAKF", "FILTER": "EAKF",
"inflation_rtps_alpha": 0.1, "inflation_rtps_alpha": 0.1,
"simulate_underdispersiveness": false, "simulate_underdispersiveness": false,
"rnseed": 181612, "rnseed": 181612,
"path": "./", "path": "./",
"filename": "experiment.pickle", "filename": "experiment.pickle",
"label": "X", "label": "X",
"nature_run_from_file": false, "nature_run_from_file": false,
"do_parameter_estimation": true "do_parameter_estimation": true
} }
``` ```
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment