anchor.diff:  call test_value(dcp,0.58784927087317074d0,'mij2dcp: dcp')
   
   ! mij2mag
  -call mij2mag(mrr,mtt,mpp,mrt,mrp,mtp,mag)
  -call test_value(mag,6.4522434199195224d0,'mij2mag: mag')
  +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
  
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: pnt(12) =  5.9647319496912486d+018
  call pnt2mom(pnt,mom)
  call test_value(mom,5.3501397941634130d+018,'pnt2mom: mom')
  call pnt2mom_nondc(pnt,mom)
  call test_value(mom,5.4550129578214554d+018,'pnt2mom_nondc: mom')
  
  ! pnt2sdr
  ! P, N, and T axes to best-fitting double couple strike, dip, and rake
  