anchor.diff:  mtp = 1.1714d+18
   call mij2dcp(mrr,mtt,mpp,mrt,mrp,mtp,dcp)
   call test_value(dcp,0.58784927087317074d0,'mij2dcp: dcp')
  +
  +! TODO: mij2mag
   
   ! mij2mom
   ! Moment tensor to scalar moment
  
positive.code: call test_value(mrp,-3.1150000000000003d+021,'pnt2mij: mrp')
  call test_value(mtp, 1.0330000000000004d+021,'pnt2mij: mtp')
  
  ! pnt2mom
  ! P, N, and T axes to scalar moment
  pnt(1) = -0.64401587537772065
  
negative.code: double precision :: fth, fss, fno
  double precision :: dcp
  double precision :: mom, mag
  double precision :: wid, len
  
  ! mij2dcp
  ! Moment tensor to double couple percentage
  