 | F2PY-generated modules treat storage schemes transparently |
| If input array has C storage, a copy is taken, calculated with, and returned as output |
| F2PY needs to know whether arguments are input, output or both |
| To monitor (hidden) array copying, turn on the flag
f2py ... -DF2PY_REPORT_ON_ARRAY_COPY=1 |
| In-place operations on NumPy arrays are possible in Fortran, but the default is to work on a copy, that is why our gridloop1 function does not work |
