1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
|
/* $Id: x19c.c,v 1.10 2002/12/17 20:48:14 airwin Exp $
Illustrates backdrop plotting of world, US maps.
Contributed by Wesley Ebisuzaki.
*/
#include "plcdemos.h"
/*--------------------------------------------------------------------------*\
* mapform19
*
* Defines specific coordinate transformation for example 19.
* Not to be confused with mapform in src/plmap.c.
* x[], y[] are the coordinates to be plotted.
\*--------------------------------------------------------------------------*/
void
mapform19(PLINT n, PLFLT *x, PLFLT *y)
{
int i;
double xp, yp, radius;
for (i = 0; i < n; i++) {
radius = 90.0 - y[i];
xp = radius * cos(x[i] * PI / 180.0);
yp = radius * sin(x[i] * PI / 180.0);
x[i] = xp;
y[i] = yp;
}
}
/*--------------------------------------------------------------------------*\
* main
*
* Shows two views of the world map.
\*--------------------------------------------------------------------------*/
int
main(int argc, char **argv)
{
PLFLT minx, maxx, miny, maxy;
int c;
/* Parse and process command line arguments */
(void) plParseOpts(&argc, argv, PL_PARSE_FULL);
/* Longitude (x) and latitude (y) */
miny = -70;
maxy = 80;
plinit();
/* Cartesian plots */
/* Most of world */
minx = 190;
maxx = 190+360;
plcol0(1);
plenv(minx, maxx, miny, maxy, 1, -1);
plmap(NULL, "usaglobe", minx, maxx, miny, maxy);
/* The Americas */
minx = 190;
maxx = 340;
plcol0(1);
plenv(minx, maxx, miny, maxy, 1, -1);
plmap(NULL, "usaglobe", minx, maxx, miny, maxy);
/* Polar, Northern hemisphere */
minx = 0;
maxx = 360;
plenv(-75., 75., -75., 75., 1, -1);
plmap(mapform19,"globe", minx, maxx, miny, maxy);
pllsty(2);
plmeridians(mapform19,10.0, 10.0, 0.0, 360.0, -10.0, 80.0);
plend();
exit(0);
}
|
