#include "udf.h"
#define INT_YI 0.03
#define UDF_FILENAME "udf_loc_velo"
real AVG_YI;


DEFINE_ADJUST(outlet_avg_yi, d)
{
FILE *fp;
real sum_yi_A=0.0;
real sum_yi=0.0;
int j=0;
face_t f;
Thread *thread_out=Lookup_Thread(d,4);
  begin_f_loop(f,thread_out)
{
  sum_yi_A+=F_YI(f,thread_out,0);
  j++;      
}
end_f_loop(f,thread_out)

  AVG_YI=sum_yi_A/j;
fp=fopen(UDF_FILENAME,"w");

fprintf(fp, "%f",sum_yi_A); /* write out kount to data file */
fprintf(fp, "%f",AVG_YI);
fprintf(fp, "%d",j);
fclose(fp);

printf("%f\n",sum_yi_A);/* write out kount to data file */
printf("%f\n",AVG_YI);
printf("%d\n",j);
}


DEFINE_PROFILE(inlet,thread_in,index)
{
    real inlet_yi=0.0;
    face_t inlet_face;
  
    real time=RP_Get_Real("flow-time");   
    begin_f_loop(inlet_face,thread_in)
    {
          if(time==0)
        {
         inlet_yi=INT_YI;
          F_PROFILE(inlet_face,thread_in,index)=inlet_yi;
        }
        
          else
        {
          inlet_yi=AVG_YI+0.2*0.5;
          F_PROFILE(inlet_face,thread_in,index)=inlet_yi;
        }
    }
        end_f_loop(inlet_face,thread_in)   

printf("%f\n",inlet_yi);                     
}