anchor.diff:          end do
       endif
   
  +    ypos = (ubound(domain%accumulated_precipitation,2)-lbound(domain%accumulated_precipitation,2))/2
  +    ypos = ypos + lbound(domain%accumulated_precipitation,2)
       do i=1,200
           ! print *,"Microphysics"
  -        ! note this should be wrapped into the domain object(?)
  +        ! note should this be wrapped into the domain object(?)
           call microphysics(domain, dt = 20.0)
           ! print *,"domain%advect(dt = 4.0)"
           call domain%advect(dt = 1.0)
   
  -
           ! print *,"domain%halo_exchange()"
           call domain%halo_exchange()
   
  -        if (this_image()==1) then
  -            print*, i, domain%accumulated_precipitation(::20,2)
  +        if (this_image()==(num_images()/2)) then
  +            print*, this_image(), i, domain%accumulated_precipitation(::10,ypos)
           endif
           ! call domain%enforce_limits()
       end do
  
positive.code: 
    block
      type(domain_t) :: domain
      integer :: i,nz
      print *,this_image(),"domain%initialize_from_file('input-parameters.txt')"
      call domain%initialize_from_file('input-parameters.txt')
  
  
negative.code:   end block
  
   if (this_image()==1) print *,"Test passed."
  
  end program
  