magma_conduits/volcanoes.lua

local depth_root = -3000 -- TODO: should add some kind of optional magma chamber down there
local depth_base = -50 -- point where the mountain root starts expanding
local depth_maxwidth = -30 -- point of maximum width

local radius_vent = 3 -- approximate minimum radius of vent - noise adds a lot to this
local radius_lining = 5 -- the difference between this and the vent radius is about how thick the layer of lining nodes is, though noise will affect it
local caldera_min = 5 -- minimum radius of caldera
local caldera_max = 20 -- maximum radius of caldera

local snow_line = 120 -- above this elevation snow is added to the dirt type
local snow_border = 15 -- transitional zone where there's dirt with snow on it

local depth_maxpeak = magma_conduits.config.volcano_max_height
local depth_minpeak = magma_conduits.config.volcano_min_height
local slope_min = magma_conduits.config.volcano_min_slope
local slope_max = magma_conduits.config.volcano_max_slope

local region_mapblocks = magma_conduits.config.volcano_region_mapblocks
local mapgen_chunksize = tonumber(minetest.get_mapgen_setting("chunksize"))
local volcano_region_size = region_mapblocks * mapgen_chunksize * 16


local p_active = magma_conduits.config.volcano_probability_active
local p_dormant = magma_conduits.config.volcano_probability_dormant
local p_extinct = magma_conduits.config.volcano_probability_extinct

if p_active + p_dormant + p_extinct > 1.0 then
	minetest.log("error", "[magma_conduits] probabilities of various volcano types adds up to more than 1")
end

local state_dormant = 1 - p_active
local state_extinct = 1 - p_active - p_dormant
local state_none = 1 - p_active - p_dormant - p_extinct

local c_air = minetest.get_content_id("air")
local c_lava = minetest.get_content_id("default:lava_source")
local c_water = minetest.get_content_id("default:water_source")

local c_lining = minetest.get_content_id("default:obsidian")
local c_hot_lining = minetest.get_content_id("default:obsidian")
local c_cone = minetest.get_content_id("default:stone")

local c_ash = minetest.get_content_id("default:gravel")
local c_soil = minetest.get_content_id("default:dirt")
local c_soil_grass = minetest.get_content_id("default:dirt_with_grass")
local c_soil_snow = minetest.get_content_id("default:dirt_with_snow")
local c_snow = minetest.get_content_id("default:snow")
local c_snow_block = minetest.get_content_id("default:snowblock")

local c_underwater_soil = minetest.get_content_id("default:sand")
local c_plug = minetest.get_content_id("default:obsidian")

if magma_conduits.config.glowing_rock then
	c_hot_lining = minetest.get_content_id("magma_conduits:glow_obsidian")
end

local water_level = tonumber(minetest.get_mapgen_setting("water_level"))
local mapgen_seed = tonumber(minetest.get_mapgen_setting("seed"))

-- Mapgen v7 has a bizzare glitch where it will sometimes cut slices out of the generated volcano
-- cone *after* mapgen is finished. The slices are taken at maxp.y or minp.y and resemble
-- the curvy form of a vein or a cave that's one node tall.
-- Some annoying hackery is needed to patch those slices back up
-- again, and I only want to do that hackery if we're actually in mapgen v7.
local mg_name = minetest.get_mapgen_setting("mg_name")

-- derived values

local radius_cone_max = (depth_maxpeak-depth_maxwidth)/(2*slope_min) + radius_lining + 20
local depth_maxwidth_dist = depth_maxwidth-depth_base
local depth_maxpeak_dist = depth_maxpeak-depth_maxwidth

local scatter_2d = function(min_xz, gridscale, border_width)
	local bordered_scale = gridscale - 2 * border_width
	local point = {}
	point.x = math.random() * bordered_scale + min_xz.x + border_width
	point.y = 0
	point.z = math.random() * bordered_scale + min_xz.z + border_width
	return point
end

local get_volcano = function(pos)
	local corner_xz = {x = math.floor(pos.x / volcano_region_size) * volcano_region_size, z = math.floor(pos.z / volcano_region_size) * volcano_region_size}

	local next_seed = math.random(1, 1000000000)
	math.randomseed(corner_xz.x + corner_xz.z * 2 ^ 8 + mapgen_seed)

	local state = math.random()
	if state < state_none then
		math.randomseed(next_seed)
		return nil
	end
	
	local location = scatter_2d(corner_xz, volcano_region_size, radius_cone_max)
	--local location = {x=corner_xz.x+volcano_region_size/2, z = corner_xz.z+volcano_region_size/2} -- For testing, puts volcanoes in a consistent grid
	local depth_peak = math.random(depth_minpeak, depth_maxpeak)
	local depth_lava = math.random(depth_peak - 25, depth_peak)
	local slope = math.random() * (slope_max - slope_min) + slope_min
	local caldera = math.random() * (caldera_max - caldera_min) + caldera_min
		
	math.randomseed(next_seed)
	return {location = location, depth_peak = depth_peak, depth_lava = depth_lava, slope = slope, state = state, caldera = caldera}
