#pragma hint clr color
#pragma label clr "color"
#pragma label multdan "multiply"
#pragma label diff "Diffuse"
#pragma label spec "Specular"
#pragma label amb "Ambient"
#pragma label rough "Roughness"
#pragma choice lmodel "0" "Lambert"
#pragma choice lmodel "1" "Phong"
#pragma choice lmodel "2" "Blinn"


float pulse(float a,b,fuzz,x)
{
    return(smooth(a-fuzz,a,x)-smooth(b-fuzz,b,x));
}

surface
danLmodel  (
        vector clr={0,1,0};
        float multdan = 1;
        float diff = 1;
        float spec = 1;
        float amb = 0.1;
    float rough = 0.1;
    string lmodel = "0";
        )
{
    vector dan_noise;
    vector Nn = normalize(N);
    vector In = normalize(I);
    vector Nf = frontface(Nn, I);    
    vector tpos = wo_space(P);

        dan_noise = snoise(tpos)*(clr);
        vector opac = dan_noise;
        vector newcolr = multdan*dan_noise;

    Cf =  newcolr;
        Of =  opac;
    Cf *= Of;

    if (lmodel == "0") Cf += amb * ambient() + diff * diffuse(Nf);
    if (lmodel == "1") Cf += amb * ambient() + diff * diffuse(Nf) + spec * phong(Nf, -In, 1.0/rough);
    if (lmodel == "2") Cf += amb * ambient() + diff * diffuse(Nf) + spec * blinn(Nf, -In, rough);

       

}