anchor.diff:  double precision :: hue(2)
   double precision :: lightness(2)
   double precision :: chroma(2)
  +double precision :: chroma_exp
   character(len=16) :: color_path_type
   character(len=8) :: convert(2)
   double precision :: gmt_cpt_lims(2)
  +character(len=8) :: output_mode
   
   integer :: verbose
   
  
positive.code: integer :: i
  integer :: uout
  double precision :: lcheck,ccheck,hcheck
  
  ! Parse command line
  call gcmdln(option,hue1,hue2,light1,light2,chroma1,chroma2,ncolors,convert1,convert2,&
  
negative.code:     call lch2rgb(l,c,h,red,grn,blu)
      call checkrgb(red,grn,blu,clip)
      call rgb2lch(red,grn,blu,lcheck,ccheck,hcheck)
      if (output_mode.eq.'RGB') then
          ! Print nicely formatted RGB coordinates
          call formatrgb(red,grn,blu,rgbstring)
      elseif (output_mode.eq.'HSV') then
          ! Print nicely formatted HSV coordinates
          call rgb2hsv(red,grn,blu,hsv1,hsv2,hsv3)
          call formathsv(hsv1,hsv2,hsv3,rgbstring)
      else
          call usage('colortool: no output_mode "'//trim(output_mode)//'"')
      endif
      if (gmtfile.eq.'none') then
          write(uout,'(A)') trim(rgbstring)
      else
  