#include"udf.h"

#define k 2.0*M_PI/0.95

#define w 2.0*M_PI*1

DEFINE_GRID_MOTION(vel,domain,dt,time,dtime)

{

Thread *tf=DT_THREAD(dt);

face_t f;

int n;

Node *v;

double theta;

int i=0;

SET_DEFORMING_THREAD_FLAG(THREAD_T0(tf));

begin_f_loop(f,tf)

{

f_node_loop(f,tf,n)

{

v=F_NODE(f,tf,n);

if(NODE_POS_NEED_UPDATE(v))

{

NODE_POS_UPDATED(v);

NODE_Y(v)=NODE_Y(v)+(0.02-0.08*NODE_X(v)+0.16*NODE_X(v)*NODE_X(v))*(sin(k*NODE_X(v)-w*time)-sin(k*NODE_X(v)-w*(time-dtime)));

}

}

}

end_f_loop(f,tf);

}