Added a function to calculate the average u for a wurtzite structure. It
works for a simple case.