In procedure fxhyp, instead of computing twice the integral of F, store it the first time: ffdx becomes an array and we do not need xxpr any longer. The storage is the same, there is less computation. In procedure grid_noro, instead of storing the non-smoothed orography in a temporary array zmea0, compute zphi earlier.