UCODE_2005 Version 1.013 08/15/2008
                 - Updated version # to 1.013
                 - Support for nearly all of the UCODE_2005 data-exchange files
                   is now available in GW_Chart, freeware available from
                   http://water.usgs.gov/nrp/gwsoftware/GW_Chart/GW_Chart.html.
                 - Corrected error for the case when the final parameters did
                   not yield the lowest sum-of-squared residuals. In that
                   situation ucode was printing information related to
                   parameters from the final iteration rather than those that
                   yielded the lowest sum-of-squared residuals. This can affect
                   UCODE_2005 output files with filename extensions #uout, _b1,
                   _b2, _dm, _linp, _mc, _mv, _nm, _os, _pc, _r, _s1, _sc, _sd,
                   _so, _sppp, _sppr, _spsp, _spsr, _spu, _su, _w, _ws, _ww, and
                   _wt (if keyword WtOSConstant>0, see page 85 of the 
                   documentation); Residual_Analysis output files with
                   extentions _rb, _rc, _rd, and _rg; and Residual_Analysis_ADV
                   output file with extention _rdadv. This correction affects
                   example ex1a in appendix A of the documentation in that the
                   values in some of the files listed above are slightly
                   different.
                 - Corrected error in which the flag for calculating nonlinear
                   intervals on a parameter in the parameter_data block
                   was ignored for mode NonlinearIntervals=yes in the
                   UCODE_CONTROL_DATA block.
                 - Added a final substitution of the template files so they
                   contain the optimal values. Prior to this, the template
                   files contained values perturbed for sensitivity
                   calculations.
                 - To allow multiple runs to be used to calculate nonlinear
                   intervals (as often is needed because of lengthy runs and
                   convergence problems), new files with date/time stamps are
                   produced for files with extentions _intconf, _intconfpar,
                   _intpred, and _intpredpar to create, for example,
                   ucode_rootname--2008_06_21--22.17._intconf. For each
                   UCODE_2005 run, files with extensions _intconf, _intconfpar,
                   _intpred, and _intpredpar are still produced. As before,
                   they contain the intervals calculated in the most recent run.
                   To obtain a full list of inteval limits, values from the
                   time-stamped files may need to be moved into these files.
                   The nonlinear intervals can now be plotted in graph form
                   using GW_Chart.
                 - Corrected an error in which an additional forward&der run
                   was executed for the trustregion option when derivatives were
                   read from process model output
                 - Corrected value of uncertainty associated with prior
                   information on "parameters for prediction" (See note for
                   UCODE_2005 version 1.009). Previously the variance on these
                   parameters was multiplied by the cev of the regression,
                   but now the variance is used directly without multiplication.
                 - Previously parameters that were omitted from the regression
                   due to parameter constraints were marked with a !. Now
                   parameters that are omitted due to insensitivity are also
                   marked with a !.
                 - Added an underscore file _scgrp that reports composite scaled
                   sensitivities and the two observations with the largest
                   sensitivities for each observation group.
                 - Modified ucode such that regressions with parameters omitted
                   due to insensitivity or parameter value constraints could be
                   used for advanced evaluations (e.g. predictive uncertainty,
                   nonlinear intervals) without rerunning the regression with
                   those parameters as adjustable=no. The uncertainty measures
                   include the uncertainty associated with the "parameters for
                   prediction".
                 - increased number of significant figures printed in _pasub to
                   match those printed in _paopt.  Format was changed from
                   1PE15.7 to 1PE25.16.
                 - increased number of significant figures printed in #uout for
                   SUM OF SQUARED, WEIGHTED RESIDUALS.  Format was changed from
                   G12.5 to G15.7.
                 - unused subroutine UCODE_DEF_WRITEITER was deleted from
                   ucode_mod.f90
                 - removed write statements that are no longer needed from
                   ucode_2005.f90

                 - Updated Corfac_Plus to version 1.005. Corrected reading of
                   prior for "parameters for prediction" (See note for
                   UCODE_2005 version 1.009) in the case where
                   parameters had been omitted from the regression due to
                   insensitivity or constraint. This problem caused corfac_plus
                   to fail in this situation.

                 - Updated Linear_Uncertainty to version 1.006. Increased
                   format for writing the numbers of items at the start of the
                   output to accomodate more observations.

                 - Updated Model_Linearity_Adv to version 1.004 to use results
                   from projects that include parameters for prediction.