Mesh.py

ChangeState( Kd = ( 0.0, 1.0, 0.0 ),
             Ka = ( 1.0, 1.0, 0.0 ),
             Ks = ( .3, .3, .3 ),
             Tr = 1.,
             Ns = 1.,
             intensity = 0.8 )

Mesh( vertices = self.borehole_vertices.flatten( ),
      indices = self.borehole_indices,
      fmt     = self.vertex_format,
      mode = 'lines' )

self.borehole_vertices = np.array( [ [ 0, 0, 0 ], [ 1, 1, 1 ], [ 2, 2, 2 ], ... ] )
self.borehole_indices = [ 0, 1, 1, 2, 2, 3, ... ]

shaders.glsl

/* simple.glsl

simple diffuse lighting based on laberts cosine law; see e.g.:
    http://en.wikipedia.org/wiki/Lambertian_reflectance
    http://en.wikipedia.org/wiki/Lambert%27s_cosine_law
*/
---VERTEX SHADER-------------------------------------------------------
#ifdef GL_ES
    precision highp float;
#endif

attribute vec3  v_pos;
attribute vec3  v_normal;

uniform mat4 modelview_mat;
uniform mat4 projection_mat;
uniform float Tr;

varying vec4 normal_vec;
varying vec4 vertex_pos;

void main (void) {
     //compute vertex position in eye_sapce and normalize normal vector
    vec4 pos = modelview_mat * vec4(v_pos,1.0);
    vertex_pos = pos;
    normal_vec = vec4(v_normal,0.0);
    gl_Position = projection_mat * pos;
}


---FRAGMENT SHADER-----------------------------------------------------
#ifdef GL_ES
    precision highp float;
#endif

varying vec4 normal_vec;
varying vec4 vertex_pos;

uniform mat4 normal_mat;
uniform vec3 Kd;
uniform vec3 Ka;
uniform vec3 Ks;
uniform float Tr;
uniform float Ns;
uniform float intensity;

void main (void){
    //correct normal, and compute light vector (assume light at the eye)
    vec4 v_normal = normalize( normal_mat * normal_vec );
    vec4 v_light = normalize( vec4(0,0,0,1) - vertex_pos );
    //reflectance based on lamberts law of cosine
    
    vec3 Ia = intensity*Kd;
    vec3 Id = intensity*Ka * max(dot(v_light, v_normal), 0.0);
    vec3 Is = intensity*Ks * pow(max(dot(v_light, v_normal), 0.0), Ns); 
    gl_FragColor = vec4(Ia + Id + Is, Tr);