[joha]

Hi,

I'm trying to get the following UDF to work:

DEFINE_SOURCE(z_cap_force, cell, thread, dS, eqn) {

    real source=0.0;
    //skip source on first iteration (no gradient available)
    if (NULL != THREAD_STORAGE(thread,SV_VOF_G)){
        source = 150*C_VOF_G(cell, thread)[2];

        printf("Fk: %f\n", Fk);
        printf("C_VOF_G[0]: %f\n", C_VOF_G(cell, thread)[0]);
        printf("C_VOF_G[1]: %f\n", C_VOF_G(cell, thread)[1]);
        printf("C_VOF_G[2]: %f\n", C_VOF_G(cell, thread)[2]);
    }

    dS[eqn] = 0.0;
    return source;
}

The line "NULL != THREAD_STORAGE(thread,SV_VOF_G)" is "stolen" from one of the forum threads here If I got it right this macro checks if the volume fraction gradients have been stored. I also read that I have to set "Keep temporary solver memory from being freed?" to yes what I did.

Nevertheless "(NULL != THREAD_STORAGE(thread,SV_VOF_G))" is never true.

Anyone got an idea what went wrong?