vsemC function

C version of the VSEM model