I am working on a Fortran code that uses a uniformly spaced grid in two directions (x,y). Which works fine, but when I need to study a certain problem with good resolution, I need to increase the number of points in one of the directions. This works fine, but it requires a lot of computational resources.

I want to implement a non-uniform grid in one of those directions (the one I need, which is the y direction), in particular, a Gaussian spacing around a mean value. In other words, I need a bunch of points around a certain y value, then points with an increasingly larger space as we go further from this point (for both lower and higher values of y).

I have looked around but can’t seem to find anything useful. Maybe I am using the wrong keywords. I found many things for Python or Matlab that use the already implemented logspace or other functions. I also found a logspace in Fortran (I did not find the source/algorithm), but I would prefer to code my own. If it has already been solved, please let me know.