Utility program to generate an ensemble of perturbed ensemble member restart files. This program can be run in parallel and used as a stand alone program.
This namelist is read from the file input.nml. Namelists start with an ampersand '&' and terminate with a slash '/'. Character strings that contain a '/' must be enclosed in quotes to prevent them from prematurely terminating the namelist.
&perturb_single_instance
ens_size = ''
input_files = ''
output_files = ''
output_file_list = ''
perturbation_amplitude = 0.0
single_restart_file_in = .false.
/| Item | Type | Description |
|---|---|---|
| ens_size | integer | Total number of ensemble members. |
| input_files | character(len=256),dimension(num_domains) | The restart file you would like to perturb from. |
| output_file_list | character(len=256) | A file containing a list of the desired output names. |
| output_files | character(len=256) | An array of filenames |
| perturbation_amplitude | real(r8) | The desired perturbation amplitude. If the model provides an interface then it will use that subroutine, otherwise it will simply add gaussian noise to the entire state, and this is the standard deviation. |
| single_restart_file_in | logical | A boolean, specifying if you have a single file restart, such as the case for lower order models. |
Below is an example of a typical namelist for the perturb_single_instance.
&perturb_single_instance_nml
ens_size = 3
input_files = 'caminput.nc'
output_files = 'cam_pert1.nc','cam_pert2.nc','cam_pert3.nc'
/- inputfile.nc (description file that will be perturbed)
- output_file_list.txt (a file containing a list of restart files) and,
- perturb_single_instance.nml
- none