#include <stdlib.h>
#include <unistd.h>
+/* This allows us to use RSL_ftype, RSL_f_list, RSL_invf_list from rsl.h. */
+#define USE_RSL_VARS
#include "rsl.h"
extern int radar_verbose_flag;
-typedef short UF_buffer[16384]; /* Some UF files are bigger than 4096
- * that the UF doc's specify.
- */
+/* Changed old buffer size (16384) for larger dualpol files. BLK 5/18/2011 */
+typedef short UF_buffer[20000]; /* Some UF files are bigger than 4096
+ * that the UF doc's specify.
+ */
#define UF_MORE 0
#define UF_DONE 1
int i;
if (volume == NULL) return NULL;
for (i=volume->h.nsweeps; i>0; i--)
- if (volume->sweep[i-1] != NULL) {
- volume->h.nsweeps = i;
- break;
- }
+ if (volume->sweep[i-1] != NULL) {
+ volume->h.nsweeps = i;
+ break;
+ }
return volume;
}
Radar *reset_nsweeps_in_all_volumes(Radar *radar)
int i;
if (radar == NULL) return NULL;
for (i=0; i<radar->h.nvolumes; i++)
- radar->v[i] = reset_nsweeps_in_volume(radar->v[i]);
+ radar->v[i] = reset_nsweeps_in_volume(radar->v[i]);
return radar;
}
new_volume->h = old_volume->h;
new_volume->h.nsweeps = nsweeps;
for (i=0; i<old_volume->h.nsweeps; i++)
- new_volume->sweep[i] = old_volume->sweep[i]; /* Just copy pointers. */
+ new_volume->sweep[i] = old_volume->sweep[i]; /* Just copy pointers. */
/* Free the old sweep array, not the pointers to sweeps. */
free(old_volume->sweep);
return new_volume;
short *end_addr;
end_addr = buf + n;
while (buf < end_addr) {
- swap_2_bytes(buf);
- buf++;
+ swap_2_bytes(buf);
+ buf++;
}
}
int current_fh_index;
float scale_factor;
- int nfields, isweep, ifield, iray, i, m;
+ int nfields, isweep, ifield, iray, i, j, m;
static int pulled_time_from_first_ray = 0;
+ static int need_scan_mode = 1;
char *field_type; /* For printing the field type upon error. */
- short proj_name[4];
+ short proj_name[4];
Ray *ray;
Sweep *sweep;
Radar *radar;
isweep = uf_ma[9] - 1;
if (rsl_qsweep != NULL) {
- if (isweep > rsl_qsweep_max) return UF_DONE;
- if (rsl_qsweep[isweep] == 0) return UF_MORE;
+ if (isweep > rsl_qsweep_max) return UF_DONE;
+ if (rsl_qsweep[isweep] == 0) return UF_MORE;
}
/* Sticky solution here. */
#else
if (radar == NULL) {
- radar = RSL_new_radar(MAX_RADAR_VOLUMES);
- *the_radar = radar;
- pulled_time_from_first_ray = 0;
- for (i=0; i<MAX_RADAR_VOLUMES; i++)
- if (rsl_qfield[i]) /* See RSL_select_fields in volume.c */
- radar->v[i] = RSL_new_volume(20);
+ radar = RSL_new_radar(MAX_RADAR_VOLUMES);
+ *the_radar = radar;
+ pulled_time_from_first_ray = 0;
+ for (i=0; i<MAX_RADAR_VOLUMES; i++)
+ if (rsl_qfield[i]) /* See RSL_select_fields in volume.c */
+ radar->v[i] = RSL_new_volume(20);
}
#endif
+
+ if (need_scan_mode) {
+ /* PPI and RHI are enum constants defined in rsl.h */
+ if (uf_ma[34] == 1) radar->h.scan_mode = PPI;
+ else if (uf_ma[34] == 3) radar->h.scan_mode = RHI;
+ else {
+ fprintf(stderr,"Warning: UF sweep mode = %d\n", uf_ma[34]);
+ fprintf(stderr," Expected 1 or 3 (PPI or RHI)\n");
+ fprintf(stderr," Setting radar->h.scan_mode to PPI\n");
+ radar->h.scan_mode = PPI;
+ }
+ need_scan_mode = 0;
+ }
+
/* For LITTLE ENDIAN:
* WE "UNSWAP" character strings. Because there are so few of them,
* it is easier to catch them here. The entire UF buffer is swapped prior
*/
for (i=0; i<nfields; i++) {
- if (little_endian()) swap_2_bytes(&uf_dh[3+2*i]); /* Unswap. */
- if (strncmp((char *)&uf_dh[3+2*i], "DZ", 2) == 0) ifield = DZ_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "VR", 2) == 0) ifield = VR_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "SW", 2) == 0) ifield = SW_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "CZ", 2) == 0) ifield = CZ_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "UZ", 2) == 0) ifield = ZT_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "ZT", 2) == 0) ifield = ZT_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "DR", 2) == 0) ifield = DR_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "ZD", 2) == 0) ifield = ZD_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "DM", 2) == 0) ifield = DM_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "RH", 2) == 0) ifield = RH_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "PH", 2) == 0) ifield = PH_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "XZ", 2) == 0) ifield = XZ_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "CD", 2) == 0) ifield = CD_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "MZ", 2) == 0) ifield = MZ_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "MD", 2) == 0) ifield = MD_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "ZE", 2) == 0) ifield = ZE_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "VE", 2) == 0) ifield = VE_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "KD", 2) == 0) ifield = KD_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "TI", 2) == 0) ifield = TI_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "DX", 2) == 0) ifield = DX_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "CH", 2) == 0) ifield = CH_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "AH", 2) == 0) ifield = AH_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "CV", 2) == 0) ifield = CV_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "AV", 2) == 0) ifield = AV_INDEX;
- else if (strncmp((char *)&uf_dh[3+2*i], "SQ", 2) == 0) ifield = SQ_INDEX;
- else { /* etc. DON'T know how to handle this yet. */
- field_type = (char *)&uf_dh[3+2*i];
- fprintf(stderr, "Unknown field type %c%c\n", (char)field_type[0], (char)field_type[1]);
- continue;
- }
- switch (ifield) {
- case DZ_INDEX: f = DZ_F; invf = DZ_INVF; break;
- case VR_INDEX: f = VR_F; invf = VR_INVF; break;
- case SW_INDEX: f = SW_F; invf = SW_INVF; break;
- case CZ_INDEX: f = CZ_F; invf = CZ_INVF; break;
- case ZT_INDEX: f = ZT_F; invf = ZT_INVF; break;
- case DR_INDEX: f = DR_F; invf = DR_INVF; break;
- case ZD_INDEX: f = ZD_F; invf = ZD_INVF; break;
- case DM_INDEX: f = DM_F; invf = DM_INVF; break;
- case RH_INDEX: f = RH_F; invf = RH_INVF; break;
- case PH_INDEX: f = PH_F; invf = PH_INVF; break;
- case XZ_INDEX: f = XZ_F; invf = XZ_INVF; break;
- case CD_INDEX: f = CD_F; invf = CD_INVF; break;
- case MZ_INDEX: f = MZ_F; invf = MZ_INVF; break;
- case MD_INDEX: f = MD_F; invf = MD_INVF; break;
- case ZE_INDEX: f = ZE_F; invf = ZE_INVF; break;
- case VE_INDEX: f = VE_F; invf = VE_INVF; break;
- case KD_INDEX: f = KD_F; invf = KD_INVF; break;
- case TI_INDEX: f = TI_F; invf = TI_INVF; break;
- case DX_INDEX: f = DX_F; invf = DX_INVF; break;
- case CH_INDEX: f = CH_F; invf = CH_INVF; break;
- case AH_INDEX: f = AH_F; invf = AH_INVF; break;
- case CV_INDEX: f = CV_F; invf = CV_INVF; break;
- case AV_INDEX: f = AV_F; invf = AV_INVF; break;
- case SQ_INDEX: f = SQ_F; invf = SQ_INVF; break;
- default: f = DZ_F; invf = DZ_INVF; break;
- }
-
- /* Do we place the data into this volume? */
- if (radar->v[ifield] == NULL) continue; /* Nope. */
-
- if (isweep >= radar->v[ifield]->h.nsweeps) { /* Exceeded sweep limit.
- * Allocate more sweeps.
- * Copy all previous sweeps.
- */
- if (radar_verbose_flag)
- fprintf(stderr,"Exceeded sweep allocation of %d. Adding 20 more.\n", isweep);
- new_volume = RSL_new_volume(radar->v[ifield]->h.nsweeps+20);
- new_volume = copy_sweeps_into_volume(new_volume, radar->v[ifield]);
- radar->v[ifield] = new_volume;
- }
-
- if (radar->v[ifield]->sweep[isweep] == NULL) {
- if (radar_verbose_flag)
- fprintf(stderr,"Allocating new sweep for field %d, isweep %d\n", ifield, isweep);
- radar->v[ifield]->sweep[isweep] = RSL_new_sweep(1000);
- radar->v[ifield]->sweep[isweep]->h.nrays = 0; /* Increment this for each
- * ray encountered.
- */
- radar->v[ifield]->h.f = f;
- radar->v[ifield]->h.invf = invf;
- radar->v[ifield]->sweep[isweep]->h.f = f;
- radar->v[ifield]->sweep[isweep]->h.invf = invf;
- radar->v[ifield]->sweep[isweep]->h.sweep_num = uf_ma[9];
- radar->v[ifield]->sweep[isweep]->h.elev = uf_ma[35] / 64.0;
- }
-
-
-
- current_fh_index = uf_dh[4+2*i];
- uf_fh = uf + current_fh_index - 1;
- sweep = radar->v[ifield]->sweep[isweep];
- iray = sweep->h.nrays;
- nbins = uf_fh[5];
- radar->v[ifield]->sweep[isweep]->ray[iray] = RSL_new_ray(nbins);
- ray = radar->v[ifield]->sweep[isweep]->ray[iray];
- sweep->h.nrays += 1;
-
-
- if (ray) {
- /*
- * ---- Beginning of MANDATORY HEADER BLOCK.
- */
- ray->h.ray_num = uf_ma[7];
- if (little_endian()) swap2(&uf_ma[10], 8);
- memcpy(radar->h.radar_name, &uf_ma[10], 8);
- if (little_endian()) swap2(&uf_ma[10], 8/2);
- memcpy(radar->h.name, &uf_ma[14], 8);
- if (little_endian()) swap2(&uf_ma[14], 8/2);
-
+ if (little_endian()) swap_2_bytes(&uf_dh[3+2*i]); /* Unswap. */
+ ifield = -1;
+ field_type = (char *)&uf_dh[3+2*i];
+ for (j=0; j<MAX_RADAR_VOLUMES; j++) {
+ if (strncmp(field_type, RSL_ftype[j], 2) == 0) {
+ ifield = j;
+ break;
+ }
+ }
+ if (ifield < 0) { /* DON'T know how to handle this yet. */
+ fprintf(stderr, "Unknown field type %c%c\n", (char)field_type[0],
+ (char)field_type[1]);
+ continue;
+ }
+
+ f = RSL_f_list[ifield];
+ invf = RSL_invf_list[ifield];
+
+ /* Do we place the data into this volume? */
+ if (radar->v[ifield] == NULL) continue; /* Nope. */
+
+ if (isweep >= radar->v[ifield]->h.nsweeps) { /* Exceeded sweep limit.
+ * Allocate more sweeps.
+ * Copy all previous sweeps.
+ */
+ if (radar_verbose_flag)
+ fprintf(stderr,"Exceeded sweep allocation of %d. Adding 20 more.\n", isweep);
+ new_volume = RSL_new_volume(radar->v[ifield]->h.nsweeps+20);
+ new_volume = copy_sweeps_into_volume(new_volume, radar->v[ifield]);
+ radar->v[ifield] = new_volume;
+ }
+
+ if (radar->v[ifield]->sweep[isweep] == NULL) {
+ if (radar_verbose_flag)
+ fprintf(stderr,"Allocating new sweep for field %d, isweep %d\n", ifield, isweep);
+ radar->v[ifield]->sweep[isweep] = RSL_new_sweep(1000);
+ radar->v[ifield]->sweep[isweep]->h.nrays = 0; /* Increment this for each
+ * ray encountered.
+ */
+ radar->v[ifield]->h.f = f;
+ radar->v[ifield]->h.invf = invf;
+ radar->v[ifield]->sweep[isweep]->h.f = f;
+ radar->v[ifield]->sweep[isweep]->h.invf = invf;
+ radar->v[ifield]->sweep[isweep]->h.sweep_num = uf_ma[9];
+ radar->v[ifield]->sweep[isweep]->h.elev = uf_ma[35] / 64.0;
+ }
+
+
+
+ current_fh_index = uf_dh[4+2*i];
+ uf_fh = uf + current_fh_index - 1;
+ sweep = radar->v[ifield]->sweep[isweep];
+ iray = sweep->h.nrays;
+ nbins = uf_fh[5];
+ radar->v[ifield]->sweep[isweep]->ray[iray] = RSL_new_ray(nbins);
+ ray = radar->v[ifield]->sweep[isweep]->ray[iray];
+ sweep->h.nrays += 1;
+
+
+ if (ray) {
+ /*
+ * ---- Beginning of MANDATORY HEADER BLOCK.
+ */
+ ray->h.ray_num = uf_ma[7];
+ if (little_endian()) swap2(&uf_ma[10], 8);
+ memcpy(radar->h.radar_name, &uf_ma[10], 8);
+ if (little_endian()) swap2(&uf_ma[10], 8/2);
+ memcpy(radar->h.name, &uf_ma[14], 8);
+ if (little_endian()) swap2(&uf_ma[14], 8/2);
+
/* All components of lat/lon are the same sign. If not, then
* what ever wrote the UF was in error. A simple RSL program
* can repair the damage, however, not here.
*/
- ray->h.lat = uf_ma[18] + uf_ma[19]/60.0 + uf_ma[20]/64.0/3600;
- ray->h.lon = uf_ma[21] + uf_ma[22]/60.0 + uf_ma[23]/64.0/3600;
- ray->h.alt = uf_ma[24];
- ray->h.year = uf_ma[25];
- if (ray->h.year < 1900) {
- ray->h.year += 1900;
- if (ray->h.year < 1980) ray->h.year += 100; /* Year >= 2000. */
- }
- ray->h.month = uf_ma[26];
- ray->h.day = uf_ma[27];
- ray->h.hour = uf_ma[28];
- ray->h.minute = uf_ma[29];
- ray->h.sec = uf_ma[30];
- ray->h.azimuth = uf_ma[32] / 64.0;
-
- /* If Local Use Header is present and contains azimuth, use that
- * azimuth for VR and SW. This is for WSR-88D, which runs separate
- * scans for DZ and VR/SW at the lower elevations, which means DZ
- * VR/SW and have different azimuths in the "same" ray.
- */
- len_lu = uf_ma[4] - uf_ma[3];
- if (len_lu == 2 && (ifield == VR_INDEX || ifield == SW_INDEX)) {
- if (strncmp((char *)uf_lu,"ZA",2) == 0 ||
- strncmp((char *)uf_lu,"AZ",2) == 0)
- ray->h.azimuth = uf_lu[1] / 64.0;
- }
- if (ray->h.azimuth < 0.) ray->h.azimuth += 360.; /* make it 0 to 360. */
- ray->h.elev = uf_ma[33] / 64.0;
- ray->h.elev_num = sweep->h.sweep_num;
- ray->h.fix_angle = sweep->h.elev = uf_ma[35] / 64.0;
- ray->h.azim_rate = uf_ma[36] / 64.0;
- ray->h.sweep_rate = ray->h.azim_rate * (60.0/360.0);
- missing_data = uf_ma[44];
-
- if (pulled_time_from_first_ray == 0) {
- radar->h.height = uf_ma[24];
- radar->h.latd = uf_ma[18];
- radar->h.latm = uf_ma[19];
- radar->h.lats = uf_ma[20] / 64.0;
- radar->h.lond = uf_ma[21];
- radar->h.lonm = uf_ma[22];
- radar->h.lons = uf_ma[23] / 64.0;
- radar->h.year = ray->h.year;
- radar->h.month = ray->h.month;
- radar->h.day = ray->h.day;
- radar->h.hour = ray->h.hour;
- radar->h.minute = ray->h.minute;
- radar->h.sec = ray->h.sec;
- strcpy(radar->h.radar_type, "uf");
- pulled_time_from_first_ray = 1;
- }
- /*
- * ---- End of MANDATORY HEADER BLOCK.
- */
-
- /* ---- Optional header used for MCTEX files. */
- /* If this is a MCTEX file, the first 4 words following the
- mandatory header contain the string 'MCTEX'. */
- memcpy(proj_name, (short *)(uf + uf_ma[2] - 1), 8);
- if (little_endian()) swap2(proj_name, 4);
-
-
- /* ---- Local Use Header (if present) was checked during Mandatory
- * Header processing above.
- */
-
- /* ---- Begining of FIELD HEADER. */
- uf_fh = uf+current_fh_index - 1;
- scale_factor = uf_fh[1];
- ray->h.range_bin1 = uf_fh[2] * 1000.0 + uf_fh[3];
- ray->h.gate_size = uf_fh[4];
-
- ray->h.nbins = uf_fh[5];
- ray->h.pulse_width = uf_fh[6]/(RSL_SPEED_OF_LIGHT/1.0e6);
-
- if (strncmp((char *)proj_name, "MCTEX", 5) == 0) /* MCTEX? */
- {
- /* The beamwidth values are not correct in Mctex UF files. */
- ray->h.beam_width = 1.0;
- sweep->h.beam_width = ray->h.beam_width;
- sweep->h.horz_half_bw = ray->h.beam_width/2.0;
- sweep->h.vert_half_bw = ray->h.beam_width/2.0;
- }
- else /* Not MCTEX */
- {
- ray->h.beam_width = uf_fh[7] / 64.0;
- sweep->h.beam_width = uf_fh[7] / 64.0;
- sweep->h.horz_half_bw = uf_fh[7] / 128.0; /* DFF 4/4/95 */
- sweep->h.vert_half_bw = uf_fh[8] / 128.0; /* DFF 4/4/95 */
- }
- /* fprintf (stderr, "uf_fh[7] = %d, [8] = %d\n", (int)uf_fh[7], (int)uf_fh[8]); */
- if((int)uf_fh[7] == -32768) {
- ray->h.beam_width = 1;
- sweep->h.beam_width = 1;
- sweep->h.horz_half_bw = .5;
- sweep->h.vert_half_bw = .5;
- }
-
- ray->h.frequency = uf_fh[9] / 64.0;
- ray->h.wavelength = uf_fh[11] / 64.0 / 100.0; /* cm to m. */
- ray->h.pulse_count = uf_fh[12];
- if (ifield == DZ_INDEX || ifield == ZT_INDEX) {
- radar->v[ifield]->h.calibr_const = uf_fh[16] / 100.0;
- /* uf value scaled by 100 */
- }
- else {
- radar->v[ifield]->h.calibr_const = 0.0;
- }
- if (uf_fh[17] == (short)UF_NO_DATA) x = 0;
- else x = uf_fh[17] / 1000000.0; /* PRT in seconds. */
- if (x != 0) {
- ray->h.prf = 1/x;
- ray->h.unam_rng = RSL_SPEED_OF_LIGHT / (2.0 * ray->h.prf * 1000.0);
- }
- else {
- ray->h.prf = 0.0;
- ray->h.unam_rng = 0.0;
- }
-
- if (VR_INDEX == ifield || VE_INDEX == ifield) {
- ray->h.nyq_vel = uf_fh[19] / scale_factor;
- }
-
- /* ---- End of FIELD HEADER. */
-
- ray->h.f = f;
- ray->h.invf = invf;
-
- /* ---- Begining of FIELD DATA. */
- uf_data = uf+uf_fh[0] - 1;
-
- len_data = ray->h.nbins; /* Known because of RSL_new_ray. */
- for (m=0; m<len_data; m++) {
- if (uf_data[m] == (short)UF_NO_DATA)
- ray->range[m] = invf(BADVAL); /* BADVAL */
- else {
- if(uf_data[m] == missing_data)
- ray->range[m] = invf(NOECHO); /* NOECHO */
- else
- ray->range[m] = invf((float)uf_data[m]/scale_factor);
- }
- }
- }
+ ray->h.lat = uf_ma[18] + uf_ma[19]/60.0 + uf_ma[20]/64.0/3600;
+ ray->h.lon = uf_ma[21] + uf_ma[22]/60.0 + uf_ma[23]/64.0/3600;
+ ray->h.alt = uf_ma[24];
+ ray->h.year = uf_ma[25];
+ if (ray->h.year < 1900) {
+ ray->h.year += 1900;
+ if (ray->h.year < 1980) ray->h.year += 100; /* Year >= 2000. */
+ }
+ ray->h.month = uf_ma[26];
+ ray->h.day = uf_ma[27];
+ ray->h.hour = uf_ma[28];
+ ray->h.minute = uf_ma[29];
+ ray->h.sec = uf_ma[30];
+ ray->h.azimuth = uf_ma[32] / 64.0;
+
+ /* If Local Use Header is present and contains azimuth, use that
+ * azimuth for VR and SW. This is for WSR-88D, which runs separate
+ * scans for DZ and VR/SW at the lower elevations, which means DZ
+ * VR/SW and have different azimuths in the "same" ray.
+ */
+ len_lu = uf_ma[4] - uf_ma[3];
+ if (len_lu == 2 && (ifield == VR_INDEX || ifield == SW_INDEX)) {
+ if (strncmp((char *)uf_lu,"ZA",2) == 0 ||
+ strncmp((char *)uf_lu,"AZ",2) == 0)
+ ray->h.azimuth = uf_lu[1] / 64.0;
+ }
+ if (ray->h.azimuth < 0.) ray->h.azimuth += 360.; /* make it 0 to 360. */
+ ray->h.elev = uf_ma[33] / 64.0;
+ ray->h.elev_num = sweep->h.sweep_num;
+ ray->h.fix_angle = sweep->h.elev = uf_ma[35] / 64.0;
+ ray->h.azim_rate = uf_ma[36] / 64.0;
+ ray->h.sweep_rate = ray->h.azim_rate * (60.0/360.0);
+ missing_data = uf_ma[44];
+
+ if (pulled_time_from_first_ray == 0) {
+ radar->h.height = uf_ma[24];
+ radar->h.latd = uf_ma[18];
+ radar->h.latm = uf_ma[19];
+ radar->h.lats = uf_ma[20] / 64.0;
+ radar->h.lond = uf_ma[21];
+ radar->h.lonm = uf_ma[22];
+ radar->h.lons = uf_ma[23] / 64.0;
+ radar->h.year = ray->h.year;
+ radar->h.month = ray->h.month;
+ radar->h.day = ray->h.day;
+ radar->h.hour = ray->h.hour;
+ radar->h.minute = ray->h.minute;
+ radar->h.sec = ray->h.sec;
+ strcpy(radar->h.radar_type, "uf");
+ pulled_time_from_first_ray = 1;
+ }
+ /*
+ * ---- End of MANDATORY HEADER BLOCK.
+ */
+
+ /* ---- Optional header used for MCTEX files. */
+ /* If this is a MCTEX file, the first 4 words following the
+ mandatory header contain the string 'MCTEX'. */
+ memcpy(proj_name, (short *)(uf + uf_ma[2] - 1), 8);
+ if (little_endian()) swap2(proj_name, 4);
+
+
+ /* ---- Local Use Header (if present) was checked during Mandatory
+ * Header processing above.
+ */
+
+ /* ---- Begining of FIELD HEADER. */
+ uf_fh = uf+current_fh_index - 1;
+ scale_factor = uf_fh[1];
+ ray->h.range_bin1 = uf_fh[2] * 1000.0 + uf_fh[3];
+ ray->h.gate_size = uf_fh[4];
+
+ ray->h.nbins = uf_fh[5];
+ ray->h.pulse_width = uf_fh[6]/(RSL_SPEED_OF_LIGHT/1.0e6);
+
+ if (strncmp((char *)proj_name, "MCTEX", 5) == 0) /* MCTEX? */
+ {
+ /* The beamwidth values are not correct in Mctex UF files. */
+ ray->h.beam_width = 1.0;
+ sweep->h.beam_width = ray->h.beam_width;
+ sweep->h.horz_half_bw = ray->h.beam_width/2.0;
+ sweep->h.vert_half_bw = ray->h.beam_width/2.0;
+ }
+ else /* Not MCTEX */
+ {
+ ray->h.beam_width = uf_fh[7] / 64.0;
+ sweep->h.beam_width = uf_fh[7] / 64.0;
+ sweep->h.horz_half_bw = uf_fh[7] / 128.0; /* DFF 4/4/95 */
+ sweep->h.vert_half_bw = uf_fh[8] / 128.0; /* DFF 4/4/95 */
+ }
+ /* fprintf (stderr, "uf_fh[7] = %d, [8] = %d\n", (int)uf_fh[7], (int)uf_fh[8]); */
+ if((int)uf_fh[7] == -32768) {
+ ray->h.beam_width = 1;
+ sweep->h.beam_width = 1;
+ sweep->h.horz_half_bw = .5;
+ sweep->h.vert_half_bw = .5;
+ }
+
+ ray->h.frequency = uf_fh[9] / 64.0;
+ ray->h.wavelength = uf_fh[11] / 64.0 / 100.0; /* cm to m. */
+ ray->h.pulse_count = uf_fh[12];
+ if (ifield == DZ_INDEX || ifield == ZT_INDEX) {
+ radar->v[ifield]->h.calibr_const = uf_fh[16] / 100.0;
+ /* uf value scaled by 100 */
+ }
+ else {
+ radar->v[ifield]->h.calibr_const = 0.0;
+ }
+ if (uf_fh[17] == (short)UF_NO_DATA) x = 0;
+ else x = uf_fh[17] / 1000000.0; /* PRT in seconds. */
+ if (x != 0) {
+ ray->h.prf = 1/x;
+ ray->h.unam_rng = RSL_SPEED_OF_LIGHT / (2.0 * ray->h.prf * 1000.0);
+ }
+ else {
+ ray->h.prf = 0.0;
+ ray->h.unam_rng = 0.0;
+ }
+
+ if (VR_INDEX == ifield || VE_INDEX == ifield) {
+ ray->h.nyq_vel = uf_fh[19] / scale_factor;
+ }
+
+ /* ---- End of FIELD HEADER. */
+
+ ray->h.f = f;
+ ray->h.invf = invf;
+
+ /* ---- Begining of FIELD DATA. */
+ uf_data = uf+uf_fh[0] - 1;
+
+ len_data = ray->h.nbins; /* Known because of RSL_new_ray. */
+ for (m=0; m<len_data; m++) {
+ if (uf_data[m] == (short)UF_NO_DATA)
+ ray->range[m] = invf(BADVAL); /* BADVAL */
+ else {
+ if(uf_data[m] == missing_data)
+ ray->range[m] = invf(NOECHO); /* NOECHO */
+ else
+ ray->range[m] = invf((float)uf_data[m]/scale_factor);
+ }
+ }
+ }
}
return UF_MORE;
}
addr_end = uf + sizeof(UF_buffer)/sizeof(short);
while (uf < addr_end)
- swap_2_bytes(uf++);
+ swap_2_bytes(uf++);
}
enum UF_type {NOT_UF, TRUE_UF, TWO_BYTE_UF, FOUR_BYTE_UF};
Radar *RSL_uf_to_radar_fp(FILE *fp)
{
union {
- char buf[6];
- short sword;
- int word;
+ char buf[6];
+ short sword;
+ int word;
} magic;
Radar *radar;
int nbytes;
radar = NULL;
- setvbuf(fp,NULL,_IOFBF,(size_t)NEW_BUFSIZ); /* Faster i/o? */
+ /* setvbuf(fp,NULL,_IOFBF,(size_t)NEW_BUFSIZ); * Faster i/o? */
if (fread(magic.buf, sizeof(char), 6, fp) <= 0) return NULL;
/*
* Check for fortran record length delimeters, NCAR kludge.
switch (uf_type) {
case FOUR_BYTE_UF:
- if (radar_verbose_flag) fprintf(stderr,"UF file with 4 byte FORTRAN record delimeters.\n");
- /* Handle first record specially, since we needed magic information. */
- nbytes = magic.word;
- if (little_endian()) swap_4_bytes(&nbytes);
- memcpy(uf, &magic.buf[4], 2);
- (void)fread(&uf[1], sizeof(char), nbytes-2, fp);
- if (little_endian()) swap_uf_buffer(uf);
- (void)fread(&nbytes, sizeof(int), 1, fp);
- if (uf_into_radar(uf, &radar) == UF_DONE) break;
- /* Now the rest of the file. */
- while(fread(&nbytes, sizeof(int), 1, fp) > 0) {
- if (little_endian()) swap_4_bytes(&nbytes);
-
- (void)fread(uf, sizeof(char), nbytes, fp);
- if (little_endian()) swap_uf_buffer(uf);
-
- (void)fread(&nbytes, sizeof(int), 1, fp);
-
- if (uf_into_radar(uf, &radar) == UF_DONE) break;
- }
- break;
+ if (radar_verbose_flag) fprintf(stderr,"UF file with 4 byte FORTRAN record delimeters.\n");
+ /* Handle first record specially, since we needed magic information. */
+ nbytes = magic.word;
+ if (little_endian()) swap_4_bytes(&nbytes);
+ memcpy(uf, &magic.buf[4], 2);
+ (void)fread(&uf[1], sizeof(char), nbytes-2, fp);
+ if (little_endian()) swap_uf_buffer(uf);
+ (void)fread(&nbytes, sizeof(int), 1, fp);
+ if (uf_into_radar(uf, &radar) == UF_DONE) break;
+ /* Now the rest of the file. */
+ while(fread(&nbytes, sizeof(int), 1, fp) > 0) {
+ if (little_endian()) swap_4_bytes(&nbytes);
+
+ (void)fread(uf, sizeof(char), nbytes, fp);
+ if (little_endian()) swap_uf_buffer(uf);
+
+ (void)fread(&nbytes, sizeof(int), 1, fp);
+
+ if (uf_into_radar(uf, &radar) == UF_DONE) break;
+ }
+ break;
case TWO_BYTE_UF:
- if (radar_verbose_flag) fprintf(stderr,"UF file with 2 byte FORTRAN record delimeters.\n");
- /* Handle first record specially, since we needed magic information. */
- sbytes = magic.sword;
- if (little_endian()) swap_2_bytes(&sbytes);
- memcpy(uf, &magic.buf[2], 4);
- (void)fread(&uf[2], sizeof(char), sbytes-4, fp);
- if (little_endian()) swap_uf_buffer(uf);
- (void)fread(&sbytes, sizeof(short), 1, fp);
- uf_into_radar(uf, &radar);
- /* Now the rest of the file. */
- while(fread(&sbytes, sizeof(short), 1, fp) > 0) {
- if (little_endian()) swap_2_bytes(&sbytes);
-
- (void)fread(uf, sizeof(char), sbytes, fp);
- if (little_endian()) swap_uf_buffer(uf);
-
- (void)fread(&sbytes, sizeof(short), 1, fp);
-
- if (uf_into_radar(uf, &radar) == UF_DONE) break;
- }
- break;
+ if (radar_verbose_flag) fprintf(stderr,"UF file with 2 byte FORTRAN record delimeters.\n");
+ /* Handle first record specially, since we needed magic information. */
+ sbytes = magic.sword;
+ if (little_endian()) swap_2_bytes(&sbytes);
+ memcpy(uf, &magic.buf[2], 4);
+ (void)fread(&uf[2], sizeof(char), sbytes-4, fp);
+ if (little_endian()) swap_uf_buffer(uf);
+ (void)fread(&sbytes, sizeof(short), 1, fp);
+ uf_into_radar(uf, &radar);
+ /* Now the rest of the file. */
+ while(fread(&sbytes, sizeof(short), 1, fp) > 0) {
+ if (little_endian()) swap_2_bytes(&sbytes);
+
+ (void)fread(uf, sizeof(char), sbytes, fp);
+ if (little_endian()) swap_uf_buffer(uf);
+
+ (void)fread(&sbytes, sizeof(short), 1, fp);
+
+ if (uf_into_radar(uf, &radar) == UF_DONE) break;
+ }
+ break;
case TRUE_UF:
- if (radar_verbose_flag) fprintf(stderr,"UF file with no FORTRAN record delimeters. Good.\n");
- /* Handle first record specially, since we needed magic information. */
- memcpy(&sbytes, &magic.buf[2], 2); /* Record length is in word #2. */
- if (little_endian()) swap_2_bytes(&sbytes); /* # of 2 byte words. */
-
- memcpy(uf, &magic.buf[0], 6);
- (void)fread(&uf[3], sizeof(short), sbytes-3, fp);
- if (little_endian()) swap_uf_buffer(uf);
- uf_into_radar(uf, &radar);
- /* Now the rest of the file. */
- while(fread(uf, sizeof(short), 2, fp) > 0) {
- memcpy(&sbytes, &uf[1], 2); /* Record length is in word #2. */
- if (little_endian()) swap_2_bytes(&sbytes);
-
- (void)fread(&uf[2], sizeof(short), sbytes-2, fp); /* Have words 1,2. */
- if (little_endian()) swap_uf_buffer(uf);
-
- if (uf_into_radar(uf, &radar) == UF_DONE) break;
- }
- break;
-
+ if (radar_verbose_flag) fprintf(stderr,"UF file with no FORTRAN record delimeters. Good.\n");
+ /* Handle first record specially, since we needed magic information. */
+ memcpy(&sbytes, &magic.buf[2], 2); /* Record length is in word #2. */
+ if (little_endian()) swap_2_bytes(&sbytes); /* # of 2 byte words. */
+
+ memcpy(uf, &magic.buf[0], 6);
+ (void)fread(&uf[3], sizeof(short), sbytes-3, fp);
+ if (little_endian()) swap_uf_buffer(uf);
+ uf_into_radar(uf, &radar);
+ /* Now the rest of the file. */
+ while(fread(uf, sizeof(short), 2, fp) > 0) {
+ memcpy(&sbytes, &uf[1], 2); /* Record length is in word #2. */
+ if (little_endian()) swap_2_bytes(&sbytes);
+
+ (void)fread(&uf[2], sizeof(short), sbytes-2, fp); /* Have words 1,2. */
+ if (little_endian()) swap_uf_buffer(uf);
+
+ if (uf_into_radar(uf, &radar) == UF_DONE) break;
+ }
+ break;
+
case NOT_UF: return NULL; break;
}
radar = reset_nsweeps_in_all_volumes(radar);
radar = NULL;
if (infile == NULL) {
- int save_fd;
- save_fd = dup(0);
- fp = fdopen(save_fd, "r");
+ int save_fd;
+ save_fd = dup(0);
+ fp = fdopen(save_fd, "r");
} else if ((fp = fopen(infile, "r")) == NULL) {
- perror(infile);
- return radar;
+ perror(infile);
+ return radar;
}
fp = uncompress_pipe(fp); /* Transparently gunzip. */
radar = RSL_uf_to_radar_fp(fp);
rsl_pclose(fp);
-
+
return radar;
}
projects / ~andy / rsl / blobdiff