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.
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.
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.
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.
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.
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.
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.