How to access the porosity set by SetPorosity in PhreeqcRM using GET_POR

[conghaoyi]

Dear all,

Is there a way to access the porosity set by the SetPorosity function in PhreeqcRM using the GET_POR BASIC function in Phreeqc? In my simulation, the volume of the system is changing as expected using the SetPorosity function. However, when trying to access the porosity using the GET_POR function, the output porosity is constant at 0.1.

The simulation does not contain a TRANSPORT block. It is only a kinetic model.

The workflow of the simulation is shown in the pseudocode below:

Code: [Select]

For steps in range(nsteps):
  SetPorosity(por) # por changes at each time step
  SetConcentrations(c)
  RunCells()
  time = time + time_step
  SetTime(time)

However, when printing the porosity using the code below, the porosity from GET_POR stays 0.1

Code: [Select]

USER_PUNCH
-headings GET_PORO
-start
10 PUNCH GET_POR(1)
-end


[dlparkhurst]

I think you need to be doing transport calculations (TRANSPORT) for that function to work.

An alternative may be to use PUT to save the porosity value and GET when you want to retrieve it.