end

local perlin_params = {
	offset = 0,
	scale = 1,
	spread = {x=30, y=30, z=30},
	seed = -40901,
	octaves = 3,
	persist = 0.67
}
local nvals_perlin_buffer = {}
local nobj_perlin = nil
local data = {}

-- Used as part of the post-mapgen hackery to fix the weird slices mapgen v7 sometimes takes out of volcano cones
local patch_func = function(patch_area, patch_content)
	local minp = patch_area.MinEdge
	local maxp = patch_area.MaxEdge
	
	local map_vm = minetest.get_voxel_manip(minp, maxp)
	local emin, emax = map_vm:get_emerged_area()
	local patch_data = {}
	map_vm:get_data(patch_data)
	
	local map_area = VoxelArea:new{MinEdge=emin, MaxEdge=emax}
	local map_iterator = map_area:iterp(minp, maxp)
	
	changes_made = false
	
	for vi in patch_area:iterp(minp, maxp) do
		local mi = map_iterator()
		if patch_content[vi] ~= nil and patch_data[mi] ~= patch_content[vi] then
			patch_data[mi] = patch_content[vi]
			changes_made = true
		end	
	end

	if changes_made then	
		--send data back to voxelmanip
		map_vm:set_data(patch_data)
		--calc lighting
		map_vm:set_lighting({day = 0, night = 0})
		map_vm:calc_lighting()
		--write it to world
		map_vm:write_to_map()
	end
end

minetest.register_on_generated(function(minp, maxp, seed)
	if minp.y > depth_maxpeak or maxp.y < depth_root then
		return
	end

	local sidelen = maxp.x - minp.x + 1 --length of a mapblock

	local volcano = get_volcano(minp)
	
	if volcano == nil then
		return -- no volcano in this map region
	end
	
	local depth_peak = volcano.depth_peak
	local base_radius = (depth_peak - depth_maxwidth) * volcano.slope + radius_lining

	-- early out if the volcano is too far away to matter
	if	vector.distance(volcano.location, {x=minp.x, y=0, z=minp.z}) > base_radius + 20 and 
		vector.distance(volcano.location, {x=maxp.x, y=0, z=minp.z}) > base_radius + 20 and 
		vector.distance(volcano.location, {x=maxp.x, y=0, z=maxp.z}) > base_radius + 20 and 
		vector.distance(volcano.location, {x=minp.x, y=0, z=maxp.z}) > base_radius + 20
	then
		return
	end
	
	local patch_area_max
	local patch_content_max	
	local patch_area_min
	local patch_content_min
	if mg_name == "v7" then
		patch_area_max = VoxelArea:new{MinEdge={x=minp.x, y=maxp.y, z=minp.z}, MaxEdge=maxp}
		patch_content_max = {}
		patch_area_min = VoxelArea:new{MinEdge=minp, MaxEdge={x=maxp.x, y=minp.y, z=maxp.z}}
		patch_content_min = {}
	end
	
	local vm, emin, emax = minetest.get_mapgen_object("voxelmanip")
	local area = VoxelArea:new{MinEdge=emin, MaxEdge=emax}
	vm:get_data(data)
	
	local chunk_lengths = {x = sidelen, y = sidelen, z = sidelen} --table of chunk edges

	nobj_perlin = nobj_perlin or minetest.get_perlin_map(perlin_params, chunk_lengths)
	local nvals_perlin = nobj_perlin:get3dMap_flat(minp, nvals_perlin_buffer) -- switch to get_3d_map_flat for minetest v0.5
	local noise_area = VoxelArea:new{MinEdge=minp, MaxEdge=maxp}
	local noise_iterator = noise_area:iterp(minp, maxp)
	
	local x_coord = volcano.location.x
	local z_coord = volcano.location.z
	local depth_lava = volcano.depth_lava
	local caldera = volcano.caldera
	
	local state = volcano.state
	
	for vi, x, y, z in area:iterp_xyz(minp, maxp) do

		local vi3d = noise_iterator()

		local distance_perturbation = (nvals_perlin[vi3d]+1)*10
		local distance = vector.distance({x=x, y=y, z=z}, {x=x_coord, y=y, z=z_coord}) - distance_perturbation

		local dirtstuff
		local replace_soil = false -- determines if the soil type should be replaced with c_soil if there's layers on top of it
		if state < state_dormant then
			if y < water_level then
				dirtstuff = c_underwater_soil
			elseif y < snow_line then
				dirtstuff = c_soil_grass
				replace_soil = true
			elseif y < snow_line + snow_border then
				dirtstuff = c_soil_snow
				replace_soil = true
			else
				dirtstuff = c_snow_block
			end
		else
			dirtstuff = c_ash
		end
		
		local pipestuff
		local liningstuff
		if y < depth_lava + math.random() * 1.1 then
			if state < state_extinct then
				pipestuff = c_plug -- extinct volcano
				liningstuff = c_lining
			else
				pipestuff = c_lava
				liningstuff = c_hot_lining
			end
		else
			if state < state_dormant then
				pipestuff = c_plug -- dormant volcano
				liningstuff = c_lining
			else
				pipestuff = c_air -- active volcano
				liningstuff = c_lining
			end
		end
		
		if y < depth_base then -- pipe
			if distance < radius_vent then
				data[vi] = pipestuff
			elseif distance < radius_lining then
				if data[vi] ~= c_air and data[vi] ~= c_lava then -- leave holes into caves and into existing lava
					data[vi] = liningstuff
				end
			end
		elseif y < depth_maxwidth then -- root
			if distance < radius_vent then
				data[vi] = pipestuff
			elseif distance < radius_lining then
				data[vi] = liningstuff
			elseif distance < radius_lining + ((y - depth_base)/depth_maxwidth_dist) * base_radius then
				data[vi] = c_cone
			end
		elseif y < depth_peak + 5 then -- cone
			local current_elevation = y - depth_maxwidth
			local peak_elevation = depth_peak - depth_maxwidth
			if current_elevation > peak_elevation - caldera and distance < current_elevation - peak_elevation + caldera then
				data[vi] = c_air -- caldera
			elseif distance < radius_vent then
				data[vi] = pipestuff
			elseif distance < radius_lining then
				data[vi] = liningstuff
			elseif distance <  current_elevation * -volcano.slope + base_radius then
				data[vi] = c_cone
			elseif distance < current_elevation * -volcano.slope + base_radius + nvals_perlin[vi3d]*-4 then
				data[vi] = dirtstuff
				if replace_soil and data[vi - area.ystride] == dirtstuff then
					data[vi - area.ystride] = c_soil -- soil underneath a layer of other soil shouldn't have grass on top
				end
				if y >= snow_line then
					if data[vi + area.ystride] == c_air then
						data[vi + area.ystride] = c_snow -- generation advances in a positive y direction so this will be overwritten if more solid stuff is placed above
					end
				end
			end
			
			-- Used as part of the post-mapgen hackery to fix the weird slices mapgen v7 sometimes takes out of volcano cones
			if mg_name == "v7" then
				if y == maxp.y then
					patch_content_max[patch_area_max:index(x,y,z)] = data[vi]
				elseif y == minp.y then
					patch_content_min[patch_area_min:index(x,y,z)] = data[vi]
				end
			end
		end
	end

	-- Used as part of the post-mapgen hackery to fix the weird slices mapgen v7 sometimes takes out of volcano cones
	if mg_name == "v7" then
		minetest.after(2, patch_func, patch_area_max, patch_content_max)
		minetest.after(2, patch_func, patch_area_min, patch_content_min)
	end
	
	--send data back to voxelmanip
	vm:set_data(data)
	--calc lighting
	vm:set_lighting({day = 0, night = 0})
	vm:calc_lighting()
	vm:update_liquids()
	--write it to world
	vm:write_to_map()
end)

