@@ -569,8 +569,7 @@ def _parse_radial_fmt(self):
569
569
az [az > 360 ] -= 360
570
570
azi = np .deg2rad (az )
571
571
# self.azi = np.deg2rad(azi)
572
- dist = np .arange (start_range // reso + 1 , end_range // reso + 2 , 1 ) * reso
573
- dist = dist [: nbins ]
572
+ dist = np .arange (start_range // reso + 1 , end_range // reso + 1 , 1 ) * reso
574
573
lon , lat = get_coordinate (
575
574
dist , azi , self .params ["elevation" ], self .stationlon , self .stationlat
576
575
)
@@ -677,8 +676,7 @@ def _parse_cappi_fmt(self):
677
676
az += azi0 [0 ]
678
677
az [az > 360 ] -= 360
679
678
azi = np .deg2rad (az )
680
- dist = np .arange (start_range // reso + 1 , end_range // reso + 2 , 1 ) * reso
681
- dist = dist [: nbins ]
679
+ dist = np .arange (start_range + reso , end_range + reso , reso )
682
680
raw = np .vstack (data ).astype (int )
683
681
raw = np .ma .masked_less (raw , 5 )
684
682
data = (raw - offset ) / scale
@@ -703,8 +701,6 @@ def _parse_cappi_fmt(self):
703
701
)
704
702
ds ["longitude" ] = (["azimuth" , "distance" ], lon )
705
703
ds ["latitude" ] = (["azimuth" , "distance" ], lat )
706
- if len (height ) == 1 :
707
- ds = ds .squeeze ("height" )
708
704
self ._dataset = ds
709
705
710
706
def _parse_vwp_fmt (self ):
0 commit comments