Const. momentum surce

TUU

I want to add an additional source of momentum. This source is a function of the position only (S=S(x,y,z)). I have S tabulated in a large ASCII table containing x,y,z,Sx,Sy,Sz. The way I wanted to use the table was by use of the user subroutine, sormom.f. This subroutine will then have to load the table-file and find the right values for Sx(x,y,z), Sy(x,y,z) and Sz(x,y,z). This does work but is very time consuming caused by the size of the file that has to be loaded for every single point in the domain during every iteration.

Question: Is it possible to load this table just once and then save it in the memory as it does not very with any values solved for?

Regards TUU