[Iman Sengupta]

I want to simulate flow of a non Newtonian power law fluid. I got this code from a thread in cfd-online.com. But This code was showing some error and I modified it. But still the code is showing parse error and not running. Please help me out of it.

#include "udf.h"

DEFINE_PROPERTY(user_power_law,c,t)

{

real visc_dyn_min=0.001;
real visc_dyn_max=100;

real index_n=0.6;

real k=1.7;

real mu_app;

real mu_temp;

real N_ITER=0;

mu_temp=k*pow(C_STRAIN_RATE_MAG(c,t),index_n-1);

if (N_ITER<1)

{

mu_app=k;

}

else if((N_ITER>=1)&&(mu_temp>visc_dyn_min)&&(mu_temp<v isc_dyn_max))
{
mu_app=mu_temp;
}
else if ((N_ITER>=1)&&(mu_temp>=visc_dyn_max))
{

mu_app=visc_dyn_max;
}
else if ((N_ITER>=1)&&(mu_temp<=visc_dyn_min))

{
mu_app=visc_dyn_min;
}

return mu_app;

}