I was hoping to speed up the code I mentioned last month implementing a line integral convolution (LIC) flow visualization, but I haven’t had time to improve it (and probably won’t for awhile). I would like to eventually write this as a DLM, but until then here’s my pure IDL draft. It contains a main-level program that is an example of using the main routine MG_LIC.

Global winds

Let me know if it helps or you make any improvements.