C.....Color mapping appropriate for two-sided data.  zero = blue.
      subroutine cmap2(h,s,v,xpnt)
      real h, s, v, xpnt
      if (xpnt .gt. 0.8) then
         h = 330. + 30*(1.0-xpnt)/0.2
         s = (1.0-xpnt)/0.2
         v = 1.0
      else
         if (xpnt .lt. 0.) xpnt = 1.5*xpnt
         h = (240.0 + 120.*xpnt)
         s = 1.0
         v = 1.0
      endif
      return
      end