void w3j(l1,s,l2,w)
int l1,s,l2;
float **w;
{
double *wk,*darray1();
int dim,m1,m2,ind;

  dim=2*l1+1;
  wk=darray1(0,dim*(2*l2+1)-1);
  w3j_(&l1,&s,&l2,wk,&dim);
  ind=0;
  for(m2=-l2;m2<=l2;m2++) 
  for(m1=-l1;m1<=l1;m1++) w[m1][m2]=wk[ind++];
  free_darray1(wk,0);
}