We extract part of function `calculate_radii_rossby` into new function
`node_to_prop`. So we now separate reading data from the shapefiles
and computing the average radius and Rossby number. Thus, we prepare
for rewriting the algorithm in the main part of the script. The
quantity of reading is still the same, for now.