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:     !
      ! Convenience print routine
      !
      subroutine print_multidomain(md, global_stats_only)
          class(multidomain_type), intent(inout) :: md ! inout allows for timer to run
          logical, optional, intent(in) :: global_stats_only
  
          integer(ip) :: i, j, k, ecw
          real(dp) :: minstage, maxstage, minspeed, maxspeed, stg1, speed_sq, depth_C, depth_E, depth_N
  
negative.code:             timing_tol_local = 0.0_dp
          end if
  
          energy_is_finite_local = .TRUE.
  
          ! All the domain times should be the same
          model_time = md%domains(1)%time
  
  