Fixed input_xmls setting of size for -yn moments (size is (N+1)^2, not N)

This commit is contained in:
Adam Nelson 2014-04-04 12:06:01 -04:00
parent 8a752e6f51
commit 26edd77b45

View file

@ -2371,8 +2371,8 @@ contains
end if
t % score_bins(j : j + n_order) = SCORE_FLUX_YN
t % moment_order(j : j + n_order) = n_order
j = j + n_order
t % moment_order(j : j + n_order) = (n_order + 1)**2
j = j + (n_order + 1)**2
case ('total')
t % score_bins(j) = SCORE_TOTAL
@ -2390,8 +2390,8 @@ contains
end if
t % score_bins(j : j + n_order) = SCORE_TOTAL_YN
t % moment_order(j : j + n_order) = n_order
j = j + n_order
t % moment_order(j : j + n_order) = (n_order + 1)**2
j = j + (n_order + 1)**2
case ('scatter')
t % score_bins(j) = SCORE_SCATTER
@ -2439,15 +2439,15 @@ contains
t % estimator = ESTIMATOR_ANALOG
! Setup P0:Pn
t % score_bins(j : j + n_order) = SCORE_SCATTER_YN
t % moment_order(j : j + n_order) = n_order
j = j + n_order
t % moment_order(j : j + n_order) = (n_order + 1)**2
j = j + (n_order + 1)**2
case ('nu-scatter-yn')
t % estimator = ESTIMATOR_ANALOG
! Setup P0:Pn
t % score_bins(j : j + n_order) = SCORE_NU_SCATTER_YN
t % moment_order(j : j + n_order) = n_order
j = j + n_order
t % moment_order(j : j + n_order) = (n_order + 1)**2
j = j + (n_order + 1)**2
case('transport')
t % score_bins(j) = SCORE_TRANSPORT