-- this function needs to be fed house x, z dimensions and rotation
-- it will then calculate the minimum point (xbmin, avsurfy, zbmin) where the house should be spawned
-- and mark a mapchunk-sized 'house area' for terrain blending

mg_villages.village_area_mark_single_house_area = function(village_area, minp, maxp, pos, pr, village_nr, village)

	local YFLATMIN = 2 -- Lowest flat area height
	local FFAPROP = 0.5 -- front flat area proportion of dimension
	local np_blend = {
		offset = 0,
		scale = 1,
		spread = {x=12, y=12, z=12},
		seed = 38393,
		octaves = 3,
		persist = 0.67
	}

	local sidelen = maxp.x - minp.x + 1
	
	local xdim, zdim -- dimensions of house plus front flat area
	if pos.brotate == 0 or pos.brotate == 2 then
		xdim = pos.bsizex
		zdim = pos.bsizez + math.floor(FFAPROP * pos.bsizez)
	else
		xdim = pos.bsizex + math.floor(FFAPROP * pos.bsizex)
		zdim = pos.bsizez
	end

	-- 2D noise perlinmap
	local chulens = {x=sidelen, y=sidelen, z=sidelen}
	local minpos = {x=minp.x, y=minp.z}
	local nvals_blend = minetest.get_perlin_map(np_blend, chulens):get2dMap_flat(minpos)

	-- mark mapchunk-sized house area
	local ni = 1
--	for z = minp.z, maxp.z do
--	for x = minp.x, maxp.x do -- for each column do
	for z = math.max( village.vz - village.vs, minp.z), math.min(village.vz + village.vs, maxp.z), 1 do
	for x = math.max( village.vx - village.vs, minp.x), math.min(village.vx + village.vs, maxp.x), 1 do -- for each column do

		local xrm = x - minp.x -- relative to mapchunk minp
		local zrm = z - minp.z
		local xr = x - village.vx -- relative to blend centre
		local zr = z - village.vz
		local xre1 = (zdim / 2) * (xr / zr)
		local zre1 = zdim / 2
		local xre2 = xdim / 2
		local zre2 = (xdim / 2) * (zr / xr)
		local rade1 = math.sqrt(xre1 ^ 2 + zre1 ^ 2)
		local rade2 = math.sqrt(xre2 ^ 2 + zre2 ^ 2)
		local flatrad = math.min(rade1, rade2) -- radius at edge of rectangular house flat area
		local n_absblend = math.abs(nvals_blend[ni])
		local blenradn = village.vs - n_absblend * 2 -- vary blend radius
		local flatradn = flatrad + n_absblend * 2 -- vary shape of house flat area
		local nodrad = math.sqrt(xr ^ 2 + zr ^ 2) -- node radius

		-- only blend the terrain if it does not already belong to another village
		if( village_area[ x ][ z ][ 2 ] == 0 ) then
			if    x >= (pos.x-1) and x <= (pos.x + pos.bsizex + 1) -- area reserved for house
			  and z >= (pos.z-1) and z <= (pos.z + pos.bsizez + 1) then
				village_area[ x ][ z ] = {village_nr, 4}
			elseif nodrad <= flatradn or (xr == 0 and zr == 0) then -- irregular flat area around house
				village_area[ x ][ z ] = {village_nr, 1}
			elseif nodrad <= blenradn then -- terrain blend area
				local blenprop = ((nodrad - flatradn) / (blenradn - flatradn))
				village_area[ x ][ z ] = {village_nr, -1 * blenprop} -- terrain blending
			else -- no change to terrain
				--village_area[xrm][zrm] = {village_nr, 0}
			end
		end
		ni = ni + 1
	end
	end
end