Writing the calculated value in DPM_SCALAR_UPDATE to file
 

Hello,

I have a multiphase (multi-species) gas-liquid eulerian flow and I want to use DPM for using massless-particles as a separate phase and record the concentrations that particles pass through (if we consider particles as microorganisms) of several species(substrate and O2) in the Eulerian liquid phase in their path.

I am trying to calculate the concentrations of 2 species in the cells where a particle exists (through macro DEFINE_DPM_SCALAR_UPDATE) and store them in a variable (P_USER_REAL(p,0) = (C_YI(c,pt[0],0)/180)*998.2; Cmole=P_USER_REAL(p,0);) and then use (DEFINE_EXECUTE_AT_END) to write these values in each time step into a file with:

fyleS = fopen("FileS.out","a");

fprintf( fyleS,"%e \n",Cmole)




for each particle, there should be two concentrations. I have 100 particles, but in each time step it only writes 3 values (which is still 3 values when I have only one particle)

Can anyone help me with this?