anchor.diff:  subroutine lch2rgb(lin,cin,hin,r,g,b)
   ! Lightness-chroma-hue to red-green-blue
   implicit none
  -real :: lin,cin,hin,r,g,b
  -real :: a1,b1
  +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
  
positive.code:     ! Calculate RGB coordinates
      call lch2rgb(l,c,h,red,grn,blu)
      call checkrgb(red,grn,blu,clip)
      ! Print nicely formatted RGB coordinates
      call formatrgb(red,grn,blu,rgbstring)
      if (gmtfile.eq.'none') then
  
negative.code: ! Check that RGB values are in range 0-255 and return whether they needed to be clipped
  use colormodule
  implicit none
  double precision :: r,g,b
  integer :: clip(3)
  clip = 0
  if (r.lt.0.0) then
  