`pkGetSources`

`pkGetSources`

-- Gets the source terms for the continuity equations and the relative
contribution of each processThe subroutine `pkGetSources`

evaluates the collisional source
and loss terms of the continuity equations for each species. The relative
contribution of each reaction can also be evaluated.

pkGetSources( LDensity, SrcC[, SrcP][, Ratios][, locus][, STAT] )

Name | Description | Data type and attributes |
---|---|---|

{ }LDensity(NnTV+NPhot) | Local density of species with variable concentration | REAL(double), INTENT(IN) |

`SrcC(NnTV)` | Collisional gain term of continuity equations | REAL(double), INTENT(OUT) |

`SrcP(NnTV)` | Collisional loss term of continuity equations divided by
`LDensity` | REAL(double), OPTIONAL, INTENT(OUT) |

`Ratios(NKgas,NnTV)` | Percent contribution of each reaction to the collisional source and loss terms of each species | REAL(double), OPTIONAL, INTENT(OUT) |

`locus` | Identifies where the reactions take place | CHARACTER(10), OPTIONAL, INTENT(IN) |

`STAT` | Execution status indicator | INTEGER, OPTIONAL, INTENT(OUT) |

The gain term is in L^{-3}s^{-1} units where
L^{-3} is the unit used to indicate the gas density in the data
file; the loss term is in s^{-1} units.

`STAT`

= `0`: no errors;
`STAT`

= `n`: error number.

If `locus`

is not present or is equal to `'gas'`, the
reactions are in the gas. Otherwise the reactions are on a surface and the value
of `locus`

identifies the surface. This allows having different
reactions on different surfaces.

The dimensions of array ` LDensity` can be

In the previous version of *PLASMAKIN* the photon density was set by an
independent routine `pkSetPhoton`

. With the increase of the size of
` LDensity` to

**Example E-8. pkGetSources: Computing the source
terms in a 1D discharge, taking into account gas temperature gradients.**

INTEGER :: i, Nxp REAL(double) :: Tgas(Nxp) REAL(double), ALLOCATABLE :: Dens(:,:), SrcC(:,:), SrcP(:,:) ... ALLOCATE( Dens(NnTV,Nxp), SrcC(NnTV,Nxp), SrcP(NnTV,Nxp) ) ... DO i = 1, Nxp CALL pkSetValue( GasTemp=Tgas(i) ) CALL pkGetSources( Dens(:,i), SrcC(:,i), SrcP(:,i) ) END DO