float mult = chramp('size', fit01(@N.y,0,1));

float density = chf('density');

float min = chf('min');

 

mult = fit(mult, 0, 1, min, 1);

@Cd = mult;

 

@pscale = fit01(rand(@ptnum), 0.01, 0.035);

@pscale *= mult;

 

if(rand(@ptnum)>mult*density){ 

    removepoint(0,@ptnum);

}

@orient = rand(@ptnum); 

vector direction  = chv('direction');