C.....Color mapping appropriate for two-sided data.  Zero = black.
      subroutine cmap2(h,s,v,xpnt)
      real h, s, v, xpnt
      if (xpnt .gt. 0.8) then
         h = 30. - 30.*(xpnt-0.8)*(xpnt-0.8)/0.04
         s = (1.0-xpnt)/0.2
         v = 1.0
      else if (xpnt .gt. .3) then
         h = 135.0 - 105.*(xpnt-.3)/.5
         s = 1.0
         v = 1.0
      else if (xpnt .gt. 0.) then
         h = 180.0 - 45.*xpnt/.3
         s = 1.0
         v = xpnt/.3
      else if (xpnt .gt. -.3) then
         h = (180. - 45.*xpnt/.3)
         s = 1.0
         v = abs(xpnt)/.3
      else
         h = 225. - 135.*(xpnt+.3)/.7
         s = 1.0
         v = 1.0
      endif
      return
      end