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:     !
      ! Convenience print routine
      !
      subroutine print_multidomain(md, global_stats_only, energy_is_finite)
          class(multidomain_type), intent(inout) :: md ! inout allows for timer to run
          logical, optional, intent(in) :: global_stats_only
          logical, optional, intent(out) :: energy_is_finite
  
          integer(ip) :: i, j, k, ecw
          real(dp) :: minstage, maxstage, minspeed, maxspeed, stg1, speed_sq, depth_C, depth_E, depth_N
  