MZ_DLLEXPORT
void c_plshades( PLFLT **a, PLINT nx, PLINT ny, PLINT (*defined) (PLFLT, PLFLT),
PLFLT xmin, PLFLT xmax, PLFLT ymin, PLFLT ymax,
PLFLT *clevel, PLINT nlevel, PLINT fill_width,
PLINT cont_color, PLINT cont_width,
void (*fill) (PLINT, PLFLT *, PLFLT *), PLINT rectangular,
void (*pltr) (PLFLT, PLFLT, PLFLT *, PLFLT *, PLPointer),
PLPointer pltr_data )
{
PLFLT shade_min, shade_max, shade_color;
PLINT i, init_color, init_width;
for (i = 0; i < nlevel-1; i++) {
shade_min = clevel[i];
shade_max = clevel[i+1];
shade_color = i / (PLFLT) (nlevel-2);
/* The constants in order mean
* (1) color map1,
* (0, 0, 0, 0) all edge effects will be done with plcont rather
* than the normal plshade drawing which gets partially blocked
* when sequential shading is done as in the present case */
plshade(a, nx, ny, defined, xmin, xmax, ymin, ymax,
shade_min, shade_max,
1, shade_color, fill_width,
0, 0, 0, 0,
fill, rectangular, pltr, pltr_data);
}
if(cont_color > 0 && cont_width > 0) {
init_color = plsc->icol0;
init_width = plsc->width;
plcol0(cont_color);
plwid(cont_width);
plcont(a, nx, ny, 1, nx, 1, ny, clevel, nlevel, pltr, pltr_data);
plcol0(init_color);
plwid(init_width);
}
}
void
shade(void)
{
int i, j;
PLFLT x, y, argx, argy, distort;
PLFLT **z, **w, zmin, zmax;
PLFLT xg1[XPTS], yg1[YPTS];
PLcGrid cgrid1;
PLcGrid2 cgrid2;
PLFLT shade_min, shade_max, sh_color;
PLINT sh_cmap = 1, sh_width;
PLINT min_color = 1, min_width = 0, max_color = 0, max_width = 0;
/* Set up function arrays */
plAlloc2dGrid(&z, XPTS, YPTS);
plAlloc2dGrid(&w, XPTS, YPTS);
/* Set up data array */
for (i = 0; i < XPTS; i++) {
x = (double) (i - (XPTS / 2)) / (double) (XPTS / 2);
for (j = 0; j < YPTS; j++) {
y = (double) (j - (YPTS / 2)) / (double) (YPTS / 2) - 1.0;
z[i][j] = - sin(7.*x) * cos(7.*y) + x*x - y*y;
w[i][j] = - cos(7.*x) * sin(7.*y) + 2 * x * y;
}
}
f2mnmx(z, XPTS, YPTS, &zmin, &zmax);
for (i = 0; i < NCONTR; i++)
clevel[i] = zmin + (zmax - zmin) * (i + 0.5) / (PLFLT) NCONTR;
/* Set up coordinate grids */
cgrid1.xg = xg1;
cgrid1.yg = yg1;
cgrid1.nx = XPTS;
cgrid1.ny = YPTS;
plAlloc2dGrid(&cgrid2.xg, XPTS, YPTS);
plAlloc2dGrid(&cgrid2.yg, XPTS, YPTS);
cgrid2.nx = XPTS;
cgrid2.ny = YPTS;
for (i = 0; i < XPTS; i++) {
for (j = 0; j < YPTS; j++) {
mypltr((PLFLT) i, (PLFLT) j, &x, &y, NULL);
argx = x * PI/2;
argy = y * PI/2;
distort = 0.4;
cgrid1.xg[i] = x + distort * cos(argx);
cgrid1.yg[j] = y - distort * cos(argy);
cgrid2.xg[i][j] = x + distort * cos(argx) * cos(argy);
cgrid2.yg[i][j] = y - distort * cos(argx) * cos(argy);
}
}
/* Plot using identity transform */
pladv(0);
plvpor(0.1, 0.9, 0.1, 0.9);
plwind(-1.0, 1.0, -1.0, 1.0);
for (i = 0; i < NCONTR; i++) {
shade_min = zmin + (zmax - zmin) * i / (PLFLT) NCONTR;
shade_max = zmin + (zmax - zmin) * (i +1) / (PLFLT) NCONTR;
sh_color = i / (PLFLT) (NCONTR-1);
sh_width = 2;
plpsty(0);
plshade(z, XPTS, YPTS, NULL, -1., 1., -1., 1.,
shade_min, shade_max,
sh_cmap, sh_color, sh_width,
min_color, min_width, max_color, max_width,
plfill, 1, NULL, NULL);
}
plcol(1);
plbox("bcnst", 0.0, 0, "bcnstv", 0.0, 0);
plcol(2);
/*
plcont(w, XPTS, YPTS, 1, XPTS, 1, YPTS, clevel, NCONTR, mypltr, NULL);
*/
pllab("distance", "altitude", "Bogon density");
/* Clean up */
plFree2dGrid(z, XPTS, YPTS);
plFree2dGrid(w, XPTS, YPTS);
plFree2dGrid(cgrid2.xg, XPTS, YPTS);
plFree2dGrid(cgrid2.yg, XPTS, YPTS);
}
