def METHOD_NAME(query):
    lats = np.arange(query.min_lat, query.max_lat + query.res, query.res)
    lons = np.arange(query.min_lon, query.max_lon + query.res, query.res)
    num_lat = lats.size
    num_lon = lons.size
    intensity = np.empty((num_lat, num_lon))
    inclination = np.empty((num_lat, num_lon))
    declination = np.empty((num_lat, num_lon))
    for i, lat in enumerate(lats):
        for j, lon in enumerate(lons):
            mag = igrf12.igrf(date, glat=lat, glon=lon, alt_km=0., isv=0, itype=1)
            intensity[i][j] = mag.total / 1e5
            inclination[i][j] = mag.incl
            declination[i][j] = mag.decl
    return Result(query=query, lats=lats, lons=lons,
        declination=declination, inclination=inclination, intensity=intensity)