anchor.diff:              form='unformatted')
           write(10) eta(1:nx,1:ny)
           close(10)
  -      endif
  +      end if
   
         if (DumpWind .eqv. .true.) then
           open(unit = 10, status='replace',file='output/wind_x.'//num, &
  
positive.code:   do j = 0,ny+1
      do i = 0,nx
        temp(i,j) = wetmask(i,j)- wetmask(i+1,j)
      enddo
    enddo
  
    do j = 0,ny+1
      do i = 0,nx
        if (temp(i,j) .ne. 0.0) then
          hfacE(i,j) = 0d0
        endif
      enddo
    enddo
  
    ! and now for all  eastern cells
    hfacE(nx+1,:) = hfacE(1,:)
  
negative.code:     array(:,ny+1,:) = array(:,1,:)
    else
      array = default
    end if
  
    return
  end subroutine read_input_fileU
  