     date.month_name().lower()[:3],
            str(date.year)[-2:],
            np.round(date.day / 7) + 1,
        )
        step = self.run_time
        steps = np.arange(0, step // 24, np.sign(step))
        fnames = set([gdas_file(date + pd.Timedelta(s, unit=&#34;d&#34;)) for s in steps])
        return list(fnames)

    def get_coordinates(self, stations, height):
        locations = []
        for station in stations:
            lat, lon = self.locations[station]
            locations.append([lat, l