February 21, 2023 at 9:55 am
Chris_g
Subscriber

I am testing on a MWE so this is the whole code:

 #include "udf.h"
 #include "surf.h"
 #include "dpm_types.h"
 #include "dpm.h"
 #include "dpm_mem.h"
 #include "mem.h"
 #include "materials.h"
 #include "threads.h"
 
 
DEFINE_DPM_DRAG (test_drag, Re, tp)
{
    cell_t c;
    Thread *t0;
    real uvel;
    
    t0 = TP_CELL_THREAD(tp);
    c = TP_CELL(tp);
    
    uvel = C_U(c,t0);
    
    real drag_force;
    drag_force = 0.1;
         
    return (drag_force);
}