anchor.diff:  ! Check that RGB values are in range 0-255 and return whether they needed to be clipped
   use colormodule
   implicit none
  -real :: r,g,b
  +double precision :: r,g,b
   integer :: clip(3)
   clip = 0
   if (r.lt.0.0) then
  
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 lch2rgb(lin,cin,hin,r,g,b)
  ! Lightness-chroma-hue to red-green-blue
  implicit none
  double precision :: lin,cin,hin,r,g,b
  double precision :: a1,b1
  call lch2lab(cin,hin,a1,b1)
  call lab2rgb(lin,a1,b1,r,g,b)
  return
  