Contact: | Jeff Anderson |
Reviewers: | |
Revision: | $Revision: 1.5 $ |
Release Name: | $Name: iceland $ |
Change Date: | $Date: 2005/06/11 14:29:32 $ |
Change history: | see CVS log |
Computes a weighting factor to reduce the impact of observations on state variables using information from groups of ensembles. Can be run using groups or using archived summary information available from previous group filter experiments.
types_mod utilities_mod time_manager_mod
use reg_factor_mod, only : | comp_reg_factor |
Optional namelist interface ®_factor_nml may be read from file input.nml.
real(r8) :: comp_reg_factor integer, intent(in) :: num_groups real(r8), intent(in), dimension(num_groups) :: regress integer, intent(in) :: time_index integer, intent(in) :: obs_index integer, intent(in) :: state_index integer, intent(in), optional :: obs_state_ind integer, intent(in), optional :: obs_state_max
Returns a weighting factor given regression factors from each group of a group filter or retrieves a factor generated by previous group filter runs from a file.
num_groups | Number of groups. Set to 1 when using information from previously run group filter from file. |
regress | Regression factor from each group for a given state variable and observation variable pair. |
time_index | Integer index of which time in observation file is being processed by filter. Not used in currently released version. |
obs_index | Integer index of the observation being processed. Not used in current implementation . |
state_index | Integer index of state variable being processed. Not used in current implementation. |
obs_state_ind | Index into file generated for Bgrid model which could be duplicated in other large models. |
obs_state_max | Maximum number of observation state variable pairs with non-zero impacts for a given model and observation sequence. Used for generating Bgrid statistic files. |
We adhere to the F90 standard of starting a namelist with an ampersand '&' and terminating with a slash '/'.
namelist / reg_factor_nml / select_regression, input_reg_file
Contents | Type | Description |
---|---|---|
select_regression | integer | Selects the method for computing regression factor. 1 = compute using sampling theory for any ensemble size. Default. 2 = low order model format. Works from archived time mean or time median regression files generated by low-order models like Lorenz-96. 3 = selects bgrid archived file. This is not currently supported in released versions. |
input_reg_file | character(len=129) | File name from which statistics are to be read for select_regression = 3. Default is "time_mean_reg". |
Routine | Message | Comment |
---|---|---|
comp_reg_factor | Illegal value for namelist parameter select_regression | Value must be 1(sampling theory), 2 (table lookup) or 3 (read from obs_state_reg_file). |
Need to implement large model archive capability in efficient fashion.