#include "udf.h"

DEFINE_PROFILE(unsteady_pressure, thread, position) 
{

  float t, pressure;
  float pi, period, pmag, two;
  face_t f;

  t = RP_Get_Real("flow-time");
  
  two = 2.0;
  pi = 3.1415926535897932384626433832795;
  period = 1.e-2;

  /*pmag = 1.e4;    /*slow case*/

   pmag = 1.e5;   /*base case*/

   /*pmag = 1.e6;    /*fast case*/

  pressure = pmag*sin(two*pi*t/period); 

  begin_f_loop(f, thread)
  {
   F_PROFILE(f, thread, position) = pressure;
  }
  end_f_loop(f, thread)

}