Skip to content

Commit

Permalink
Factorize slicing in save_gac
Browse files Browse the repository at this point in the history
  • Loading branch information
@sfinkens
sfinkens committed Oct 14, 2019
1 parent 4012701 commit 2345910
Showing 1 changed file with 12 additions and 73 deletions.
85 changes: 12 additions & 73 deletions pygac/gac_io.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -73,6 +73,14 @@
MISSING_DATA_LATLON = -999999


def slice_channel_inplace(arr, start_line, end_line, first_valid_lat,
last_valid_lat):
arr[:], _, _ = slice_channel(arr, start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)


def save_gac(satellite_name,
xutcs,
lats, lons,
Expand Down Expand Up @@ -111,82 +119,13 @@ def save_gac(satellite_name,
miss_lines=miss_lines,
midnight_scanline=midnight_scanline)

ref1, _, _ = slice_channel(ref1,
start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)
ref2, _, _ = slice_channel(ref2,
start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)
ref3, _, _ = slice_channel(ref3,
start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)
bt3, _, _ = slice_channel(bt3,
start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)
bt4, _, _ = slice_channel(bt4,
start_line=start_line,
for arr in (ref1, ref2, ref3, bt3, bt4, bt5, sun_zen, sun_azi, sat_zen,
sat_azi, rel_azi, lons, lats, qual_flags, xutcs):
slice_channel_inplace(arr, start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)
bt5, _, _ = slice_channel(bt5,
start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)
sun_zen, _, _ = slice_channel(sun_zen,
start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)
sun_azi, _, _ = slice_channel(sun_azi,
start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)
sat_zen, _, _ = slice_channel(sat_zen,
start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)
sat_azi, _, _ = slice_channel(sat_azi,
start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)
rel_azi, _, _ = slice_channel(rel_azi,
start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)
lons, _, _ = slice_channel(lons,
start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)
lats, _, _ = slice_channel(lats,
start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)
qual_flags, _, _ = slice_channel(qual_flags,
start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)
xutcs, _, _ = slice_channel(xutcs,
start_line=start_line,
end_line=end_line,
first_valid_lat=first_valid_lat,
last_valid_lat=last_valid_lat)


total_number_of_scan_lines = lats.shape[0]

# Reading time from the body of the gac file
Expand Down

0 comments on commit 2345910

Please sign in to comment.