home *** CD-ROM | disk | FTP | other *** search
- /* Draws a contour plot from data in z(nx,ny), using the subarray */
- /* from kx to lx in the x direction and from ky to ly in the y */
- /* direction. The array of contour levels is clevel(nlevel), and */
- /* "tr" is the name of a subroutine which transforms array indicies */
- /* into world coordinates */
-
- #include <stdio.h>
- #include <stdlib.h>
- #include "plplot.h"
- #include "declare.h"
-
- void plcont(z,nx,ny,kx,lx,ky,ly,clevel,nlevel,tr)
- int nx, ny, kx, lx, ky, ly, nlevel;
- float *z, clevel[];
- void (*tr)();
- {
- int i, mx, my, nstor;
-
- mx = lx - kx + 1;
- my = ly - ky + 1;
-
- if (kx < 1 || lx > nx || kx >= lx || ky < 1 || ky > ny || ky >= ly)
- fatal("Argument error in PLCONT");
-
- nstor = mx*my/5;
-
- if( heapc != NULL)
- free(heapc);
- if(( heapc = (int *)malloc((mx+2*nstor)*sizeof(int))) == NULL)
- fatal("Not enough heap space in PLCONT.");
-
- for (i=0; i<nlevel; i++)
- plcntr(z,nx,ny,kx,lx,ky,ly,clevel[i],&heapc[0],
- &heapc[nx],&heapc[nx+nstor],nstor,tr);
-
- }
-