From fca407bdeff2fba296d8e08a8a098c139d7b342f Mon Sep 17 00:00:00 2001 From: Jason Christopherson Date: Tue, 18 Feb 2025 07:55:58 -0600 Subject: [PATCH] Code clean up --- src/fstats_mcmc.f90 | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/src/fstats_mcmc.f90 b/src/fstats_mcmc.f90 index fcaef9f..1ba0e39 100644 --- a/src/fstats_mcmc.f90 +++ b/src/fstats_mcmc.f90 @@ -202,7 +202,7 @@ function mh_get_chain(this, bin, err) result(rst) !! The resulting chain with each parameter represented by a column. ! Local Variables - integer(int32) :: npts, nvar, flag, nstart + integer(int32) :: npts, nvar, flag, nstart, n class(errors), pointer :: errmgr type(errors), target :: deferr @@ -213,6 +213,7 @@ function mh_get_chain(this, bin, err) result(rst) errmgr => deferr end if npts = this%get_chain_length() + n = npts nvar = this%get_state_variable_count() if (present(bin)) then nstart = floor(bin * npts) @@ -223,7 +224,7 @@ function mh_get_chain(this, bin, err) result(rst) ! Process allocate(rst(npts, nvar), stat = flag, & - source = this%m_buffer(nstart:,1:nvar)) + source = this%m_buffer(nstart:n,1:nvar)) if (flag /= 0) then call report_memory_error(errmgr, "mh_get_chain", flag) return