Welcome to Fractal Forums

uniform vec3 eye;   //the camera position
uniform float time; //animation time in seconds
const float maxDepth=500.0;
#define ang 0.6
const mat2 rmx=mat2(cos(ang),sin(ang),-sin(ang),cos(ang));

float wvTrn(vec2 z, float amp, float freq, float phase){//simple sine generated terrain
	float h=0.0;
	vec2 offset=vec2(2.5+phase*0.1,1.0+phase*0.07); //phase animates the waves
	for(int i=0;i<9;i++){
		z=z*rmx*freq+offset; //rotate, scale and translate
		h+=sin(z.x+phase)*sin(z.y+phase*1.33)*sin(length(z)*2.33+phase)*amp; //add modulated waves
		amp=pow(amp,1.25); //a fiesta of magic numbers, enjoy!
	}
	return h;
}

float DE(vec3 z) 
{// returns the distance to water, land or sky
	float dMax=maxDepth*0.25-25.0-length((z-eye)*vec3(0.25,1.0,0.25)); //sky radius squished
	float dCld=dMax-min(wvTrn(z.xz*0.005,0.25,2.0,time*0.0625)*25.0,0.0); //clouds
	float dMnt=wvTrn(z.xz*0.01,0.25,2.0,0.0)*25.0+wvTrn(z.xz,0.37,1.75,0.0)*0.05+z.y+1.0; //land (bad detail)
	float dWtr=wvTrn(z.xz,0.2,1.75,time*3.0)*0.25+z.y+1.0; //water
	float d=min(dMnt,min(dWtr,dCld)); //find the closest

	//this coloring section sets global vars for glow, specularity, specular exponent and diffuse color
	//the diffuse color has an alpha value that i use to determine reflection (-) and refraction (+)
	if(bColoring){ //for fragmentarium this is a seperate function: color()
		cGlow=vec3(0.3);spec=0.25;specExp=32.0; 
		diffuse=vec4(0.5+wvTrn(z.xz*0.25,0.25,2.0,0.0),0.7,0.4,0.0)*(0.2+dMnt*4.0);
		if(d==dWtr){ //if we hit water mix with land color if shallow and then reflect
			spec=1.0;cGlow=mix(cGlow,vec3(0.8),min((dMnt-dWtr),1.0));
			diffuse=mix(diffuse,vec4(0.007,0.01,0.015,-0.5),min((dMnt-dWtr)*2.0,1.0));
		}
		if(d==dCld){ //we hit a cloud (not very good coloring yet)
			spec=0.5;cGlow=vec3(0.6,0.7,8.0);specExp=8.0;
			diffuse=mix(vec4(0.5,0.6,1.0,0.0),vec4(1.0,1.0,1.0,0.0),(dCld-dMax)*0.5-0.5);
		}else{ //we didn't hit a cloud so add fog
			float f=distance(z,eye);
			diffuse=mix(diffuse,vec4(0.45,0.55,0.95,0.0),min(f/maxDepth,1.0));
			cGlow=mix(cGlow,vec3(0.45,0.55,0.95),min(f/maxDepth,1.0));
		}
	}
	return d;
}

float rand( vec3 co ){//2d version from lumina extended
	return fract(sin(dot(vec3(co)*0.123,vec3(12.9898,78.233,112.166))) * 43758.5453);
}
float noyz(vec3 co){
	vec3 d=smoothstep(0.0,1.0,fract(co));
	co=floor(co);const vec2 v=vec2(1.0,0.0);
	return mix(mix(mix(rand(co),rand(co+v.xyy),d.x),
		mix(rand(co+v.yxy),rand(co+v.xxy),d.x),
		d.y),mix(mix(rand(co+v.yyx),rand(co+v.xyx),d.x),
		mix(rand(co+v.yxx),rand(co+v.xxx),d.x),d.y),d.z);
}
float HeightMap(vec3 z){
	float w=0.5,f=0.0;
	for(int i=0;i<12;i++){
		f+=w*noyz(z);
		w*=0.5;
		z*=2.0;
	}
	return f;
}

vec4 dnoise3d(vec3 z){
	vec3 iv=floor(z),uv=fract(z);
	vec2 iV=vec2(0.0,1.0);
	vec3 du=30.0*uv*uv*(uv*(uv-2.0)+1.0);
	uv=uv*uv*uv*(uv*(uv*6.0-15.0)+10.0);
	float a = rand( iv );
	vec3 b=vec3(rand( iv+iV.yxx ),rand( iv+iV.xyx ),rand( iv+iV.xxy ));
	vec3 d=vec3(rand( iv+iV.yyx ),rand( iv+iV.xyy ),rand( iv+iV.yxy ));
	vec3 k1=b-vec3(a),k2=vec3(a)-b-b.yzx+d;
	float k3 = - a + b.x + b.y +b.z - d.x - d.y - d.z + rand( iv+iV.yyy );
	return vec4(du*(k1+k2*uv.yzx+k2.zxy*uv.zxy+k3*uv.yzx*uv.zxy),
		a+dot(k1,uv)+dot(k2,uv*uv.yzx)+k3*uv.x*uv.y*uv.z);
}

float iqFBM(vec3 z){
	float w=0.5,f=0.0;
	vec3 dv=vec3(0.0);
	for(int i=0;i<12;i++){
		vec4 n=dnoise3d(z);
		dv+=n.xyz;
		f+=w*n.w/(1.0+dot(dv,dv));
		w*=0.5;
		z*=2.0;
	}
	return f;
}

amp=pow(amp,1.25);

amp=amp*sqrt(sqrt(amp));