displace
d_disp  (
        float height=1; 
        )
{

        P = wo_space(P);

        // Code goes here...

        P.y += 1-pow(abs(s-.5)*2,2);
        P.y *= height;

        P=ow_space(P);
        N=computenormal(P);


}