2 * Ingest NEXRAD (wsr88d) data and print the azimuth hash table created.
4 * This example is the most minimum of coding that you need to do
5 * to achieve good results from using the RSL code.
7 * This is short and sweet to demonstrate the simplicity of use for
14 void print_link_list(Azimuth_hash *list)
20 printf("ray# %d azim %f, hi# %d lo# %d|", list->ray->h.ray_num, list->ray->h.azimuth, list->ray_high->ray->h.ray_num, list->ray_low->ray->h.ray_num);
21 print_link_list(list->next);
26 * Cannot compile if the hash table is static in volume.c. For
27 * testing we make it globally known.
35 extern int RSL_max_sweeps; /* Initial allocation for sweep_list.
36 * RSL_new_sweep will allocate the space first
39 extern int RSL_nsweep_addr; /* A count of sweeps in the table. */
40 extern Sweep_list *RSL_sweep_list;
41 extern int RSL_nextents;
43 void print_hash_table (Sweep *s)
51 if (s == NULL) return;
52 sweep_index = SWEEP_INDEX(s);
53 res = 360.0/RSL_sweep_list[sweep_index].hash->nindexes;
54 printf("Azimuth resolution = %f for %d bins.\n", res, RSL_sweep_list[sweep_index].hash->nindexes);
55 for (i=0; i<RSL_sweep_list[sweep_index].hash->nindexes; i++) {
56 index = RSL_sweep_list[sweep_index].hash->indexes[i];
58 printf("RSL_sweep_list[%d].hash->indexes[%d] = ", sweep_index, i);
63 print_link_list(index);
67 void poke_about_sweep(Sweep *s)
69 /* This routine demonstrates that the azimuth we want is the azimuth
75 ray = RSL_get_first_ray_of_sweep(s);
76 res = 360.0/s->h.nrays;
77 for (azim=0; azim<360; azim+=res) {
78 ray = RSL_get_ray_from_sweep(s, azim);
80 printf("Azimuth %f matched in ray # %d, h.azimuth= %f, diff=%f\n",
81 azim, ray->h.ray_num, ray->h.azimuth,
84 printf("Azimuth %f NOT FOUND within 1/2 beamwidth; 1/2beam=%f\n",
85 azim, s->h.horz_half_bw);
90 void main(int argc, char **argv)
95 if (argc != 3) {fprintf(stderr, "Usage: %s infile callid_or_firstfile\n", argv[0]); exit(-1);}
96 RSL_radar_verbose_on();
97 radar = RSL_wsr88d_to_radar(argv[1], argv[2]);
98 if (radar == NULL) exit(-1);
99 /***********************************************************************/
101 /* You now have a pointer to Radar. */
102 /* Now use *radar all you like. */
104 /***********************************************************************/
106 /* Use radar->v[DZ_INDEX] for REFELECTIVITY
107 * radar->v[VR_INDEX] for VELOCITY
108 * radar->v[SW_INDEX] for SPECTRUM_WIDTH
111 sweep = RSL_get_first_sweep_of_volume(radar->v[DZ_INDEX]);
112 print_hash_table(sweep);
114 poke_about_sweep(sweep);