static void MakeColorMap (ColorMap& cmap)
{
  cmap.Load("copper-256"); 			       // load a default colormap
  if (! cmap) Matpack.Error("Can't load colormap");
  cmap.Reverse();					      // reverse ordering
  int i, n = cmap.Size();				   // change colormap now
  for (i = 0;  i < 20; i++) cmap[i] = ColorB(0,100+2*i,160+2*i);         // water
  for (i = 20; i < 80; i++) cmap[i] = ColorB(2*i+20,2*i+20,i/2+10);   // mangrove
  for (i = 1;  i < 60; i++) cmap[n-i] = ColorB(255-i,250-i,250-2*i);      // snow
}