anchor.diff:          ! ideally this will just become timer_start, save_file()...
           ! the output object needs a pointer to model_time and will have to know how to create output file names
           call output_timer%start()
  -        if (domain%model_time >= next_output) then
  +        if ((domain%model_time + small_time_delta) >= next_output) then
               if (this_image()==1) write(*,*) "Writing output file"
               if (i>24) then
                   write(file_name, '(A,I6.6,"_",A,".nc")')    &
  
positive.code:         !
          ! -----------------------------------------------------
          call input_timer%start()
          if (boundary%current_time <= domain%model_time ) then
              if (this_image()==1) write(*,*) ""
              if (this_image()==1) write(*,*) " ----------------------------------------------------------------------"
              if (this_image()==1) write(*,*) "Updating Boundary conditions"
  
negative.code:     use initialization,     only : init_model
      use timer_interface,    only : timer_t
      use time_object,        only : Time_type
      use time_delta_object,  only : time_delta_t
      use wind,               only : update_winds
      use restart_interface,  only : restart_model
      use icar_constants,     only : kVARS
  