anchor.diff:          
           if (.not.data_allocated) then
               if (debug) then
  -                write(*,*) "Allocating Linear Wind Data and FFTW plans"
  +                write(*,"(A,A)") char(13),"Allocating Linear Wind Data and FFTW plans"
               endif
               ! using fftw_alloc routines to ensure better allignment for vectorization
               n_elements=nx*ny*nz
  
positive.code:         if (.not.options%ideal) then
              call add_buffer_topo(domain%terrain,complex_terrain_firstpass,5)
              buffer=2
              call add_buffer_topo(real(real(complex_terrain_firstpass)),complex_terrain,0)
              buffer=buffer+original_buffer
          else
              call add_buffer_topo(domain%terrain,complex_terrain,0)
          endif
          
          nx=size(complex_terrain,1)
  
negative.code:                 end do
              end do
          endif
          
          if (present(debug)) then
              if (debug) then
                  allocate(real_terrain(nx,ny))
                  real_terrain=buffer_topo
                  call io_write2d("buffered_terrain.nc","data",real_terrain)
                  deallocate(real_terrain)
              endif
          endif
          
      end subroutine add_buffer_topo
      
  