ifield = HC_INDEX;
f = HC_F;
invf = HC_INVF;
+ case NSIG_DTB_DBZ2:
+ ifield = CZ_INDEX;
+ f = CZ_F;
+ invf = CZ_INVF;
+ case NSIG_DTB_ZDRC2:
+ ifield = ZD_INDEX;
+ f = ZD_F;
+ invf = ZD_INVF;
break;
default:
fprintf(stderr,"Unknown field type: %d Skipping it.\n", data_type);
ray->h.wavelength = wave/100.0; /* meters */
ray->h.nyq_vel = max_vel; /* m/s */
if (elev == 0.) elev = sweep->h.elev;
- ray->h.elev = elev;
+ ray->h.elev = (nsig_from_bang(ray_p->h.end_elev)+nsig_from_bang(ray_p->h.beg_elev))/2.0;
/* Compute mean azimuth angle for ray. */
az1 = nsig_from_bang(ray_p->h.beg_azm);
az2 = nsig_from_bang(ray_p->h.end_azm);