== MS_thom

   CIJ_THOM - generate elastic constants from Thomsen parameters

   [C]=CIJ_thom(vp,vs,rh,eps,gam,del)

   Inputs: 
       rh  : Density (kg/m2 or g/cm3) (automatic conversion < 50) 
       vp  : m/s or km/s (input is automatically to m/s converted if < 50)
       vs  : m/s or km/s (input is automatically to m/s converted if < 50)
       eps, gam, del : Dimensionless
