{
  dipvec[0] = cos(dip*.01745 ) * sin (azi*.01745) ;
  dipvec[1] = cos(dip*.01745) * cos(azi*.01745);
  dipvec[2] = -sin (dip*.01745); 
   
  strikevec[0] =  sin (azi*.01745 + 1.5708);
  strikevec[1] =  cos(azi*.01745 + 1.5708);
  strikevec[2] =  0; 
 
  normalvec[0] = dipvec[1] * strikevec[2] - dipvec[2] * strikevec[1] ;
  normalvec[1] = dipvec[2] * strikevec[0] - dipvec[0] * strikevec[2] ;
  normalvec[2] = dipvec[0] * strikevec[1] - dipvec[1] * strikevec[0] ;
 }