anchor.diff:  real(kind=sgl)                          :: res(4)
   
   real(kind=sgl)                          :: t, omega, qq(4), o(3,3)
  -real(kind=sgl)                          :: VL(3,3), VR(3,3), Wr(3), Wi(3), WORK(10)
  +real(kind=sgl)                          :: VL(3,3), VR(3,3), Wr(3), Wi(3), WORK(20)
   complex(kind=sgl)                       :: ev
   complex(kind=sgl),parameter             :: cone = cmplx(1.0,0.0)
   real(kind=sgl),parameter                :: thr = 1.0E-7
  
positive.code: ! 
  !> @date 8/12/13  MDG 1.0 original
  !> @date 07/08/14 MDG 2.0 replaced by direct solution
  !--------------------------------------------------------------------------
  recursive function om2ax(om) result(res)
  
  
negative.code: !> @date 8/12/13  MDG 1.0 original
  !> @date 07/08/14 MDG 2.0 replaced by direct solution
  !> @date 08/20/14 MDG 3.0 replaced by eigenvalue-based method
  !> @date 09/23/19 MDG 3.1 make size(WORK,1) and LWORK equal (comment by M. Diehl)
  !--------------------------------------------------------------------------
  recursive function om2ax_d(om) result(res)
  
  