/* Color mapping appropriate for one-sided data. */
#include <math.h>
#ifndef M_PI
#define M_PI 3.14159265358979323846
#endif

void cmap1(float *h,float *s,float *v,float xpnt) {
  float t;
  if (xpnt > 0.8) {
    *h = 330. + 30*(1.-xpnt)/0.2;
    *s = (1.-xpnt)/0.2;
    *v = 1.0;
  } else if (xpnt > 0.2) {
    t = 2.0*M_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;
  }
  return;
}