== MS_rot3C4

 Script to rotate a set of elastic constants (Cij) by 3 angles:

 alp = clockwise about 1-axis (looking at origin, == yaw)
 bet = clockwise about 2-axis (looking at origin, == -dip)
 gam = clockwise about 3-axis (looking at origin, == azimuth)

 [CR] = CIJ_rot3(C,alp,bet,gam)

  Inputs: the Cij elastic stiffness tensor, 3 angles (degrees)

  Output is the rotated Cij elastic stiffness tensor                  

  NOTE: The rotations are applied in order, ie: alpha, then beta then gamma


