diff --git a/tunnel.glsl b/tunnel.glsl
new file mode 100644
index 0000000..e1e6774
--- /dev/null
+++ b/tunnel.glsl
@@ -0,0 +1,65 @@
+#ifdef GL_ES
+precision mediump float;
+#extension GL_OES_standard_derivatives : enable
+#endif
+
+uniform float time;
+uniform vec2 resolution;
+
+const float PI = 3.14159;
+
+vec3 hsv(float h, float s, float v) {
+	float c = s * v;
+	float _ = mod(h * 6.0, 6.0);
+	vec3 C = vec3(c, c*(1.0 - abs(mod(_, 2.0) - 1.0)), 0.0);
+	if (_ < 1.0) {
+		C = vec3(C.x, C.y, C.z);
+	} else if (_ < 2.0) {
+		C = vec3(C.y, C.x, C.z);
+	} else if (_ < 3.0) {
+		C = vec3(C.z, C.x, C.y);
+	} else if (_ < 4.0) {
+		C = vec3(C.z, C.y, C.x);
+	} else if (_ < 5.0) {
+		C = vec3(C.y, C.z, C.x);
+	} else {
+		C = vec3(C.x, C.z, C.y);
+	}
+	return C + (v - c);
+}
+
+float map(vec3 p) {
+	return 2.0 - length(p.xz);
+}
+
+float noise(vec2 co){
+    return fract(sin(dot(co.xy ,vec2(12.9898,78.233))) * 43758.5453);
+}
+
+void main( void ) {
+	vec2 p = (2.0 * gl_FragCoord.xy - resolution) / resolution.y;
+	vec3 cp = vec3(cos(time * 0.2), 0.0, sin(time * 0.45)) * 0.5;
+	vec3 cl = vec3(-sin(time), 10.0, cos(time));
+	vec3 cf = normalize(cl - cp);
+	vec3 cs = normalize(cross(cf, vec3(sin(time * 0.1), 0.0, cos(time * 0.1))));
+	vec3 cu = normalize(cross(cs, cf));
+	float focus = 0.5;
+	vec3 rd = normalize(cs * p.x + cu * p.y + cf * focus);
+	vec3 rp = cp;
+	for (int i = 0; i < 64; ++i) {
+		float d = map(rp);
+		if (d < 0.001)
+			break;
+		rp += rd * d;
+	}
+	float a = (atan(rp.z, rp.x)) * 16.0 / PI;
+	float div = 1./pow(2.,3.+floor(sin(time*3.)*2.));
+	float ai = floor(a*div);
+	float af = fract(a*div);
+	float d = (rp.y + 2.5 * time) * 10.0;
+	float di = floor(d*div);
+	float df = fract(d*div);
+	float v = 32.0 * af * (1.0 - af) * df * (1.0 - df) * exp(-rp.y * 0.8);
+	gl_FragColor = vec4(hsv(noise(vec2(ai, di) * 0.01), 1.0, v), 1.0);
+}
+