[Pyrokrates]

Hey, simple question: I would like to calculate a field magUp and save it with UDM for face and cell.


When I use the following code, everything works perfectly expect the last part for the F_UDMI. The program crashes, when I try to use the "C1" commands. Can you explain this? What am I doing wrong?


Side question: How can I write the part with UgofU... and Up... with the NV_D, NV_DS etc. commands? I run lots of errors trying to do this.


Thanks in advance


Chris

#include "udf.h"
#include "math.h"

DEFINE_EXECUTE_ON_LOADING(Load_UDM_Up, libname)
{
    Set_User_Memory_Name(0,"Up");
    Message("UDM id 0 reserved for Up.\n");
}

DEFINE_ADJUST(ADJUST_Up,d)
{
    Thread *t;
    cell_t c;
    face_t f;

    real K = 5; // (unit: s)
    real Up_mag;

    thread_loop_c(t,d)
    {
        begin_c_loop(c,t)
        {
            real UgofU = C_U(c,t) * C_U_G(c,t)[0] + C_V(c,t) * C_U_G(c,t)[1] + C_W(c,t) * C_U_G(c,t)[2];
            real VgofV = C_U(c,t) * C_V_G(c,t)[0] + C_V(c,t) * C_V_G(c,t)[1] + C_W(c,t) * C_V_G(c,t)[2];
            real WgofW = C_U(c,t) * C_W_G(c,t)[0] + C_V(c,t) * C_W_G(c,t)[1] + C_W(c,t) * C_W_G(c,t)[2];

            real Up = C_U(c,t) - UgofU * K;
            real Vp = C_V(c,t) - VgofV * K;
            real Wp = C_W(c,t) - WgofW * K;

            Up_mag = sqrt(pow(Up,2) + pow(Vp,2) + pow(Wp,2));

            C_UDMI(c,t,1) = Up_mag;
        }
        end_c_loop(c,t)

        begin_f_loop(f,t)
        {
            if(BOUNDARY_FACE_THREAD_P(t)) // Face is part of boundary field
            {
                F_UDMI(f,t,1) = C_UDMI(F_C0(f,t),F_C0_THREAD(f,t),1);
            }
            else // Face is part of internal field
            {
                F_UDMI(f,t,1) = (C_UDMI(F_C0(f,t),F_C0_THREAD(f,t),1) + C_UDMI(F_C1(f,t),F_C1_THREAD(f,t),1)) / 2;
            }
        }
        end_f_loop(f,t)
    }
}

[lolno]

connectivity macro of C1 doesnt exist in boundary face thread of a wall