Package: gplume
## Gaussian Plume Model for Atmospheric DispersionThis repository contains a Python implementation of the Gaussian Plume Model for atmospheric dispersion of contaminants (PM2.5), along with an inverse modeling approach to estimate emission rates from receptor measurements. The model calculates the contaminant concentration at receptor locations based on emission rates and wind conditions. It assumes 4 sources, which can be customized according to the code that generates an input box. This code works for multiple sources and multiple receptors.### Code Description#### Forward Atmospheric Dispersion Modeling (gpm.py)The main code filegpm.pyincludes the following components:Contaminant Parameters: Set the parameters related to the contaminant being modeled, such as gravitational acceleration, dynamic viscosity of air, density of the contaminant, diameter of particles, deposition velocity, and molar mass of the contaminant.Source and Receptor Data: Define the locations and characteristics of the emission source and receptors where deposition measurements are made. This includes the number of sources, x-y-z coordinates, labels, and emission rates.gplumeFunction: This function computes the contaminant concentration (in kg/m^3) at a given set of receptor locations using the standard Gaussian plume solution. It takes into account the source characteristics, receptor locations, and wind speed.forward_atmospheric_dispersionFunction: This function calculates and plots the ground-level contaminant concentration contours based on the Gaussian Plume Model. It takes the wind speed as an input and calls thegplumefunction to calculate the concentrations. The resulting contours are displayed using thematplotliblibrary.#### Inverse Modeling (inverse.py)Theinverse.pyfile introduces an inverse modeling approach to estimate emission rates from observed contaminant concentrations at receptor locations. The key components in this file are:ermakFunction: This function computes the contaminant concentration at receptor locations using the Ermak dispersion model. It takes into account the emission rates, wind speed, and other parameters.Objective Function: Theobjective_functioncalculates the difference between predicted and observed contaminant concentrations based on the Ermak model. It sets up an optimization problem to find the optimal emission rates that minimize this difference.Optimization: Theminimizefunction from thescipy.optimizemodule is used to find the optimal emission rates that best fit the observed receptor measurements.### Output and UsageFor forward modeling, the code generates plots showing the contours of ground-level contaminant concentration in mg/m^3. The maximum concentration value is displayed in the plot title. (Use “from gplume import gpm”)For inverse modeling, the code outputs the optimal emission rates that best fit the observed receptor measurements. (Use “from gplume import inverse”)### ContributingContributions to this project are welcome. If you find any issues or have suggestions for improvements, feel free to create a pull request or submit an issue on the GitHub repository.https://github.com/VaibhavVasdev/Gaussian-Plume_Model### ContactFor any inquiries or questions, please contact Vaibhav Vasdev atvaibhavvasdev63@gmail.com.