pavankonchada123
Subscriber

I have modified according to your suggestion but no use. I am trying to check for the if condition as mentioned below in the entire domain and trying to inject accordingly

if(4.15e34*val1*saturatedvapordensity(C_T(cell,cthread))*saturatedvapordensity(C_T(cell,cthread)) >10e12)
 {
     PP_RHO(p) = 998;  
      PP_FLOW_RATE(p) = nucl*(PP_RHO(p)*M_PI/6*pow((2*criticalradius),3))/0.001;
      PP_INIT_DIAM(p) = 2*criticalradius;
 PP_DIAM(p) = 2*criticalradius;     
     
 PP_VEL(p)[0] = C_U(cell,cthread);
PP_VEL(p)[1] =C_V(cell,cthread);
PP_VEL(p)[2] = C_W(cell,cthread);//C_VOLUME(cell,c_thread)*C_YI(cell,cthread,0)
      //Message("Volume %f\n",k); 
       PP_T(p) = C_T(cell,cthread);
      
   PP_MASS(p) = PP_RHO(p)*M_PI/6*pow(PP_DIAM(p),3);
   PP_INIT_MASS(p) = PP_RHO(p)*M_PI/6*pow(PP_DIAM(p),3);
 }