----------------------------------------------------------------------------------------------
-- Debugging and sightseeing commands

minetest.register_privilege("findvolcano", { description = "Allows players to use a console command to find volcanoes", give_to_singleplayer = false})

local send_volcano_state = function(pos, name)

	local xmin = math.floor(pos.x / volcano_region_size) * volcano_region_size
	local zmin = math.floor(pos.z / volcano_region_size) * volcano_region_size
	
	local text = "In region (" .. tostring(xmin) .. ", 0, " .. tostring(zmin) ..") to ("
		.. tostring(xmin+volcano_region_size) .. ", 0, " .. tostring(zmin+volcano_region_size) ..")\n"

	volcano = get_volcano(pos)
	if volcano == nil then
		minetest.chat_send_player(name, text.."No volcano present")
		return
	end
	text = text .. "Nearest volcano is at " .. minetest.pos_to_string(volcano.location, 0)
		.. "\nHeight: " .. tostring(volcano.depth_peak) .. " Slope: " .. tostring(volcano.slope)
		.. "\nState: "
	if volcano.state < state_extinct then
		text = text .. "Extinct"
	elseif volcano.state < state_dormant then
		text = text .. "Dormant"
	else
		text = text .. "Active"
	end
	minetest.chat_send_player(name, text)
end

minetest.register_chatcommand("findvolcano", {
    params = "pos", -- Short parameter description
    description = "find the volcano in the player's map region, or in the map region containing pos if provided",
    func = function(name, param)
		if minetest.check_player_privs(name, {findvolcano = true}) then
			local pos = {}
			pos.x, pos.y, pos.z = string.match(param, "^([%d.-]+)[, ] *([%d.-]+)[, ] *([%d.-]+)$")
			pos.x = tonumber(pos.x)
			pos.y = tonumber(pos.y)
			pos.z = tonumber(pos.z)
			if pos.x and pos.y and pos.z then
				send_volcano_state(pos, name)
				return true
			else
				playerobj = minetest.get_player_by_name(name)
				send_volcano_state(playerobj:get_pos(), name)
				return true
			end
		else
			return false, "You need the findvolcano privilege to use this command."
		end
	end,
})