[Dart-dev] DART/branches Revision: 12226
dart at ucar.edu
dart at ucar.edu
Tue Dec 19 16:11:58 MST 2017
nancy at ucar.edu
2017-12-19 16:11:58 -0700 (Tue, 19 Dec 2017)
177
fix the comments in the ANDERSON 2007 section, and use
the accessor function instead of magic numbers in the test
to see if we're using the enhanced inflation or the original.
Modified: DART/branches/enhanced_inflation/assimilation_code/modules/assimilation/adaptive_inflate_mod.f90
===================================================================
--- DART/branches/enhanced_inflation/assimilation_code/modules/assimilation/adaptive_inflate_mod.f90 2017-12-19 22:41:52 UTC (rev 12225)
+++ DART/branches/enhanced_inflation/assimilation_code/modules/assimilation/adaptive_inflate_mod.f90 2017-12-19 23:11:58 UTC (rev 12226)
@@ -538,11 +538,13 @@
! on the inflation itself are provided in the inflate_handle.
! select which method to update with
-if (inflate_handle%inflation_flavor == 5) then
+if (do_enhanced_ss_inflate(inflate_handle)) then
inf_type = GHA2017
else
+ ! the only other option is the original AND2007 and no one should
+ ! be trying to use that without talking to the dart team first.
+ ! (talk to jeff anderson about the details if you're interested.)
inf_type = AND2009
- ! change here if you want to compare to anderson 2007
endif
! Use bayes theorem to update
@@ -598,13 +600,9 @@
! Use ONLY the linear approximation, cubic solution below can be numerically
! unstable for extreme cases. Should look at this later.
+ !if(gamma_corr > 1.01_r8) then ! use this test to disable the entire first section
if(gamma_corr > 0.99_r8) then
- !if(gamma_corr > 1.01_r8) then
- ! FIXME: rectify the comment in the following line (about gamma
- ! and 1.0) and the test above, which clearly used to be true for
- ! 1.0 but is not anymore.
-
! The solution of the cubic below only works if gamma is 1.0
! Can analytically find the maximum of the product: d/dlambda is a
! cubic polynomial in lambda**2; solve using cubic formula for real root
More information about the Dart-dev
mailing list