anchor.diff:      ! User-defined variables
       real :: hue1,hue2,light1,light2,chroma1,chroma2,lim1,lim2
       character(len=200) :: option,convert1,convert2,gmt
  -    integer :: ncolors
  +    integer :: ncolors,saturate
   
       ! Parse command line
       call gcmdln(option,hue1,hue2,light1,light2,chroma1,chroma2,ncolors,convert1,convert2,&
  -                gmt,lim1,lim2)
  +                gmt,lim1,lim2,saturate)
       if (vrb.gt.0) then
           write(0,*) 'Color mapping option: ',trim(option)
           write(0,*) 'Hue range:            ',hue1,hue2
  
positive.code:             read(tag,*) lim1,lim2,dz
          elseif (tag(1:2).eq.'-v') then
              vrb = 1
          elseif (tag(1:2).eq.'-d') then
              call usage('long')
          else
  
negative.code: !~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~!
  
  subroutine gcmdln(option,hue1,hue2,light1,light2,chroma1,chroma2,ncolors,convert1,convert2,&
                    gmt,lim1,lim2,saturate)
      use colormodule
      implicit none
      integer :: i,j,narg
  