[Dart-dev] DART/branches Revision: 11870

[dart at ucar.edu]

nancy at ucar.edu
2017-08-03 10:38:47 -0600 (Thu, 03 Aug 2017)
38
merge in changes from the rma_trunk.




Index: DART/branches/gigg
===================================================================
--- DART/branches/gigg	2017-08-03 15:00:42 UTC (rev 11869)
+++ DART/branches/gigg	2017-08-03 16:38:47 UTC (rev 11870)

Property changes on: DART/branches/gigg
___________________________________________________________________
Modified: svn:mergeinfo
## -33,6 +33,7 ##
 /DART/branches/rma_sprint:8095-8166
 /DART/branches/rma_state_structure_test_dir:9472-9676
 /DART/branches/rma_state_structure_unit_test:9691-10268
+/DART/branches/rma_trunk:11851-11864
 /DART/branches/rma_trunk_clamping:8786-8801
 /DART/branches/rma_trunk_filename:8011-8157
 /DART/branches/rma_trunk_get_state_fun:8568-8596
Modified: DART/branches/gigg/CHANGELOG
===================================================================
--- DART/branches/gigg/CHANGELOG	2017-08-03 15:00:42 UTC (rev 11869)
+++ DART/branches/gigg/CHANGELOG	2017-08-03 16:38:47 UTC (rev 11870)
@@ -491,7 +491,7 @@
   names changed from 'restart' to 'state' for input and output files.
 
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
-+ July 31 2017 :: single filenames, random distributions, bug fixes.  $Revision$
++ Aug 2 2017 :: single filenames, random distributions, bug fixes.  $Revision$
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
 - added code to support listing input and output filenames directly in the
@@ -510,6 +510,11 @@
 - various updates to the CAM scripts to work more smoothly with the most
   recent CIME changes and DART Manhattan updates.
 
+- added QTY_CWP_PATH and QTY_CWP_PATH_ZERO to the default quantities list for
+  the obs_def_cwp_mod.f90 forward operator.
+
+- improved some error messages in the diagnostic matlab scripts
+
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 + MMM DD 2017 :: next revision summary.  Revision: 
 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

Deleted: DART/branches/gigg/assimilation_code/modules/assimilation/README_reg_factor_mod
===================================================================
--- DART/branches/gigg/assimilation_code/modules/assimilation/README_reg_factor_mod	2017-08-03 15:00:42 UTC (rev 11869)
+++ DART/branches/gigg/assimilation_code/modules/assimilation/README_reg_factor_mod	2017-08-03 16:38:47 UTC (rev 11870)
@@ -1,77 +0,0 @@
-# DART software - Copyright UCAR. This open source software is provided
-# by UCAR, "as is", without charge, subject to all terms of use at
-# http://www.image.ucar.edu/DAReS/DART/DART_download
-#
-# DART $Id$
-
-The module reg_factor_mod is an attempt to deal with sampling errors
-in ensemble filters using a hierarchical Monte Carlo approach in which
-a group of ensembles is used and the group's sample of the regression
-factor between an observation and a state variable is used to compute
-an error estimate. The presence of error generally means that the regression
-factor should be reduced by some factor in order to give the optimal 
-result. 
-
-The regression factor used here is computed as a function of the 
-group size and the ratio of the sample standard deviation to the sample
-mean of the regression coefficient. Files are generated offline for each
-desired group size and contain a list of ratios, Q, and corresponding 
-regression confidence factors. At present, files have been generated for
-groups of 2, 4, 8 and 16. 
-
-The first step in generating these files is to generate a file for an
-infinite group size using the program sys_sim401.f90. In this case, the
-only input parameter is the number of samples that are used to generate
-the regression confidence factors for a series of Q's. In the results used 
-to date, 1000000 samples are used to generate this. Results are generated
-for every 0.01 from Q = 0.0 up to 6.00. The output from this procedure can
-be found in the file undamped_base.
-
-Because the regression confidence factor as a function of Q has a very long
-but small valued tail, for computational efficiency (and error tolerance), 
-the program sys_sim401.f90 has been modified to linearly damp the tail of 
-the distribution for values of Q > 3.0 so that they approach zero for 6.0.
-The output of this procedure can be found in the file damped_base.
-
-Because of sampling error in the Monte Carlo algorithm used to compute the
-regression confidence factors as a function of Q, there is some noise in the
-results from the two steps mentioned so far, especially for values of Q 
-larger than 3.0. A simple smoother is used to reduce this noise by doing 
-7 point centered averages of the result for the damped case. These can be 
-found in the files smoothed_damped_base. 
-
-Next, an attempt is made, rather feebly, to account for additional sampling 
-error that comes from using small groups. The program sys_sim402.f90 does
-this operation. It takes as input the group size, and a sample size that is
-used to construct Monte Carlo statistics as above. Again, results used here
-make use of 1000000 samples. 
-
-sys_sim402.f90 attempts to account for the uncertainty in the sample statistics
-of Q that are obtained from a group of ensembles. It generates group size 
-samples of a distribution with mean 1.0 and standard deviation Q. It then 
-computes Q from each of these samples and finally computes the 1000000 sample
-mean value of Q. In general, for all but the smallest input Q's, the mean
-sample Q is larger. The value of the regression confidence factor for the input