The current surface routines (RHOU, RHOQ, etc.) are created by calling bdref.f. All the code is hidden there. It would be much easier to create the surface arrays in numpy and would remove all fortran dependencies in the pre-processing. It would also be testable (the current code isn't easily testable).
Feature request
The current surface routines (RHOU, RHOQ, etc.) are created by calling bdref.f. All the code is hidden there. It would be much easier to create the surface arrays in numpy and would remove all fortran dependencies in the pre-processing. It would also be testable (the current code isn't easily testable).
Priority
Low