Conceptually, this is just adding additional columns to an index entry, with each column corresponding to an environmental variable value at the geographic location coincident with the geographic coordinate of the entry.

The algorithm is something like:

for each entry:
  coord = entry.coordinate
  for each env_var:
    v = env_var.value(coord)
    entry[varname] = v