anchor.diff:                      floor((model_time - md%last_write_time)/approx_writeout_freq)*approx_writeout_freq
               end if
           end if
  +
  +        if(present(energy_is_finite)) energy_is_finite = energy_is_finite_local
   
       end subroutine
   
  
positive.code:         write(log_output_unit, "(A)") '-----------'
          call md%report_mass_conservation_statistics()
  
  
          TIMER_STOP('printing_stats')
  
  
negative.code:     ! @param print_less_often optional integer (default 1) -- only print every 'nth' time we would otherwise writeout
      ! @param timing_tol real (default 0) if the time since the last write out is > "approximate_writeout_frequency - timing_tol" 
      ! , then do the write. This is an attempt to avoid round-off causing a shift in the write-out times
      ! @param energy_is_finite optional logical variable -- if present it will be passed to the statistics printing routine. In that
      ! case it will be .TRUE. if the global_energy_on_rho is finite, and .FALSE. otherwise. If the statistics are not printed (e.g.
      ! print_less_often > 1) then we do not compute energy, and it will be set to .TRUE.
      subroutine write_outputs_and_print_statistics(md, &
          approximate_writeout_frequency, &
          write_grids_less_often, &
          write_gauges_less_often, &
          print_less_often,&
          timing_tol, &
          energy_is_finite)
  
          class(multidomain_type), intent(inout) :: md
          real(dp), optional, intent(in) :: approximate_writeout_frequency, timing_tol
          integer(ip), optional, intent(in) :: write_grids_less_often, write_gauges_less_often, print_less_often
          logical, optional, intent(out) :: energy_is_finite
  
          real(dp) :: approx_writeout_freq, model_time, timing_tol_local
          integer(ip) :: write_grids_n, write_gauges_n, print_n, j
          logical :: energy_is_finite_local
  
          if(present(approximate_writeout_frequency)) then
              approx_writeout_freq = approximate_writeout_frequency
  