anchor.diff:  ! Local
   integer :: i
   character(len=3) :: typ1, typ2
  -real :: x1, y1, z1
  -real :: x2, y2, z2
  +double precision :: x1, y1, z1
  +double precision :: x2, y2, z2
   integer :: clip(3)
   
   ! Check formatting of input strings
  
positive.code: subroutine xyz2rgb(xin,yin,zin,r,g,b)
  ! X-y-z tristimulus to red-green-blue
  implicit none
  real :: xin,yin,zin,x,y,z,r,g,b
  x = xin/100.0
  y = yin/100.0
  z = zin/100.0
  
negative.code: 
  subroutine rgb2lch(red,grn,blu,l,c,h)
  implicit none
  double precision :: red,grn,blu,l,c,h
  double precision :: a,b
  call rgb2lab(red,grn,blu,l,a,b)
  call lab2lch(a,b,c,h)
  return
  