anchor.diff:  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
  
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: mrt = -3.5502d+18
  mrp = -3.843d+18
  mtp = 1.1714d+18
  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
  