Changes between Version 4 and Version 5 of DataAssimilationDiagnostics


Ignore:
Timestamp:
Feb 22, 2023, 2:05:27 PM (14 months ago)
Author:
lnerger
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • DataAssimilationDiagnostics

    v4 v5  
    1919Rank histograms are frequently used to assess the distribution of an ensemble around an observation or, in twin experiments, the true state. The histograms use bins computed form the ensemble distribution and count how frequent e.g. the observation falls into which bin. A flat histogram typically indicates a good ensemble. A concave (U-shaped) histogram indicates too little ensemble spread, while a convex histogram is obtained when the ensemble spread is too large. Further, a sloped histogram indicates bias. (A discussion on the interpretation of rank histograms can be found in Hamill, Monthly Weather Review, 129 (2001) 550-560)
    2020
    21 For a documention on `PDAF_diag_histogram` see the [wiki:PDAF_diag_histogram detail page on PDAF_diag_histogram]. The routine is used in the Lorenz-96 model example in `testsuite/src/lorenz96/compute_truerms.F90`.
     21For a documention on `PDAF_diag_histogram` see the [wiki:PDAF_diag_histogram detail page on PDAF_diag_histogram]. The routine is used in the Lorenz-96 model example in `models/lorenz96/compute_truerms.F90`.
    2222
    2323
     
    2626Ensemble Kalman filters assume that the ensemble is Gaussian distributed. In this case the distribution is symmetric and only the first and second moments of the distribution (the mean and standard deviation) are non-zero. The routine `PDAF_diag_ensstats` allows a data assimilation program to check the values of the third (skewness) and fourth (kurtosis) moment of the distribution. As there are different definition of the kurtosis, please note that PDAF uses the definition used by Lawson and Hansen, Mon. Wea. Rev. 132 (2004) 1966.
    2727
    28 For a documention on `PDAF_diag_ensstats` see the [wiki:PDAF_diag_ensstats detail page on PDAF_diag_ensstats]. The routine is used in the Lorenz-96 model example in `testsuite/src/lorenz96/compute_truerms.F90`.
     28For a documention on `PDAF_diag_ensstats` see the [wiki:PDAF_diag_ensstats detail page on PDAF_diag_ensstats]. The routine is used in the Lorenz-96 model example in `models/lorenz96/compute_truerms.F90`.
    2929
    3030
     
    3434This routine was added in PDAF V2.0. It computes the continuous ranked probability score (CRPS) and its decomposition in reliability and resoliution (Hersbach, Weather and Forecasting, 2000).
    3535
    36 For a documention on `PDAF_diag_CRPS` see the [wiki:PDAF_diag_crps detail page on PDAF_diag_CRPS]. The routine is used in the Lorenz-63 model example in `testsuite/src/lorenz63/compute_truerms.F90`.
     36For a documention on `PDAF_diag_CRPS` see the [wiki:PDAF_diag_crps detail page on PDAF_diag_CRPS]. The routine is used in the Lorenz-63 model example in `models/lorenz63/compute_truerms.F90`.