 |
MPTRAC
|
|
MPTRAC
|
MPTRAC library declarations. More...
#include <ctype.h>#include <gsl/gsl_fft_complex.h>#include <gsl/gsl_math.h>#include <gsl/gsl_randist.h>#include <gsl/gsl_rng.h>#include <gsl/gsl_spline.h>#include <gsl/gsl_statistics.h>#include <math.h>#include <netcdf.h>#include <omp.h>#include <stdint.h>#include <stdio.h>#include <stdlib.h>#include <string.h>#include <time.h>#include <sys/time.h>Go to the source code of this file.
Data Structures | |
| struct | ctl_t |
| Control parameters. More... | |
| struct | atm_t |
| Air parcel data. More... | |
| struct | cache_t |
| Cache data structure. More... | |
| struct | clim_photo_t |
| Climatological data in the form of photolysis rates. More... | |
| struct | clim_ts_t |
| Climatological data in the form of time series. More... | |
| struct | clim_zm_t |
| Climatological data in the form of zonal means. More... | |
| struct | clim_t |
| Climatological data. More... | |
| struct | met_t |
| Meteo data structure. More... | |
Macros | |
| #define | AVO 6.02214076e23 |
| Avogadro constant [1/mol]. More... | |
| #define | CPD 1003.5 |
| Specific heat of dry air at constant pressure [J/(kg K)]. More... | |
| #define | EPS (MH2O / MA) |
| Ratio of the specific gas constant of dry air and water vapor [1]. More... | |
| #define | G0 9.80665 |
| Standard gravity [m/s^2]. More... | |
| #define | H0 7.0 |
| Scale height [km]. More... | |
| #define | LV 2501000. |
| Latent heat of vaporization of water [J/kg]. More... | |
| #define | KB 1.3806504e-23 |
| Boltzmann constant [kg m^2/(K s^2)]. More... | |
| #define | MA 28.9644 |
| Molar mass of dry air [g/mol]. More... | |
| #define | MH2O 18.01528 |
| Molar mass of water vapor [g/mol]. More... | |
| #define | MO3 48.00 |
| Molar mass of ozone [g/mol]. More... | |
| #define | P0 1013.25 |
| Standard pressure [hPa]. More... | |
| #define | RA (1e3 * RI / MA) |
| Specific gas constant of dry air [J/(kg K)]. More... | |
| #define | RE 6367.421 |
| Mean radius of Earth [km]. More... | |
| #define | RI 8.3144598 |
| Ideal gas constant [J/(mol K)]. More... | |
| #define | T0 273.15 |
| Standard temperature [K]. More... | |
| #define | LEN 5000 |
| Maximum length of ASCII data lines. More... | |
| #define | NP 10000000 |
| Maximum number of atmospheric data points. More... | |
| #define | NQ 15 |
| Maximum number of quantities per data point. More... | |
| #define | NCSI 1000000 |
| Maximum number of data points for CSI calculation. More... | |
| #define | EP 140 |
| Maximum number of pressure levels for meteo data. More... | |
| #define | EX 1202 |
| Maximum number of longitudes for meteo data. More... | |
| #define | EY 602 |
| Maximum number of latitudes for meteo data. More... | |
| #define | NENS 2000 |
| Maximum number of data points for ensemble analysis. More... | |
| #define | NOBS 10000000 |
| Maximum number of observation data points. More... | |
| #define | NTHREADS 512 |
| Maximum number of OpenMP threads. More... | |
| #define | CY 250 |
| Maximum number of latitudes for climatological data. More... | |
| #define | CO3 30 |
| Maximum number of total column ozone data for climatological data. More... | |
| #define | CP 70 |
| Maximum number of pressure levels for climatological data. More... | |
| #define | CSZA 50 |
| Maximum number of solar zenith angles for climatological data. More... | |
| #define | CT 12 |
| Maximum number of time steps for climatological data. More... | |
| #define | CTS 1000 |
| Maximum number of data points of climatological time series. More... | |
| #define | ALLOC(ptr, type, n) |
| Allocate memory for a pointer with error handling. More... | |
| #define | ARRAY_2D(ix, iy, ny) ((ix) * (ny) + (iy)) |
| Macro for computing the linear index of a 2D array element. More... | |
| #define | ARRAY_3D(ix, iy, ny, iz, nz) (((ix)*(ny) + (iy)) * (nz) + (iz)) |
| Compute the linear index of a 3D array element. More... | |
| #define | ARRHENIUS(a, b, t) ((a) * exp( -(b) / (t))) |
| Calculate the Arrhenius rate constant. More... | |
| #define | DEG2DX(dlon, lat) ((dlon) * M_PI * RE / 180. * cos((lat) / 180. * M_PI)) |
| Convert a longitude difference to a distance in the x-direction (east-west) at a specific latitude. More... | |
| #define | DEG2DY(dlat) ((dlat) * M_PI * RE / 180.) |
| Convert a latitude difference to a distance in the y-direction (north-south). More... | |
| #define | DP2DZ(dp, p) (- (dp) * H0 / (p)) |
| Convert a pressure difference to a height difference in the vertical direction. More... | |
| #define | DX2DEG(dx, lat) |
| Convert a distance in kilometers to degrees longitude at a given latitude. More... | |
| #define | DY2DEG(dy) ((dy) * 180. / (M_PI * RE)) |
| Convert a distance in kilometers to degrees latitude. More... | |
| #define | DZ2DP(dz, p) (-(dz) * (p) / H0) |
| Convert a change in altitude to a change in pressure. More... | |
| #define | DIST(a, b) sqrt(DIST2(a, b)) |
| Calculate the distance between two points in Cartesian coordinates. More... | |
| #define | DIST2(a, b) ((a[0]-b[0])*(a[0]-b[0])+(a[1]-b[1])*(a[1]-b[1])+(a[2]-b[2])*(a[2]-b[2])) |
| Calculate the squared Euclidean distance between two points in Cartesian coordinates. More... | |
| #define | DOTP(a, b) (a[0]*b[0]+a[1]*b[1]+a[2]*b[2]) |
| Calculate the dot product of two vectors. More... | |
| #define | FMOD(x, y) ((x) - (int) ((x) / (y)) * (y)) |
| Calculate the floating-point remainder of dividing x by y. More... | |
| #define | FREAD(ptr, type, size, in) |
| Read data from a file stream and store it in memory. More... | |
| #define | FWRITE(ptr, type, size, out) |
| Write data from memory to a file stream. More... | |
| #define | INTPOL_INIT double cw[4] = {0.0, 0.0, 0.0, 0.0}; int ci[3] = {0, 0, 0}; |
| Initialize arrays for interpolation. More... | |
| #define | INTPOL_2D(var, init) |
| Perform 2D interpolation for a meteorological variable. More... | |
| #define | INTPOL_3D(var, init) |
| Perform 3D interpolation for a meteorological variable. More... | |
| #define | INTPOL_SPACE_ALL(p, lon, lat) |
| Interpolate multiple meteorological variables in space. More... | |
| #define | INTPOL_TIME_ALL(time, p, lon, lat) |
| Interpolate multiple meteorological variables in time. More... | |
| #define | LAPSE(p1, t1, p2, t2) |
| Calculate lapse rate. More... | |
| #define | LIN(x0, y0, x1, y1, x) ((y0)+((y1)-(y0))/((x1)-(x0))*((x)-(x0))) |
| Linear interpolation. More... | |
| #define | MAX(a, b) (((a)>(b))?(a):(b)) |
| Macro to determine the maximum of two values. More... | |
| #define | MET_HEADER |
| Write header for meteorological data file. More... | |
| #define | MIN(a, b) (((a)<(b))?(a):(b)) |
| Macro to determine the minimum of two values. More... | |
| #define | MOLEC_DENS(p, t) (AVO * 1e-6 * ((p) * 100) / (RI * (t))) |
| Calculate the density of a gas molecule. More... | |
| #define | NC(cmd) |
| Execute a NetCDF command and check for errors. More... | |
| #define | NC_DEF_VAR(varname, type, ndims, dims, long_name, units) |
| Define a NetCDF variable with attributes. More... | |
| #define | NC_GET_DOUBLE(varname, ptr, force) |
| Retrieve a double-precision variable from a NetCDF file. More... | |
| #define | NC_INQ_DIM(dimname, ptr, min, max) |
| Inquire the length of a dimension in a NetCDF file. More... | |
| #define | NC_PUT_DOUBLE(varname, ptr, hyperslab) |
| Write double precision data to a NetCDF variable. More... | |
| #define | NC_PUT_FLOAT(varname, ptr, hyperslab) |
| Write a float array to a NetCDF file. More... | |
| #define | NC_PUT_INT(varname, ptr, hyperslab) |
| Write integer data to a NetCDF variable. More... | |
| #define | NC_PUT_ATT(varname, attname, text) |
| Add a text attribute to a NetCDF variable. More... | |
| #define | NC_PUT_ATT_GLOBAL(attname, text) NC(nc_put_att_text(ncid, NC_GLOBAL, attname, strnlen(text, LEN), text)); |
| Add a global text attribute to a NetCDF file. More... | |
| #define | NN(x0, y0, x1, y1, x) (fabs((x) - (x0)) <= fabs((x) - (x1)) ? (y0) : (y1)) |
| Perform nearest-neighbor interpolation. More... | |
| #define | NORM(a) sqrt(DOTP(a, a)) |
| Compute the norm (magnitude) of a vector. More... | |
| #define | PARTICLE_LOOP(ip0, ip1, check_dt, ...) |
| Loop over particle indices with OpenACC acceleration. More... | |
| #define | P(z) (P0 * exp(-(z) / H0)) |
| Compute pressure at given altitude. More... | |
| #define | PSAT(t) (6.112 * exp(17.62 * ((t) - T0) / (243.12 + (t) - T0))) |
| Compute saturation pressure over water. More... | |
| #define | PSICE(t) (6.112 * exp(22.46 * ((t) - T0) / (272.62 + (t) - T0))) |
| Compute saturation pressure over ice (WMO, 2018). More... | |
| #define | PW(p, h2o) ((p) * MAX((h2o), 0.1e-6) / (1. + (1. - EPS) * MAX((h2o), 0.1e-6))) |
| Calculate partial water vapor pressure. More... | |
| #define | RH(p, t, h2o) (PW(p, h2o) / PSAT(t) * 100.) |
| Compute relative humidity over water. More... | |
| #define | RHICE(p, t, h2o) (PW(p, h2o) / PSICE(t) * 100.) |
| Compute relative humidity over ice. More... | |
| #define | RHO(p, t) (100. * (p) / (RA * (t))) |
| Compute density of air. More... | |
| #define | ROETH_PHOTOL(a, b, c, sza) ((c)*(sza) < M_PI/2. ? (a) * exp((b) * (1 - 1/cos((c) * (sza)))) : 0) |
| Roeth approximation formula for photolysis reactions. More... | |
| #define | SET_ATM(qnt, val) |
| Set atmospheric quantity value. More... | |
| #define | SET_QNT(qnt, name, longname, unit) |
| Set atmospheric quantity index. More... | |
| #define | SH(h2o) (EPS * MAX((h2o), 0.1e-6)) |
| Compute specific humidity from water vapor volume mixing ratio. More... | |
| #define | SQR(x) ((x)*(x)) |
| Compute the square of a value. More... | |
| #define | SWAP(x, y, type) do {type tmp = x; x = y; y = tmp;} while(0); |
| Swap two values. More... | |
| #define | TDEW(p, h2o) |
| Calculate dew point temperature. More... | |
| #define | TICE(p, h2o) |
| Calculate frost point temperature (WMO, 2018). More... | |
| #define | THETA(p, t) ((t) * pow(1000. / (p), 0.286)) |
| Compute potential temperature. More... | |
| #define | THETAVIRT(p, t, h2o) (TVIRT(THETA((p), (t)), MAX((h2o), 0.1e-6))) |
| Compute virtual potential temperature. More... | |
| #define | TOK(line, tok, format, var) |
| Get string tokens. More... | |
| #define | TVIRT(t, h2o) ((t) * (1. + (1. - EPS) * MAX((h2o), 0.1e-6))) |
| Compute virtual temperature. More... | |
| #define | Z(p) (H0 * log(P0 / (p))) |
| Convert pressure to altitude. More... | |
| #define | ZDIFF(lnp0, t0, h2o0, lnp1, t1, h2o1) |
| Calculate geopotential height difference. More... | |
| #define | ZETA(ps, p, t) |
| Calculate potential vorticity using the Zeta approximation. More... | |
| #define | LOGLEV 2 |
| Level of log messages (0=none, 1=basic, 2=detailed, 3=debug). More... | |
| #define | LOG(level, ...) |
| Print a log message with a specified logging level. More... | |
| #define | WARN(...) |
| Print a warning message with contextual information. More... | |
| #define | ERRMSG(...) |
| Print an error message with contextual information and terminate the program. More... | |
| #define | PRINT(format, var) |
| Print the value of a variable with contextual information. More... | |
| #define | NTIMER 100 |
| Maximum number of timers. More... | |
| #define | PRINT_TIMERS timer("END", "END", 1); |
| Print the current state of all timers. More... | |
| #define | SELECT_TIMER(id, group, color) |
| Select and start a timer with specific attributes. More... | |
| #define | START_TIMERS NVTX_PUSH("START", NVTX_CPU); |
| Starts a timer for tracking. More... | |
| #define | STOP_TIMERS NVTX_POP; |
| Stop the current timer. More... | |
| #define | NVTX_PUSH(range_title, range_color) {} |
| #define | NVTX_POP {} |
Functions | |
| void | thrustSortWrapper (double *__restrict__ c, int n, int *__restrict__ index) |
| Wrapper to Thrust sorting function. More... | |
| void | broadcast_large_data (void *data, size_t N) |
| Broadcasts large data across all processes in an MPI communicator. More... | |
| void | cart2geo (const double *x, double *z, double *lon, double *lat) |
| Converts Cartesian coordinates to geographic coordinates. More... | |
| double | clim_oh (const ctl_t *ctl, const clim_t *clim, const double t, const double lon, const double lat, const double p) |
| Calculates the hydroxyl radical (OH) concentration from climatology data, with an optional diurnal correction based on solar zenith angle. More... | |
| void | clim_oh_diurnal_correction (ctl_t *ctl, clim_t *clim) |
| Applies a diurnal correction to the hydroxyl radical (OH) concentration in climatology data. More... | |
| double | clim_photo (double rate[CP][CSZA][CO3], clim_photo_t *photo, double p, double sza, double o3c) |
| Calculates the photolysis rate for a given set of atmospheric conditions. More... | |
| double | clim_tropo (const clim_t *clim, const double t, const double lat) |
| Calculates the tropopause pressure based on climatological data. More... | |
| void | clim_tropo_init (clim_t *clim) |
| Initializes the tropopause data in the climatology structure. More... | |
| double | clim_ts (const clim_ts_t *ts, const double t) |
| Interpolates a time series of climatological variables. More... | |
| double | clim_zm (const clim_zm_t *zm, const double t, const double lat, const double p) |
| Interpolates monthly mean zonal mean climatological variables. More... | |
| void | compress_cms (ctl_t *ctl, char *varname, float *array, size_t nx, size_t ny, size_t np, int decompress, FILE *inout) |
| Compresses or decompresses a 3D array of floats using a custom multiscale compression algorithm. More... | |
| void | compress_pck (char *varname, float *array, size_t nxy, size_t nz, int decompress, FILE *inout) |
| Compresses or decompresses a 3D array of floats. More... | |
| void | compress_zfp (char *varname, float *array, int nx, int ny, int nz, int precision, double tolerance, int decompress, FILE *inout) |
| Compresses or decompresses a 3D array of floats using the ZFP library. More... | |
| void | compress_zstd (char *varname, float *array, size_t n, int decompress, FILE *inout) |
| Compresses or decompresses an array of floats using the Zstandard (ZSTD) library. More... | |
| void | day2doy (const int year, const int mon, const int day, int *doy) |
| Get day of year from date. More... | |
| void | doy2day (const int year, const int doy, int *mon, int *day) |
| Converts a given day of the year (DOY) to a date (month and day). More... | |
| void | fft_help (double *fcReal, double *fcImag, int n) |
| Computes the Fast Fourier Transform (FFT) of a complex sequence. More... | |
| void | geo2cart (const double z, const double lon, const double lat, double *x) |
| Converts geographic coordinates (longitude, latitude, altitude) to Cartesian coordinates. More... | |
| void | get_met (ctl_t *ctl, clim_t *clim, double t, met_t **met0, met_t **met1) |
| Retrieves meteorological data for the specified time. More... | |
| void | get_met_help (ctl_t *ctl, double t, int direct, char *metbase, double dt_met, char *filename) |
| Helper function to generate the filename for meteorological data. More... | |
| void | get_met_replace (char *orig, char *search, char *repl) |
| Replaces occurrences of a substring in a string with another substring. More... | |
| void | get_tropo (int met_tropo, ctl_t *ctl, clim_t *clim, met_t *met, double *lons, int nx, double *lats, int ny, double *pt, double *zt, double *tt, double *qt, double *o3t, double *ps, double *zs) |
| Calculate tropopause data. More... | |
| void | intpol_met_4d_coord (met_t *met0, float height0[EX][EY][EP], float array0[EX][EY][EP], met_t *met1, float height1[EX][EY][EP], float array1[EX][EY][EP], double ts, double height, double lon, double lat, double *var, int *ci, double *cw, int init) |
| Interpolates meteorological variables to a given position and time. More... | |
| void | intpol_met_space_3d (met_t *met, float array[EX][EY][EP], double p, double lon, double lat, double *var, int *ci, double *cw, int init) |
| Interpolates meteorological variables in 3D space. More... | |
| void | intpol_met_space_3d_ml (met_t *met, float array[EX][EY][EP], double p, double lon, double lat, double *var) |
| Interpolates a meteorological variable in 3D space (longitude, latitude, pressure). More... | |
| void | intpol_met_space_2d (met_t *met, float array[EX][EY], double lon, double lat, double *var, int *ci, double *cw, int init) |
| Interpolates meteorological variables in 2D space. More... | |
| void | intpol_met_time_3d (met_t *met0, float array0[EX][EY][EP], met_t *met1, float array1[EX][EY][EP], double ts, double p, double lon, double lat, double *var, int *ci, double *cw, int init) |
| Interpolates meteorological data in 3D space and time. More... | |
| void | intpol_met_time_3d_ml (met_t *met0, float array0[EX][EY][EP], met_t *met1, float array1[EX][EY][EP], double ts, double p, double lon, double lat, double *var) |
| Interpolates a meteorological variable in time and 3D space (longitude, latitude, pressure). More... | |
| void | intpol_met_time_2d (met_t *met0, float array0[EX][EY], met_t *met1, float array1[EX][EY], double ts, double lon, double lat, double *var, int *ci, double *cw, int init) |
| Interpolates meteorological data in 2D space and time. More... | |
| void | intpol_tropo_3d (double time0, float array0[EX][EY], double time1, float array1[EX][EY], double lons[EX], double lats[EY], int nlon, int nlat, double time, double lon, double lat, int method, double *var, double *sigma) |
| Interpolates tropopause data in 3D (latitude, longitude, and time). More... | |
| void | jsec2time (const double jsec, int *year, int *mon, int *day, int *hour, int *min, int *sec, double *remain) |
| Converts Julian seconds to calendar date and time components. More... | |
| double | kernel_weight (const double kz[EP], const double kw[EP], const int nk, const double p) |
| Calculates the kernel weight based on altitude and given kernel data. More... | |
| double | lapse_rate (const double t, const double h2o) |
| Calculates the moist adiabatic lapse rate in Kelvin per kilometer. More... | |
| void | level_definitions (ctl_t *ctl) |
| Defines pressure levels for meteorological data. More... | |
| int | locate_irr (const double *xx, const int n, const double x) |
| Locate the index of the interval containing a given value in a sorted array. More... | |
| int | locate_irr_float (const float *xx, const int n, const double x, const int ig) |
| Locate the index of the interval containing a given value in an irregularly spaced array. More... | |
| int | locate_irr_3d (float profiles[EX][EY][EP], int np, int ind_lon, int ind_lat, double x) |
| Locate the index of the interval containing a given value in a 3D irregular grid. More... | |
| int | locate_reg (const double *xx, const int n, const double x) |
| Locate the index of the interval containing a given value in a regular grid. More... | |
| void | locate_vert (float profiles[EX][EY][EP], int np, int lon_ap_ind, int lat_ap_ind, double alt_ap, int *ind) |
| Locate the four vertical indizes of a box for a given height value. More... | |
| void | module_advect (ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, double *dt) |
| Performs the advection of atmospheric particles using meteorological data. More... | |
| void | module_advect_init (ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm) |
| Initializes the advection module by setting up pressure fields. More... | |
| void | module_bound_cond (ctl_t *ctl, clim_t *clim, met_t *met0, met_t *met1, atm_t *atm, double *dt) |
| Apply boundary conditions to particles based on meteorological and climatological data. More... | |
| void | module_chemgrid (ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, double t) |
| Calculate grid data for chemistry modules. More... | |
| void | module_chem_init (ctl_t *ctl, clim_t *clim, met_t *met0, met_t *met1, atm_t *atm) |
| Initializes the chemistry modules by setting atmospheric composition. More... | |
| void | module_convection (ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, double *dt, double *rs) |
| Simulate convective processes for atmospheric particles. More... | |
| void | module_decay (ctl_t *ctl, clim_t *clim, atm_t *atm, double *dt) |
| Simulate exponential decay processes for atmospheric particles. More... | |
| void | module_diffusion_meso (ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, cache_t *cache, double *dt, double *rs) |
| Simulate mesoscale diffusion for atmospheric particles. More... | |
| void | module_diffusion_turb (ctl_t *ctl, clim_t *clim, atm_t *atm, double *dt, double *rs) |
| Simulate turbulent diffusion for atmospheric particles. More... | |
| void | module_dry_deposition (ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, double *dt) |
| Simulate dry deposition of atmospheric particles. More... | |
| void | module_h2o2_chem (ctl_t *ctl, clim_t *clim, met_t *met0, met_t *met1, atm_t *atm, double *dt) |
| Perform chemical reactions involving H2O2 within cloud particles. More... | |
| void | module_isosurf_init (ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, cache_t *cache) |
| Initialize the isosurface module based on atmospheric data. More... | |
| void | module_isosurf (ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, cache_t *cache, double *dt) |
| Apply the isosurface module to adjust atmospheric properties. More... | |
| void | module_kpp_chem (ctl_t *ctl, clim_t *clim, met_t *met0, met_t *met1, atm_t *atm, double *dt) |
| KPP chemistry module. More... | |
| void | module_meteo (ctl_t *ctl, clim_t *clim, met_t *met0, met_t *met1, atm_t *atm, double *dt) |
| Update atmospheric properties using meteorological data. More... | |
| void | module_mixing (ctl_t *ctl, clim_t *clim, atm_t *atm, double t) |
| Update atmospheric properties through interparcel mixing. More... | |
| void | module_mixing_help (ctl_t *ctl, clim_t *clim, atm_t *atm, const int *ixs, const int *iys, const int *izs, int qnt_idx) |
| Perform interparcel mixing for a specific quantity. More... | |
| void | module_oh_chem (ctl_t *ctl, clim_t *clim, met_t *met0, met_t *met1, atm_t *atm, double *dt) |
| Perform hydroxyl chemistry calculations for atmospheric particles. More... | |
| void | module_position (ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, double *dt) |
| Update the positions and pressure levels of atmospheric particles. More... | |
| void | module_rng_init (int ntask) |
| Initialize random number generators for parallel tasks. More... | |
| void | module_rng (ctl_t *ctl, double *rs, size_t n, int method) |
| Generate random numbers using various methods and distributions. More... | |
| void | module_sedi (ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, double *dt) |
| Simulate sedimentation of particles in the atmosphere. More... | |
| void | module_sort (ctl_t *ctl, met_t *met0, atm_t *atm) |
| Sort particles according to box index. More... | |
| void | module_sort_help (double *a, const int *p, const int np) |
| Reorder an array based on a given permutation. More... | |
| void | module_timesteps (ctl_t *ctl, met_t *met0, atm_t *atm, double *dt, double t) |
| Calculate time steps for air parcels based on specified conditions. More... | |
| void | module_timesteps_init (ctl_t *ctl, atm_t *atm) |
| Initialize start time and time interval for time-stepping. More... | |
| void | module_tracer_chem (ctl_t *ctl, clim_t *clim, met_t *met0, met_t *met1, atm_t *atm, double *dt) |
| Simulate chemical reactions involving long-lived atmospheric tracers. More... | |
| void | module_wet_deposition (ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, double *dt) |
| Perform wet deposition calculations for air parcels. More... | |
| double | nat_temperature (const double p, const double h2o, const double hno3) |
| Calculates the nitric acid trihydrate (NAT) temperature. More... | |
| int | read_atm (const char *filename, ctl_t *ctl, atm_t *atm) |
| Reads air parcel data from a specified file into the given atmospheric structure. More... | |
| int | read_atm_asc (const char *filename, ctl_t *ctl, atm_t *atm) |
| Reads air parcel data from an ASCII file and populates the given atmospheric structure. More... | |
| int | read_atm_bin (const char *filename, ctl_t *ctl, atm_t *atm) |
| Reads air parcel data from a binary file and populates the given atmospheric structure. More... | |
| int | read_atm_clams (const char *filename, ctl_t *ctl, atm_t *atm) |
| Reads air parcel data from a CLaMS netCDF file and populates the given atmospheric structure. More... | |
| int | read_atm_nc (const char *filename, ctl_t *ctl, atm_t *atm) |
| Reads air parcel data from a generic netCDF file and populates the given atmospheric structure. More... | |
| void | read_clim (ctl_t *ctl, clim_t *clim) |
| Reads various climatological data and populates the given climatology structure. More... | |
| void | read_clim_photo (const char *filename, clim_photo_t *photo) |
| Reads photolysis rates from a NetCDF file and populates the given photolysis structure. More... | |
| int | read_clim_ts (const char *filename, clim_ts_t *ts) |
| Reads a climatological time series from a file and populates the given time series structure. More... | |
| void | read_clim_zm (const char *filename, char *varname, clim_zm_t *zm) |
| Reads zonally averaged climatological data from a netCDF file and populates the given structure. More... | |
| void | read_ctl (const char *filename, int argc, char *argv[], ctl_t *ctl) |
| Reads control parameters from a configuration file and populates the given structure. More... | |
| void | read_kernel (const char *filename, double kz[EP], double kw[EP], int *nk) |
| Reads kernel function data from a file and populates the provided arrays. More... | |
| int | read_met (const char *filename, ctl_t *ctl, clim_t *clim, met_t *met) |
| Reads meteorological data from a file and populates the provided structures. More... | |
| void | read_met_bin_2d (FILE *in, met_t *met, float var[EX][EY], char *varname) |
| Reads a 2-dimensional meteorological variable from a binary file and stores it in the provided array. More... | |
| void | read_met_bin_3d (FILE *in, ctl_t *ctl, met_t *met, float var[EX][EY][EP], char *varname, float bound_min, float bound_max) |
| Reads 3D meteorological data from a binary file, potentially using different compression methods. More... | |
| void | read_met_cape (clim_t *clim, met_t *met) |
| Calculates Convective Available Potential Energy (CAPE) for each grid point. More... | |
| void | read_met_cloud (met_t *met) |
| Calculates cloud-related variables for each grid point. More... | |
| void | read_met_detrend (ctl_t *ctl, met_t *met) |
| Detrends meteorological data. More... | |
| void | read_met_extrapolate (met_t *met) |
| Extrapolates meteorological data. More... | |
| void | read_met_geopot (ctl_t *ctl, met_t *met) |
| Calculates geopotential heights from meteorological data. More... | |
| void | read_met_grid (const char *filename, int ncid, ctl_t *ctl, met_t *met) |
| Reads meteorological grid information from a NetCDF file. More... | |
| void | read_met_levels (int ncid, ctl_t *ctl, met_t *met) |
| Reads meteorological variables at different vertical levels from a NetCDF file. More... | |
| void | read_met_ml2pl (ctl_t *ctl, met_t *met, float var[EX][EY][EP], char *varname) |
| Interpolates meteorological data to specified pressure levels. More... | |
| void | read_met_monotonize (met_t *met) |
| Makes zeta and pressure profiles monotone. More... | |
| int | read_met_nc_2d (int ncid, char *varname, char *varname2, char *varname3, char *varname4, ctl_t *ctl, met_t *met, float dest[EX][EY], float scl, int init) |
| Reads a 2-dimensional meteorological variable from a NetCDF file. More... | |
| int | read_met_nc_3d (int ncid, char *varname, char *varname2, char *varname3, char *varname4, ctl_t *ctl, met_t *met, float dest[EX][EY][EP], float scl) |
| Reads a 3-dimensional meteorological variable from a NetCDF file. More... | |
| void | read_met_pbl (met_t *met) |
| Calculates the planetary boundary layer (PBL) height for each grid point. More... | |
| void | read_met_periodic (met_t *met) |
| Applies periodic boundary conditions to meteorological data along longitudinal axis. More... | |
| void | read_met_polar_winds (met_t *met) |
| Applies a fix for polar winds in meteorological data. More... | |
| void | read_met_pv (met_t *met) |
| Calculates potential vorticity (PV) from meteorological data. More... | |
| void | read_met_ozone (met_t *met) |
| Calculates the total column ozone from meteorological ozone data. More... | |
| void | read_met_sample (ctl_t *ctl, met_t *met) |
| Downsamples meteorological data based on specified parameters. More... | |
| void | read_met_surface (int ncid, met_t *met, ctl_t *ctl) |
| Reads surface meteorological data from a netCDF file and stores it in the meteorological data structure. More... | |
| void | read_met_tropo (ctl_t *ctl, clim_t *clim, met_t *met) |
| Calculates the tropopause and related meteorological variables based on various methods and stores the results in the meteorological data structure. More... | |
| void | read_obs (const char *filename, ctl_t *ctl, double *rt, double *rz, double *rlon, double *rlat, double *robs, int *nobs) |
| Reads observation data from a file and stores it in arrays. More... | |
| void | read_obs_asc (const char *filename, double *rt, double *rz, double *rlon, double *rlat, double *robs, int *nobs) |
| Reads observation data from an ASCII file. More... | |
| void | read_obs_nc (const char *filename, double *rt, double *rz, double *rlon, double *rlat, double *robs, int *nobs) |
| Reads observation data from a NetCDF file. More... | |
| double | scan_ctl (const char *filename, int argc, char *argv[], const char *varname, int arridx, const char *defvalue, char *value) |
| Scans a control file or command-line arguments for a specified variable. More... | |
| double | sedi (const double p, const double T, const double rp, const double rhop) |
| Calculates the sedimentation velocity of a particle in air. More... | |
| void | spline (const double *x, const double *y, const int n, const double *x2, double *y2, const int n2, const int method) |
| Performs spline interpolation or linear interpolation. More... | |
| float | stddev (const float *data, const int n) |
| Calculates the standard deviation of a set of data. More... | |
| double | sza_calc (const double sec, const double lon, const double lat) |
| Calculates the solar zenith angle. More... | |
| void | time2jsec (const int year, const int mon, const int day, const int hour, const int min, const int sec, const double remain, double *jsec) |
| Converts time components to seconds since January 1, 2000, 12:00:00 UTC. More... | |
| void | timer (const char *name, const char *group, int output) |
| Measures and reports elapsed time for named and grouped timers. More... | |
| double | time_from_filename (const char *filename, int offset) |
| Extracts and converts a timestamp from a filename to Julian seconds. More... | |
| double | tropo_weight (const clim_t *clim, const double t, const double lat, const double p) |
| Computes the weighting factor for a given pressure with respect to the tropopause. More... | |
| void | write_atm (const char *filename, ctl_t *ctl, atm_t *atm, double t) |
| Writes air parcel data to a file in various formats. More... | |
| void | write_atm_asc (const char *filename, ctl_t *ctl, atm_t *atm, double t) |
| Writes air parcel data to an ASCII file or gnuplot. More... | |
| void | write_atm_bin (const char *filename, ctl_t *ctl, atm_t *atm) |
| Writes air parcel data to a binary file. More... | |
| void | write_atm_clams (const char *filename, ctl_t *ctl, atm_t *atm) |
| Writes air parcel data to a NetCDF file in the CLaMS format. More... | |
| void | write_atm_clams_traj (const char *dirname, ctl_t *ctl, atm_t *atm, double t) |
| Writes CLaMS trajectory data to a NetCDF file. More... | |
| void | write_atm_nc (const char *filename, ctl_t *ctl, atm_t *atm) |
| Writes air parcel data to a NetCDF file. More... | |
| void | write_csi (const char *filename, ctl_t *ctl, atm_t *atm, double t) |
| Writes Critical Success Index (CSI) data to a file. More... | |
| void | write_ens (const char *filename, ctl_t *ctl, atm_t *atm, double t) |
| Writes ensemble data to a file. More... | |
| void | write_grid (const char *filename, ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, double t) |
| Writes grid data to a file in ASCII or netCDF format. More... | |
| void | write_grid_asc (const char *filename, ctl_t *ctl, double *cd, double *mean[NQ], double *sigma[NQ], double *vmr_impl, double t, double *z, double *lon, double *lat, double *area, double dz, int *np) |
| Writes grid data to an ASCII file. More... | |
| void | write_grid_nc (const char *filename, ctl_t *ctl, double *cd, double *mean[NQ], double *sigma[NQ], double *vmr_impl, double t, double *z, double *lon, double *lat, double *area, double dz, int *np) |
| Writes grid data to a NetCDF file. More... | |
| int | write_met (const char *filename, ctl_t *ctl, met_t *met) |
| Writes meteorological data to a binary file. More... | |
| void | write_met_bin_2d (FILE *out, met_t *met, float var[EX][EY], char *varname) |
| Writes a 2-dimensional meteorological variable to a binary file. More... | |
| void | write_met_bin_3d (FILE *out, ctl_t *ctl, met_t *met, float var[EX][EY][EP], char *varname, int precision, double tolerance) |
| Writes a 3-dimensional meteorological variable to a binary file. More... | |
| void | write_output (const char *dirname, ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, double t) |
| Writes various types of output data to files in a specified directory. More... | |
| void | write_prof (const char *filename, ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, double t) |
| Writes profile data to a specified file. More... | |
| void | write_sample (const char *filename, ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, double t) |
| Writes sample data to a specified file. More... | |
| void | write_station (const char *filename, ctl_t *ctl, atm_t *atm, double t) |
| Writes station data to a specified file. More... | |
| void | write_vtk (const char *filename, ctl_t *ctl, atm_t *atm, double t) |
| Writes VTK (Visualization Toolkit) data to a specified file. More... | |
MPTRAC library declarations.
Definition in file mptrac.h.
| #define CPD 1003.5 |
| #define LV 2501000. |
| #define KB 1.3806504e-23 |
| #define NP 10000000 |
| #define NCSI 1000000 |
| #define EP 140 |
| #define EX 1202 |
| #define NENS 2000 |
| #define NOBS 10000000 |
| #define CY 250 |
| #define CO3 30 |
| #define CP 70 |
| #define CSZA 50 |
| #define CT 12 |
| #define CTS 1000 |
| #define ALLOC | ( | ptr, | |
| type, | |||
| n | |||
| ) |
Allocate memory for a pointer with error handling.
This macro allocates memory for a pointer of a given type and size using the calloc function. It includes error handling to check if memory allocation was successful. If the code is being compiled with OpenACC support (_OPENACC macro defined), it additionally checks if the code is running on a GPU device, and if not, it raises an error.
| ptr | Pointer variable to be allocated. |
| type | Data type of the pointer. |
| n | Number of elements to allocate memory for. |
| #define ARRAY_2D | ( | ix, | |
| iy, | |||
| ny | |||
| ) | ((ix) * (ny) + (iy)) |
Macro for computing the linear index of a 2D array element.
The ARRAY_2D macro computes the linear index of a 2D array element based on the specified row index (ix), column index (iy), and number of columns (ny).
| ix | Integer representing the row index of the 2D array element. |
| iy | Integer representing the column index of the 2D array element. |
| ny | Integer representing the number of columns in the 2D array. |
The macro computes the linear index using the formula: (ix) * (ny) + (iy). This formula assumes row-major storage, where elements of each row are stored sequentially in memory.
| #define ARRAY_3D | ( | ix, | |
| iy, | |||
| ny, | |||
| iz, | |||
| nz | |||
| ) | (((ix)*(ny) + (iy)) * (nz) + (iz)) |
Compute the linear index of a 3D array element.
This macro computes the linear index of a 3D array element based on the specified row index (ix), column index (iy), depth index (iz), number of columns (ny), and number of depths (nz).
| ix | Row index of the 3D array element. |
| iy | Column index of the 3D array element. |
| ny | Number of columns in the 3D array. |
| iz | Depth index of the 3D array element. |
| nz | Number of depths in the 3D array. |
| #define ARRHENIUS | ( | a, | |
| b, | |||
| t | |||
| ) | ((a) * exp( -(b) / (t))) |
Calculate the Arrhenius rate constant.
The Arrhenius equation is commonly used in chemical kinetics to describe the temperature dependence of reaction rates. This macro calculates the rate constant (k) based on the Arrhenius equation:
\[ k = a \times \exp( -b / T ), \]
where:
| a | Pre-exponential factor or frequency factor. |
| b | Activation energy. |
| t | Temperature in Kelvin. |
| #define DEG2DX | ( | dlon, | |
| lat | |||
| ) | ((dlon) * M_PI * RE / 180. * cos((lat) / 180. * M_PI)) |
Convert a longitude difference to a distance in the x-direction (east-west) at a specific latitude.
This macro calculates the distance in the x-direction (east-west) corresponding to a given longitude difference at a specific latitude using the formula:
\[ dx = dlon \times \pi \times RE / 180 times \cos(lat), \]
where:
| dlon | Difference in longitudes in degrees. |
| lat | Latitude in degrees. |
| #define DEG2DY | ( | dlat | ) | ((dlat) * M_PI * RE / 180.) |
Convert a latitude difference to a distance in the y-direction (north-south).
This macro calculates the distance in the y-direction (north-south) corresponding to a given latitude difference using the formula:
\[ dy = dlat \times \pi \times RE / 180, \]
where:
| dlat | Difference in latitudes in degrees. |
| #define DP2DZ | ( | dp, | |
| p | |||
| ) | (- (dp) * H0 / (p)) |
Convert a pressure difference to a height difference in the vertical direction.
This macro calculates the change in height (altitude) corresponding to a given pressure difference using the formula:
\[ dz = - (dp) \times H_0 / p \]
where:
| dp | Pressure difference in hPa. |
| p | Reference pressure in hPa. |
| #define DX2DEG | ( | dx, | |
| lat | |||
| ) |
Convert a distance in kilometers to degrees longitude at a given latitude.
This macro calculates the change in longitude in degrees corresponding to a given distance in kilometers at a specified latitude on the Earth's surface. It uses the formula:
\[ dlon = \frac{dx \times 180}{\pi \times RE \times \cos(lat)} \]
| dx | Distance in kilometers. |
| lat | Latitude in degrees. |
| #define DY2DEG | ( | dy | ) | ((dy) * 180. / (M_PI * RE)) |
Convert a distance in kilometers to degrees latitude.
This macro calculates the change in latitude in degrees corresponding to a given distance in kilometers on the Earth's surface. It uses the formula:
\[ dlat = \frac{dy \times 180}{\pi \times RE} \]
| dy | Distance in kilometers. |
| #define DZ2DP | ( | dz, | |
| p | |||
| ) | (-(dz) * (p) / H0) |
Convert a change in altitude to a change in pressure.
This macro calculates the change in pressure corresponding to a given change in altitude. It uses the hydrostatic equation:
\[ dp = -\left(dz \times \frac{p}{H_0}\right) \]
| dz | Change in altitude in kilometers. |
| p | Current pressure in hPa. |
| #define DIST | ( | a, | |
| b | |||
| ) | sqrt(DIST2(a, b)) |
Calculate the distance between two points in Cartesian coordinates.
This macro calculates the Euclidean distance between two points in Cartesian coordinates. It uses the square root of the square of the distance obtained from the DIST2 macro.
| a | Coordinates of the first point as an array of doubles. |
| b | Coordinates of the second point as an array of doubles. |
| #define DIST2 | ( | a, | |
| b | |||
| ) | ((a[0]-b[0])*(a[0]-b[0])+(a[1]-b[1])*(a[1]-b[1])+(a[2]-b[2])*(a[2]-b[2])) |
Calculate the squared Euclidean distance between two points in Cartesian coordinates.
This macro calculates the squared Euclidean distance between two points in Cartesian coordinates. It computes the sum of the squares of the differences of corresponding coordinates.
| a | Coordinates of the first point as an array of doubles. |
| b | Coordinates of the second point as an array of doubles. |
| #define DOTP | ( | a, | |
| b | |||
| ) | (a[0]*b[0]+a[1]*b[1]+a[2]*b[2]) |
Calculate the dot product of two vectors.
This macro computes the dot product of two vectors represented as arrays of doubles. It multiplies corresponding components of the vectors and sums the results.
| a | The first vector as an array of doubles. |
| b | The second vector as an array of doubles. |
| #define FMOD | ( | x, | |
| y | |||
| ) | ((x) - (int) ((x) / (y)) * (y)) |
Calculate the floating-point remainder of dividing x by y.
This macro computes the floating-point remainder of dividing x by y. It calculates this remainder as x minus the integer part of (x / y) times y.
| x | The dividend. |
| y | The divisor. |
| #define FREAD | ( | ptr, | |
| type, | |||
| size, | |||
| in | |||
| ) |
Read data from a file stream and store it in memory.
This macro reads data of a specified type from the given input file stream and stores it in the specified memory location. It ensures that the correct amount of data is read from the file stream, and if not, it raises an error.
| ptr | Pointer to the memory location where the data will be stored. |
| type | Type of the data elements to be read. |
| size | Number of elements to read. |
| in | File stream from which to read the data. |
| #define FWRITE | ( | ptr, | |
| type, | |||
| size, | |||
| out | |||
| ) |
Write data from memory to a file stream.
This macro writes data of a specified type from the specified memory location to the given output file stream. It ensures that the correct amount of data is written to the file stream, and if not, it raises an error.
| ptr | Pointer to the memory location containing the data to be written. |
| type | Type of the data elements to be written. |
| size | Number of elements to write. |
| out | File stream to which the data will be written. |
| #define INTPOL_INIT double cw[4] = {0.0, 0.0, 0.0, 0.0}; int ci[3] = {0, 0, 0}; |
| #define INTPOL_2D | ( | var, | |
| init | |||
| ) |
Perform 2D interpolation for a meteorological variable.
This macro performs 2D interpolation for a given meteorological variable at a specific time and location.
| var | The variable to interpolate. |
| init | A flag indicating whether to initialize the interpolation arrays (cw and ci). Set to 1 for initialization, 0 otherwise. |
var.| #define INTPOL_3D | ( | var, | |
| init | |||
| ) |
Perform 3D interpolation for a meteorological variable.
This macro performs 3D interpolation for a given meteorological variable at a specific time, pressure level, and location.
| var | The variable to interpolate. |
| init | A flag indicating whether to initialize the interpolation arrays (cw and ci). Set to 1 for initialization, 0 otherwise. |
var.| #define INTPOL_SPACE_ALL | ( | p, | |
| lon, | |||
| lat | |||
| ) |
Interpolate multiple meteorological variables in space.
This macro performs spatial interpolation for multiple meteorological variables at a given pressure level, longitude, and latitude.
| p | The pressure level at which to interpolate the variables. |
| lon | The longitude at which to interpolate the variables. |
| lat | The latitude at which to interpolate the variables. |
| #define INTPOL_TIME_ALL | ( | time, | |
| p, | |||
| lon, | |||
| lat | |||
| ) |
Interpolate multiple meteorological variables in time.
This macro performs temporal interpolation for multiple meteorological variables at a given time, pressure level, longitude, and latitude.
| time | The time at which to interpolate the variables. |
| p | The pressure level at which to interpolate the variables. |
| lon | The longitude at which to interpolate the variables. |
| lat | The latitude at which to interpolate the variables. |
| #define LAPSE | ( | p1, | |
| t1, | |||
| p2, | |||
| t2 | |||
| ) |
Calculate lapse rate.
This macro calculates the lapse rate between two pressure levels given their temperatures and pressures.
| p1 | Pressure at the first level (in hPa). |
| t1 | Temperature at the first level (in K). |
| p2 | Pressure at the second level (in hPa). |
| t2 | Temperature at the second level (in K). |
| #define LIN | ( | x0, | |
| y0, | |||
| x1, | |||
| y1, | |||
| x | |||
| ) | ((y0)+((y1)-(y0))/((x1)-(x0))*((x)-(x0))) |
Linear interpolation.
This macro performs linear interpolation to estimate the value of y at a given x based on two points (x0, y0) and (x1, y1).
| x0 | X-coordinate of the first point. |
| y0 | Y-coordinate of the first point. |
| x1 | X-coordinate of the second point. |
| y1 | Y-coordinate of the second point. |
| x | The x-coordinate at which to estimate the y-value. |
| #define MAX | ( | a, | |
| b | |||
| ) | (((a)>(b))?(a):(b)) |
Macro to determine the maximum of two values.
This macro evaluates to the larger of its two arguments, a and b. It uses a ternary conditional operator to compare the values of a and b and returns a if a is greater than b; otherwise, it returns b.
| a | The first value to compare. Can be of any type that supports comparison. |
| b | The second value to compare. Can be of any type that supports comparison. |
a or b.a and b are evaluated twice. If a or b have side effects (e.g., increment operators, function calls), the side effects will occur more than once. This can lead to unexpected behavior.a and b should be of compatible types to avoid potential issues with comparison and return value.| #define MET_HEADER |
Write header for meteorological data file.
This macro writes a header to a meteorological data file, providing information about the variables stored in the file and their corresponding columns.
| out | Pointer to the file stream where the header will be written. |
| #define MIN | ( | a, | |
| b | |||
| ) | (((a)<(b))?(a):(b)) |
Macro to determine the minimum of two values.
This macro evaluates to the smaller of its two arguments, a and b. It uses a ternary conditional operator to compare the values of a and b and returns a if a is less than b; otherwise, it returns b.
| a | The first value to compare. Can be of any type that supports comparison. |
| b | The second value to compare. Can be of any type that supports comparison. |
a or b.a and b are evaluated twice. If a or b have side effects (e.g., increment operators, function calls), the side effects will occur more than once. This can lead to unexpected behavior.a and b should be of compatible types to avoid potential issues with comparison and return value.Calculate the density of a gas molecule.
This macro calculates the density of a gas molecule using the provided pressure and temperature values.
| p | Pressure of the gas in Pascals. |
| t | Temperature of the gas in Kelvin. |
| #define NC | ( | cmd | ) |
Execute a NetCDF command and check for errors.
This macro executes a NetCDF command and checks the result. If the result indicates an error, it prints the error message using ERRMSG.
| cmd | NetCDF command to execute. |
| #define NC_DEF_VAR | ( | varname, | |
| type, | |||
| ndims, | |||
| dims, | |||
| long_name, | |||
| units | |||
| ) |
Define a NetCDF variable with attributes.
This macro defines a NetCDF variable with the specified name, data type, dimensions, long name, and units. It also sets the "long_name" and "units" attributes for the variable.
| varname | Name of the variable. |
| type | Data type of the variable. |
| ndims | Number of dimensions for the variable. |
| dims | Array of dimension IDs. |
| long_name | Long name of the variable. |
| units | Units of the variable. |
| #define NC_GET_DOUBLE | ( | varname, | |
| ptr, | |||
| force | |||
| ) |
Retrieve a double-precision variable from a NetCDF file.
This macro retrieves a double-precision variable from a NetCDF file. It first checks if the variable exists in the file and then reads its data into the specified pointer. If the force parameter is set to true, it forces the retrieval of the variable, raising an error if the variable does not exist. If force is false, it retrieves the variable if it exists and issues a warning if it does not.
| varname | Name of the variable to retrieve. |
| ptr | Pointer to the memory location where the data will be stored. |
| force | Boolean flag indicating whether to force retrieval (true) or not (false). |
| #define NC_INQ_DIM | ( | dimname, | |
| ptr, | |||
| min, | |||
| max | |||
| ) |
Inquire the length of a dimension in a NetCDF file.
This macro retrieves the length of a specified dimension from a NetCDF file. It checks if the length of the dimension is within a specified range and assigns the length to the provided pointer. If the length is outside the specified range, an error message is raised.
| dimname | Name of the dimension to inquire. |
| ptr | Pointer to an integer where the dimension length will be stored. |
| min | Minimum acceptable length for the dimension. |
| max | Maximum acceptable length for the dimension. |
| #define NC_PUT_DOUBLE | ( | varname, | |
| ptr, | |||
| hyperslab | |||
| ) |
Write double precision data to a NetCDF variable.
This macro writes data to a specified NetCDF variable. It can handle both full variable writes and hyperslab writes depending on the hyperslab parameter. If hyperslab is true, the data is written as a hyperslab; otherwise, the entire variable is written.
| varname | Name of the NetCDF variable to write to. |
| ptr | Pointer to the data to be written. |
| hyperslab | Boolean indicating whether to write the data as a hyperslab. |
| #define NC_PUT_FLOAT | ( | varname, | |
| ptr, | |||
| hyperslab | |||
| ) |
Write a float array to a NetCDF file.
This macro writes a float array to a specified variable in a NetCDF file. Depending on the value of the hyperslab parameter, the data can be written as a hyperslab or as a whole variable.
| varname | Name of the variable to which the float array will be written. |
| ptr | Pointer to the float array to be written. |
| hyperslab | Boolean flag indicating if the data should be written as a hyperslab.
|
| #define NC_PUT_INT | ( | varname, | |
| ptr, | |||
| hyperslab | |||
| ) |
Write integer data to a NetCDF variable.
This macro writes data to a specified NetCDF variable. It can handle both full variable writes and hyperslab writes depending on the hyperslab parameter. If hyperslab is true, the data is written as a hyperslab; otherwise, the entire variable is written.
| varname | Name of the NetCDF variable to write to. |
| ptr | Pointer to the data to be written. |
| hyperslab | Boolean indicating whether to write the data as a hyperslab. |
| #define NC_PUT_ATT | ( | varname, | |
| attname, | |||
| text | |||
| ) |
Add a text attribute to a NetCDF variable.
This macro adds a text attribute to a specified NetCDF variable. It first retrieves the variable ID using its name, then it attaches the text attribute to the variable.
| varname | Name of the NetCDF variable to which the attribute will be added. |
| attname | Name of the attribute to be added. |
| text | Text of the attribute to be added. |
| #define NC_PUT_ATT_GLOBAL | ( | attname, | |
| text | |||
| ) | NC(nc_put_att_text(ncid, NC_GLOBAL, attname, strnlen(text, LEN), text)); |
Add a global text attribute to a NetCDF file.
This macro adds a text attribute to the global attributes of a NetCDF file. It directly attaches the attribute to the file, rather than to a specific variable.
| attname | Name of the global attribute to be added. |
| text | Text of the attribute to be added. |
| #define NN | ( | x0, | |
| y0, | |||
| x1, | |||
| y1, | |||
| x | |||
| ) | (fabs((x) - (x0)) <= fabs((x) - (x1)) ? (y0) : (y1)) |
Perform nearest-neighbor interpolation.
This macro returns the value of the nearest neighbor (y0 or y1) for a given x value. It compares the distances between x and x0, and between x and x1, and returns the y value corresponding to the closer x value.
| x0 | The x-coordinate of the first point. |
| y0 | The y-coordinate of the first point. |
| x1 | The x-coordinate of the second point. |
| y1 | The y-coordinate of the second point. |
| x | The x-coordinate for which the nearest neighbor is to be found. |
| #define NORM | ( | a | ) | sqrt(DOTP(a, a)) |
Compute the norm (magnitude) of a vector.
This macro computes the Euclidean norm (magnitude) of a vector a using the dot product of the vector with itself. The vector is assumed to have three components: a[0], a[1], and a[2].
| a | The vector for which the norm is to be computed. |
| #define PARTICLE_LOOP | ( | ip0, | |
| ip1, | |||
| check_dt, | |||
| ... | |||
| ) |
Loop over particle indices with OpenACC acceleration.
This macro defines a loop over particle indices from ip0 to ip1 with optional checking of dt. If _OPENACC is defined, the loop is accelerated using OpenACC directives. Otherwise, OpenMP parallelization is used.
| ip0 | The starting index of the loop (inclusive). |
| ip1 | The ending index of the loop (exclusive). |
| check_dt | Flag indicating whether to check the array dt for non-zero values. |
| ... | Optional pragma directives to be applied. |
Compute pressure at given altitude.
This macro calculates the pressure at a given altitude using the barometric formula.
| z | The altitude in kilometers. |
The barometric formula used for this calculation is:
\[ P(z) = P_0 \times e^{-(z / H_0)}, \]
where:
Note: The constants \( P_0 \) and \( H_0 \) must be defined before using this macro.
Compute saturation pressure over water.
This macro calculates the saturation pressure over water based on the WMO (2018) formula.
| t | The temperature in degrees Celsius. |
The saturation pressure over water is calculated using the formula:
\[ P_{\textrm{sat}}(t) = 6.112 \times e^{17.62 \times \frac{(t - T_0)}{243.12 + (t - T_0)}}, \]
where:
Note: The constants \( T_0 \) must be defined before using this macro.
Compute saturation pressure over ice (WMO, 2018).
This macro calculates the saturation pressure over ice based on the WMO (2018) formula.
| t | The temperature in K. |
The saturation pressure over ice is calculated using the formula:
\[ P_{\textrm{ice}}(t) = 6.112 \times e^{22.46 \times \frac{(t - T_0)}{272.62 + (t - T_0)}}, \]
where:
Note: The constant \( T_0 \) must be defined before using this macro.
Calculate partial water vapor pressure.
This macro calculates the partial water vapor pressure using the given total pressure and water vapor mixing ratio.
| p | The total pressure in hPa (hectopascals). |
| h2o | The water vapor mixing ratio in ppv (parts per volume). |
The partial water vapor pressure is calculated using the formula:
\[ P_{\textrm{w}}(p, h_2o) = \frac{p \times \max(h_2o, 0.1 \times 10^{-6})}{1 + (1 - \epsilon) \times \max(h_2o, 0.1 \times 10^{-6})}, \]
where:
Note: The constant \( \epsilon \) must be defined before using this macro.
Compute relative humidity over water.
This macro calculates the relative humidity over water using the given total pressure, temperature, and water vapor mixing ratio.
| p | The total pressure in hPa. |
| t | The temperature in K. |
| h2o | The water vapor mixing ratio in ppv (parts per volume). |
The relative humidity over water is calculated using the formula:
\[ RH_{\textrm{w}}(p, t, h_2o) = \frac{P_{\textrm{w}}(p, h_2o)}{P_{\textrm{sat}}(t)} \times 100, \]
where:
Note: The macros PW() and PSAT() must be defined before using this macro.
Compute relative humidity over ice.
This macro calculates the relative humidity over ice using the given total pressure, temperature, and water vapor mixing ratio.
| p | The total pressure in hPa. |
| t | The temperature in K. |
| h2o | The water vapor mixing ratio in ppv (parts per volume). |
The relative humidity over ice is calculated using the formula:
\[ RH_{\textrm{ice}}(p, t, h_2o) = \frac{P_{\textrm{w}}(p, h_2o)}{P_{\textrm{ice}}(t)} \times 100, \]
where:
Note: The macros PW() and PSICE() must be defined before using this macro.
| #define RHO | ( | p, | |
| t | |||
| ) | (100. * (p) / (RA * (t))) |
Compute density of air.
This macro calculates the density of air using the given total pressure and temperature.
| p | The total pressure in hPa. |
| t | The temperature in K. |
The density of air is calculated using the formula:
\[ \rho(p, t) = \frac{100 \times p}{R_a \times t}, \]
where:
| #define ROETH_PHOTOL | ( | a, | |
| b, | |||
| c, | |||
| sza | |||
| ) | ((c)*(sza) < M_PI/2. ? (a) * exp((b) * (1 - 1/cos((c) * (sza)))) : 0) |
Roeth approximation formula for photolysis reactions.
This macro calculates the rate of a photolysis reaction using the Roeth approximation formula, which takes into account the solar zenith angle (SZA).
| a | Coefficient 'a' in the Roeth formula. |
| b | Coefficient 'b' in the Roeth formula. |
| c | Coefficient 'c' in the Roeth formula. |
| sza | The solar zenith angle in radians. |
The Roeth approximation formula for photolysis reactions is given by:
\[ ROETH\_PHOTOL(a, b, c, sza) = \begin{cases} a \times \exp\left(b \times \left(1 - \frac{1}{\cos(c \times sza)}\right)\right), \textrm{if } c \times sza < \frac{\pi}{2} \\ 0, \textrm{otherwise} \end{cases} \]
where:
| #define SET_ATM | ( | qnt, | |
| val | |||
| ) |
Set atmospheric quantity value.
This macro sets the value of a specific atmospheric quantity at a given index 'ip'. The macro first checks if the control index 'ctl->qnt' is non-negative before assigning the value, ensuring that the quantity index is valid.
| qnt | The index representing the atmospheric quantity to set. |
| val | The value to set for the atmospheric quantity. |
| #define SET_QNT | ( | qnt, | |
| name, | |||
| longname, | |||
| unit | |||
| ) |
Set atmospheric quantity index.
This macro sets the index, long name, and unit of a specific atmospheric quantity based on its name. It compares the name parameter with the name of the atmospheric quantity stored in 'ctl->qnt_name'. If a match is found, it assigns the index to 'ctl->qnt', updates the long name, and updates the unit.
| qnt | The index representing the atmospheric quantity. |
| name | The name of the atmospheric quantity. |
| longname | The long name of the atmospheric quantity. |
| unit | The unit of the atmospheric quantity. |
Compute specific humidity from water vapor volume mixing ratio.
This macro calculates the specific humidity from the water vapor volume mixing ratio. Specific humidity represents the ratio of the mass of water vapor to the total mass of air and is dimensionless.
| h2o | The water vapor volume mixing ratio. |
| #define SQR | ( | x | ) | ((x)*(x)) |
| #define SWAP | ( | x, | |
| y, | |||
| type | |||
| ) | do {type tmp = x; x = y; y = tmp;} while(0); |
| #define TDEW | ( | p, | |
| h2o | |||
| ) |
Calculate dew point temperature.
This macro computes the dew point temperature using the formula provided by the World Meteorological Organization (WMO, 2018).
| p | The atmospheric pressure in hPa. |
| h2o | The water vapor volume mixing ratio. |
Formula:
\[ T_{\textrm{dew}} = T_0 + \frac{243.12 \times \ln\left(\frac{{P_W(p, h_{2}O)}}{6.112}\right)}{17.62 - \ln\left(\frac{{P_W(p, h_{2}O)}}{6.112}\right)} \]
where:
| #define TICE | ( | p, | |
| h2o | |||
| ) |
Calculate frost point temperature (WMO, 2018).
This macro computes the frost point temperature using the formula provided by the World Meteorological Organization (WMO, 2018).
| p | The atmospheric pressure in hPa. |
| h2o | The water vapor volume mixing ratio. |
Formula:
\[ T_{\textrm{ice}} = T_0 + \frac{272.62 \times \ln\left(\frac{{P_W(p, h_{2}O)}}{6.112}\right)}{22.46 - \ln\left(\frac{{P_W(p, h_{2}O)}}{6.112}\right)} \]
where:
| #define THETA | ( | p, | |
| t | |||
| ) | ((t) * pow(1000. / (p), 0.286)) |
Compute potential temperature.
This macro calculates the potential temperature of the atmosphere.
| p | The atmospheric pressure in hPa. |
| t | The temperature in Kelvin. |
Formula:
\[ \theta = T \left( \frac{1000}{P} \right)^{0.286} \]
where:
Compute virtual potential temperature.
This macro calculates the virtual potential temperature of the atmosphere, which takes into account the effect of water vapor on the atmosphere's buoyancy.
| p | The atmospheric pressure in hPa. |
| t | The temperature in Kelvin. |
| h2o | The water vapor volume mixing ratio (ppv). |
Formula:
The virtual potential temperature ( \( \theta_v \)) is computed as
\[ \theta_v = \theta \left( 1 + \frac{0.61 \times q}{\epsilon} \right), \]
where:
| #define TOK | ( | line, | |
| tok, | |||
| format, | |||
| var | |||
| ) |
Get string tokens.
This macro extracts tokens from a given string, typically used for parsing input lines.
| line | The input string containing tokens. |
| tok | A pointer to the token string. |
| format | The format string specifying the expected format of the token. |
| var | The variable to store the parsed token value. |
The macro tokenizes the input line using space and tab characters as delimiters. It then parses each token according to the specified format string and stores the parsed value in the provided variable.
Compute virtual temperature.
This macro calculates the virtual temperature of air given its temperature and water vapor volume mixing ratio.
| t | The temperature of the air in Kelvin. |
| h2o | The water vapor volume mixing ratio. |
The virtual temperature (T_v) is computed as the temperature (t) multiplied by (1 + (1 - EPS) * max(h2o, 0.1e-6)), where EPS is the ratio of the molar mass of water vapor to the molar mass of dry air.
Convert pressure to altitude.
This macro calculates the altitude from the given pressure using the barometric formula.
| p | The pressure in hPa (hectopascal). |
Formula:
The altitude (z) is computed as H0 times the natural logarithm of the ratio of the reference pressure (P0) to the given pressure (p), where H0 is the scale height and P0 is the reference pressure at sea level.
| #define ZDIFF | ( | lnp0, | |
| t0, | |||
| h2o0, | |||
| lnp1, | |||
| t1, | |||
| h2o1 | |||
| ) |
Calculate geopotential height difference.
This macro calculates the geopotential height difference between two pressure levels using the hypsometric equation.
| lnp0 | The natural logarithm of the pressure at the first level. |
| t0 | The temperature at the first level in Kelvin (K). |
| h2o0 | The water vapor volume mixing ratio at the first level. |
| lnp1 | The natural logarithm of the pressure at the second level. |
| t1 | The temperature at the second level in Kelvin (K). |
| h2o1 | The water vapor volume mixing ratio at the second level. |
Formula: The geopotential height difference (dz) is computed as a function of the difference in natural logarithm of pressure (lnp) between the two levels, the average virtual temperature (ThetaVirt) of the two levels, the specific gas constant for dry air (RI), and the acceleration due to gravity at the surface of the Earth (G0).
| #define ZETA | ( | ps, | |
| p, | |||
| t | |||
| ) |
Calculate potential vorticity using the Zeta approximation.
This macro calculates the potential vorticity using the Zeta approximation, which is a function of pressure, surface pressure, and temperature.
| ps | Surface pressure in hPa. |
| p | Pressure at the given level in hPa. |
| t | Temperature at the given level in Kelvin (K). |
Formula: The potential vorticity (Zeta) is computed based on the given pressure (p), surface pressure (ps), and temperature (t). It involves the potential temperature (Theta) at the given pressure level.
| #define LOGLEV 2 |
| #define LOG | ( | level, | |
| ... | |||
| ) |
Print a log message with a specified logging level.
This macro prints a formatted log message to the standard output if the specified logging level meets certain conditions. The message will be indented if the logging level is greater than or equal to 2.
| level | The logging level of the message. This should be an integer value. |
| ... | The formatted message string and its arguments, similar to printf. |
The LOG macro provides a simple way to log messages with different levels of importance. The message is only printed if the specified level is less than or equal to the pre-defined LOGLEV macro. If the level is greater than or equal to 2, the message is preceded by two spaces for indentation.
The macro expands to a block of code that:
level is greater than or equal to 2, and if so, prints two spaces.level is less than or equal to LOGLEV, and if so, prints the formatted message followed by a newline.LOGLEV macro must be defined with an appropriate logging level before using the LOG macro.| #define WARN | ( | ... | ) |
Print a warning message with contextual information.
This macro prints a formatted warning message to the standard output, including the file name, function name, and line number where the warning occurred. The message is then passed to the LOG macro with a logging level of 0.
| ... | The formatted warning message string and its arguments, similar to printf. |
The WARN macro is used to print warning messages with additional context about where the warning was triggered. The message includes the following contextual information:
__FILE__).__func__).__LINE__).After printing this contextual information, the macro uses the LOG macro with a logging level of 0 to print the actual warning message. This ensures that warning messages are always logged, regardless of the value of LOGLEV.
LOG macro must be defined before using the WARN macro.| #define ERRMSG | ( | ... | ) |
Print an error message with contextual information and terminate the program.
This macro prints a formatted error message to the standard output, including the file name, function name, and line number where the error occurred. After printing the message, the program is terminated with an exit status indicating failure.
| ... | The formatted error message string and its arguments, similar to printf. |
The ERRMSG macro is used to report critical errors that require the program to terminate immediately. The message includes the following contextual information:
__FILE__).__func__).__LINE__).After printing this contextual information, the macro uses the LOG macro with a logging level of 0 to print the actual error message. Finally, the program exits with a failure status (EXIT_FAILURE).
LOG macro must be defined before using the ERRMSG macro.| #define PRINT | ( | format, | |
| var | |||
| ) |
Print the value of a variable with contextual information.
This macro prints the value of a variable to the standard output, including the file name, function name, and line number where the macro is called. The output also includes the variable's name and value in a formatted string.
| format | The format string used to print the variable's value, similar to printf. |
| var | The variable to be printed. |
The PRINT macro is used to output the value of a variable along with additional context about where the macro is called. The message includes:
__FILE__).__func__).__LINE__).#var).format).This macro is particularly useful for debugging purposes, providing a convenient way to trace variable values and their locations in the code.
| #define PRINT_TIMERS timer("END", "END", 1); |
Print the current state of all timers.
This macro calls the timer function with predefined arguments to signify the end of the timer logging process. It is used to print the results of all the timers that have been tracked.
timer function must be defined elsewhere in the codebase for this macro to function correctly.| #define SELECT_TIMER | ( | id, | |
| group, | |||
| color | |||
| ) |
Select and start a timer with specific attributes.
This macro stops the current timer (if any) and starts a new timer with the specified ID, group, and color. It uses the NVTX_POP and NVTX_PUSH macros for managing timer events and the timer function to log the timer start event.
| id | The identifier for the timer. |
| group | The group name associated with the timer. |
| color | The color code associated with the timer for NVTX visualization. |
NVTX_POP, NVTX_PUSH, and timer functions/macros must be defined elsewhere in the codebase for this macro to function correctly.| #define START_TIMERS NVTX_PUSH("START", NVTX_CPU); |
Starts a timer for tracking.
This macro initializes the timer tracking process by pushing a start event onto the stack using the NVTX_PUSH macro with a predefined ID ("START") and color (NVTX_CPU).
NVTX_PUSH macro must be defined elsewhere in the codebase for this macro to function correctly.| #define STOP_TIMERS NVTX_POP; |
| void thrustSortWrapper | ( | double *__restrict__ | c, |
| int | n, | ||
| int *__restrict__ | index | ||
| ) |
Wrapper to Thrust sorting function.
This function serves as a wrapper for a Thrust sorting operation, sorting the array c and maintaining the correspondence with the index array.
| c | Pointer to the array of double values to be sorted. |
| n | The number of elements in the array c. |
| index | Pointer to the array of indices, which will be updated to reflect the sorted order. |
The thrustSortWrapper function uses the Thrust library to sort the array c in ascending order. The index array is updated to reflect the new order of elements in c after sorting.
This function is particularly useful when the sorted order of elements needs to be tracked by indices.
c and index arrays must be of the same length n.| void broadcast_large_data | ( | void * | data, |
| size_t | N | ||
| ) |
Broadcasts large data across all processes in an MPI communicator.
This function divides the data into manageable chunks and broadcasts each chunk sequentially. This approach is necessary because the data size may exceed the maximum allowable message size for a single MPI_Bcast operation.
| data | Pointer to the data to be broadcasted. |
| N | Size of the data in bytes. |
The function first broadcasts the total size of the data to all processes. Then, it calculates the number of chunks needed to broadcast the entire data. Each chunk is broadcasted in sequence until the entire data has been sent.
The maximum chunk size is defined as CHUNK_SIZE (2147483647 bytes).
| void cart2geo | ( | const double * | x, |
| double * | z, | ||
| double * | lon, | ||
| double * | lat | ||
| ) |
Converts Cartesian coordinates to geographic coordinates.
This function converts a point from Cartesian coordinates (x, y, z) to geographic coordinates (longitude, latitude, and altitude). It uses the spherical Earth approximation for the conversion.
| x | Pointer to an array containing the Cartesian coordinates (x, y, z) in kilometers. |
| z | Pointer to a double where the computed altitude (above the reference ellipsoid) will be stored, in kilometers. |
| lon | Pointer to a double where the computed longitude (in degrees) will be stored. |
| lat | Pointer to a double where the computed latitude (in degrees) will be stored. |
Converts Cartesian coordinates to geographic coordinates.
Definition at line 74 of file mptrac.c.
| double clim_oh | ( | const ctl_t * | ctl, |
| const clim_t * | clim, | ||
| const double | t, | ||
| const double | lon, | ||
| const double | lat, | ||
| const double | p | ||
| ) |
Calculates the hydroxyl radical (OH) concentration from climatology data, with an optional diurnal correction based on solar zenith angle.
This function retrieves OH data from a given climatology and applies a diurnal correction if the correction factor (oh_chem_beta) is greater than zero. The diurnal correction accounts for the variation in OH concentration due to changes in the solar zenith angle.
| ctl | Pointer to the control structure containing configuration parameters. |
| clim | Pointer to the climatology structure containing OH data. |
| t | Time at which the OH concentration is to be calculated. |
| lon | Longitude at which the OH concentration is to be calculated. |
| lat | Latitude at which the OH concentration is to be calculated. |
| p | Pressure level at which the OH concentration is to be calculated. |
Definition at line 88 of file mptrac.c.
Applies a diurnal correction to the hydroxyl radical (OH) concentration in climatology data.
This function iterates over the climatology data points for OH concentration and integrates the day/night correction factor over longitude. The correction factor is based on the solar zenith angle, and it adjusts the OH data to account for diurnal variations. The corrected OH data is scaled accordingly.
| ctl | Pointer to the control structure containing configuration parameters, including the correction factor (oh_chem_beta). |
| clim | Pointer to the climatology structure containing OH data that will be corrected. |
Definition at line 112 of file mptrac.c.
| double clim_photo | ( | double | rate[CP][CSZA][CO3], |
| clim_photo_t * | photo, | ||
| double | p, | ||
| double | sza, | ||
| double | o3c | ||
| ) |
Calculates the photolysis rate for a given set of atmospheric conditions.
This function computes the photolysis rate based on provided climatology data and input parameters such as pressure, solar zenith angle (SZA), and ozone column. It ensures that the input parameters are within the valid range of the climatology data and interpolates the photolysis rate accordingly.
| rate | 3D array containing the photolysis rates for different combinations of pressure, SZA, and ozone column. |
| photo | Pointer to the climatology data structure containing arrays of valid pressure levels, SZAs, and ozone columns. |
| p | Pressure at which the photolysis rate is to be calculated. |
| sza | Solar zenith angle at which the photolysis rate is to be calculated. |
| o3c | Ozone column at which the photolysis rate is to be calculated. |
This function performs the following steps:
Definition at line 142 of file mptrac.c.
| double clim_tropo | ( | const clim_t * | clim, |
| const double | t, | ||
| const double | lat | ||
| ) |
Calculates the tropopause pressure based on climatological data.
This function computes the tropopause pressure using climatological data for different times and latitudes. It interpolates the tropopause pressure based on the input time and latitude parameters.
| clim | Pointer to the climatology structure containing tropopause pressure data. |
| t | Time for which the tropopause pressure is to be calculated, in seconds since the beginning of the year. |
| lat | Latitude at which the tropopause pressure is to be calculated. |
This function performs the following steps:
Definition at line 193 of file mptrac.c.
| void clim_tropo_init | ( | clim_t * | clim | ) |
Initializes the tropopause data in the climatology structure.
This function initializes the tropopause data in the climatology structure. It sets the time steps, latitudes, and tropopause pressure values based on predefined arrays.
| clim | Pointer to the climatology structure to be initialized. |
This function performs the following steps:
Definition at line 221 of file mptrac.c.
| double clim_ts | ( | const clim_ts_t * | ts, |
| const double | t | ||
| ) |
Interpolates a time series of climatological variables.
This function interpolates a time series of climatological variables based on the input time and the provided data points.
| ts | Pointer to the time series structure containing data points. |
| t | Time at which to interpolate the climatological variable (in seconds). |
This function performs linear interpolation between the closest data points to the input time t. If t is outside the range of the provided time series, the value at the nearest boundary is returned.
Definition at line 376 of file mptrac.c.
| double clim_zm | ( | const clim_zm_t * | zm, |
| const double | t, | ||
| const double | lat, | ||
| const double | p | ||
| ) |
Interpolates monthly mean zonal mean climatological variables.
This function interpolates climatological variables based on pressure, latitude, and time. The climatological data is provided in the form of monthly mean zonal mean values.
| zm | Pointer to the climatological zonal mean structure containing data points. |
| t | Time at which to interpolate the climatological variable (in seconds since the beginning of the year). |
| lat | Latitude at which to interpolate the climatological variable (in degrees). |
| p | Pressure at which to interpolate the climatological variable (in hPa). |
This function performs trilinear interpolation between the nearest data points to the input time t, latitude lat, and pressure or altitude p. If the input values are outside the range of the provided data, the function extrapolates by using the nearest boundary values.
Definition at line 394 of file mptrac.c.
| void compress_cms | ( | ctl_t * | ctl, |
| char * | varname, | ||
| float * | array, | ||
| size_t | nx, | ||
| size_t | ny, | ||
| size_t | np, | ||
| int | decompress, | ||
| FILE * | inout | ||
| ) |
Compresses or decompresses a 3D array of floats using a custom multiscale compression algorithm.
This function either compresses or decompresses a 3D array of floats based on the value of the decompress parameter. The compression and decompression are performed using a custom multiscale module.
| varname | The name of the variable being processed. |
| array | Pointer to the 3D array of floats to be compressed or decompressed. |
| nx | The number of elements in the x-dimension of the array. |
| ny | The number of elements in the y-dimension of the array. |
| np | The number of elements in the p-dimension of the array. |
| decompress | If non-zero, the function will decompress the data; otherwise, it will compress the data. |
| inout | File pointer for input or output operations. It is used for reading compressed data during decompression and writing compressed data during compression. |
The function performs the following steps:
array.The function logs the compression or decompression details and frees allocated resources before returning.
array is already allocated and can hold the decompressed data.inout is correctly opened for reading or writing as required.| void compress_pck | ( | char * | varname, |
| float * | array, | ||
| size_t | nxy, | ||
| size_t | nz, | ||
| int | decompress, | ||
| FILE * | inout | ||
| ) |
Compresses or decompresses a 3D array of floats.
This function either compresses or decompresses a 3D array of floats based on the value of the decompress parameter. Compression reduces the storage size by converting float values (4 bytes) to unsigned short values (2 bytes) with scaling and offset. Decompression restores the original float values from the compressed unsigned short representation.
| varname | The name of the variable being processed. |
| array | Pointer to the 3D array of floats to be compressed or decompressed. |
| nxy | The number of elements in the first two dimensions of the array. |
| nz | The number of elements in the third dimension of the array. |
| decompress | If non-zero, the function will decompress the data; otherwise, it will compress the data. |
| inout | File pointer for input or output operations. It is used for reading compressed data during decompression and writing compressed data during compression. |
The function performs the following steps:
array.The function allocates memory for the compressed data array and frees it before returning.
Definition at line 575 of file mptrac.c.
| void compress_zfp | ( | char * | varname, |
| float * | array, | ||
| int | nx, | ||
| int | ny, | ||
| int | nz, | ||
| int | precision, | ||
| double | tolerance, | ||
| int | decompress, | ||
| FILE * | inout | ||
| ) |
Compresses or decompresses a 3D array of floats using the ZFP library.
This function either compresses or decompresses a 3D array of floats based on the value of the decompress parameter. Compression reduces the storage size using the ZFP compression algorithm, which supports fixed-precision or fixed-accuracy modes. Decompression restores the original float values from the compressed representation.
| varname | The name of the variable being processed. |
| array | Pointer to the 3D array of floats to be compressed or decompressed. |
| nx | The number of elements in the x-dimension of the array. |
| ny | The number of elements in the y-dimension of the array. |
| nz | The number of elements in the z-dimension of the array. |
| precision | The precision parameter for ZFP compression. If greater than 0, it sets the fixed precision mode. |
| tolerance | The tolerance parameter for ZFP compression. If greater than 0 and precision is 0, it sets the fixed accuracy mode. |
| decompress | If non-zero, the function will decompress the data; otherwise, it will compress the data. |
| inout | File pointer for input or output operations. It is used for reading compressed data during decompression and writing compressed data during compression. |
The function performs the following steps:
array.The function logs the compression or decompression details and frees allocated resources before returning.
| void compress_zstd | ( | char * | varname, |
| float * | array, | ||
| size_t | n, | ||
| int | decompress, | ||
| FILE * | inout | ||
| ) |
Compresses or decompresses an array of floats using the Zstandard (ZSTD) library.
This function either compresses or decompresses an array of floats based on the value of the decompress parameter. Compression reduces the storage size using the ZSTD compression algorithm. Decompression restores the original float values from the compressed representation.
| varname | The name of the variable being processed. |
| array | Pointer to the array of floats to be compressed or decompressed. |
| n | The number of elements in the array. |
| decompress | If non-zero, the function will decompress the data; otherwise, it will compress the data. |
| inout | File pointer for input or output operations. It is used for reading compressed data during decompression and writing compressed data during compression. |
The function performs the following steps:
array.The function logs the compression or decompression details and frees allocated resources before returning.
array is already allocated and can hold the decompressed data.inout is correctly opened for reading or writing as required.| void day2doy | ( | const int | year, |
| const int | mon, | ||
| const int | day, | ||
| int * | doy | ||
| ) |
Get day of year from date.
Converts a given date to the day of the year (DOY).
This function computes the day of the year (DOY) for a given date specified by the year, month, and day. It takes into account whether the given year is a leap year or not.
| year | The year of the date. |
| mon | The month of the date (1-12). |
| day | The day of the month (1-31). |
| doy | Pointer to an integer where the computed day of the year will be stored. |
The function uses two arrays, d0 and d0l, which contain the cumulative number of days at the start of each month for non-leap years and leap years respectively. It checks if the year is a leap year and calculates the day of the year accordingly.
Definition at line 816 of file mptrac.c.
| void doy2day | ( | const int | year, |
| const int | doy, | ||
| int * | mon, | ||
| int * | day | ||
| ) |
Converts a given day of the year (DOY) to a date (month and day).
This function computes the month and day for a given day of the year (DOY) and year. It accounts for whether the given year is a leap year or not.
| year | The year corresponding to the DOY. |
| doy | The day of the year (1-365 or 1-366). |
| mon | Pointer to an integer where the computed month will be stored. |
| day | Pointer to an integer where the computed day of the month will be stored. |
The function uses two arrays, d0 and d0l, which contain the cumulative number of days at the start of each month for non-leap years and leap years respectively. It checks if the year is a leap year and calculates the month and day of the month accordingly.
Definition at line 835 of file mptrac.c.
| void fft_help | ( | double * | fcReal, |
| double * | fcImag, | ||
| int | n | ||
| ) |
Computes the Fast Fourier Transform (FFT) of a complex sequence.
This function calculates the FFT of a complex sequence represented by separate arrays for the real and imaginary parts. The input arrays fcReal and fcImag are modified in place to contain the transformed data.
| fcReal | Pointer to an array of doubles representing the real part of the input sequence. The array should have at least n elements. |
| fcImag | Pointer to an array of doubles representing the imaginary part of the input sequence. The array should have at least n elements. |
| n | The number of complex data points in the input sequence. This value should not exceed PMAX. |
fcReal and fcImag must point to arrays of at least n elements. n must be less than or equal to PMAX.fcReal and fcImag will contain the real and imaginary parts of the FFT result, respectively.n exceeds PMAX, the function will trigger an error message and terminate.Definition at line 865 of file mptrac.c.
| void geo2cart | ( | const double | z, |
| const double | lon, | ||
| const double | lat, | ||
| double * | x | ||
| ) |
Converts geographic coordinates (longitude, latitude, altitude) to Cartesian coordinates.
This function converts geographic coordinates specified by longitude, latitude, and altitude into Cartesian coordinates. The Earth is approximated as a sphere with radius defined by the constant RE.
| z | The altitude above the Earth's surface in kilometers. |
| lon | The longitude in degrees. |
| lat | The latitude in degrees. |
| x | Pointer to an array of three doubles where the computed Cartesian coordinates (x, y, z) will be stored. |
The function computes the Cartesian coordinates using the given altitude, longitude, and latitude. It assumes the Earth is a perfect sphere and uses the following formulas:
RE is defined as the Earth's radius in kilometers. Definition at line 905 of file mptrac.c.
Retrieves meteorological data for the specified time.
This function retrieves meteorological data for the given time t and updates the provided pointers to the met0 and met1 structures accordingly. It handles both the initialization and subsequent updates of the meteorological data based on the direction of time integration.
| ctl | Pointer to the control structure containing configuration settings. |
| clim | Pointer to the climate structure. |
| t | The current time for which meteorological data is to be retrieved. |
| met0 | Pointer to the pointer of the first meteorological data structure. |
| met1 | Pointer to the pointer of the second meteorological data structure. |
The function performs the following steps:
ctl, clim, met0, and met1 are properly initialized before calling this function.Definition at line 919 of file mptrac.c.
| void get_met_help | ( | ctl_t * | ctl, |
| double | t, | ||
| int | direct, | ||
| char * | metbase, | ||
| double | dt_met, | ||
| char * | filename | ||
| ) |
Helper function to generate the filename for meteorological data.
This function generates the appropriate filename for the meteorological data file based on the provided time t, direction direct, and the base filename metbase. The filename is formatted according to the specified meteorological data type.
| ctl | Pointer to the control structure containing configuration settings. |
| t | The time for which the meteorological data filename is to be generated. |
| direct | The direction of time integration (-1 for backward, 1 for forward). |
| metbase | The base string for the meteorological data filenames. |
| dt_met | The time interval between meteorological data files. |
| filename | The generated filename for the meteorological data file. |
The function performs the following steps:
dt_met based on the direction.ctl.ctl and filename are properly initialized before calling this function.Definition at line 1065 of file mptrac.c.
| void get_met_replace | ( | char * | orig, |
| char * | search, | ||
| char * | repl | ||
| ) |
Replaces occurrences of a substring in a string with another substring.
This function replaces occurrences of the substring search in the string orig with the substring repl. The replacement is performed in-place.
| orig | The original string where replacements are to be made. |
| search | The substring to be replaced. |
| repl | The substring to replace occurrences of search. |
The function iterates over the original string orig and replaces each occurrence of the substring search with the substring repl. It performs the replacement operation up to three times to ensure multiple occurrences are replaced.
YYYY, MM, and DD by year, month, and day in filenames. orig, search, and repl are properly initialized and have sufficient memory allocated before calling this function.Definition at line 1130 of file mptrac.c.
| void get_tropo | ( | int | met_tropo, |
| ctl_t * | ctl, | ||
| clim_t * | clim, | ||
| met_t * | met, | ||
| double * | lons, | ||
| int | nx, | ||
| double * | lats, | ||
| int | ny, | ||
| double * | pt, | ||
| double * | zt, | ||
| double * | tt, | ||
| double * | qt, | ||
| double * | o3t, | ||
| double * | ps, | ||
| double * | zs | ||
| ) |
Calculate tropopause data.
This function reads and interpolates various meteorological parameters such as tropopause pressure, temperature, and ozone concentration at specified latitudes and longitudes. The interpolated data is stored in the provided arrays.
| met_tropo | An integer specifying the type of meteorological data to use. |
| ctl | Pointer to a ctl_t structure that controls the meteorological data processing. |
| clim | Pointer to a clim_t structure containing climatological data. |
| met | Pointer to a met_t structure containing meteorological data. |
| lons | Array of longitudes at which to interpolate data. The array should have nx elements. |
| nx | Number of longitude points. |
| lats | Array of latitudes at which to interpolate data. The array should have ny elements. |
| ny | Number of latitude points. |
| pt | Pointer to an array where the interpolated pressure values will be stored. The array should have nx * ny elements. |
| zt | Pointer to an array where the interpolated height values will be stored. The array should have nx * ny elements. |
| tt | Pointer to an array where the interpolated temperature values will be stored. The array should have nx * ny elements. |
| qt | Pointer to an array where the interpolated specific humidity values will be stored. The array should have nx * ny elements. |
| o3t | Pointer to an array where the interpolated ozone concentration values will be stored. The array should have nx * ny elements. |
| ps | Pointer to an array where the interpolated surface pressure values will be stored. The array should have nx * ny elements. |
| zs | Pointer to an array where the interpolated surface height values will be stored. The array should have nx * ny elements. |
lons must have at least nx elements. lats must have at least ny elements. pt, zt, tt, qt, o3t, ps, and zs must have at least nx * ny elements.pt, zt, tt, qt, o3t, ps, and zs will contain the interpolated meteorological data.read_met_tropo, intpol_met_space_2d, and intpol_met_space_3d for reading and interpolating the tropopause data.Definition at line 1154 of file mptrac.c.
| void intpol_met_4d_coord | ( | met_t * | met0, |
| float | height0[EX][EY][EP], | ||
| float | array0[EX][EY][EP], | ||
| met_t * | met1, | ||
| float | height1[EX][EY][EP], | ||
| float | array1[EX][EY][EP], | ||
| double | ts, | ||
| double | height, | ||
| double | lon, | ||
| double | lat, | ||
| double * | var, | ||
| int * | ci, | ||
| double * | cw, | ||
| int | init | ||
| ) |
Interpolates meteorological variables to a given position and time.
This function interpolates meteorological variables to a specified position and time. It calculates the interpolated value based on the values provided at two time steps and performs interpolation in time, longitude, latitude, and altitude dimensions.
| met0 | Pointer to the meteorological data at the first time step. |
| height0 | Array containing heights at the first time step. |
| array0 | Array containing meteorological variable values at the first time step. |
| met1 | Pointer to the meteorological data at the second time step. |
| height1 | Array containing heights at the second time step. |
| array1 | Array containing meteorological variable values at the second time step. |
| ts | Interpolation time (fractional time between met0 and met1). |
| height | Altitude at which to interpolate. |
| lon | Longitude at which to interpolate. |
| lat | Latitude at which to interpolate. |
| var | Pointer to store the interpolated variable value. |
| ci | Array to store the calculated indices. |
| cw | Array to store the weighting factors. |
| init | Flag indicating if it's the first call (1) or not (0). |
The function first restricts the longitude within the range [0, 360) degrees. It then calculates the horizontal indices (ci[0] and ci[1]) based on the provided longitude and latitude. Next, it locates the vertical indices for each edge of the column based on the provided height.
The function then calculates the weighting factors for time, longitude, latitude, and altitude. It iterates over the interpolation process to determine the altitude weighting factor. After initializing the interpolation parameters, it calculates the interpolated variable value and stores it in the memory location pointed to by var.
height0, array0, height1, array1, ci, cw) have sufficient memory allocated before calling this function.Definition at line 1197 of file mptrac.c.
| void intpol_met_space_3d | ( | met_t * | met, |
| float | array[EX][EY][EP], | ||
| double | p, | ||
| double | lon, | ||
| double | lat, | ||
| double * | var, | ||
| int * | ci, | ||
| double * | cw, | ||
| int | init | ||
| ) |
Interpolates meteorological variables in 3D space.
This function interpolates meteorological variables at a specified pressure level and geographic position. It calculates the interpolated value based on the values provided at neighboring grid points and performs interpolation in pressure, longitude, and latitude dimensions.
| met | Pointer to the meteorological data. |
| array | Array containing meteorological variable values. |
| p | Pressure level at which to interpolate. |
| lon | Longitude at which to interpolate. |
| lat | Latitude at which to interpolate. |
| var | Pointer to store the interpolated variable value. |
| ci | Array to store the calculated indices. |
| cw | Array to store the weighting factors. |
| init | Flag indicating if it's the first call (1) or not (0). |
The function first checks the longitude and adjusts it if necessary to ensure it falls within the valid range. It then calculates the interpolation indices based on the provided pressure level, longitude, and latitude. Next, it computes the interpolation weights for pressure, longitude, and latitude.
The function interpolates vertically first and then horizontally. The interpolated value is stored in the memory location pointed to by var.
array, ci, and cw arrays have sufficient memory allocated before calling this function.Definition at line 1373 of file mptrac.c.
| void intpol_met_space_3d_ml | ( | met_t * | met, |
| float | array[EX][EY][EP], | ||
| double | p, | ||
| double | lon, | ||
| double | lat, | ||
| double * | var | ||
| ) |
Interpolates a meteorological variable in 3D space (longitude, latitude, pressure).
This function performs trilinear interpolation of a meteorological variable based on the provided longitude, latitude, and pressure coordinates. The meteorological data is given in a 3D array, and the function calculates the interpolated value and stores it in the variable pointed to by var. The function operates on model level data.
| [in] | met | Pointer to a met_t structure containing the meteorological data. |
| [in] | array | 3D array of meteorological data with dimensions [EX][EY][EP]. |
| [in] | p | Pressure coordinate at which to interpolate. |
| [in] | lon | Longitude coordinate at which to interpolate. |
| [in] | lat | Latitude coordinate at which to interpolate. |
| [out] | var | Pointer to a double where the interpolated value will be stored. |
Definition at line 1430 of file mptrac.c.
| void intpol_met_space_2d | ( | met_t * | met, |
| float | array[EX][EY], | ||
| double | lon, | ||
| double | lat, | ||
| double * | var, | ||
| int * | ci, | ||
| double * | cw, | ||
| int | init | ||
| ) |
Interpolates meteorological variables in 2D space.
This function interpolates meteorological variables at a specified geographic position. It calculates the interpolated value based on the values provided at neighboring grid points and performs interpolation in longitude and latitude dimensions.
| met | Pointer to the meteorological data. |
| array | Array containing meteorological variable values. |
| lon | Longitude at which to interpolate. |
| lat | Latitude at which to interpolate. |
| var | Pointer to store the interpolated variable value. |
| ci | Array to store the calculated indices. |
| cw | Array to store the weighting factors. |
| init | Flag indicating if it's the first call (1) or not (0). |
The function first checks the longitude and adjusts it if necessary to ensure it falls within the valid range. It then calculates the interpolation indices based on the provided longitude and latitude. Next, it computes the interpolation weights for longitude and latitude.
The function interpolates horizontally and stores the interpolated value in the memory location pointed to by var. If any of the data values used in interpolation are not finite, the function handles this situation by choosing a valid value or performing a simple interpolation.
array, ci, and cw arrays have sufficient memory allocated before calling this function.Definition at line 1500 of file mptrac.c.
| void intpol_met_time_3d | ( | met_t * | met0, |
| float | array0[EX][EY][EP], | ||
| met_t * | met1, | ||
| float | array1[EX][EY][EP], | ||
| double | ts, | ||
| double | p, | ||
| double | lon, | ||
| double | lat, | ||
| double * | var, | ||
| int * | ci, | ||
| double * | cw, | ||
| int | init | ||
| ) |
Interpolates meteorological data in 3D space and time.
This function interpolates meteorological data in three dimensions (longitude, latitude, and pressure) and time. It calculates the interpolated value based on the values provided at neighboring grid points and performs interpolation both spatially and temporally.
| met0 | Pointer to the meteorological data at time t0. |
| array0 | 3D array of meteorological data at time t0. |
| met1 | Pointer to the meteorological data at time t1. |
| array1 | 3D array of meteorological data at time t1. |
| ts | Time stamp at which to interpolate. |
| p | Pressure level at which to interpolate. |
| lon | Longitude at which to interpolate. |
| lat | Latitude at which to interpolate. |
| var | Pointer to store the interpolated value. |
| ci | Array to store the calculated indices. |
| cw | Array to store the weighting factors. |
| init | Flag indicating if it's the first call (1) or not (0). |
The function first performs spatial interpolation for both time instances (t0 and t1) using the intpol_met_space_3d function. It then calculates the weighting factor wt based on the time stamp ts. Finally, it performs temporal interpolation using the interpolated values at t0 and t1 along with the weighting factor to compute the final interpolated value stored in var.
ci and cw arrays have sufficient memory allocated before calling this function.Definition at line 1557 of file mptrac.c.
| void intpol_met_time_3d_ml | ( | met_t * | met0, |
| float | array0[EX][EY][EP], | ||
| met_t * | met1, | ||
| float | array1[EX][EY][EP], | ||
| double | ts, | ||
| double | p, | ||
| double | lon, | ||
| double | lat, | ||
| double * | var | ||
| ) |
Interpolates a meteorological variable in time and 3D space (longitude, latitude, pressure).
This function performs spatiotemporal interpolation of a meteorological variable based on the provided longitude, latitude, pressure, and timestamp. The meteorological data is given in two 3D arrays corresponding to two different time steps, and the function calculates the interpolated value and stores it in the variable pointed to by var. The function operates on model level data.
| [in] | met0 | Pointer to a met_t structure containing the meteorological data for the first time step. |
| [in] | array0 | 3D array of meteorological data for the first time step with dimensions [EX][EY][EP]. |
| [in] | met1 | Pointer to a met_t structure containing the meteorological data for the second time step. |
| [in] | array1 | 3D array of meteorological data for the second time step with dimensions [EX][EY][EP]. |
| [in] | ts | Timestamp at which to interpolate. |
| [in] | p | Pressure coordinate at which to interpolate. |
| [in] | lon | Longitude coordinate at which to interpolate. |
| [in] | lat | Latitude coordinate at which to interpolate. |
| [out] | var | Pointer to a double where the interpolated value will be stored. |
Definition at line 1586 of file mptrac.c.
| void intpol_met_time_2d | ( | met_t * | met0, |
| float | array0[EX][EY], | ||
| met_t * | met1, | ||
| float | array1[EX][EY], | ||
| double | ts, | ||
| double | lon, | ||
| double | lat, | ||
| double * | var, | ||
| int * | ci, | ||
| double * | cw, | ||
| int | init | ||
| ) |
Interpolates meteorological data in 2D space and time.
This function interpolates meteorological data in two dimensions (longitude and latitude) and time. It calculates the interpolated value based on the values provided at neighboring grid points and performs interpolation both spatially and temporally.
| met0 | Pointer to the meteorological data at time t0. |
| array0 | 2D array of meteorological data at time t0. |
| met1 | Pointer to the meteorological data at time t1. |
| array1 | 2D array of meteorological data at time t1. |
| ts | Time stamp at which to interpolate. |
| lon | Longitude at which to interpolate. |
| lat | Latitude at which to interpolate. |
| var | Pointer to store the interpolated value. |
| ci | Array to store the calculated indices. |
| cw | Array to store the weighting factors. |
| init | Flag indicating if it's the first call (1) or not (0). |
The function first performs spatial interpolation for both time instances (t0 and t1) using the intpol_met_space_2d function. It then calculates the weighting factor wt based on the time stamp ts. Finally, it performs temporal interpolation using the interpolated values at t0 and t1 along with the weighting factor to compute the final interpolated value stored in var. If one of the interpolated values is not finite, it selects the valid value based on the weighting factor wt.
ci and cw arrays have sufficient memory allocated before calling this function.Definition at line 1609 of file mptrac.c.
| void intpol_tropo_3d | ( | double | time0, |
| float | array0[EX][EY], | ||
| double | time1, | ||
| float | array1[EX][EY], | ||
| double | lons[EX], | ||
| double | lats[EY], | ||
| int | nlon, | ||
| int | nlat, | ||
| double | time, | ||
| double | lon, | ||
| double | lat, | ||
| int | method, | ||
| double * | var, | ||
| double * | sigma | ||
| ) |
Interpolates tropopause data in 3D (latitude, longitude, and time).
This function performs interpolation of tropopause data at a given latitude, longitude, and time. The interpolation can be performed using either linear interpolation or nearest neighbor interpolation. The standard deviation of the data points used in the interpolation is also computed.
| time0 | Time corresponding to the first data array array0. |
| array0 | A 2D array of tropopause data at time0. The dimensions are EX by EY. |
| time1 | Time corresponding to the second data array array1. |
| array1 | A 2D array of tropopause data at time1. The dimensions are EX by EY. |
| lons | Array of longitudes with EX elements. |
| lats | Array of latitudes with EY elements. |
| nlon | Number of longitudes. |
| nlat | Number of latitudes. |
| time | The specific time at which to interpolate the data. |
| lon | The specific longitude at which to interpolate the data. |
| lat | The specific latitude at which to interpolate the data. |
| method | Interpolation method: 1 for linear interpolation, otherwise nearest neighbor interpolation is used. |
| var | Pointer to the variable where the interpolated value will be stored. |
| sigma | Pointer to the variable where the standard deviation of the data points will be stored. |
array0 and array1 must be 2D arrays of size EX by EY. lons must have at least nlon elements and lats must have at least nlat elements.var will contain the interpolated value. sigma will contain the standard deviation of the data points used in the interpolation.lons. locate_reg, LIN, and NN for locating indices and performing interpolation.EX and EY are defined appropriately to match the dimensions of array0 and array1.Definition at line 1642 of file mptrac.c.
| void jsec2time | ( | const double | jsec, |
| int * | year, | ||
| int * | mon, | ||
| int * | day, | ||
| int * | hour, | ||
| int * | min, | ||
| int * | sec, | ||
| double * | remain | ||
| ) |
Converts Julian seconds to calendar date and time components.
This function converts Julian seconds to calendar date and time components, including year, month, day, hour, minute, and second. It also calculates the fractional part of the seconds.
| jsec | Julian seconds to convert. |
| year | Pointer to store the year. |
| mon | Pointer to store the month. |
| day | Pointer to store the day. |
| hour | Pointer to store the hour. |
| min | Pointer to store the minute. |
| sec | Pointer to store the second. |
| remain | Pointer to store the fractional part of seconds. |
The function initializes a time structure t0 with a fixed starting date and time. It then converts the Julian seconds to a time_t type by adding the seconds to the epoch time. Next, it converts the time_t value to a UTC time structure t1. Finally, it extracts the year, month, day, hour, minute, and second components from t1 and calculates the fractional part of seconds, which is stored in remain.
Definition at line 1735 of file mptrac.c.
| double kernel_weight | ( | const double | kz[EP], |
| const double | kw[EP], | ||
| const int | nk, | ||
| const double | p | ||
| ) |
Calculates the kernel weight based on altitude and given kernel data.
This function calculates the kernel weight based on altitude and given kernel data. It takes arrays of altitudes (kz) and corresponding weights (kw), the number of data points (nk), and the current altitude (p) as input.
| kz | Array of altitudes. |
| kw | Array of corresponding weights. |
| nk | Number of data points. |
| p | Current altitude. |
If the number of data points is less than 2 (nk < 2), the function returns a default weight of 1.0.
The function first computes the altitude z based on the current altitude p. Then it checks whether z is outside the range of altitudes in the kernel data. If so, it returns the corresponding weight at the nearest altitude boundary. Otherwise, it interpolates linearly between the two closest altitudes in the kernel data to determine the weight at altitude z.
Definition at line 1768 of file mptrac.c.
| double lapse_rate | ( | const double | t, |
| const double | h2o | ||
| ) |
Calculates the moist adiabatic lapse rate in Kelvin per kilometer.
This function calculates the moist adiabatic lapse rate in Kelvin per kilometer from the given temperature (t) in Kelvin and water vapor volume mixing ratio (h2o).
| t | Temperature in Kelvin. |
| h2o | Water vapor volume mixing ratio. |
The moist adiabatic lapse rate is calculated using the formula:
\[ \Gamma = \frac{{1000 \times g \times \left(a + L_v \times r \times T\right)}} {{C_{pd} \times a + L_v^2 \times r \times \epsilon}} \]
where:
The constants used in the calculation are defined externally:
Definition at line 1794 of file mptrac.c.
| void level_definitions | ( | ctl_t * | ctl | ) |
Defines pressure levels for meteorological data.
This function defines pressure levels for meteorological data based on the given control structure (ctl). Pressure levels are defined differently based on the value of press_level_def in ctl.
| ctl | Control structure containing information about pressure level definitions. |
The function determines the number of pressure levels (met_np) and the corresponding pressure values (met_p) based on the value of press_level_def in the control structure ctl. It initializes the met_np and met_p fields accordingly.
press_level_def are:press_level_def will result in an error message.Definition at line 1812 of file mptrac.c.
| int locate_irr | ( | const double * | xx, |
| const int | n, | ||
| const double | x | ||
| ) |
Locate the index of the interval containing a given value in a sorted array.
This function locates the index of the interval containing a given value in a sorted array. It uses a binary search algorithm to efficiently find the interval.
| xx | Pointer to the sorted array. |
| n | Size of the array. |
| x | Value to be located. |
x.The function assumes that the array xx is sorted in ascending order. It returns the index of the interval where the value x is located. If the value x is outside the range of the array, the function returns the index of the closest interval.
Definition at line 2019 of file mptrac.c.
| int locate_irr_float | ( | const float * | xx, |
| const int | n, | ||
| const double | x, | ||
| const int | ig | ||
| ) |
Locate the index of the interval containing a given value in an irregularly spaced array.
This function performs a binary search to locate the interval in the array xx such that xx[ig] <= x < xx[ig + 1]. If the value x lies within the interval specified by the initial guess index ig, the function returns ig. Otherwise, it searches the array to find the correct interval.
| xx | Pointer to the array of floats representing the irregularly spaced intervals. The array must be of size n. |
| n | The number of elements in the array xx. |
| x | The value to locate within the intervals of the array xx. |
| ig | The initial guess index. If the interval [xx[ig], xx[ig+1]) contains x, the function returns ig directly. |
i such that xx[i] <= x < xx[i + 1]. If x is out of bounds, it returns the index of the closest interval.xx contains at least two elements. xx.xx is not sorted in either increasing or decreasing order, or if it contains less than two elements.Definition at line 2049 of file mptrac.c.
| int locate_irr_3d | ( | float | profiles[EX][EY][EP], |
| int | np, | ||
| int | ind_lon, | ||
| int | ind_lat, | ||
| double | x | ||
| ) |
Locate the index of the interval containing a given value in a 3D irregular grid.
This function locates the index of the interval containing a given value in a 3D irregular grid. It searches for the interval in the specified longitude and latitude indices of the grid.
| profiles | 3D array representing the irregular grid. |
| np | Size of the profile (number of data points). |
| ind_lon | Index of the longitude. |
| ind_lat | Index of the latitude. |
| x | Value to be located. |
x.The function assumes that the array profiles represents a 3D irregular grid. It calculates the index of the interval where the value x is located based on the data in the specified longitude and latitude indices. If the value x is outside the range of the profile, the function returns the index of the closest interval.
Definition at line 2083 of file mptrac.c.
| int locate_reg | ( | const double * | xx, |
| const int | n, | ||
| const double | x | ||
| ) |
Locate the index of the interval containing a given value in a regular grid.
This function locates the index of the interval containing a given value in a regular grid. It calculates the index based on the spacing between grid points and the value to be located.
| xx | Pointer to the array representing the regular grid. |
| n | Size of the grid (number of grid points). |
| x | Value to be located. |
x.The function assumes that the array xx represents a regular grid with equally spaced points. It calculates the index of the interval where the value x is located based on the spacing between grid points. If the value x is outside the range of the grid, the function returns the index of the closest interval.
Definition at line 2117 of file mptrac.c.
| void locate_vert | ( | float | profiles[EX][EY][EP], |
| int | np, | ||
| int | lon_ap_ind, | ||
| int | lat_ap_ind, | ||
| double | alt_ap, | ||
| int * | ind | ||
| ) |
Locate the four vertical indizes of a box for a given height value.
This function locates the vertical indices corresponding to a given height in a 3D irregular grid. It calculates the indices based on the specified longitude and latitude indices of the grid.
| profiles | 3D array representing the irregular grid. |
| np | Size of the profile (number of data points). |
| lon_ap_ind | Index of the longitude. |
| lat_ap_ind | Index of the latitude. |
| alt_ap | Height value. |
| ind | Pointer to an array to store the resulting indices. |
The function calculates the indices corresponding to the specified height in the 3D irregular grid. It stores the resulting indices in the array pointed to by ind. The indices are calculated based on the specified longitude and latitude indices of the grid.
Definition at line 2136 of file mptrac.c.
Performs the advection of atmospheric particles using meteorological data.
This function advects particles in the atmosphere using meteorological data from two time steps, updating their positions and interpolating necessary data. It supports both pressure and zeta vertical coordinate systems.
| ctl | Pointer to the control structure containing configuration flags. |
| met0 | Pointer to the initial meteorological data structure. |
| met1 | Pointer to the final meteorological data structure. |
| atm | Pointer to the air parcel data structure. |
| dt | Array of time step values for each particle. |
The function performs the following operations:
Definition at line 2156 of file mptrac.c.
Initializes the advection module by setting up pressure fields.
This function initializes the advection module, setting up the air parcel pressure to be consistent with the given zeta vertical coordinate. It utilizes meteorological data from two time steps and interpolates the pressure values accordingly.
| ctl | Pointer to the control structure containing configuration flags. |
| met0 | Pointer to the initial meteorological data structure. |
| met1 | Pointer to the final meteorological data structure. |
| atm | Pointer to the air parcel data structure. |
The function performs the following operations:
Definition at line 2314 of file mptrac.c.
| void module_bound_cond | ( | ctl_t * | ctl, |
| clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| double * | dt | ||
| ) |
Apply boundary conditions to particles based on meteorological and climatological data.
This function applies boundary conditions to particles based on specified criteria, including latitude, pressure, surface layer parameters, and climatological data. It loops over each particle and checks whether it satisfies the specified boundary conditions. If a particle satisfies the conditions, its properties such as mass, volume mixing ratio, and age of air are updated accordingly.
It checks for quantity flags to determine which properties need to be updated. If the latitude or pressure of a particle falls outside the specified ranges, it skips the particle. It also considers surface layer parameters such as surface pressure, height, zeta range, and planetary boundary layer. If a particle is within the specified surface layer boundaries, its properties are updated accordingly.
The function updates properties such as mass and volume mixing ratio if the corresponding flags are set. It retrieves volume mixing ratio values for various trace gases (e.g., CFC-10, CFC-11, N2O, SF6) from climatological time series data and updates the particle properties accordingly. Additionally, it updates the age of air for each particle based on the current simulation time.
| ctl | Pointer to the control structure containing simulation parameters. |
| clim | Pointer to the climatological data structure containing time series data. |
| met0 | Pointer to the meteorological data structure at the initial time step. |
| met1 | Pointer to the meteorological data structure at the next time step. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| dt | Pointer to the time step value. |
Definition at line 2334 of file mptrac.c.
Calculate grid data for chemistry modules.
This function initializes and updates chemical grid quantities based on atmospheric data and interpolation of meteorological variables.
| ctl | Pointer to the control structure containing simulation parameters. |
| clim | Pointer to the climate data structure containing climatological data. |
| met0 | Pointer to the first meteorological data structure. |
| met1 | Pointer to the second meteorological data structure. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| t | Time for which chemical grid is updated. |
Definition at line 2429 of file mptrac.c.
Initializes the chemistry modules by setting atmospheric composition.
This function initializes various chemical components of the atmosphere using meteorological data and climatological information. It interpolates and sets values for water vapor (H2O), ozone (O3), and several radical species such as OH, HO2, H2O2, and O1D for each air parcel.
| ctl | Pointer to the control structure containing quantity flags. |
| clim | Pointer to the climatology structure containing climatological data. |
| met0 | Pointer to the initial meteorological data structure. |
| met1 | Pointer to the final meteorological data structure. |
| atm | Pointer to the air parcel data structure. |
The function uses OpenMP for parallel processing, iterating over each point in the atmosphere (atm->np) to initialize chemical species concentrations. It performs the following steps:
Definition at line 2579 of file mptrac.c.
| void module_convection | ( | ctl_t * | ctl, |
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| double * | dt, | ||
| double * | rs | ||
| ) |
Simulate convective processes for atmospheric particles.
This function simulates convective processes for atmospheric particles based on convective available potential energy (CAPE) and convective inhibition (CIN). It loops over each particle and checks whether the CAPE exceeds a specified threshold. If CAPE is above the threshold and meets the conditions, it performs convective mixing to update the particle's pressure based on random numbers generated for each particle.
The function interpolates CAPE and CIN from meteorological data and checks whether the particle is above the cloud top level. If the conditions are met, it determines the pressure range for vertical mixing and calculates the density range based on temperature at different pressure levels. It then calculates the new density based on random numbers and updates the particle's pressure accordingly.
| ctl | Pointer to the control structure containing simulation parameters. |
| met0 | Pointer to the meteorological data structure at the initial time step. |
| met1 | Pointer to the meteorological data structure at the next time step. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| dt | Pointer to the time step value. |
| rs | Pointer to an array of random numbers generated for each particle. |
Definition at line 2616 of file mptrac.c.
Simulate exponential decay processes for atmospheric particles.
This function simulates decay processes for atmospheric particles based on their mass or volume mixing ratio. It loops over each particle and calculates the decay rate using weighting factors for tropospheric and stratospheric lifetimes. Exponential decay is then calculated, and the mass or volume mixing ratio of particles is updated accordingly. Loss rates can also be calculated and updated based on the decay process.
The function checks for quantity flags to ensure that mass or volume mixing ratio data is available. It then calculates the weighting factor based on the particle's location in the atmosphere and sets the lifetime accordingly. Exponential decay is calculated using the time step and the lifetime, and the particle's mass or volume mixing ratio is updated. Loss rates can also be updated based on the decay process.
| ctl | Pointer to the control structure containing simulation parameters. |
| clim | Pointer to the climate data structure containing atmospheric data. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| dt | Pointer to the time step value. |
Definition at line 2688 of file mptrac.c.
| void module_diffusion_meso | ( | ctl_t * | ctl, |
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| cache_t * | cache, | ||
| double * | dt, | ||
| double * | rs | ||
| ) |
Simulate mesoscale diffusion for atmospheric particles.
This function simulates mesoscale diffusion for atmospheric particles, including horizontal and vertical wind fluctuations. It calculates standard deviations of local wind data and temporal correlations for mesoscale fluctuations. Mesoscale wind fluctuations are then calculated based on the provided random numbers and turbulence parameters. The particle positions are updated accordingly.
The function loops over each particle and calculates indices for interpolation of wind data. It then computes standard deviations of local wind data and temporal correlations for mesoscale fluctuations. Based on the turbulence parameters and provided random numbers, it calculates horizontal and vertical mesoscale wind fluctuations. Finally, it updates the particle positions based on the calculated wind fluctuations.
| ctl | Pointer to the control structure containing simulation parameters. |
| met0 | Pointer to the meteorological data structure at the current time step. |
| met1 | Pointer to the meteorological data structure at the next time step. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| cache | Pointer to the cache structure for temporary storage. |
| dt | Pointer to the time step value. |
| rs | Pointer to the array of random numbers. |
Definition at line 2727 of file mptrac.c.
Simulate turbulent diffusion for atmospheric particles.
This function simulates turbulent diffusion for atmospheric particles, including horizontal and vertical diffusion. It calculates diffusivity based on the provided weighting factor and turbulence parameters. The diffusion coefficients are then used to calculate horizontal and vertical displacements of particles based on the provided random numbers and time step.
The function loops over each particle and calculates the weighting factor based on the atmospheric properties. It then computes diffusivity for horizontal and vertical diffusion. Based on the diffusivity values and provided random numbers, it calculates horizontal and vertical displacements of particles and updates their positions accordingly.
| ctl | Pointer to the control structure containing simulation parameters. |
| clim | Pointer to the climatological data structure containing atmospheric properties. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| dt | Pointer to the time step value. |
| rs | Pointer to the array of random numbers. |
Definition at line 2807 of file mptrac.c.
Simulate dry deposition of atmospheric particles.
This function simulates the dry deposition of atmospheric particles, including both particulate matter and gases. It calculates the sedimentation velocity for particles based on the atmospheric properties and applies it to determine the loss of mass or volume mixing ratio due to deposition. The function loops over each particle and calculates the loss of mass or volume mixing ratio based on the deposition velocity and time step.
| ctl | Pointer to the control structure containing simulation parameters. |
| met0 | Pointer to the meteorological data structure at the current time step. |
| met1 | Pointer to the meteorological data structure at the next time step. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| dt | Pointer to the time step value. |
Definition at line 2847 of file mptrac.c.
| void module_h2o2_chem | ( | ctl_t * | ctl, |
| clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| double * | dt | ||
| ) |
Perform chemical reactions involving H2O2 within cloud particles.
This function simulates chemical reactions involving hydrogen peroxide (H2O2) within cloud particles. It calculates the change in H2O2 concentration over time due to chemical reactions. The reaction rates are determined based on temperature and cloud properties such as liquid water content.
| ctl | Pointer to the control structure containing simulation parameters. |
| clim | Pointer to the climatological data structure. |
| met0 | Pointer to the first meteorological data structure. |
| met1 | Pointer to the second meteorological data structure. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| dt | Pointer to an array containing the time step for each particle. |
Definition at line 2910 of file mptrac.c.
Initialize the isosurface module based on atmospheric data.
This function initializes the isosurface module based on the atmospheric data provided. It calculates the necessary variables required for generating the isosurface, such as pressure, density, or potential temperature. Additionally, it can read balloon pressure data from a file if specified in the control structure. The initialized data is stored in the cache for later use.
| ctl | Pointer to the control structure containing simulation parameters. |
| met0 | Pointer to the meteorological data structure at the current time step. |
| met1 | Pointer to the meteorological data structure at the next time step. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| cache | Pointer to the cache structure for storing initialized data. |
Definition at line 2990 of file mptrac.c.
| void module_isosurf | ( | ctl_t * | ctl, |
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| cache_t * | cache, | ||
| double * | dt | ||
| ) |
Apply the isosurface module to adjust atmospheric properties.
This function applies the isosurface module to adjust atmospheric properties based on the initialized data stored in the cache. It interpolates and restores atmospheric pressure, density, or potential temperature according to the specified isosurface mode in the control structure.
| ctl | Pointer to the control structure containing simulation parameters. |
| met0 | Pointer to the meteorological data structure at the current time step. |
| met1 | Pointer to the meteorological data structure at the next time step. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| cache | Pointer to the cache structure containing initialized data. |
| dt | Array of time step values for each particle. |
Definition at line 3054 of file mptrac.c.
| void module_kpp_chem | ( | ctl_t * | ctl, |
| clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| double * | dt | ||
| ) |
KPP chemistry module.
Simulate chemical reactions using the Kinetic PreProcessor (KPP) integration scheme.
This function simulates chemical reactions using the Kinetic PreProcessor (KPP) integration scheme for atmospheric particles. It loops over each particle in the atmospheric data structure and integrates chemical reactions over a specified time step using the KPP algorithm.
| ctl | Pointer to the control structure containing simulation parameters. |
| clim | Pointer to the climatological data structure. |
| met0 | Pointer to the first meteorological data structure. |
| met1 | Pointer to the second meteorological data structure. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| dt | Pointer to an array containing the time step for each particle. |
| void module_meteo | ( | ctl_t * | ctl, |
| clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| double * | dt | ||
| ) |
Update atmospheric properties using meteorological data.
This function updates atmospheric properties based on meteorological data interpolated between two time steps. It calculates various atmospheric quantities such as pressure, temperature, wind speed, humidity, etc., and updates the corresponding fields in the atmospheric data structure.
| ctl | Pointer to the control structure containing simulation parameters. |
| clim | Pointer to the climate data structure containing climatological data. |
| met0 | Pointer to the meteorological data structure at the current time step. |
| met1 | Pointer to the meteorological data structure at the next time step. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| dt | Array of time step values for each particle. |
Definition at line 3166 of file mptrac.c.
Update atmospheric properties through interparcel mixing.
This function updates atmospheric properties by performing interparcel mixing based on the given meteorological and climatological data. It calculates the indices of grid boxes and performs mixing for various quantities such as mass, volume mixing ratio, and other chemical species concentrations.
| ctl | Pointer to the control structure containing simulation parameters. |
| clim | Pointer to the climate data structure containing climatological data. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| t | Time at which mixing is performed. |
Definition at line 3266 of file mptrac.c.
| void module_mixing_help | ( | ctl_t * | ctl, |
| clim_t * | clim, | ||
| atm_t * | atm, | ||
| const int * | ixs, | ||
| const int * | iys, | ||
| const int * | izs, | ||
| int | qnt_idx | ||
| ) |
Perform interparcel mixing for a specific quantity.
This function performs interparcel mixing for a specific quantity based on the given indices of grid boxes. It calculates the mean concentration within each grid box and adjusts the quantity for each particle accordingly.
| ctl | Pointer to the control structure containing simulation parameters. |
| clim | Pointer to the climate data structure containing climatological data. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| ixs | Pointer to the array of grid box indices along the longitude direction. |
| iys | Pointer to the array of grid box indices along the latitude direction. |
| izs | Pointer to the array of grid box indices along the vertical direction. |
| qnt_idx | Index of the quantity for which mixing is performed. |
Definition at line 3356 of file mptrac.c.
| void module_oh_chem | ( | ctl_t * | ctl, |
| clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| double * | dt | ||
| ) |
Perform hydroxyl chemistry calculations for atmospheric particles.
This function calculates the OH chemistry for each atmospheric particle based on the specified reaction mechanism and updates the particle quantities accordingly. The OH chemistry includes bimolecular and termolecular reactions, and the reaction rates are calculated based on the provided climatological data and atmospheric conditions. The function supports both mass and volume mixing ratio quantities for the particles.
| ctl | Pointer to the control structure containing simulation parameters. |
| clim | Pointer to the climate data structure containing climatological data. |
| met0 | Pointer to the first meteorological data structure. |
| met1 | Pointer to the second meteorological data structure. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| dt | Array of time steps for each particle. |
Definition at line 3451 of file mptrac.c.
Update the positions and pressure levels of atmospheric particles.
This function updates the positions and pressure levels of atmospheric particles based on the meteorological data and the specified time step. It loops over each particle in the atmospheric data structure and performs the following operations:
| ctl | Pointer to the control structure containing simulation parameters. |
| met0 | Pointer to the first meteorological data structure. |
| met1 | Pointer to the second meteorological data structure. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| dt | Pointer to an array containing the time step for each particle. |
Definition at line 3533 of file mptrac.c.
| void module_rng_init | ( | int | ntask | ) |
Initialize random number generators for parallel tasks.
This function initializes random number generators for parallel tasks using both GSL (GNU Scientific Library) and cuRAND (NVIDIA CUDA Random Number Generation Library) if available. It sets up GSL random number generators for each OpenMP thread and initializes them with unique seeds. For cuRAND, it creates a pseudo-random number generator and sets its seed. The initialization ensures that each task or thread has its own independent random number generator to prevent interference between parallel executions.
| ntask | The number of tasks or parallel threads for which random number generators are initialized. |
Definition at line 3592 of file mptrac.c.
| void module_rng | ( | ctl_t * | ctl, |
| double * | rs, | ||
| size_t | n, | ||
| int | method | ||
| ) |
Generate random numbers using various methods and distributions.
This function generates random numbers using different methods and distributions based on the specified method and random number generator type. It supports uniform and normal distributions and can utilize GSL, Squares (Widynski, 2022), or cuRAND random number generators.
| ctl | Pointer to the control structure containing parameters and settings. |
| rs | Pointer to the array where the generated random numbers will be stored. |
| n | The number of random numbers to generate. |
| method | The method for generating random numbers:
|
Definition at line 3623 of file mptrac.c.
Simulate sedimentation of particles in the atmosphere.
This function calculates the sedimentation velocity of particles based on atmospheric pressure, temperature, and particle properties such as radius and density. It then updates the pressure of each particle based on the sedimentation velocity and the specified time step.
| ctl | Pointer to the control structure containing parameters and settings. |
| met0 | Pointer to the meteorological data at the current time step. |
| met1 | Pointer to the meteorological data at the next time step. |
| atm | Pointer to the atmospheric data containing particle information. |
| dt | Pointer to the array of time steps for each particle. |
sedi function, which takes atmospheric pressure, temperature, particle radius, and particle density as inputs. DZ2DP function.Definition at line 3727 of file mptrac.c.
Sort particles according to box index.
This function sorts particles within the atmosphere data structure based on their geographical coordinates (longitude and latitude) and pressure level. It allocates temporary arrays to store indices and auxiliary data for sorting, then performs the sorting operation. After sorting, it updates the order of particles in the atmosphere data structure.
| ctl | Pointer to the control structure containing parameters and settings. |
| met0 | Pointer to the meteorological data at the current time step. |
| atm | Pointer to the atmospheric data containing particle information. |
locate_reg and locate_irr functions to determine the appropriate index for sorting particles based on their longitude, latitude, and pressure level. Definition at line 3756 of file mptrac.c.
| void module_sort_help | ( | double * | a, |
| const int * | p, | ||
| const int | np | ||
| ) |
Reorder an array based on a given permutation.
This function reorders the elements of a given array based on a specified permutation array. It allocates temporary memory to store the reordered elements, performs the reordering operation, and then updates the original array with the reordered elements.
| a | Pointer to the array to be reordered. |
| p | Pointer to the permutation array defining the order of elements. |
| np | The number of elements in the array. |
p, which defines the new order of elements in the array a.Definition at line 3816 of file mptrac.c.
Calculate time steps for air parcels based on specified conditions.
This function calculates the time steps for air parcels based on specified conditions, including the direction of simulation, start and stop times, and a given target time. It adjusts the time step for each air parcel accordingly and checks for horizontal boundary conditions of local meteorological data.
| ctl | Pointer to the control structure containing simulation parameters. |
| met0 | Pointer to the initial meteorological data structure. |
| atm | Pointer to the atmospheric data structure containing air parcel information. |
| dt | Pointer to the array storing the calculated time steps for air parcels. |
| t | The target time for which time steps are calculated. |
t. Definition at line 3852 of file mptrac.c.
Initialize start time and time interval for time-stepping.
This function initializes the start time and time interval for time-stepping based on the direction of simulation and the provided atmospheric data. It sets the start time according to the minimum or maximum time in the atmospheric data, depending on the simulation direction. Additionally, it checks the time interval and adjusts the start time accordingly for rounding purposes.
| ctl | Pointer to the control structure containing simulation parameters. |
| atm | Pointer to the atmospheric data structure containing air parcel information. |
Definition at line 3889 of file mptrac.c.
| void module_tracer_chem | ( | ctl_t * | ctl, |
| clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| double * | dt | ||
| ) |
Simulate chemical reactions involving long-lived atmospheric tracers.
This function simulates chemical reactions involving atmospheric tracers, such as CFC-10, CFC-11, CFC-12, and N2O. It calculates the change in tracer concentrations over time based on reaction rates and environmental factors such as temperature, ozone concentration, solar zenith angle, and O(1D) volume mixing ratio.
| ctl | Pointer to the control structure containing simulation parameters. |
| clim | Pointer to the climatological data structure. |
| met0 | Pointer to the first meteorological data structure. |
| met1 | Pointer to the second meteorological data structure. |
| atm | Pointer to the atmospheric data structure containing particle information. |
| dt | Pointer to an array containing the time step for each particle. |
Definition at line 3920 of file mptrac.c.
Perform wet deposition calculations for air parcels.
This function calculates the wet deposition process for each air parcel based on provided atmospheric and meteorological data. It estimates the precipitation rate and scavenging coefficients for particles and gases inside and below cloud layers. The scavenging coefficients are used to calculate the exponential decay of mass or volume mixing ratio over time due to wet deposition.
| ctl | Pointer to the control structure containing simulation parameters. |
| met0 | Pointer to the initial meteorological data structure. |
| met1 | Pointer to the updated meteorological data structure. |
| atm | Pointer to the atmospheric data structure containing air parcel information. |
| dt | Array containing the time step for each air parcel. |
Definition at line 3988 of file mptrac.c.
| double nat_temperature | ( | const double | p, |
| const double | h2o, | ||
| const double | hno3 | ||
| ) |
Calculates the nitric acid trihydrate (NAT) temperature.
This function computes the temperature at which nitric acid trihydrate (NAT) can form given the partial pressures of water vapor and nitric acid in the atmosphere.
| p | The total atmospheric pressure (in hPa). |
| h2o | The volume mixing ratio of water vapor (H2O). |
| hno3 | The volume mixing ratio of nitric acid (HNO3). |
This function follows these steps:
The calculations are based on empirical relationships involving logarithms of the partial pressures of H2O and HNO3.
Definition at line 4123 of file mptrac.c.
Reads air parcel data from a specified file into the given atmospheric structure.
This function reads air parcel data from a file and populates the provided atm_t structure based on the type of data specified in the ctl_t control structure. It supports various data formats including ASCII, binary, netCDF, and CLaMS.
| filename | The name of the file containing the atmospheric data. |
| ctl | A pointer to the control structure (ctl_t) that specifies the type of data. |
| atm | A pointer to the atmospheric structure (atm_t) that will be populated with the data. |
This function performs the following steps:
ctl->atm_type):0 for ASCII data1 for binary data2 for netCDF data3 or 4 for CLaMS dataThe function utilizes several helper functions and macros:
SELECT_TIMER for setting the timer.LOG for logging information.ERRMSG for handling error messages.gsl_stats_minmax for calculating minimum and maximum values.Z for converting altitude.Definition at line 4147 of file mptrac.c.
Reads air parcel data from an ASCII file and populates the given atmospheric structure.
This function reads air parcel data from an ASCII file and stores the data in the provided atm_t structure. It reads each line of the file, extracts the necessary data fields, and converts the altitude to pressure.
| filename | The name of the ASCII file containing the atmospheric data. |
| ctl | A pointer to the control structure (ctl_t) that specifies the number of quantities. |
| atm | A pointer to the atmospheric structure (atm_t) that will be populated with the data. |
This function performs the following steps:
NP) and logs an error message if so.The function utilizes several macros and helper functions:
WARN for logging warnings.ERRMSG for handling error messages.TOK for tokenizing and reading values from the line.P for converting altitude to pressure.Definition at line 4218 of file mptrac.c.
Reads air parcel data from a binary file and populates the given atmospheric structure.
This function reads air parcel data from a binary file and stores the data in the provided atm_t structure. It checks the version of the binary data, reads the data values, and verifies the integrity of the data read.
| filename | The name of the binary file containing the atmospheric data. |
| ctl | A pointer to the control structure (ctl_t) that specifies the number of quantities. |
| atm | A pointer to the atmospheric structure (atm_t) that will be populated with the data. |
This function performs the following steps:
np).The function utilizes several macros and helper functions:
ERRMSG for handling error messages.FREAD for reading data from the binary file.Definition at line 4260 of file mptrac.c.
Reads air parcel data from a CLaMS netCDF file and populates the given atmospheric structure.
This function reads air parcel data from a CLaMS (Chemical Lagrangian Model of the Stratosphere) netCDF file and stores the data in the provided atm_t structure. It handles various coordinate systems and ensures the necessary data is correctly read and stored.
| filename | The name of the netCDF file containing the atmospheric data. |
| ctl | A pointer to the control structure (ctl_t) that specifies the vertical coordinate system and quantities. |
| atm | A pointer to the atmospheric structure (atm_t) that will be populated with the data. |
This function performs the following steps:
np) from the "NPARTS" dimension.TIME_INIT) if available, otherwise uses the "time" variable.vert_coord_ap):vert_coord_ap is 1, reads "ZETA" and optionally "PRESS".The function utilizes several macros and helper functions:
NC_INQ_DIM for inquiring about dimensions in the netCDF file.NC_GET_DOUBLE for reading double values from the netCDF file.NC for checking netCDF function return values.WARN for logging warnings.Definition at line 4316 of file mptrac.c.
Reads air parcel data from a generic netCDF file and populates the given atmospheric structure.
This function reads air parcel data from a netCDF file and stores the data in the provided atm_t structure. It retrieves the dimensions, geolocations (time, pressure, longitude, latitude), and specified variables from the file.
| filename | The name of the netCDF file containing the atmospheric data. |
| ctl | A pointer to the control structure (ctl_t) that specifies the number of quantities and their names. |
| atm | A pointer to the atmospheric structure (atm_t) that will be populated with the data. |
This function performs the following steps:
np) from the "obs" dimension.atm_t structure.The function utilizes several macros and helper functions:
NC_INQ_DIM for inquiring about dimensions in the netCDF file.NC_GET_DOUBLE for reading double values from the netCDF file.NC for checking netCDF function return values.Definition at line 4372 of file mptrac.c.
Reads various climatological data and populates the given climatology structure.
This function reads a range of climatological datasets based on the specified control settings and stores the data in the provided clim_t structure. It handles initialization of tropopause climatology, photolysis rates, and multiple gas species' climatologies and time series.
| ctl | A pointer to the control structure (ctl_t) that specifies file names and parameters for climatology data. |
| clim | A pointer to the climatology structure (clim_t) that will be populated with the data. |
This function performs the following steps:
ctl.ctl.ctl and applies a diurnal correction if specified.ctl.ctl.The function utilizes several helper functions:
clim_tropo_init for initializing tropopause climatology.read_clim_photo for reading photolysis rates.read_clim_zm for reading zonal mean climatologies.clim_oh_diurnal_correction for applying diurnal correction to OH climatology.read_clim_ts for reading time series data.Definition at line 4405 of file mptrac.c.
| void read_clim_photo | ( | const char * | filename, |
| clim_photo_t * | photo | ||
| ) |
Reads photolysis rates from a NetCDF file and populates the given photolysis structure.
This function opens a NetCDF file specified by the filename, reads various dimensions and data related to photolysis rates, and stores this data in the provided clim_photo_t structure. It includes checks for data consistency and logs detailed information about the loaded data.
| filename | A string containing the path to the NetCDF file containing photolysis rate data. |
| photo | A pointer to the photolysis structure (clim_photo_t) that will be populated with the data. |
The function performs the following steps:
clim_photo_t structure.Definition at line 4465 of file mptrac.c.
| int read_clim_ts | ( | const char * | filename, |
| clim_ts_t * | ts | ||
| ) |
Reads a climatological time series from a file and populates the given time series structure.
This function reads time and volume mixing ratio (VMR) data from a specified file, processes the data, and stores it in the provided clim_ts_t structure. It also includes checks for data consistency and logs detailed information about the loaded data.
| filename | A string containing the path to the file containing the climatological time series data. |
| ts | A pointer to the time series structure (clim_ts_t) that will be populated with the data. |
The function performs the following steps:
Definition at line 4600 of file mptrac.c.
| void read_clim_zm | ( | const char * | filename, |
| char * | varname, | ||
| clim_zm_t * | zm | ||
| ) |
Reads zonally averaged climatological data from a netCDF file and populates the given structure.
This function reads data from a specified netCDF file, including pressure levels, latitudes, and volume mixing ratios (VMR) for a specified variable. It performs necessary checks and logs detailed information about the loaded data.
| filename | A string containing the path to the netCDF file. |
| varname | A string containing the name of the variable to be read from the netCDF file. |
| zm | A pointer to the structure (clim_zm_t) that will be populated with the data. |
The function performs the following steps:
Definition at line 4654 of file mptrac.c.
| void read_ctl | ( | const char * | filename, |
| int | argc, | ||
| char * | argv[], | ||
| ctl_t * | ctl | ||
| ) |
Reads control parameters from a configuration file and populates the given structure.
This function reads control parameters from a specified configuration file and command line arguments, populating the provided ctl_t structure with the parsed data. It handles a wide range of parameters, performing necessary checks and providing default values where applicable.
| filename | A string containing the path to the configuration file. |
| argc | An integer representing the number of command line arguments. |
| argv | An array of strings containing the command line arguments. |
| ctl | A pointer to the structure (ctl_t) that will be populated with the control parameters. |
The function performs the following steps:
Definition at line 4753 of file mptrac.c.
| void read_kernel | ( | const char * | filename, |
| double | kz[EP], | ||
| double | kw[EP], | ||
| int * | nk | ||
| ) |
Reads kernel function data from a file and populates the provided arrays.
This function reads kernel function data from a specified file, populating the provided arrays kz and kw with the parsed data. It also updates the variable pointed to by nk with the number of data points read. The function ensures that the height levels are in ascending order and performs checks for the number of height levels read.
| filename | A string containing the path to the file containing kernel function data. |
| kz | A double array to store the height levels of the kernel function. |
| kw | A double array to store the weights corresponding to the height levels. |
| nk | A pointer to an integer variable representing the number of data points read. |
The function performs the following steps:
nk with the number of data points read.Definition at line 5534 of file mptrac.c.
Reads meteorological data from a file and populates the provided structures.
This function reads meteorological data from a specified file, populating the provided structures ctl, clim, and met with the parsed data. The function handles both netCDF and binary data formats. For netCDF data, it performs various data processing steps such as reading coordinates, surface data, atmospheric levels, and calculating derived quantities. For binary data, it reads the data directly into the appropriate arrays and structures.
| filename | A string containing the path to the file containing meteorological data. |
| ctl | A pointer to a structure containing control parameters. |
| clim | A pointer to a structure containing climatological data. |
| met | A pointer to a structure to store the meteorological data. |
The function performs the following steps:
Definition at line 5575 of file mptrac.c.
| void read_met_bin_2d | ( | FILE * | in, |
| met_t * | met, | ||
| float | var[EX][EY], | ||
| char * | varname | ||
| ) |
Reads a 2-dimensional meteorological variable from a binary file and stores it in the provided array.
This function reads a 2-dimensional meteorological variable from a binary file, which is assumed to be uncompressed, and stores it in the provided 2-dimensional array var. The variable name is used for logging purposes to identify the data being read.
| in | A pointer to the FILE structure representing the binary file to read from. |
| met | A pointer to a structure containing meteorological data. |
| var | A 2-dimensional array to store the read variable. |
| varname | A string containing the name of the variable being read. |
The function performs the following steps:
Definition at line 5822 of file mptrac.c.
| void read_met_bin_3d | ( | FILE * | in, |
| ctl_t * | ctl, | ||
| met_t * | met, | ||
| float | var[EX][EY][EP], | ||
| char * | varname, | ||
| float | bound_min, | ||
| float | bound_max | ||
| ) |
Reads 3D meteorological data from a binary file, potentially using different compression methods.
This function reads 3-dimensional meteorological data from a binary file into a specified variable array. The data can be read in uncompressed form or using one of several supported compression methods. The data is then clamped to specified minimum and maximum bounds.
| [in] | in | Pointer to the input file from which to read the data. |
| [in] | ctl | Pointer to the control structure that contains metadata about the type of data and how it is stored. |
| [in] | met | Pointer to the meteorological structure that contains the dimensions of the data. |
| [out] | var | 3D array to store the read data, with dimensions [EX][EY][EP]. |
| [in] | varname | Name of the variable being read, used for logging and debugging. |
| [in] | bound_min | Minimum bound to which data values should be clamped. |
| [in] | bound_max | Maximum bound to which data values should be clamped. |
The function supports the following types of data:
Depending on the compression type specified in the control structure, the appropriate reading and decompression function is used. The data is read into a temporary buffer, then copied into the output array, applying the specified bounds to each value.
met structure. Definition at line 5851 of file mptrac.c.
Calculates Convective Available Potential Energy (CAPE) for each grid point.
This function calculates the Convective Available Potential Energy (CAPE) at each grid point based on the provided meteorological data. CAPE is a measure of the energy available for deep convection, which is essential for severe weather development.
| clim | A pointer to a structure containing climatological data. |
| met | A pointer to a structure containing meteorological data. |
The function performs the following steps:
Definition at line 5937 of file mptrac.c.
| void read_met_cloud | ( | met_t * | met | ) |
Calculates cloud-related variables for each grid point.
This function calculates cloud-related variables, such as cloud cover, cloud top pressure, cloud bottom pressure, and total cloud water content, based on the provided meteorological data.
| met | A pointer to a structure containing meteorological data. |
The function performs the following steps:
Definition at line 6047 of file mptrac.c.
Detrends meteorological data.
This function detrends meteorological data by removing spatially varying backgrounds from each grid point. Detrending helps in removing systematic biases and trends from the data, enabling better analysis and modeling.
| ctl | A pointer to a structure containing control parameters. |
| met | A pointer to a structure containing meteorological data. |
The function performs the following steps:
Definition at line 6098 of file mptrac.c.
| void read_met_extrapolate | ( | met_t * | met | ) |
Extrapolates meteorological data.
This function extrapolates meteorological data by filling missing or invalid data points with values from the nearest valid point above. Extrapolation is performed column-wise, ensuring that missing data points are replaced with valid values.
| met | A pointer to a structure containing meteorological data. |
The function performs the following steps:
Definition at line 6202 of file mptrac.c.
Calculates geopotential heights from meteorological data.
This function calculates geopotential heights from provided meteorological data using the hydrostatic equation. Geopotential heights are computed column-wise for each grid point based on the temperature, pressure, and surface height information. Optionally, the calculated geopotential heights can be smoothed horizontally using a weighted averaging scheme.
| ctl | A pointer to a structure containing control parameters. |
| met | A pointer to a structure containing meteorological data. |
The function performs the following steps:
Definition at line 6242 of file mptrac.c.
Reads meteorological grid information from a NetCDF file.
This function reads meteorological grid information from a NetCDF file, including time, spatial dimensions, and pressure levels. It also extracts longitudes, latitudes, and pressure levels from the NetCDF file based on the specified control parameters. The function determines the time information either from the filename or from the data file, depending on the file type.
| filename | The filename of the NetCDF file. |
| ncid | The NetCDF file identifier. |
| ctl | A pointer to a structure containing control parameters. |
| met | A pointer to a structure to store meteorological data. |
The function performs the following steps:
Definition at line 6369 of file mptrac.c.
Reads meteorological variables at different vertical levels from a NetCDF file.
This function reads meteorological variables such as temperature, wind components, specific humidity, ozone data, cloud parameters, and cloud cover at various vertical levels from a NetCDF file. The function supports reading meteorological data from both MPTRAC and CLaMS formats. Depending on the file format, it reads specific variables and performs necessary conversions or interpolations.
| ncid | The NetCDF file identifier. |
| ctl | A pointer to a structure containing control parameters. |
| met | A pointer to a structure to store meteorological data. |
The function performs the following steps:
Definition at line 6487 of file mptrac.c.
Interpolates meteorological data to specified pressure levels.
This function interpolates meteorological data from model levels to pressure levels. The interpolation is performed in parallel over the spatial grid defined in the meteorological data structure.
| [in] | ctl | A pointer to a control structure containing the number of pressure levels (met_np) and the pressure levels themselves (met_p). |
| [in] | met | A pointer to a meteorological data structure containing the grid dimensions (nx, ny) and the pressure profile (pl). |
| [in,out] | var | A 3D array containing the meteorological variable to be interpolated. On output, this array will contain the interpolated values at the specified pressure levels. |
| [in] | varname | A string representing the name of the meteorological variable being interpolated. |
This function performs the following steps:
var array.Definition at line 6630 of file mptrac.c.
| void read_met_monotonize | ( | met_t * | met | ) |
Makes zeta and pressure profiles monotone.
This function ensures that zeta and pressure profiles are monotone increasing and decreasing with altitude. It iterates over each grid point and each level to identify inversions and linearly interpolate between them to maintain monotonicity. The interpolation is performed for both zeta and pressure profiles.
| met | A pointer to a structure containing meteorological data. |
The function performs the following steps:
Definition at line 6672 of file mptrac.c.
| int read_met_nc_2d | ( | int | ncid, |
| char * | varname, | ||
| char * | varname2, | ||
| char * | varname3, | ||
| char * | varname4, | ||
| ctl_t * | ctl, | ||
| met_t * | met, | ||
| float | dest[EX][EY], | ||
| float | scl, | ||
| int | init | ||
| ) |
Reads a 2-dimensional meteorological variable from a NetCDF file.
This function reads a 2-dimensional meteorological variable from a NetCDF file and stores it in a specified destination array. It supports both packed and unpacked data formats and handles missing values and scaling factors accordingly. The function also checks the meteorological data layout to ensure correct data copying.
| ncid | The NetCDF file ID. |
| varname | The name of the variable to read. |
| varname2 | An alternative name of the variable to read (in case varname is not found). |
| varname3 | An alternative name of the variable to read (in case varname2 is not found). |
| varname4 | An alternative name of the variable to read (in case varname3 is not found). |
| ctl | A pointer to a structure containing control parameters. |
| met | A pointer to a structure containing meteorological data. |
| dest | The destination array to store the read data. |
| scl | A scaling factor to apply to the read data. |
| init | Flag indicating whether to initialize the destination array before reading. |
The function performs the following steps:
Definition at line 6751 of file mptrac.c.
| int read_met_nc_3d | ( | int | ncid, |
| char * | varname, | ||
| char * | varname2, | ||
| char * | varname3, | ||
| char * | varname4, | ||
| ctl_t * | ctl, | ||
| met_t * | met, | ||
| float | dest[EX][EY][EP], | ||
| float | scl | ||
| ) |
Reads a 3-dimensional meteorological variable from a NetCDF file.
This function reads a 3-dimensional meteorological variable from a NetCDF file and stores it in a specified destination array. It supports both packed and unpacked data formats and handles missing values and scaling factors accordingly. The function also checks the meteorological data layout to ensure correct data copying.
| ncid | The NetCDF file ID. |
| varname | The name of the variable to read. |
| varname2 | An alternative name of the variable to read (in case varname is not found). |
| varname3 | An alternative name of the variable to read (in case varname2 is not found). |
| varname4 | An alternative name of the variable to read (in case varname3 is not found). |
| ctl | A pointer to a structure containing control parameters. |
| met | A pointer to a structure containing meteorological data. |
| dest | The destination array to store the read data. |
| scl | A scaling factor to apply to the read data. |
The function performs the following steps:
Definition at line 6901 of file mptrac.c.
| void read_met_pbl | ( | met_t * | met | ) |
Calculates the planetary boundary layer (PBL) height for each grid point.
This function estimates the height of the planetary boundary layer (PBL) based on various meteorological parameters. The method used is based on empirical relationships, such as those proposed by Vogelezang and Holtslag (1996) or Seidel et al. (2012). It computes the PBL height by analyzing the vertical profiles of temperature, wind speed, humidity, and pressure.
| met | A pointer to a structure containing meteorological data. |
The function performs the following steps:
Definition at line 7047 of file mptrac.c.
| void read_met_periodic | ( | met_t * | met | ) |
Applies periodic boundary conditions to meteorological data along longitudinal axis.
This function applies periodic boundary conditions to meteorological data along the longitudinal axis. It checks if the difference between the last and first longitudes and the difference between the second and first longitudes are approximately equal to 360 degrees, indicating periodicity. If the condition is met, the function increases the longitude counter, sets the longitude value for the new grid point, and copies meteorological data from the first grid point to the last grid point to ensure periodicity.
| met | A pointer to a structure containing meteorological data. |
The function performs the following steps:
Definition at line 7118 of file mptrac.c.
| void read_met_polar_winds | ( | met_t * | met | ) |
Applies a fix for polar winds in meteorological data.
This function applies a fix for polar winds in meteorological data, particularly focusing on the u and v wind components. It checks if the latitudes at the top and bottom of the grid are close to the poles. If so, it transforms the winds at 89-degree latitude into Cartesian coordinates, takes their mean, and replaces the winds at 90-degree latitude with this mean, effectively fixing the unrealistic behavior of winds at the poles.
| met | A pointer to a structure containing meteorological data. |
The function performs the following steps:
Definition at line 7173 of file mptrac.c.
| void read_met_pv | ( | met_t * | met | ) |
Calculates potential vorticity (PV) from meteorological data.
This function calculates the potential vorticity (PV) from the provided meteorological data. It employs finite difference methods to estimate gradients of temperature, wind components, and pressure in longitude, latitude, and pressure dimensions. These gradients are then used to compute PV at each grid point. Additionally, a fix for polar regions is applied to ensure smoothness of PV values in these regions.
| met | A pointer to a structure containing meteorological data. |
The function performs the following steps:
Definition at line 7232 of file mptrac.c.
| void read_met_ozone | ( | met_t * | met | ) |
Calculates the total column ozone from meteorological ozone data.
This function calculates the total column ozone from the provided meteorological ozone data. It integrates ozone concentrations over altitude to obtain the column ozone density. The result is then converted to Dobson units, which represent the thickness of the ozone layer if compressed into one layer at standard temperature and pressure.
| met | A pointer to a structure containing meteorological ozone data. |
The function performs the following steps:
Definition at line 7338 of file mptrac.c.
Downsamples meteorological data based on specified parameters.
This function downsamples meteorological data based on the provided control parameters. It reduces the resolution of meteorological data by averaging over specified intervals in longitude, latitude, and altitude.
| ctl | A pointer to a structure containing control parameters for downsampling. |
| met | A pointer to a structure containing meteorological data to be downsampled. |
The function performs the following steps:
Definition at line 7366 of file mptrac.c.
Reads surface meteorological data from a netCDF file and stores it in the meteorological data structure.
This function reads various surface meteorological variables from a netCDF file and stores them in the provided meteorological data structure. Depending on the configuration, it may read data for surface pressure, geopotential height, temperature, zonal and meridional wind components, land-sea mask, and sea surface temperature.
| ncid | NetCDF file identifier. |
| met | A pointer to the meteorological data structure to store the read data. |
| ctl | A pointer to a structure containing control parameters. |
The function performs the following steps:
Definition at line 7514 of file mptrac.c.
Calculates the tropopause and related meteorological variables based on various methods and stores the results in the meteorological data structure.
This function calculates the tropopause and related meteorological variables using different methods specified by the control parameters. The calculated tropopause pressure is stored in the provided meteorological data structure.
| ctl | A pointer to a structure containing control parameters. |
| clim | A pointer to the climatological data structure. |
| met | A pointer to the meteorological data structure to store the calculated tropopause pressure and related variables. |
The function performs the following steps:
Definition at line 7634 of file mptrac.c.
| void read_obs | ( | const char * | filename, |
| ctl_t * | ctl, | ||
| double * | rt, | ||
| double * | rz, | ||
| double * | rlon, | ||
| double * | rlat, | ||
| double * | robs, | ||
| int * | nobs | ||
| ) |
Reads observation data from a file and stores it in arrays.
This function reads observation data from a specified file in either ASCII or NetCDF format, depending on the value of the OBS_TYPE control parameter. It stores the time, altitude, longitude, latitude, and observation values in the provided arrays.
| filename | The path to the observation data file. |
| ctl | A pointer to a structure containing control parameters. |
| rt | An array to store the time values of the observations. |
| rz | An array to store the altitude values of the observations. |
| rlon | An array to store the longitude values of the observations. |
| rlat | An array to store the latitude values of the observations. |
| robs | An array to store the observation values. |
| nobs | A pointer to an integer variable to store the number of observations read. |
The function performs the following steps:
Definition at line 7806 of file mptrac.c.
| void read_obs_asc | ( | const char * | filename, |
| double * | rt, | ||
| double * | rz, | ||
| double * | rlon, | ||
| double * | rlat, | ||
| double * | robs, | ||
| int * | nobs | ||
| ) |
Reads observation data from an ASCII file.
This function reads observation data from a specified ASCII file. It extracts time, altitude, longitude, latitude, and observation values from each line of the file and stores them in the provided arrays.
| filename | The path to the ASCII file containing the observation data. |
| rt | An array to store the time values of the observations. |
| rz | An array to store the altitude values of the observations. |
| rlon | An array to store the longitude values of the observations. |
| rlat | An array to store the latitude values of the observations. |
| robs | An array to store the observation values. |
| nobs | A pointer to an integer variable to store the number of observations read. |
The function performs the following steps:
Definition at line 7850 of file mptrac.c.
| void read_obs_nc | ( | const char * | filename, |
| double * | rt, | ||
| double * | rz, | ||
| double * | rlon, | ||
| double * | rlat, | ||
| double * | robs, | ||
| int * | nobs | ||
| ) |
Reads observation data from a NetCDF file.
This function reads observation data from a specified NetCDF file. It extracts time, altitude, longitude, latitude, and observation values from the variables in the NetCDF file and stores them in the provided arrays.
| filename | The path to the NetCDF file containing the observation data. |
| rt | An array to store the time values of the observations. |
| rz | An array to store the altitude values of the observations. |
| rlon | An array to store the longitude values of the observations. |
| rlat | An array to store the latitude values of the observations. |
| robs | An array to store the observation values. |
| nobs | A pointer to an integer variable to store the number of observations read. |
The function performs the following steps:
Definition at line 7878 of file mptrac.c.
| double scan_ctl | ( | const char * | filename, |
| int | argc, | ||
| char * | argv[], | ||
| const char * | varname, | ||
| int | arridx, | ||
| const char * | defvalue, | ||
| char * | value | ||
| ) |
Scans a control file or command-line arguments for a specified variable.
This function scans either a control file or command-line arguments for a specified variable name and retrieves its value. It searches for the variable name in the control file or command-line arguments and returns its corresponding value. If the variable is not found, it returns a default value specified by the user.
| filename | The name of the control file to be scanned. If NULL, only command-line arguments will be scanned. |
| argc | The number of command-line arguments. |
| argv | An array of command-line arguments. |
| varname | The name of the variable to be searched. |
| arridx | The index of the variable array, if applicable. Set to -1 if not an array. |
| defvalue | The default value to be returned if the variable is not found. |
| value | A pointer to a character array to store the retrieved value. |
The function performs the following steps:
Definition at line 7907 of file mptrac.c.
| double sedi | ( | const double | p, |
| const double | T, | ||
| const double | rp, | ||
| const double | rhop | ||
| ) |
Calculates the sedimentation velocity of a particle in air.
This function calculates the sedimentation velocity of a particle in air using the given parameters.
| p | The atmospheric pressure [hPa]. |
| T | The temperature [K]. |
| rp | The radius of the particle [microns]. |
| rhop | The density of the particle [kg/m^3]. |
The function performs the following steps:
Definition at line 7979 of file mptrac.c.
| void spline | ( | const double * | x, |
| const double * | y, | ||
| const int | n, | ||
| const double * | x2, | ||
| double * | y2, | ||
| const int | n2, | ||
| const int | method | ||
| ) |
Performs spline interpolation or linear interpolation.
This function interpolates a set of data points using either cubic spline interpolation or linear interpolation, depending on the specified method.
| x | The array of x-coordinates of the data points. |
| y | The array of y-coordinates of the data points. |
| n | The number of data points. |
| x2 | The array of x-coordinates where interpolation is required. |
| y2 | The array to store the interpolated y-values. |
| n2 | The number of points to interpolate. |
| method | The interpolation method: 1 for cubic spline, 0 for linear interpolation. |
If the method is set to 1 (cubic spline interpolation):
If the method is set to 0 (linear interpolation):
Definition at line 8012 of file mptrac.c.
| float stddev | ( | const float * | data, |
| const int | n | ||
| ) |
Calculates the standard deviation of a set of data.
This function calculates the standard deviation of a set of floating-point data values.
| data | Pointer to the array of data values. |
| n | Number of data values in the array. |
The standard deviation is calculated using the formula:
\[ \sigma = \sqrt{\frac{\sum_{i=1}^{n} (x_i - \bar{x})^2}{n}} \]
where:
| double sza_calc | ( | const double | sec, |
| const double | lon, | ||
| const double | lat | ||
| ) |
Calculates the solar zenith angle.
This function calculates the solar zenith angle, which is the angle between the zenith (straight up) and the line connecting the observer to the center of the sun.
| sec | Seconds elapsed since 2000-01-01T12:00Z. |
| lon | Observer's longitude in degrees. |
| lat | Observer's latitude in degrees. |
The solar zenith angle is calculated based on the observer's position (longitude and latitude) and the time specified in seconds elapsed since 2000-01-01T12:00Z.
Definition at line 8082 of file mptrac.c.
| void time2jsec | ( | const int | year, |
| const int | mon, | ||
| const int | day, | ||
| const int | hour, | ||
| const int | min, | ||
| const int | sec, | ||
| const double | remain, | ||
| double * | jsec | ||
| ) |
Converts time components to seconds since January 1, 2000, 12:00:00 UTC.
This function calculates the number of seconds elapsed since January 1, 2000, 12:00:00 UTC, based on the provided year, month, day, hour, minute, and second. It also includes a fractional part to represent the remaining seconds.
| year | The year. |
| mon | The month (1-12). |
| day | The day of the month (1-31). |
| hour | The hour of the day (0-23). |
| min | The minute (0-59). |
| sec | The second (0-59). |
| remain | The fractional part of seconds. |
| jsec | Pointer to store the calculated number of seconds since January 1, 2000, 12:00:00 UTC. |
The function calculates the time elapsed since January 1, 2000, 12:00:00 UTC, up to the specified time and includes any fractional seconds indicated by the "remain" parameter.
Definition at line 8123 of file mptrac.c.
| void timer | ( | const char * | name, |
| const char * | group, | ||
| int | output | ||
| ) |
Measures and reports elapsed time for named and grouped timers.
The timer function measures elapsed time for a specified named timer and an optional group of timers, accumulating time statistics such as minimum, maximum, and mean elapsed times. It also provides an option to log the timing statistics to an output.
| name | A string representing the name of the timer. |
| group | A string representing the group to which the timer belongs. |
| output | An integer flag indicating whether to report the timing statistics (non-zero to report). |
The function keeps track of multiple timers and groups. When called, it:
output parameter is non-zero.name and group.name and group are new, and if so, initializes them.omp_get_wtime() to get the current wall time. NTIMER macro.NTIMER, the function will trigger an error message.Definition at line 8154 of file mptrac.c.
| double time_from_filename | ( | const char * | filename, |
| int | offset | ||
| ) |
Extracts and converts a timestamp from a filename to Julian seconds.
The time_from_filename function parses a given filename to extract a timestamp and converts it to Julian seconds. The timestamp is expected to follow a specific format and position within the filename, defined by the offset parameter.
| filename | A string representing the filename containing the timestamp. |
| offset | An integer indicating the position from the end of the filename where the timestamp starts. |
The function performs the following steps:
time2jsec function.YYYY-MM-DD_HH-MM (e.g., "2023-05-27_14-45").Definition at line 8222 of file mptrac.c.
| double tropo_weight | ( | const clim_t * | clim, |
| const double | t, | ||
| const double | lat, | ||
| const double | p | ||
| ) |
Computes the weighting factor for a given pressure with respect to the tropopause.
The tropo_weight function calculates a weighting factor that indicates how much a given pressure p is influenced by the tropopause based on climatological data.
| clim | A pointer to a clim_t structure containing climatological tropopause data. |
| t | The time parameter, typically representing the time of year or specific temporal context. |
| lat | The latitude for which the weighting factor is being calculated. |
| p | The pressure for which the weighting factor is to be computed. |
The function performs the following steps:
pt using the clim_tropo function.p1 (lower bound) and p0 (upper bound).p:1 if p is greater than p0 (troposphere).0 if p is less than p1 (stratosphere).1 and 0 for pressures between p0 and p1.clim_tropo function is assumed to provide the tropopause pressure based on climatological data. 0.866877899 and its reciprocal) are specific to the function's logic.Definition at line 8257 of file mptrac.c.
Writes air parcel data to a file in various formats.
The write_atm function writes the air parcel data stored in the atm structure to a file specified by filename. The format of the output file is determined by the atm_type_out field in the ctl control structure.
| filename | A string representing the name of the file to write the data to. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| atm | A pointer to an atm_t structure containing atmospheric data. |
| t | The current time, used for certain output formats. |
The function performs the following steps:
SELECT_TIMER macro.atm_type_out value in the ctl structure, writes the data in one of the following formats:atm_type_out == 0): Calls write_atm_asc.atm_type_out == 1): Calls write_atm_bin.atm_type_out == 2): Calls write_atm_nc.atm_type_out == 3): Calls write_atm_clams_traj.atm_type_out == 4): Calls write_atm_clams.atm_type_out value is not supported, triggers an error message.ctl structure.Definition at line 8281 of file mptrac.c.
Writes air parcel data to an ASCII file or gnuplot.
The write_atm_asc function writes the atmospheric data stored in the atm structure to an ASCII file specified by filename or to pipe to gnuplot if requested.
| filename | A string representing the name of the file to write the data to. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| atm | A pointer to an atm_t structure containing atmospheric data. |
| t | The current time used for filtering and timestamping. |
The function performs the following steps:
atm structure, filtering by time if specified, and writes the data to the output file.Definition at line 8341 of file mptrac.c.
Writes air parcel data to a binary file.
The write_atm_bin function writes the air parcel data stored in the atm structure to a binary file specified by filename. The function includes versioning information and ensures that all relevant data arrays are written in a consistent binary format.
| filename | A string representing the name of the file to write the data to. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| atm | A pointer to an atm_t structure containing atmospheric data. |
The function performs the following steps:
ctl structure and writes each quantity array to the file.Definition at line 8423 of file mptrac.c.
Writes air parcel data to a NetCDF file in the CLaMS format.
The write_atm_clams function creates a NetCDF file and writes air parcel data into it. The data includes time, latitude, longitude, pressure, and other specified quantities. The function defines the dimensions and variables, sets global attributes, and writes the data to the file.
| filename | A string representing the name of the output file. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| atm | A pointer to an atm_t structure containing atmospheric data. |
The function performs the following steps:
Definition at line 8473 of file mptrac.c.
Writes CLaMS trajectory data to a NetCDF file.
The write_atm_clams_traj function writes trajectory data for the CLaMS model to a NetCDF file. The file is created and populated with data including time, latitude, longitude, pressure, and other quantities. The function also handles the creation of a final initialization file at the last time step.
| dirname | A string representing the directory name where the file will be created. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| atm | A pointer to an atm_t structure containing atmospheric data. |
| t | The current time in seconds since a reference epoch. |
The function performs the following steps:
Definition at line 8523 of file mptrac.c.
Writes air parcel data to a NetCDF file.
The write_atm_nc function creates a NetCDF file and writes air parcel data into it. The data includes time, pressure, longitude, latitude, and other specified quantities. The function defines the dimensions and variables, sets global attributes, and writes the data to the file.
| filename | A string representing the name of the output file. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| atm | A pointer to an atm_t structure containing atmospheric data. |
The function performs the following steps:
Definition at line 8672 of file mptrac.c.
Writes Critical Success Index (CSI) data to a file.
The write_csi function processes air parcel and observation data to calculate and write various verification statistics, including the Critical Success Index (CSI), to a specified output file at regular intervals. The statistics include measures such as the number of hits, misses, and false alarms, bias, probability of detection, false alarm rate, equitable threat score, and correlation coefficients.
| filename | A string representing the name of the output file. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| atm | A pointer to an atm_t structure containing atmospheric data. |
| t | A double representing the current time. |
The function performs the following steps:
Definition at line 8717 of file mptrac.c.
Writes ensemble data to a file.
The write_ens function processes air parcel data to calculate ensemble means and standard deviations for various quantities and writes them to a specified output file. It handles ensemble members and calculates statistics such as means and standard deviations for each ensemble, along with latitude, longitude, altitude, and time information.
| filename | A string representing the name of the output file. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| atm | A pointer to an atm_t structure containing atmospheric data. |
| t | A double representing the current time. |
The function performs the following steps:
Definition at line 8977 of file mptrac.c.
| void write_grid | ( | const char * | filename, |
| ctl_t * | ctl, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| double | t | ||
| ) |
Writes grid data to a file in ASCII or netCDF format.
The write_grid function processes air parcel data to calculate various grid-based statistics such as column density, mean, and standard deviation for specified quantities. It then writes this data to a specified output file either in ASCII or netCDF format based on the configuration parameters provided in the ctl structure.
| filename | A string representing the name of the output file. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| met0 | A pointer to a met_t structure containing meteorological data for the initial time step. |
| met1 | A pointer to a met_t structure containing meteorological data for the final time step. |
| atm | A pointer to an atm_t structure containing atmospheric data. |
| t | A double representing the current time. |
The function performs the following steps:
grid_type in the control parameters.Definition at line 9074 of file mptrac.c.
| void write_grid_asc | ( | const char * | filename, |
| ctl_t * | ctl, | ||
| double * | cd, | ||
| double * | mean[NQ], | ||
| double * | sigma[NQ], | ||
| double * | vmr_impl, | ||
| double | t, | ||
| double * | z, | ||
| double * | lon, | ||
| double * | lat, | ||
| double * | area, | ||
| double | dz, | ||
| int * | np | ||
| ) |
Writes grid data to an ASCII file.
The write_grid_asc function writes gridded air parcel data, including column density, mean and standard deviation for specified quantities, and volume mixing ratio (if available), to an ASCII file. The function also supports writing gnuplot commands to generate plots if requested in the control parameters.
| filename | A string representing the name of the output file. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| cd | An array of doubles representing column density values. |
| mean | An array of arrays of doubles representing the mean values for specified quantities. |
| sigma | An array of arrays of doubles representing the standard deviation values for specified quantities. |
| vmr_impl | An array of doubles representing the volume mixing ratio (implicit) values. |
| t | A double representing the current time. |
| z | An array of doubles representing vertical coordinates (altitude). |
| lon | An array of doubles representing longitudinal coordinates. |
| lat | An array of doubles representing latitudinal coordinates. |
| area | An array of doubles representing surface area values. |
| dz | A double representing the layer depth. |
| np | An array of integers representing the number of particles. |
The function performs the following steps:
Definition at line 9266 of file mptrac.c.
| void write_grid_nc | ( | const char * | filename, |
| ctl_t * | ctl, | ||
| double * | cd, | ||
| double * | mean[NQ], | ||
| double * | sigma[NQ], | ||
| double * | vmr_impl, | ||
| double | t, | ||
| double * | z, | ||
| double * | lon, | ||
| double * | lat, | ||
| double * | area, | ||
| double | dz, | ||
| int * | np | ||
| ) |
Writes grid data to a NetCDF file.
The write_grid_nc function writes gridded air parcel data, including column density, mean and standard deviation for specified quantities, and volume mixing ratio (if available), to a NetCDF file. NetCDF is a self-describing, machine-independent data format for storing scientific data.
| filename | A string representing the name of the output file. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| cd | An array of doubles representing column density values. |
| mean | An array of arrays of doubles representing the mean values for specified quantities. |
| sigma | An array of arrays of doubles representing the standard deviation values for specified quantities. |
| vmr_impl | An array of doubles representing the volume mixing ratio (implicit) values. |
| t | A double representing the current time. |
| z | An array of doubles representing vertical coordinates (altitude). |
| lon | An array of doubles representing longitudinal coordinates. |
| lat | An array of doubles representing latitudinal coordinates. |
| area | An array of doubles representing surface area values. |
| dz | A double representing the layer depth. |
| np | An array of integers representing the number of particles. |
The function performs the following steps:
Definition at line 9370 of file mptrac.c.
Writes meteorological data to a binary file.
The write_met function writes meteorological data, including surface data and level data, to a binary file specified by the filename parameter. The data is written in binary format for efficient storage and transfer.
| filename | A string representing the name of the output file. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| met | A pointer to a met_t structure containing meteorological data. |
The function performs the following steps:
ctl_t structure.Definition at line 9494 of file mptrac.c.
| void write_met_bin_2d | ( | FILE * | out, |
| met_t * | met, | ||
| float | var[EX][EY], | ||
| char * | varname | ||
| ) |
Writes a 2-dimensional meteorological variable to a binary file.
The write_met_bin_2d function writes a 2-dimensional meteorological variable to a binary file specified by the out parameter. The variable data is provided in a 2-dimensional array var with maximum dimensions EX by EY. The variable name is provided as a string in the varname parameter.
| out | A pointer to a FILE structure representing the output file. |
| met | A pointer to a met_t structure containing meteorological data. |
| var | An array of floats representing the 2-dimensional variable data. |
| varname | A string containing the name of the variable being written. |
The function performs the following steps:
var to the temporary buffer help.out.Definition at line 9617 of file mptrac.c.
| void write_met_bin_3d | ( | FILE * | out, |
| ctl_t * | ctl, | ||
| met_t * | met, | ||
| float | var[EX][EY][EP], | ||
| char * | varname, | ||
| int | precision, | ||
| double | tolerance | ||
| ) |
Writes a 3-dimensional meteorological variable to a binary file.
The write_met_bin_3d function writes a 3-dimensional meteorological variable to a binary file specified by the out parameter. The variable data is provided in a 3-dimensional array var with maximum dimensions EX by EY by EP. The variable name is provided as a string in the varname parameter. Additionally, the function takes parameters for specifying the compression precision and tolerance.
| out | A pointer to a FILE structure representing the output file. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| met | A pointer to a met_t structure containing meteorological data. |
| var | An array of floats representing the 3-dimensional variable data. |
| varname | A string containing the name of the variable being written. |
| precision | An integer specifying the precision of compression (for certain compression methods). |
| tolerance | A double specifying the tolerance for compression (for certain compression methods). |
The function performs the following steps:
var to the temporary buffer help.out using the specified compression method (uncompressed, packed, zfp, zstd, cmultiscale).ctl->met_type, the function writes the variable data using different compression methods. If ctl->met_type is not supported, an error message is logged.Definition at line 9646 of file mptrac.c.
| void write_output | ( | const char * | dirname, |
| ctl_t * | ctl, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| double | t | ||
| ) |
Writes various types of output data to files in a specified directory.
The write_output function writes various types of output data to files in the directory specified by the dirname parameter. The function takes control parameters (ctl), two meteorological data structures (met0 and met1), an atmospheric data structure (atm), and a time value (t) as input.
| dirname | A string representing the directory path where output files will be written. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| met0 | A pointer to a met_t structure representing the first set of meteorological data. |
| met1 | A pointer to a met_t structure representing the second set of meteorological data. |
| atm | A pointer to an atm_t structure representing atmospheric data. |
| t | A double value representing the time at which the output is being written. |
The function performs the following steps:
t) to extract year, month, day, hour, minute, and second.Definition at line 9722 of file mptrac.c.
| void write_prof | ( | const char * | filename, |
| ctl_t * | ctl, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| double | t | ||
| ) |
Writes profile data to a specified file.
The write_prof function writes profile data to a file specified by the filename parameter. It takes control parameters (ctl), two meteorological data structures (met0 and met1), an atmospheric data structure (atm), and a time value (t) as input.
| filename | A string representing the filename where the profile data will be written. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| met0 | A pointer to a met_t structure representing the first set of meteorological data. |
| met1 | A pointer to a met_t structure representing the second set of meteorological data. |
| atm | A pointer to an atm_t structure representing atmospheric data. |
| t | A double value representing the time at which the profile data is being written. |
The function performs the following steps:
Definition at line 9817 of file mptrac.c.
| void write_sample | ( | const char * | filename, |
| ctl_t * | ctl, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| double | t | ||
| ) |
Writes sample data to a specified file.
The write_sample function writes sample data to a file specified by the filename parameter. It takes control parameters (ctl), two meteorological data structures (met0 and met1), an atmospheric data structure (atm), and a time value (t) as input.
| filename | A string representing the filename where the sample data will be written. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| met0 | A pointer to a met_t structure representing the first set of meteorological data. |
| met1 | A pointer to a met_t structure representing the second set of meteorological data. |
| atm | A pointer to an atm_t structure representing atmospheric data. |
| t | A double value representing the time at which the sample data is being written. |
The function performs the following steps:
Definition at line 10045 of file mptrac.c.
Writes station data to a specified file.
The write_station function writes station data to a file specified by the filename parameter. It takes control parameters (ctl), an atmospheric data structure (atm), and a time value (t) as input.
| filename | A string representing the filename where the station data will be written. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| atm | A pointer to an atm_t structure representing atmospheric data. |
| t | A double value representing the time at which the station data is being written. |
The function performs the following steps:
Definition at line 10208 of file mptrac.c.
Writes VTK (Visualization Toolkit) data to a specified file.
The write_vtk function writes VTK data to a file specified by the filename parameter. It takes control parameters (ctl), an atmospheric data structure (atm), and a time value (t) as input.
| filename | A string representing the filename where the VTK data will be written. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| atm | A pointer to an atm_t structure representing atmospheric data. |
| t | A double value representing the time at which the VTK data is being written. |
The function performs the following steps:
Definition at line 10294 of file mptrac.c.
1.9.4