28char des_equipot_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Graphics/des_equipot.C,v 1.5 2014/10/13 08:53:22 j_novak Exp $" ;
74void des_equipot(
float* uutab,
int nx,
int ny,
float xmin,
float xmax,
75 float ymin,
float ymax,
int ncour,
const char* nomx,
const char* nomy,
76 const char* title,
const char* device,
int newgraph,
int nxpage,
82 float uumin = uutab[0] ;
83 float uumax = uutab[0] ;
84 for (
int i=1; i<nx*ny; i++) {
85 uumin = (uutab[i] < uumin) ? uutab[i] : uumin ;
86 uumax = (uutab[i] > uumax) ? uutab[i] : uumax ;
89 cout <<
" " << title <<
" : min, max : " << uumin <<
" " << uumax
95 float* isopot =
new float [ncour] ;
96 float hh = (uumax-uumin) /
float(ncour) ;
97 for (
int i=0; i<ncour; i++) {
98 isopot[i] = uumin + hh * float(i) ;
103 float hx = (xmax - xmin)/
float(nx-1) ;
104 float hy = (ymax - ymin)/
float(ny-1) ;
117 if ( (newgraph == 1) || (newgraph == 3) ) {
119 if (device == 0x0) device =
"?" ;
121 int ier = cpgbeg(0, device, nxpage, nypage) ;
123 cout <<
"des_equipot: problem in opening PGPLOT display !" << endl ;
129 float size = float(1.3) ;
140 cpgenv(xmin, xmax, ymin, ymax, 1, 0 ) ;
141 cpglab(nomx,nomy,title) ;
145 float dynamique = float(fabs(uumax - uumin)) ;
147 if (dynamique > 1.e-14) {
149 cpgcont(uutab, nx, ny, 1, nx, 1, ny, isopot, ncour, tr) ;
156 if ( (newgraph == 2) || (newgraph == 3) ) {
void des_equipot(float *uutab, int nx, int ny, float xmin, float xmax, float ymin, float ymax, int ncour, const char *nomx, const char *nomy, const char *title, const char *device=0x0, int newgraph=3, int nxpage=1, int nypage=1)
Basic routine for drawing isocontours.