2 * For RSL version 0.28 and higher.
4 * v1.0 Began 2/16/94 by John Merritt.
6 * Demonstrates reading NEXRAD files and loading the Radar structure.
13 usage(int argc, char **argv)
15 fprintf(stderr,"Usage: %s infile [firstfile | callid]\n", argv[0]);
19 process_args(int argc, char **argv,
20 char **in_file, char **out_file)
24 *in_file = strdup(argv[1]);
26 } else if (argc == 3) {
27 *in_file = strdup(argv[1]);
28 *out_file = strdup(argv[2]);
35 main(int argc, char **argv)
37 char *infile, *firstfile;
45 /* 1. Process the arguments. */
46 process_args(argc, argv, &infile, &firstfile); /* malloc for in/outfile */
51 RSL_radar_verbose_on();
52 if ((radar = RSL_anyformat_to_radar(infile, firstfile)) == NULL) exit(-1);
54 printf("Radar date: %2.2d/%2.2d/%2.2d\n", radar->h.month, radar->h.day, radar->h.year);
55 printf("Radar time: %2.2d:%2.2d:%f\n", radar->h.hour, radar->h.minute, radar->h.sec);
56 printf("Radar file: %s\n", infile);
57 printf("Radar site: %c%c%c%c\n",
62 printf("Radar date: %2.2d/%2.2d/%2.2d\n", radar->h.month, radar->h.day, radar->h.year);
63 printf("Radar time: %2.2d:%2.2d:%f\n", radar->h.hour, radar->h.minute, radar->h.sec);
66 if (i == DZ_INDEX) RSL_load_refl_color_table();
67 if (i == VR_INDEX) RSL_load_vel_color_table();
68 if (i == SW_INDEX) RSL_load_sw_color_table();
69 if (i == VR_INDEX) RSL_rebin_velocity_volume(radar->v[i]); /* Modifies v[i]. */
71 for (j=1; j<=5; j++) {
72 cappi = RSL_cappi_at_h( radar -> v[i], (float)j/2.0, 200.0);
73 sprintf(outfile, "cappi_%2.2d.gif", j);
74 printf("Writing %s.\n", outfile);
75 RSL_sweep_to_gif(cappi->sweep, outfile, 400, 400, 200.0);
76 RSL_free_cappi(cappi);