int wandc(double Hc,double Mr,double d,double delta,
    double SS,double g,double *a1,double *a2,double *ad,int option);