**Trajectory Ensemble Receptor Models:**

Conditional Probability Function (CPF)

**CPF**

Calculate CPF value.

**Syntax**

cpf = CPF (conc,direction,speed,c,bin,discard,plot)

**Description**

cpf = CPF (conc,direction,speed,c,bin,discard,plot) returns the cpf vector of given [conc,direction,speed] data (N by 1 vectors). c stands for the criterion value of CPF model; bin stands for the wind direction numbers; discard stands for the lowest acceptable wind speed, for example if discard equals 5, the data with speed less than 5 will not be included in CPF analysis; and the figure of CPF result will be plotted if plot equals 1.

Concentration Field Analysis (CFA)

**CFA**

Calculate CFA value.

**Syntax**

cfa = CFA (conc,trajectory,grid,plot)

**Description**

cfa = CFA (conc,trajectory,grid,plot) returns a grid matrix (latitude by longitude) cfa of calculated CFA value. conc is a N by S (S is site number) matrix of pollutant measurements; trajectory is a N by S by T (T is trajectory time) matrix of the corresponding trajectory endpoints; grid is the size of the grid cell, the unit is degree; and a world map will be plotted to display the CFA result if plot equals 1.

Concentration Weighted Trajectory (CWT)

**CWT**

Calculate CWT value.

**Syntax**

cwt = CWT (conc,trajectory,grid,plot)

**Description**

cwt = CWT (conc,trajectory,grid,plot) returns a grid matrix (latitude by longitude) cwt of calculated CWT value. conc is a N by S (S is site number) matrix of pollutant measurements; trajectory is a N by S by T (T is trajectory time) matrix of the corresponding trajectory endpoints; grid is the size of the grid cell, the unit is degree; and a world map will be plotted to display the CWT result if plot equals 1.

Residence Time Weighted Concentration (RTWC)

**RTWC**

Calculate RTWC value.

**Syntax**

rtwc= RTWC (conc,trajectory,grid,plot,iteration)

**Description**

rtwc= RTWC (conc,trajectory,grid,plot,iteration) returns a grid matrix (latitude by longitude) rtwc of calculated RTWC value. conc is a N by S (S is site number) matrix of pollutant measurements; trajectory is a N by S by T (T is trajectory time) matrix of the corresponding trajectory endpoints; grid is the size of the grid cell, the unit is degree; iteration stands for the criterion value to stop iteration, for example if iteration equals 0.5, iteration will stop if the average difference between two successive iterations is below 0.5%; and a world map will be plotted to display the RTWC result if plot equals 1.

Potential Source Contribution Function (PSCF)

**PSCF**

Calculate PSCF value.

**Syntax**

pscf = PSCF (conc,trajectory,grid,plot,c)

**Description**

pscf = PSCF (conc,trajectory,grid,plot,c) returns a grid matrix (latitude by longitude) pscf of calculated PSCF value. conc is a N by S (S is site number) matrix of pollutant measurements; trajectory is a N by S by T (T is trajectory time) matrix of the corresponding trajectory endpoints; grid is the size of the grid cell, the unit is degree; c stands for the criterion value of PSCF model; and a world map will be plotted to display the PSCF result if plot equals 1. Multiple-site PSCF method will automatically be applied if S is larger than 1.

Simplified Quantitative Transport Bias Analysis (SQTBA)

**SQTBA**

Calculate SQTBA value.

**Syntax**

sqtba = SQTBA (conc,trajectory,grid,plot)

**Description**

sqtba = SQTBA (conc,trajectory,grid,plot) returns a grid matrix (latitude by longitude) sqtba of calculated SQTBA value. conc is a N by S (S is site number) matrix of pollutant measurements; trajectory is a N by S by T (T is trajectory time) matrix of the corresponding trajectory endpoints; grid is the size of the grid cell, the unit is degree; and a world map will be plotted to display the SQTBA result if plot equals 1.