[Dart-dev] DART/branches Revision: 12691
dart at ucar.edu
dart at ucar.edu
Wed Jun 27 16:29:37 MDT 2018
jla at ucar.edu
2018-06-27 16:29:37 -0600 (Wed, 27 Jun 2018)
22
Cleaning up comments.
Modified: DART/branches/gigg/assimilation_code/programs/perfect_model_obs/perfect_model_obs.f90
===================================================================
--- DART/branches/gigg/assimilation_code/programs/perfect_model_obs/perfect_model_obs.f90 2018-06-27 22:22:30 UTC (rev 12690)
+++ DART/branches/gigg/assimilation_code/programs/perfect_model_obs/perfect_model_obs.f90 2018-06-27 22:29:37 UTC (rev 12691)
@@ -563,7 +563,7 @@
if( qc_ens_handle%vars(i, 1) == 0 ) then
! For GIGG filters, need to define different distributions for observation uncertainty
- ! The obs_err_var in the obs_seq file has different meanings for different optiona
+ ! The obs_err_var in the obs_seq file has different meanings for different option
if(gigg_obs_error_pdf == 1) then
! Standard gaussian
obs_value(1) = random_gaussian(random_seq, true_obs(1), &
@@ -570,7 +570,6 @@
sqrt(get_obs_def_error_variance(obs_def)))
elseif(gigg_obs_error_pdf == 2) then
! Observation pdf given truth is inverse gamma and value in file is type 1 relative error variance
- ! Get alpha and beta for the random number generation
T1Rr = get_obs_def_error_variance(obs_def)
shape = 1.0_r8 / T1Rr + 2.0_r8
scale = true_obs(1) * (shape - 1.0_r8)
@@ -577,7 +576,6 @@
obs_value(1) = random_inverse_gamma(random_seq, shape, scale)
elseif(gigg_obs_error_pdf == 3) then
! Observation pdf given truth is gamma and value in file is type 1 relative error variance
- ! Get alpha and beta for the random number generation
T1Rr = get_obs_def_error_variance(obs_def)
shape = 1.0_r8 / T1Rr
scale = true_obs(1) / shape
More information about the Dart-dev
mailing list