Corresponding relation between RM_GetSpeciesName and RM_GetSpeciesConcentrations

[Yongqiang]

Hi Phreeqc developpers,

I wonder if the output array can correspond to each other between  RM_GetSpeciesName and RM_GetSpeciesConcentrations.
Does the position of species name (e.g., Na+) corresponds to the concentration in the same position in the output array from RM_GetSpeciesConcentrations?

Thanks

[dlparkhurst]

Yes. The order of the species in RM_GetSpeciesConcentrations and RM_SetSpeciesConcentrations is given by the list retrieved with calls to RM_GetSpeciesName.

Code: [Select]

status = RM_SetSpeciesSaveOn(id, 1)
 ncomps = RM_FindComponents(id)
 nspecies = RM_GetSpeciesCount(id)
 do i = 1, nspecies
   status = RM_GetSpeciesName(id, i, name)
   write(*,*) name
 enddo

For transport modeling with the same dispersion coefficient for all elements, it is only necessary to transport the total concentrations (RM_GetConcentrations). If you are modeling multicomponent diffusion where each species can diffuse with a different diffusion coefficient (usually subject to a zero-net-current constraint), then you need to get the concentrations of each individual species for the transport calculations.

[Yongqiang]

Thanks, David.

Just to confirm again: The order of the species CONCENTRATION in RM_GetSpeciesConcentrations is the same with order of species NAME in RM_GetSpeciesName()?

Thanks
Yongqiang

[dlparkhurst]

Yes.

[Yongqiang]

Thanks for the confirmation, David!!