void euler(double rmap[3][3], double phi, double rho, double g);
void orient(double rmap[3][3]);
