C.....Color mapping appropriate for one-sided data. subroutine cmap1(h,s,v,xpnt) real h, s, v, xpnt, pi parameter (pi=3.141593) if (xpnt .gt. 0.8) then h = 330. + 30*(1.-xpnt)/0.2 s = (1.-xpnt)/0.2 v = 1.0 elseif (xpnt .gt. 0.2) then t = 2.0*pi*((xpnt-0.2)/.6) h = 240.0*(0.8-xpnt)/0.6 + 5.0*sin(t) s = 1.0 v = 1.0 else h = 270.0*(.2-xpnt)+240.0 s = 1.0 v = xpnt/0.2 endif return end