[Dart-dev] DART/branches Revision: 12687

[dart at ucar.edu]

jla at ucar.edu
2018-06-27 16:02:00 -0600 (Wed, 27 Jun 2018)
158
Corrected confusion about how scale factors should be related in call from inverse gamma to gamma.
These are confirmed by the test harness test_inverse_gamma



Modified: DART/branches/gigg/assimilation_code/modules/utilities/random_seq_mod.f90
===================================================================
--- DART/branches/gigg/assimilation_code/modules/utilities/random_seq_mod.f90	2018-06-27 20:55:57 UTC (rev 12686)
+++ DART/branches/gigg/assimilation_code/modules/utilities/random_seq_mod.f90	2018-06-27 22:02:00 UTC (rev 12687)
@@ -257,6 +257,7 @@
    call error_handler(E_ERR, 'random_inverse_gamma', errstring, source, revision, revdate)
 endif
 
+! Note that the reciprocal of the scale from the inverse gamma is required as the scale input for the gamma here
 g = ran_gamma(r, shape, 1.0_r8 / scale)
 
 ! ran_gamma can't return 0 so its safe to divide by it