EN | 中文
TraPSA
Trajectory-based Potential Sourece Apportionment Project
.

Trajectory Ensemble Receptor Models:
Conditional Probability Function (CPF)
CPF
Calculate CPF value.

Syntax

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.

© 2016 by Chuanlong Zhou. The Center for Air Resource Engineering & Science (CARES), Clarkson University.