I recommend using shape here:
```python
latitude = np.ones(conductivity.shape) * lat
```

I still don't know if this is the best way to make an array with all the same values though.