anchor.diff:  mrt = -3.5502d+18
   mrp = -3.843d+18
   mtp = 1.1714d+18
  -call mij2mom(mrr,mtt,mpp,mrt,mrp,mtp,mom)
  -call test_value(mom,5.3501397941634130d+018,'mij2mom: mom')
  +call mij2mom(mrr,mtt,mpp,mrt,mrp,mtp,mom,'dc')
  +call test_value(mom,5.3501397941634130d+018,'mij2mom: mom (dc)')
  +call mij2mom(mrr,mtt,mpp,mrt,mrp,mtp,mom,'nondc')
  +call test_value(mom,5.4550129578214554d+018,'mij2mom: mom (non-dc)')
   
   ! mij2pnt
   ! Moment tensor to P, N, and T axes
  
positive.code: ! pnt2mag
  call pnt2mag(pnt,mag)
  call test_value(mag,8.4494036748523378d0,'pnt2mag: mag')
  
  ! pnt2mom
  ! P, N, and T axes to scalar moment
  
negative.code: call test_value(dcp,0.58784927087317074d0,'mij2dcp: dcp')
  
  ! mij2mag
  call mij2mag(mrr,mtt,mpp,mrt,mrp,mtp,mag,'dc')
  call test_value(mag,6.4522434199195224d0,'mij2mag: mag (dc)')
  call mij2mag(mrr,mtt,mpp,mrt,mrp,mtp,mag,'nondc')
  call test_value(mag,6.4578638576982641d0,'mij2mag: mag (non-dc)')
  
  ! mij2mom
  ! Moment tensor to scalar moment
  