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.
22 fprintf(stderr,"Usage: bscan infile [callid]\n");
26 process_args(int argc, char **argv,
27 char **in_file, char **callid)
30 if (argc < 2) usage();
31 else if (argc == 2) *in_file = strdup(argv[1]);
33 *in_file = strdup(argv[1]);
34 *callid = strdup(argv[2]);
41 main(int argc, char **argv)
43 char *infile, *callid;
47 char *index_str[] = {"DZ", "VR", "SW"};
51 /* 1. Process the arguments. */
53 process_args(argc, argv, &infile, &callid); /* malloc for in/outfile */
57 * Pass bitwise or of DZ_MASK, VR_MASK, SW_MASK
61 RSL_radar_verbose_on();
62 if ((radar = RSL_anyformat_to_radar(infile, callid)) == NULL) exit(-1);
64 printf("Radar date: %2.2d/%2.2d/%2.2d\n", radar->h.month, radar->h.day, radar->h.year);
65 printf("Radar time: %2.2d:%2.2d:%f\n", radar->h.hour, radar->h.minute, radar->h.sec);
67 printf("Radar file: %s\n", infile);
68 printf("Radar site: %c%c%c%c\n",
73 printf("Radar date: %2.2d/%2.2d/%2.2d\n", radar->h.month, radar->h.day, radar->h.year);
74 printf("Radar time: %2.2d:%2.2d:%f\n", radar->h.hour, radar->h.minute, radar->h.sec);
77 if (i == DZ_INDEX) RSL_load_refl_color_table();
78 if (i == VR_INDEX) RSL_load_vel_color_table();
79 if (i == SW_INDEX) RSL_load_sw_color_table();
80 if (i == VR_INDEX) RSL_rebin_velocity_volume(radar->v[i]); /* Modifies v[i]. */
82 printf("Generating bscan ppm images of %s\n", index_str[i]);
83 RSL_bscan_volume((Volume *) radar->v[i], "bscan.ppm");
84 printf("----> BSCAN complete.\n");