 |
MPTRAC
|
|
MPTRAC
|
MPTRAC library declarations. More...
#include <ctype.h>#include <gsl/gsl_fft_complex.h>#include <gsl/gsl_randist.h>#include <gsl/gsl_rng.h>#include <gsl/gsl_sort.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 | particle_t |
| Particle data. More... | |
| struct | cache_t |
| Cache data structure. More... | |
| struct | depo_t |
| Ground inventories of deposited radionuclides. 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... | |
| struct | dd_t |
| Domain decomposition data structure. More... | |
Macros | |
| #define | codes_handle void* |
| Placeholder when ECCODES is not available. More... | |
| #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 | KARMAN 0.40 |
| Karman's constant. 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 | KAPPA 0.286 |
| Exponent used for potential-temperature calculations [1]. More... | |
| #define | OMEGA_EARTH 7.2921e-5 |
| Angular velocity of Earth [s^-1]. More... | |
| #define | DOBSON_UNIT 2.1415e-5 |
| Ozone column mass corresponding to one Dobson unit [kg/m^2]. More... | |
| #define | M_AIR_MOLECULE 4.8096e-26 |
| Mean mass of an air molecule [kg]. More... | |
| #define | CHEM_REF_TEMP 298.15 |
| Reference temperature for chemical equilibrium constants [K]. More... | |
| #define | H2O2_SO2_RATE_REF 9.1e7 |
| Reference rate for aqueous H2O2-SO2 chemistry [L^2/(mol^2 s)]. More... | |
| #define | H2O2_SO2_RATE_TEMP 29700. |
| Temperature parameter for the aqueous H2O2-SO2 rate [K]. More... | |
| #define | SO2_HENRY_REF 1.3e-2 |
| Henry constant of SO2 at CHEM_REF_TEMP [mol/(L atm)]. More... | |
| #define | SO2_HENRY_TEMP 2900. |
| Temperature dependence of the SO2 Henry constant [K]. More... | |
| #define | H2O2_HENRY_REF 8.3e2 |
| Henry constant of H2O2 at CHEM_REF_TEMP [mol/(L atm)]. More... | |
| #define | H2O2_HENRY_TEMP 7600. |
| Temperature dependence of the H2O2 Henry constant [K]. More... | |
| #define | SO2_CORR_A 3.12541941e-06 |
| Scale factor of the high-SO2 correction [1]. More... | |
| #define | SO2_CORR_B -5.72532259e-01 |
| Exponent of the high-SO2 correction [1]. More... | |
| #define | O1D_RATE_CCL4_A 3.30e-10 |
| O(1D) reaction pre-factor for CCl4 [cm^3/s]. More... | |
| #define | O1D_RATE_CCL4_B 0. |
| O(1D) reaction temperature parameter for CCl4 [K]. More... | |
| #define | O1D_RATE_CFC11_A 2.30e-10 |
| O(1D) reaction pre-factor for CFC-11 [cm^3/s]. More... | |
| #define | O1D_RATE_CFC11_B 0. |
| O(1D) reaction temperature parameter for CFC-11 [K]. More... | |
| #define | O1D_RATE_CFC12_A 1.40e-10 |
| O(1D) reaction pre-factor for CFC-12 [cm^3/s]. More... | |
| #define | O1D_RATE_CFC12_B -25. |
| O(1D) reaction temperature parameter for CFC-12 [K]. More... | |
| #define | O1D_RATE_N2O_A 1.19e-10 |
| O(1D) reaction pre-factor for N2O [cm^3/s]. More... | |
| #define | O1D_RATE_N2O_B -20. |
| O(1D) reaction temperature parameter for N2O [K]. More... | |
| #define | SO2_DISS_K1_REF 1.23e-2 |
| First SO2 dissociation constant at CHEM_REF_TEMP [mol/L]. More... | |
| #define | SO2_DISS_K1_TEMP 2.01e3 |
| Temperature dependence of the first SO2 dissociation constant [K]. More... | |
| #define | SO2_DISS_K2_REF 6e-8 |
| Second SO2 dissociation constant at CHEM_REF_TEMP [mol/L]. More... | |
| #define | SO2_DISS_K2_TEMP 1.12e3 |
| Temperature dependence of the second SO2 dissociation constant [K]. More... | |
| #define | WET_DEPO_T_LIQUID T0 |
| Upper temperature of the ice-cloud retention transition [K]. More... | |
| #define | WET_DEPO_T_ICE 238.15 |
| Lower temperature of the ice-cloud retention transition [K]. More... | |
| #define | WET_DEPO_T_LIQUID_BC 270. |
| Temperature threshold for below-cloud retention [K]. More... | |
| #define | RADIO_HALF_LIFE_RN222 (3.8235 * 86400.0) |
| Half-life of Rn-222 [s]. More... | |
| #define | RADIO_HALF_LIFE_PB210 (22.3 * 365.25 * 86400.0) |
| Half-life of Pb-210 [s]. More... | |
| #define | RADIO_HALF_LIFE_BE7 (53.22 * 86400.0) |
| Half-life of Be-7 [s]. More... | |
| #define | RADIO_HALF_LIFE_CS137 (30.05 * 365.25 * 86400.0) |
| Half-life of Cs-137 [s]. More... | |
| #define | RADIO_HALF_LIFE_I131 (8.02 * 86400.0) |
| Half-life of I-131 [s]. More... | |
| #define | RADIO_HALF_LIFE_XE133 (5.2474 * 86400.0) |
| Half-life of Xe-133 [s]. More... | |
| #define | RADIO_DRY_VDEP_PB210 0.002 |
| Dry deposition velocity of Pb-210 [m/s]. More... | |
| #define | RADIO_DRY_VDEP_BE7 0.001 |
| Dry deposition velocity of Be-7 [m/s]. More... | |
| #define | RADIO_DRY_VDEP_CS137 0.002 |
| Dry deposition velocity of Cs-137 [m/s]. More... | |
| #define | RADIO_DRY_VDEP_I131 0.0005 |
| Dry deposition velocity of aerosol-bound I-131 [m/s]. More... | |
| #define | RADIO_WET_COEFF_PB210 1e-4 |
| Wet deposition coefficient of Pb-210 [s^-1]. More... | |
| #define | RADIO_WET_COEFF_BE7 5e-5 |
| Wet deposition coefficient of Be-7 [s^-1]. More... | |
| #define | RADIO_WET_COEFF_CS137 1e-4 |
| Wet deposition coefficient of Cs-137 [s^-1]. More... | |
| #define | RADIO_WET_COEFF_I131 5e-5 |
| Wet deposition coefficient of aerosol-bound I-131 [s^-1]. More... | |
| #define | T0 273.15 |
| Standard temperature [K]. More... | |
| #define | EP 140 |
| Maximum number of pressure levels for meteo data. More... | |
| #define | EX 1444 |
| Maximum number of longitudes for meteo data. More... | |
| #define | EY 724 |
| Maximum number of latitudes for meteo data. More... | |
| #define | LEN 5000 |
| Maximum length of ASCII data lines. More... | |
| #define | METVAR 13 |
| Number of 3-D meteorological variables. 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 | 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 | DD_EX_GLOB (EX * 16) |
| Maximum number of longitudes of global meteo data. More... | |
| #define | DD_EY_GLOB (EY * 16) |
| Maximum number of latitudes of global meteo data. 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 | CLAMP(v, lo, hi) (((v) < (lo)) ? (lo) : (((v) > (hi)) ? (hi) : (v))) |
| Clamp a value to a specified range. More... | |
| #define | COMPRESS_SPEED(nbytes, dt) ((dt) > 0 ? ((double)(nbytes)) / ((dt) * 1024. * 1024.) : NAN) |
| Calculate compression throughput in MiB/s. More... | |
| #define | COMPRESS_RATIO(raw_size, stored_size) ((double) (raw_size) / (double) (stored_size)) |
| Calculate the compression ratio from raw and stored byte counts. More... | |
| #define | COMPRESS_BPV(n, stored_size) ((8.0 * (double) (stored_size)) / (double) (n)) |
| Calculate bits per value from stored size and element count. More... | |
| #define | DEG2DX(dlon, lat) (RE * DEG2RAD(dlon) * cos(DEG2RAD(lat))) |
| Convert a longitude difference to a distance in the x-direction (east-west) at a specific latitude. More... | |
| #define | DEG2DY(dlat) (RE * DEG2RAD(dlat)) |
| Convert a latitude difference to a distance in the y-direction (north-south). More... | |
| #define | DEG2RAD(deg) ((deg) * (M_PI / 180.0)) |
| Converts degrees to radians. 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 | DX2COORD(met, dx, lat) (((met)->coord_type == 0)? DX2DEG((dx) / 1000.0, (lat)) : ((dx) )) |
| Convert a distance in meters to a coordinate value based on grid type. More... | |
| #define | DY2COORD(met, dy) (((met)->coord_type == 0)? DY2DEG((dy) / 1000.0) : ((dy) )) |
| Convert a distance to coordinate value based on grid type. 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 | ECC(cmd) |
| Execute an ECCODES command and check for errors. More... | |
| #define | ECC_READ_2D(variable, target, scaling_factor, found_flag) |
| Writes 2-D data from a grib message into the meteo struct. More... | |
| #define | ECC_READ_3D(variable, level, target, scaling_factor, found_flag) |
| Writes 3D data from a grib message into the meteo struct. 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, level, quant) |
| Define a NetCDF variable with attributes. More... | |
| #define | NC_DEF_VAR_QUANTIZE(quant) |
| #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, check) |
| 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 | 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 | RAD2DEG(rad) ((rad) * (180.0 / M_PI)) |
| Converts radians to degrees. 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 | 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), KAPPA)) |
| 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 | USAGE |
Print usage information on -h or --help. 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) |
| Computes the value of the zeta vertical coordinate. 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) timer(id, group, 0); |
| Select and start a timer with specific attributes. More... | |
Functions | |
| 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 (const ctl_t *ctl, clim_t *clim) |
| Applies a diurnal correction to the hydroxyl radical (OH) concentration in climatology data. More... | |
| double | clim_photo (const double rate[CP][CSZA][CO3], const clim_photo_t *photo, const double p, const double sza, const 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_log_level (FILE *out, const char *codec, const char *varname, const size_t lev, const double plev, const double ratio, const double bpv, const double t_comp, const double t_decomp, const size_t n, const size_t nbytes, const float *org, const float *cmp) |
| Write one row of per-level compression diagnostics. More... | |
| void | compress_log_levels_3d (FILE *out, const char *codec, const char *varname, const met_t *met, const float *org_all, const float *cmp_all, const size_t nxy, const size_t nz, const double ratio, const double bpv, const double t_comp, const double t_decomp, const size_t nbytes) |
| Write per-level compression diagnostics for a full 3-D field. More... | |
| void | compress_cms (const ctl_t *ctl, const met_t *met, const char *varname, float *array, const int decompress, FILE *level_log, FILE *inout) |
| Compresses or decompresses a 3-D meteorological field using cmultiscale. More... | |
| void | compress_scale_to_unit (float *array, const size_t nxy, const size_t nz, double *off, double *scl) |
| Scales each vertical level of a 3-D field independently to the interval [0,1]. More... | |
| void | compress_unscale_from_unit (float *array, const size_t nxy, const size_t nz, const double *off, const double *scl) |
| Restores a levelwise [0,1]-scaled 3-D field to physical units. More... | |
| int | compress_read_lossy_scale (FILE *in, const size_t nz, double **off, double **scl) |
| Read optional lossyscaling metadata for a 3-D field. More... | |
| void | compress_write_lossy_scale (FILE *out, const int enabled, float *array, const size_t nxy, const size_t nz, double **off, double **scl) |
| Write optional lossyscaling metadata for a 3-D field. More... | |
| int | compress_metvar_index (const char *varname) |
| Maps a meteorological variable name to its internal MPTRAC variable index. More... | |
| void | compress_pck (const ctl_t *ctl, const met_t *met, const char *varname, float *array, const int decompress, FILE *level_log, FILE *inout) |
| Compresses or decompresses a 3‑D float array using the PCK format. More... | |
| void | compress_sz3 (const ctl_t *ctl, const met_t *met, const char *varname, float *array, const int decompress, FILE *level_log, FILE *inout) |
| Compresses or decompresses a 3-D float array using the SZ3 library. More... | |
| void | compress_zfp (const ctl_t *ctl, const met_t *met, const char *varname, float *array, const int decompress, FILE *level_log, FILE *inout) |
| Compresses or decompresses a 3D array of floats using the ZFP library. More... | |
| void | compress_zstd (const ctl_t *ctl, const met_t *met, const char *varname, float *array, const int decompress, FILE *level_log, FILE *inout) |
| Compresses or decompresses a float array using ZSTD. More... | |
| void | compress_lz4 (const ctl_t *ctl, const met_t *met, const char *varname, float *array, const int decompress, FILE *level_log, FILE *inout) |
| Compresses or decompresses a float array using LZ4. More... | |
| double | cos_sza (const double sec, const double lon, const double lat) |
| Calculates the cosine of the solar zenith angle. More... | |
| void | day2doy (const int year, const int mon, const int day, int *doy) |
| Get day of year from date. More... | |
| void | dd_assign_subdomains (const ctl_t *ctl, const dd_t *dd, atm_t *atm, const int init) |
| Assign or update particle subdomain ownership. More... | |
| void | dd_atm2particles (const ctl_t *ctl, cache_t *cache, atm_t *atm, particle_t *particles, const int npart) |
| Copy migratable atmospheric particles from the ATM state into a particle buffer. More... | |
| int | dd_calc_subdomain_from_coords (const ctl_t *ctl, const dd_t *dd, const double lon, const double lat) |
| Determine MPI subdomain from particle coordinates. More... | |
| void | dd_communicate_particles (const ctl_t *ctl, const dd_t *dd, particle_t **particles, int *npart, int *capacity) |
| Exchange particles between MPI ranks according to their destination rank. More... | |
| void | dd_init (const ctl_t *ctl, dd_t *dd, atm_t *atm) |
| Initialize the domain decomposition infrastructure. More... | |
| void | dd_normalize_lon_lat (const dd_t *dd, double *lon, double *lat) |
| Normalize geographic coordinates to the global grid convention. More... | |
| void | dd_particles2atm (const ctl_t *ctl, cache_t *cache, const particle_t *particles, const int npart, atm_t *atm) |
| Copy received particles from the communication buffer into the atmospheric state. More... | |
| void | dd_push (const ctl_t *ctl, atm_t *atm, cache_t *cache, int *npart) |
| Compact and stage non-local particles for domain decomposition transfer. More... | |
| void | dd_read_met_nc_grid (dd_t *dd, const ctl_t *ctl, met_t *met, const int ncid) |
| Read meteorological grid information and construct the domain-decomposed grid with halo regions. More... | |
| void | dd_sort (const ctl_t *ctl, const met_t *met0, atm_t *atm, dd_t *dd, int *npart) |
| Sort local atmospheric particles and determine export counts for domain decomposition. More... | |
| void | dd_sort_help (double *a, dd_t *dd, const int np) |
| Apply the sorting permutation to a particle data array. 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, const 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_filename (const ctl_t *ctl, const double t, const int direct, const char *metbase, const double dt_met, char *filename) |
| Generates a formatted filename for meteorological data files based on the input parameters. More... | |
| void | get_met_replace (char *orig, const char *search, const char *repl) |
| Replaces occurrences of a substring in a string with another substring. More... | |
| void | get_tropo (const int met_tropo, ctl_t *ctl, const clim_t *clim, met_t *met, const double *lons, const int nx, const double *lats, const int ny, double *pt, double *zt, double *tt, double *qt, double *o3t, double *ps, double *zs) |
| Calculate tropopause data. More... | |
| void | intpol_check_lon_lat (const double *lons, const int nlon, const double *lats, const int nlat, const double lon, const double lat, double *lon2, double *lat2) |
| Adjusts longitude and latitude to ensure they fall within valid bounds. More... | |
| void | intpol_check_cartesian (const double *lons, const int nlon, const double *lats, const int nlat, const double lon, const double lat, double *lon2, double *lat2) |
| Clamps UTM coordinates to the valid bounds. More... | |
| void | intpol_met_4d_zeta (const met_t *met0, float height0[EX][EY][EP], float array0[EX][EY][EP], const met_t *met1, float height1[EX][EY][EP], float array1[EX][EY][EP], const double ts, const double height, const double lon, const double lat, double *var, int *ci, double *cw, const int init) |
| Interpolates meteorological variables to a given position and time. More... | |
| void | intpol_met_space_3d (const met_t *met, float array[EX][EY][EP], const double p, const double lon, const double lat, double *var, int *ci, double *cw, const int init) |
| Interpolates meteorological variables in 3D space. More... | |
| void | intpol_met_space_2d (const met_t *met, float array[EX][EY], const double lon, const double lat, double *var, int *ci, double *cw, const int init) |
| Interpolates meteorological variables in 2D space. More... | |
| void | intpol_met_time_3d (const met_t *met0, float array0[EX][EY][EP], const met_t *met1, float array1[EX][EY][EP], const double ts, const double p, const double lon, const double lat, double *var, int *ci, double *cw, const int init) |
| Interpolates meteorological data in 3D space and time. More... | |
| void | intpol_met_time_2d (const met_t *met0, float array0[EX][EY], const met_t *met1, float array1[EX][EY], const double ts, const double lon, const double lat, double *var, int *ci, double *cw, const int init) |
| Interpolates meteorological data in 2D space and time. More... | |
| void | intpol_tropo_3d (const double time0, float array0[EX][EY], const double time1, float array1[EX][EY], const double lons[EX], const double lats[EY], const int nlon, const int nlat, const double time, const double lon, const double lat, const 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_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], const int np, const int lon_ap_ind, const int lat_ap_ind, const double alt_ap, int *ind) |
| Locate the four vertical indizes of a box for a given height value. More... | |
| void | module_advect (const ctl_t *ctl, const cache_t *cache, met_t *met0, met_t *met1, atm_t *atm) |
| Advances particle positions using different advection schemes. More... | |
| void | module_advect_init (const ctl_t *ctl, const cache_t *cache, met_t *met0, met_t *met1, atm_t *atm) |
| Initializes the advection module by setting up pressure fields. More... | |
| void | module_bound_cond (const ctl_t *ctl, const cache_t *cache, const clim_t *clim, met_t *met0, met_t *met1, atm_t *atm) |
| Apply boundary conditions to particles based on meteorological and climatological data. More... | |
| void | module_chem_grid (const ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, const double t) |
| Computes gridded chemical tracer concentrations (volume mixing ratio) from individual air parcel mass data and assigns them back to the parcels. More... | |
| void | module_chem_init (const ctl_t *ctl, const cache_t *cache, const clim_t *clim, met_t *met0, met_t *met1, atm_t *atm) |
| Initializes the chemistry modules by setting atmospheric composition. More... | |
| void | module_convection (const ctl_t *ctl, cache_t *cache, met_t *met0, met_t *met1, atm_t *atm) |
| Performs convective mixing of atmospheric particles. More... | |
| void | module_dd (double t, const ctl_t *ctl, cache_t *cache, dd_t *dd, atm_t *atm, met_t **met) |
| Perform domain decomposition and exchange particles between MPI ranks. More... | |
| void | module_decay (const ctl_t *ctl, const cache_t *cache, const clim_t *clim, atm_t *atm) |
| Simulate exponential decay processes for atmospheric particles. More... | |
| void | module_diff_meso (const ctl_t *ctl, cache_t *cache, const met_t *met0, const met_t *met1, atm_t *atm) |
| Simulate mesoscale diffusion for atmospheric particles. More... | |
| void | module_diff_pbl (const ctl_t *ctl, cache_t *cache, met_t *met0, met_t *met1, atm_t *atm) |
| Computes particle diffusion within the planetary boundary layer (PBL). More... | |
| void | module_diff_turb (const ctl_t *ctl, cache_t *cache, const clim_t *clim, met_t *met0, met_t *met1, atm_t *atm) |
| Applies turbulent diffusion processes to atmospheric particles. More... | |
| void | module_dry_depo (const ctl_t *ctl, const cache_t *cache, met_t *met0, met_t *met1, atm_t *atm) |
| Simulate dry deposition of atmospheric particles. More... | |
| void | module_h2o2_chem (const ctl_t *ctl, const cache_t *cache, const clim_t *clim, met_t *met0, met_t *met1, atm_t *atm) |
| Perform chemical reactions involving H2O2 within cloud particles. More... | |
| void | module_isosurf_init (const ctl_t *ctl, cache_t *cache, met_t *met0, met_t *met1, atm_t *atm) |
| Initialize the isosurface module based on atmospheric data. More... | |
| void | module_isosurf (const ctl_t *ctl, const cache_t *cache, met_t *met0, met_t *met1, atm_t *atm) |
| Apply the isosurface module to adjust atmospheric properties. More... | |
| void | module_kpp_chem (ctl_t *ctl, cache_t *cache, clim_t *clim, met_t *met0, met_t *met1, atm_t *atm) |
| KPP chemistry module. More... | |
| void | module_meteo (const ctl_t *ctl, const cache_t *cache, const clim_t *clim, met_t *met0, met_t *met1, atm_t *atm) |
| Update atmospheric properties using meteorological data. More... | |
| void | module_mixing (const ctl_t *ctl, const clim_t *clim, atm_t *atm, const double t) |
| Update atmospheric properties through interparcel mixing. More... | |
| void | module_mixing_help (const ctl_t *ctl, atm_t *atm, const int *ixs, const int *iys, const int *izs, const double *mixparam, const int qnt_idx, const int use_ensemble) |
| Perform subgrid-scale interparcel mixing of a given quantity. More... | |
| void | module_oh_chem (const ctl_t *ctl, const cache_t *cache, const clim_t *clim, met_t *met0, met_t *met1, atm_t *atm) |
| Perform hydroxyl chemistry calculations for atmospheric particles. More... | |
| void | module_position (const cache_t *cache, met_t *met0, met_t *met1, atm_t *atm, const int reflect) |
| Update the positions and pressure levels of atmospheric particles. More... | |
| void | module_rng_init (const int ntask) |
| Initialize random number generators for parallel tasks. More... | |
| void | module_rng (const ctl_t *ctl, double *rs, const size_t n, const int method) |
| Generate random numbers using various methods and distributions. More... | |
| void | module_radio_decay (const ctl_t *ctl, const cache_t *cache, atm_t *atm) |
| Apply radioactive decay to atmospheric tracer species. More... | |
| void | module_radio_depo (const ctl_t *ctl, const cache_t *cache, met_t *met0, met_t *met1, atm_t *atm, depo_t *depo) |
| Deposit supported radionuclides from air parcels onto the ground grid. More... | |
| void | module_sedi (const ctl_t *ctl, const cache_t *cache, met_t *met0, met_t *met1, atm_t *atm) |
| Simulate sedimentation of particles in the atmosphere. More... | |
| void | module_sort (const ctl_t *ctl, const 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 (const ctl_t *ctl, cache_t *cache, const met_t *met0, atm_t *atm, const double t) |
| Calculate time steps for air parcels based on specified conditions. More... | |
| void | module_timesteps_init (ctl_t *ctl, const atm_t *atm) |
| Initialize start time and time interval for time-stepping. More... | |
| void | module_tracer_chem (const ctl_t *ctl, const cache_t *cache, const clim_t *clim, met_t *met0, met_t *met1, atm_t *atm) |
| Simulate chemical reactions involving long-lived atmospheric tracers. More... | |
| void | module_wet_depo (const ctl_t *ctl, const cache_t *cache, met_t *met0, met_t *met1, atm_t *atm) |
| Perform wet deposition calculations for air parcels. More... | |
| void | mptrac_alloc (ctl_t **ctl, cache_t **cache, clim_t **clim, met_t **met0, met_t **met1, atm_t **atm, depo_t **depo, dd_t **dd) |
| Allocates and initializes memory resources for MPTRAC. More... | |
| void | mptrac_free (ctl_t *ctl, cache_t *cache, clim_t *clim, met_t *met0, met_t *met1, atm_t *atm, depo_t *depo, dd_t *dd) |
| Frees memory resources allocated for MPTRAC. More... | |
| void | mptrac_get_met (ctl_t *ctl, clim_t *clim, const double t, met_t **met0, met_t **met1, dd_t *dd) |
| Retrieves meteorological data for the specified time. More... | |
| void | mptrac_init (ctl_t *ctl, cache_t *cache, clim_t *clim, atm_t *atm, depo_t *depo, const int ntask) |
| Initializes the MPTRAC model and its associated components. More... | |
| int | mptrac_read_atm (const char *filename, const ctl_t *ctl, atm_t *atm) |
| Reads air parcel data from a specified file into the given atmospheric structure. More... | |
| void | mptrac_read_clim (const ctl_t *ctl, clim_t *clim) |
| Reads various climatological data and populates the given climatology structure. More... | |
| void | mptrac_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... | |
| int | mptrac_read_met (const char *filename, const ctl_t *ctl, const clim_t *clim, met_t *met, dd_t *dd) |
| Reads meteorological data from a file, supporting multiple formats and MPI broadcasting. More... | |
| void | mptrac_run_timestep (ctl_t *ctl, cache_t *cache, clim_t *clim, met_t **met0, met_t **met1, atm_t *atm, depo_t *depo, double t, dd_t *dd) |
| Executes a single timestep of the MPTRAC model simulation. More... | |
| void | mptrac_write_atm (const char *filename, const ctl_t *ctl, const atm_t *atm, const double t) |
| Writes air parcel data to a file in various formats. More... | |
| void | mptrac_write_met (const char *filename, const ctl_t *ctl, met_t *met) |
| Writes meteorological data to a file, supporting multiple formats and compression options. More... | |
| void | mptrac_write_output (const char *dirname, const ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, depo_t *depo, const double t) |
| Writes various types of output data to files in a specified directory. More... | |
| void | mptrac_update_device (const ctl_t *ctl, const cache_t *cache, const clim_t *clim, met_t **met0, met_t **met1, const atm_t *atm) |
| Updates device memory for specified data structures. More... | |
| void | mptrac_update_host (const ctl_t *ctl, const cache_t *cache, const clim_t *clim, met_t **met0, met_t **met1, const atm_t *atm) |
| Updates host memory for specified data structures. More... | |
| double | nat_temperature (const double p, const double h2o, const double hno3) |
| Calculates the nitric acid trihydrate (NAT) temperature. More... | |
| double | pbl_weight (const ctl_t *ctl, const atm_t *atm, const int ip, const double pbl, const double ps) |
| Computes a weighting factor based on planetary boundary layer pressure. More... | |
| int | read_atm_asc (const char *filename, const 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, const 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, const ctl_t *ctl, atm_t *atm) |
| Reads atmospheric data from a CLAMS NetCDF file. More... | |
| int | read_atm_nc (const char *filename, const 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_photo (const char *filename, clim_photo_t *photo) |
| Reads photolysis rates from a NetCDF file and populates the given photolysis structure. More... | |
| void | read_clim_photo_help (const int ncid, const char *varname, const clim_photo_t *photo, double var[CP][CSZA][CO3]) |
| Reads a 3D climatological photochemistry variable from a NetCDF file. 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, const char *varname, clim_zm_t *zm) |
| Reads zonally averaged climatological data from a netCDF 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_bin (const char *filename, const ctl_t *ctl, met_t *met) |
| Reads meteorological data from a binary file. More... | |
| void | read_met_bin_2d (FILE *in, const met_t *met, float var[EX][EY], const 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, const ctl_t *ctl, const met_t *met, float var[EX][EY][EP], const char *varname, const float bound_min, const float bound_max) |
| Reads 3D meteorological data from a binary file, potentially using different compression methods. More... | |
| void | read_met_cape (const ctl_t *ctl, const 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 (const 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 (const ctl_t *ctl, met_t *met) |
| Calculates geopotential heights from meteorological data. More... | |
| int | read_met_grib (const char *filename, const ctl_t *ctl, met_t *met) |
| Reads meteorological data from a grib file and processes it. More... | |
| void | read_met_grib_grid (codes_handle **handles, int count_handles, met_t *met) |
| Reads global meteorological information from a grib file. More... | |
| void | read_met_grib_levels (codes_handle **handles, const int num_messages, const ctl_t *ctl, met_t *met) |
| Reads meteorological variables at different vertical levels from a grib file. More... | |
| void | read_met_grib_surface (codes_handle **handles, const int num_messages, const ctl_t *ctl, met_t *met) |
| Reads surface meteorological data from a grib file and stores it in the meteorological data structure. More... | |
| void | read_met_ml2pl (const ctl_t *ctl, const met_t *met, float var[EX][EY][EP], const char *varname) |
| Interpolates meteorological data to specified pressure levels. More... | |
| void | read_met_monotonize (const ctl_t *ctl, met_t *met) |
| Makes zeta and pressure profiles monotone. More... | |
| int | read_met_nc (const char *filename, const ctl_t *ctl, met_t *met, dd_t *dd) |
| Reads meteorological data from a NetCDF file and processes it. More... | |
| void | read_met_nc_grid (const char *filename, const int ncid, const ctl_t *ctl, met_t *met, dd_t *dd) |
| Reads meteorological grid data from NetCDF files with domain decomposition. More... | |
| void | read_met_nc_levels (const int ncid, const ctl_t *ctl, met_t *met, dd_t *dd) |
| Reads and processes meteorological level data from NetCDF files with domain decomposition. More... | |
| void | read_met_nc_surface (const int ncid, const ctl_t *ctl, met_t *met, dd_t *dd) |
| Reads and processes surface meteorological data from NetCDF files with domain decomposition. More... | |
| int | read_met_nc_2d (const int ncid, const char *const varnames[], const size_t nvarnames, const ctl_t *ctl, const met_t *met, dd_t *dd, float dest[EX][EY], const float scl, const int init) |
| Reads a 2-dimensional meteorological variable from a NetCDF file. More... | |
| int | read_met_nc_3d (const int ncid, const char *const varnames[], const size_t nvarnames, const ctl_t *ctl, const met_t *met, dd_t *dd, float dest[EX][EY][EP], const float scl) |
| Reads a 3-dimensional meteorological variable from a NetCDF file. More... | |
| void | read_met_pbl (const ctl_t *ctl, met_t *met) |
| Computes the planetary boundary layer (PBL) pressure based on meteorological data. 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 (const ctl_t *ctl, met_t *met) |
| Downsamples meteorological data based on specified parameters. More... | |
| void | read_met_tropo (const ctl_t *ctl, const 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, const 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, const 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... | |
| 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, const int output) |
| Measures and reports elapsed time for named and grouped timers. More... | |
| double | time_from_filename (const char *filename, const int offset, const int with_seconds) |
| Extracts and converts a timestamp from a filename to Julian seconds. More... | |
| double | tropo_weight (const ctl_t *ctl, const clim_t *clim, const atm_t *atm, const int ip) |
| Computes a weighting factor based on tropopause pressure. More... | |
| void | write_atm_asc (const char *filename, const ctl_t *ctl, const atm_t *atm, const double t) |
| Writes air parcel data to an ASCII file or gnuplot. More... | |
| void | write_atm_bin (const char *filename, const ctl_t *ctl, const atm_t *atm) |
| Writes air parcel data to a binary file. More... | |
| void | write_atm_clams (const char *filename, const ctl_t *ctl, const atm_t *atm) |
| Writes air parcel data to a NetCDF file in the CLaMS format. More... | |
| void | write_atm_clams_traj (const char *dirname, const ctl_t *ctl, const atm_t *atm, const double t) |
| Writes CLaMS trajectory data to a NetCDF file. More... | |
| void | write_atm_nc (const char *filename, const ctl_t *ctl, const atm_t *atm) |
| Writes air parcel data to a NetCDF file. More... | |
| void | write_budget (const char *filename, const ctl_t *ctl, const atm_t *atm, const double t) |
| Writes a time series of the atmospheric tracer mass budget. More... | |
| void | write_csi (const char *filename, const ctl_t *ctl, const atm_t *atm, const double t) |
| Writes Critical Success Index (CSI) data to a file. More... | |
| void | write_csi_ens (const char *filename, const ctl_t *ctl, const atm_t *atm, const double t) |
| Writes ensemble-based Critical Success Index (CSI) and other verification statistics to an output file. More... | |
| void | write_ens (const char *filename, const ctl_t *ctl, const atm_t *atm, const double t) |
| Writes ensemble data to a file. More... | |
| void | write_depo (const char *filename, const ctl_t *ctl, const depo_t *depo, const double t) |
| Convert cumulative ground inventories to Bq m^-2 and write them. More... | |
| void | write_depo_asc (const char *filename, const ctl_t *ctl, const double *data, const double t, const double *lon, const double *lat, const double *area) |
| Write radioactive deposition densities as a gnuplot-compatible table. More... | |
| void | write_depo_nc (const char *filename, const ctl_t *ctl, const double *data, const double t, const double *lon, const double *lat, const double *area) |
| Write radioactive deposition densities as a CF-style netCDF file. More... | |
| void | write_grid (const char *filename, const ctl_t *ctl, met_t *met0, met_t *met1, const atm_t *atm, const double t) |
| Writes grid data to a file in ASCII or netCDF format. More... | |
| void | write_grid_asc (const char *filename, const ctl_t *ctl, const double *cd, double *mean[NQ], double *sigma[NQ], const double *vmr_impl, const double t, const double *z, const double *lon, const double *lat, const double *area, const double dz, const int *np) |
| Writes grid data to an ASCII file. More... | |
| void | write_grid_nc (const char *filename, const ctl_t *ctl, const double *cd, double *mean[NQ], double *sigma[NQ], const double *vmr_impl, const double t, const double *z, const double *lon, const double *lat, const double *area, const double dz, const int *np) |
| Writes grid data to a NetCDF file. More... | |
| void | write_met_bin (const char *filename, const ctl_t *ctl, met_t *met) |
| Writes meteorological data in binary format to a specified file. More... | |
| void | write_met_bin_2d (FILE *out, met_t *met, float var[EX][EY], const char *varname) |
| Writes a 2-dimensional meteorological variable to a binary file. More... | |
| void | write_met_bin_3d (FILE *out, const ctl_t *ctl, met_t *met, float var[EX][EY][EP], const char *varname, const int metvar, FILE *level_log) |
| Writes a 3-dimensional meteorological variable to a binary file. More... | |
| void | write_met_nc (const char *filename, const ctl_t *ctl, met_t *met) |
| Writes meteorological data to a NetCDF file. More... | |
| void | write_met_nc_2d (const int ncid, const char *varname, met_t *met, float var[EX][EY], const float scl) |
| Writes a 2D meteorological variable to a NetCDF file. More... | |
| void | write_met_nc_3d (const int ncid, const char *varname, met_t *met, float var[EX][EY][EP], const float scl) |
| Writes a 3D meteorological variable to a NetCDF file. More... | |
| void | write_prof (const char *filename, const ctl_t *ctl, met_t *met0, met_t *met1, const atm_t *atm, const double t) |
| Writes profile data to a specified file. More... | |
| void | write_sample (const char *filename, const ctl_t *ctl, met_t *met0, met_t *met1, const atm_t *atm, const double t) |
| Writes sample data to a specified file. More... | |
| void | write_station (const char *filename, const ctl_t *ctl, atm_t *atm, const double t) |
| Writes station data to a specified file. More... | |
| void | write_vtk (const char *filename, const ctl_t *ctl, const atm_t *atm, const double t) |
| Writes VTK (Visualization Toolkit) data to a specified file. More... | |
MPTRAC library declarations.
Definition in file mptrac.h.
| #define codes_handle void* |
| #define CPD 1003.5 |
| #define LV 2501000. |
| #define KB 1.3806504e-23 |
| #define KAPPA 0.286 |
| #define OMEGA_EARTH 7.2921e-5 |
| #define DOBSON_UNIT 2.1415e-5 |
| #define M_AIR_MOLECULE 4.8096e-26 |
| #define CHEM_REF_TEMP 298.15 |
| #define H2O2_SO2_RATE_REF 9.1e7 |
| #define H2O2_SO2_RATE_TEMP 29700. |
| #define SO2_HENRY_REF 1.3e-2 |
| #define SO2_HENRY_TEMP 2900. |
| #define H2O2_HENRY_REF 8.3e2 |
| #define H2O2_HENRY_TEMP 7600. |
| #define SO2_CORR_A 3.12541941e-06 |
| #define SO2_CORR_B -5.72532259e-01 |
| #define O1D_RATE_CCL4_A 3.30e-10 |
| #define O1D_RATE_CCL4_B 0. |
| #define O1D_RATE_CFC11_A 2.30e-10 |
| #define O1D_RATE_CFC11_B 0. |
| #define O1D_RATE_CFC12_A 1.40e-10 |
| #define O1D_RATE_CFC12_B -25. |
| #define O1D_RATE_N2O_A 1.19e-10 |
| #define O1D_RATE_N2O_B -20. |
| #define SO2_DISS_K1_REF 1.23e-2 |
| #define SO2_DISS_K1_TEMP 2.01e3 |
| #define SO2_DISS_K2_REF 6e-8 |
| #define SO2_DISS_K2_TEMP 1.12e3 |
| #define WET_DEPO_T_LIQUID T0 |
| #define WET_DEPO_T_ICE 238.15 |
| #define WET_DEPO_T_LIQUID_BC 270. |
| #define RADIO_HALF_LIFE_RN222 (3.8235 * 86400.0) |
| #define RADIO_HALF_LIFE_PB210 (22.3 * 365.25 * 86400.0) |
| #define RADIO_HALF_LIFE_BE7 (53.22 * 86400.0) |
| #define RADIO_HALF_LIFE_CS137 (30.05 * 365.25 * 86400.0) |
| #define RADIO_HALF_LIFE_I131 (8.02 * 86400.0) |
| #define RADIO_HALF_LIFE_XE133 (5.2474 * 86400.0) |
| #define RADIO_DRY_VDEP_PB210 0.002 |
| #define RADIO_DRY_VDEP_BE7 0.001 |
| #define RADIO_DRY_VDEP_CS137 0.002 |
| #define RADIO_DRY_VDEP_I131 0.0005 |
| #define RADIO_WET_COEFF_PB210 1e-4 |
| #define RADIO_WET_COEFF_BE7 5e-5 |
| #define RADIO_WET_COEFF_CS137 1e-4 |
| #define RADIO_WET_COEFF_I131 5e-5 |
| #define EP 140 |
| #define EX 1444 |
| #define NP 10000000 |
| #define NCSI 1000000 |
| #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 DD_EX_GLOB (EX * 16) |
| #define DD_EY_GLOB (EY * 16) |
| #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 CLAMP | ( | v, | |
| lo, | |||
| hi | |||
| ) | (((v) < (lo)) ? (lo) : (((v) > (hi)) ? (hi) : (v))) |
Clamp a value to a specified range.
Ensures that v lies between lo and hi. If v < lo, returns lo. If v > hi, returns hi. Otherwise, returns v unchanged.
This macro works with any numeric type (e.g., int, float, double). All arguments are evaluated exactly once — avoid passing expressions with side effects (e.g., ++ operators or function calls).
| v | Input value to clamp. |
| lo | Lower bound. |
| hi | Upper bound. |
lo and hi.| #define COMPRESS_SPEED | ( | nbytes, | |
| dt | |||
| ) | ((dt) > 0 ? ((double)(nbytes)) / ((dt) * 1024. * 1024.) : NAN) |
Calculate compression throughput in MiB/s.
This helper converts a processed byte count and elapsed wall-clock time into a throughput value expressed in MiB/s.
| nbytes | Number of processed bytes. |
| dt | Elapsed time in seconds. |
NAN if dt is not positive.| #define COMPRESS_RATIO | ( | raw_size, | |
| stored_size | |||
| ) | ((double) (raw_size) / (double) (stored_size)) |
| #define COMPRESS_BPV | ( | n, | |
| stored_size | |||
| ) | ((8.0 * (double) (stored_size)) / (double) (n)) |
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. |
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 DEG2RAD | ( | deg | ) | ((deg) * (M_PI / 180.0)) |
Converts degrees to radians.
This macro converts an angle from degrees to radians using the formula: radians = degrees * (Ï€ / 180)
| deg | The angle in degrees to be converted. |
| #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 DX2COORD | ( | met, | |
| dx, | |||
| lat | |||
| ) | (((met)->coord_type == 0)? DX2DEG((dx) / 1000.0, (lat)) : ((dx) )) |
Convert a distance in meters to a coordinate value based on grid type.
This macro converts a distance in meters to either degrees longitude (for latitude/longitude grids) or keeps it as metres (for UTM grids), depending on the coordinate system type specified in the met_t structure.
For coord_type = 0 (latitude/longitude grid):
For coord_type = 1 (UTM grid):
| met | Pointer to met_t structure containing coordinate system information. |
| dx | Distance in meters to be converted. |
| lat | Latitude in degrees (used for degree conversion in lat/lon grids). |
| #define DY2COORD | ( | met, | |
| dy | |||
| ) | (((met)->coord_type == 0)? DY2DEG((dy) / 1000.0) : ((dy) )) |
Convert a distance to coordinate value based on grid type.
This macro converts a distance in meters to either degrees latitude (for latitude/longitude grids) or keeps it as metres (for UTM grids), depending on the coordinate system type specified in the met_t structure.
For coord_type = 0 (latitude/longitude grid):
For coord_type = 1 (UTM grid):
| met | Pointer to met_t structure containing coordinate system information. |
| dy | Distance in meters to be converted. |
| #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 ECC | ( | cmd | ) |
Execute an ECCODES command and check for errors.
This macro executes an ECCODES command and checks the result. If the result indicates an error, it prints the error message using ERRMSG.
| cmd | ECCODES command to execute. |
| #define ECC_READ_2D | ( | variable, | |
| target, | |||
| scaling_factor, | |||
| found_flag | |||
| ) |
Writes 2-D data from a grib message into the meteo struct.
This macro writes 2-D data from a one-dimensional grib message into the corresponding 2-D variable in the meteo struct.
| variable | Name of the current meteorological variable |
| target | Pointer to the 2-D array in the meteo struct where the data will be stored. |
| scaling_factor | Scaling factor to apply to the data. |
| found_flag | Flag to store, that the variable was found in the grib message. |
| #define ECC_READ_3D | ( | variable, | |
| level, | |||
| target, | |||
| scaling_factor, | |||
| found_flag | |||
| ) |
Writes 3D data from a grib message into the meteo struct.
This macro writes 3D data from a one-dimensional grib message into the corresponding 3D variable in the meteo struct.
| variable | Name of the current meteorological variable. |
| level | Index of the vertical level in the 3D array where the data should be stored. |
| target | Pointer to the 3D array in the meteo struct where the data will be stored. |
| scaling_factor | Scaling factor to apply to the data. |
| found_flag | Counter to store, how many messages containing data for this variable have been read. |
| #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, | |||
| level, | |||
| quant | |||
| ) |
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. It enables compression and quantizatio of the data.
| 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. |
| level | zlib compression level (0 = off). |
| quant | Number of digits for quantization (0 = off). |
nc_def_var_deflate() by nc_def_var_filter() below. Use dynamic linking, static linking does not work. Set environment variable HDF5_PLUGIN_PATH to ./libs/build/share/netcdf-plugins/.| #define NC_DEF_VAR_QUANTIZE | ( | quant | ) |
| #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, | |||
| check | |||
| ) |
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. |
| check | Flag to check bounds. Set to 1 for bounds check. |
| #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 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:
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:
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:
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:
| #define RAD2DEG | ( | rad | ) | ((rad) * (180.0 / M_PI)) |
Converts radians to degrees.
This macro converts an angle from radians to degrees using the formula: degrees = radians * (180 / π)
| rad | The angle in radians to be converted. |
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:
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:
| #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 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), KAPPA)) |
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)^{\kappa} \]
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.
| #define USAGE |
Print usage information on -h or --help.
Calls the local usage() function and exits successfully when any command argument is -h or --help.
argc and argv variables and a local usage() function.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 | |||
| ) |
Computes the value of the zeta vertical coordinate.
This macro calculates the zeta vertical coordinate based on the given surface pressure (ps), pressure (p), and temperature (t). The calculation depends on the ratio p/ps:
p/ps <= 0.3, the function returns 1.0 multiplied by THETA(p, t).THETA(p, t).| ps | Surface pressure. |
| p | Pressure at the given level. |
| t | Temperature at the given 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 | |||
| ) | timer(id, group, 0); |
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 and group. It uses the timer function to log the timer start event.
| id | The identifier for the timer. |
| group | The group name associated with the timer. |
| 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.
| 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 89 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 124 of file mptrac.c.
| double clim_photo | ( | const double | rate[CP][CSZA][CO3], |
| const clim_photo_t * | photo, | ||
| const double | p, | ||
| const double | sza, | ||
| const 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 156 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 213 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 241 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 396 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 414 of file mptrac.c.
| void compress_log_level | ( | FILE * | out, |
| const char * | codec, | ||
| const char * | varname, | ||
| const size_t | lev, | ||
| const double | plev, | ||
| const double | ratio, | ||
| const double | bpv, | ||
| const double | t_comp, | ||
| const double | t_decomp, | ||
| const size_t | n, | ||
| const size_t | nbytes, | ||
| const float * | org, | ||
| const float * | cmp | ||
| ) |
Write one row of per-level compression diagnostics.
This helper computes the error statistics for one level (from org and cmp) and writes one ASCII table row to out. If rho is NAN, the correlation coefficient is derived from org and cmp using GSL.
| [in,out] | out | Output stream. |
| [in] | codec | Compression codec name. |
| [in] | varname | Variable name. |
| [in] | lev | Level index. |
| [in] | plev | Pressure level in hPa. |
| [in] | ratio | Compression ratio. |
| [in] | bpv | Bits per value. |
| [in] | t_comp | Compression time in seconds. |
| [in] | t_decomp | Decompression time in seconds. |
| [in] | n | Number of data values in the level slice. |
| [in] | nbytes | Byte count used for throughput reporting. |
| [in] | org | Original level data. |
| [in] | cmp | Reconstructed level data. |
Definition at line 666 of file mptrac.c.
| void compress_log_levels_3d | ( | FILE * | out, |
| const char * | codec, | ||
| const char * | varname, | ||
| const met_t * | met, | ||
| const float * | org_all, | ||
| const float * | cmp_all, | ||
| const size_t | nxy, | ||
| const size_t | nz, | ||
| const double | ratio, | ||
| const double | bpv, | ||
| const double | t_comp, | ||
| const double | t_decomp, | ||
| const size_t | nbytes | ||
| ) |
Write per-level compression diagnostics for a full 3-D field.
This helper extracts each vertical level from the original and reconstructed 3-D fields and forwards the diagnostics to compress_log_level().
| [in,out] | out | Output stream, or NULL to disable logging. |
| [in] | codec | Compression codec name. |
| [in] | varname | Variable name. |
| [in] | met | Meteorological meta-data with pressure levels. |
| [in] | org_all | Original 3-D field in horizontal-major order. |
| [in] | cmp_all | Reconstructed 3-D field in horizontal-major order. |
| [in] | nxy | Number of horizontal points per level. |
| [in] | nz | Number of vertical levels. |
| [in] | ratio | Compression ratio. |
| [in] | bpv | Bits per value. |
| [in] | t_comp | Compression time in seconds. |
| [in] | t_decomp | Decompression time in seconds. |
| [in] | nbytes | Byte count used for throughput reporting. |
Definition at line 765 of file mptrac.c.
| void compress_cms | ( | const ctl_t * | ctl, |
| const met_t * | met, | ||
| const char * | varname, | ||
| float * | array, | ||
| const int | decompress, | ||
| FILE * | level_log, | ||
| FILE * | inout | ||
| ) |
Compresses or decompresses a 3-D meteorological field using cmultiscale.
This routine compresses or decompresses a meteorological field described by the metadata in met. Compression thresholds are selected from ctl based on varname.
| [in] | ctl | Control/settings structure providing CMS parameters. |
| [in] | met | Meteorological metadata describing the field layout. |
| [in] | varname | Variable name used for selecting the CMS epsilon threshold. |
| [in,out] | array | Field data buffer (decompress: output, compress: input). |
| [in] | decompress | Non-zero to decompress; zero to compress. |
| [in,out] | level_log | Optional output stream for per-level compression diagnostics, or NULL to disable. |
| [in,out] | inout | Binary stream used for CMS I/O (decompress: read, compress: write). |
| void compress_scale_to_unit | ( | float * | array, |
| const size_t | nxy, | ||
| const size_t | nz, | ||
| double * | off, | ||
| double * | scl | ||
| ) |
Scales each vertical level of a 3-D field independently to the interval [0,1].
The routine computes per-level offsets and scales from the minimum and maximum values of each vertical slice and transforms the field in place. Constant or nearly constant levels are mapped to zero and can be reconstructed exactly with the stored offsets.
| [in,out] | array | Field values stored in horizontal-major order. |
| [in] | nxy | Number of horizontal points per level. |
| [in] | nz | Number of vertical levels. |
| [out] | off | Per-level offsets (minimum values). |
| [out] | scl | Per-level scales (ranges). |
Definition at line 811 of file mptrac.c.
| void compress_unscale_from_unit | ( | float * | array, |
| const size_t | nxy, | ||
| const size_t | nz, | ||
| const double * | off, | ||
| const double * | scl | ||
| ) |
Restores a levelwise [0,1]-scaled 3-D field to physical units.
This routine applies the inverse transformation of compress_scale_to_unit() using the supplied per-level offsets and scales. Constant levels are restored directly from their offsets.
| [in,out] | array | Field values stored in horizontal-major order. |
| [in] | nxy | Number of horizontal points per level. |
| [in] | nz | Number of vertical levels. |
| [in] | off | Per-level offsets (minimum values). |
| [in] | scl | Per-level scales (ranges). |
Definition at line 854 of file mptrac.c.
| int compress_read_lossy_scale | ( | FILE * | in, |
| const size_t | nz, | ||
| double ** | off, | ||
| double ** | scl | ||
| ) |
Read optional lossyscaling metadata for a 3-D field.
This helper reads the stored MET_LOSSY_SCALE flag and, if enabled, allocates and fills the per-level offset and scale arrays used by compress_unscale_from_unit(). The caller owns the returned arrays and must free them.
| [in,out] | in | Input stream. |
| [in] | nz | Number of vertical levels. |
| [out] | off | Returned per-level offsets, or NULL if disabled. |
| [out] | scl | Returned per-level scales, or NULL if disabled. |
0 or 1).Definition at line 874 of file mptrac.c.
| void compress_write_lossy_scale | ( | FILE * | out, |
| const int | enabled, | ||
| float * | array, | ||
| const size_t | nxy, | ||
| const size_t | nz, | ||
| double ** | off, | ||
| double ** | scl | ||
| ) |
Write optional lossyscaling metadata for a 3-D field.
This helper writes the MET_LOSSY_SCALE flag and, if enabled, allocates the per-level offset and scale arrays, scales the field in place to [0,1], and writes the resulting metadata to out. The caller owns the returned arrays and must free them.
| [in,out] | out | Output stream. |
| [in] | enabled | Lossyscaling flag (0 or 1). |
| [in,out] | array | Field values stored in horizontal-major order. |
| [in] | nxy | Number of horizontal points per level. |
| [in] | nz | Number of vertical levels. |
| [out] | off | Returned per-level offsets, or NULL if disabled. |
| [out] | scl | Returned per-level scales, or NULL if disabled. |
Definition at line 909 of file mptrac.c.
| int compress_metvar_index | ( | const char * | varname | ) |
Maps a meteorological variable name to its internal MPTRAC variable index.
The returned index can be used to access variable-specific control arrays such as MET_LOSSY_SCALE, MET_ZFP_PREC, and MET_SZ3_PREC.
| varname | Meteorological variable name as stored in the meteorological file. |
Definition at line 969 of file mptrac.c.
| void compress_pck | ( | const ctl_t * | ctl, |
| const met_t * | met, | ||
| const char * | varname, | ||
| float * | array, | ||
| const int | decompress, | ||
| FILE * | level_log, | ||
| FILE * | inout | ||
| ) |
Compresses or decompresses a 3‑D float array using the PCK format.
The routine works on the regular longitude/latitude grid described by the supplied met structure (nx, ny, np and pressure levels). In compress mode the data are scaled to unsigned short values; in decompress mode the values are reconstructed.
| [in] | ctl | Control structure with compression parameters. |
| [in] | met | Meteorological meta‑data providing grid dimensions and pressure levels. |
| [in] | varname | Variable name for logging and error‑bound selection. |
| [in,out] | array | Float array to compress or decompress. |
| [in] | decompress | Non‑zero to decompress; zero to compress. |
| [in,out] | level_log | Optional per‑level diagnostics stream, or NULL. |
| [in,out] | inout | File stream for reading/writing the PCK payload. |
Definition at line 986 of file mptrac.c.
| void compress_sz3 | ( | const ctl_t * | ctl, |
| const met_t * | met, | ||
| const char * | varname, | ||
| float * | array, | ||
| const int | decompress, | ||
| FILE * | level_log, | ||
| FILE * | inout | ||
| ) |
Compresses or decompresses a 3-D float array using the SZ3 library.
This function either compresses a 3-D floating-point array and writes it to a file stream, or reads compressed SZ3 data from a file stream and decompresses it into the provided array. The SZ3 error bound can be specified either by relative precision (bits) or absolute tolerance.
| [in] | ctl | Control structure with compression parameters. |
| [in] | met | Meteorological data providing dimensions. |
| [in] | varname | Variable name for logging. |
| [in,out] | array | Float array to compress or decompress. |
| [in] | decompress | Non-zero to decompress; zero to compress. |
| [in,out] | level_log | Optional per-level diagnostics stream, or NULL. |
| [in,out] | inout | File stream for reading/writing compressed data. |
ctl->met_sz3_prec[] and ctl->met_sz3_tol[] to determine the error bound. ctl->met_sz3_prec[] or ctl->met_sz3_tol[] must be set to a positive value. SZ_FLOAT for this function.| ERRMSG | if input parameters are invalid, memory allocation fails, compression/decompression fails, or file I/O errors occur. |
| void compress_zfp | ( | const ctl_t * | ctl, |
| const met_t * | met, | ||
| const char * | varname, | ||
| float * | array, | ||
| const int | decompress, | ||
| FILE * | level_log, | ||
| 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.
| [in] | ctl | Control structure with compression parameters. |
| [in] | met | Meteorological data providing dimensions. |
| [in] | varname | Variable name for logging. |
| [in,out] | array | Float array to compress or decompress. |
| [in] | decompress | Non-zero to decompress; zero to compress. |
| [in,out] | level_log | Optional per-level diagnostics stream, or NULL. |
| [in,out] | inout | File stream for reading/writing compressed data. |
The function performs the following steps:
array.The function logs the compression or decompression details and frees allocated resources before returning.
ctl->met_zfp_prec[] and ctl->met_zfp_tol[] to determine the error bound. | void compress_zstd | ( | const ctl_t * | ctl, |
| const met_t * | met, | ||
| const char * | varname, | ||
| float * | array, | ||
| const int | decompress, | ||
| FILE * | level_log, | ||
| FILE * | inout | ||
| ) |
Compresses or decompresses a float array using ZSTD.
The routine compresses or decompresses one full 3-D meteorological field as a raw byte stream. It does not apply level-wise scaling or quantization; the reconstructed field is the direct result of the ZSTD frame decode.
| [in] | ctl | Control structure with compression parameters. |
| [in] | met | Meteorological data providing grid dimensions and pressure levels. |
| [in] | varname | Variable name for logging. |
| [in,out] | array | Float array to compress or decompress. |
| [in] | decompress | Non-zero to decompress; zero to compress. |
| [in,out] | level_log | Optional per-level diagnostics stream, or NULL. |
| [in,out] | inout | File stream for reading/writing the ZSTD payload. |
ctl.| void compress_lz4 | ( | const ctl_t * | ctl, |
| const met_t * | met, | ||
| const char * | varname, | ||
| float * | array, | ||
| const int | decompress, | ||
| FILE * | level_log, | ||
| FILE * | inout | ||
| ) |
Compresses or decompresses a float array using LZ4.
The routine compresses or decompresses one full 3-D meteorological field as a raw byte stream. It does not apply level-wise scaling or quantization; the reconstructed field is the direct result of the LZ4 block decode.
| [in] | ctl | Control structure with compression parameters. |
| [in] | met | Meteorological data providing grid dimensions and pressure levels. |
| [in] | varname | Variable name for logging. |
| [in,out] | array | Float array to compress or decompress. |
| [in] | decompress | Non-zero to decompress; zero to compress. |
| [in,out] | level_log | Optional per-level diagnostics stream, or NULL. |
| [in,out] | inout | File pointer for input/output. Used for reading or writing compressed data. |
LZ4_compress_fast / LZ4_decompress_safe). ctl->met_lz4_accel is not positive.| double cos_sza | ( | const double | sec, |
| const double | lon, | ||
| const double | lat | ||
| ) |
Calculates the cosine of the solar zenith angle.
This function computes the cosine of the solar zenith angle (SZA), which describes the angle between the local zenith (straight up) and the line connecting the observer to the center of the Sun. The cosine of the SZA is often used directly in radiative transfer and photochemical calculations to avoid unnecessary use of trigonometric inverse functions.
| sec | Seconds elapsed since 2000-01-01T12:00Z. |
| lon | Observer's longitude in degrees. |
| lat | Observer's latitude in degrees. |
The cosine of the solar zenith angle is computed based on the observer's position (longitude and latitude) and the specified time in seconds elapsed since 2000-01-01T12:00Z.
Definition at line 1857 of file mptrac.c.
| 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 1898 of file mptrac.c.
Assign or update particle subdomain ownership.
This routine determines whether each atmospheric particle belongs to the local MPI subdomain and updates the domain-decomposition bookkeeping quantities stored in the atmospheric state.
The particle longitude and latitude are first normalized using dd_normalize_lon_lat() to ensure consistency with the longitude convention and pole handling of the global meteorological grid stored in the domain decomposition structure.
The particle position is then mapped to global grid indices using locate_reg() for longitude and locate_irr() for latitude. These indices are compared with the index range of the local subdomain owned by the current MPI rank.
If init is nonzero, the routine performs the initial ownership assignment. Particles inside the local subdomain are assigned to the current MPI rank by setting both ctl->qnt_current_subdomain and ctl->qnt_target_subdomain to that rank. Particles outside the local subdomain are marked invalid by setting both quantities to -1.
If init is zero, the routine updates only the destination rank. Particles already marked invalid (qnt_current_subdomain == -1) are skipped. For particles that remain inside the local subdomain, ctl->qnt_target_subdomain is set to the current MPI rank. For particles that have left the local subdomain, the destination rank is recomputed with dd_calc_subdomain_from_coords().
| [in] | ctl | Control structure containing the indices of the domain-decomposition quantities qnt_current_subdomain and qnt_target_subdomain. |
| [in] | dd | Domain decomposition structure containing the global grid and subdomain layout information. |
| [in,out] | atm | Atmospheric state. Particle positions are read from this structure and the subdomain and/or destination quantities are updated in place. |
| [in] | init | Initialization flag. If nonzero, perform the initial subdomain assignment; otherwise only update destination ranks for particle exchange. |
MPI_COMM_WORLD._OPENACC, the particle loop may execute on an accelerator.| void dd_atm2particles | ( | const ctl_t * | ctl, |
| cache_t * | cache, | ||
| atm_t * | atm, | ||
| particle_t * | particles, | ||
| const int | npart | ||
| ) |
Copy migratable atmospheric particles from the ATM state into a particle buffer.
This routine scans the domain-decomposition particle region of the atmospheric state and copies particles that must be transferred to a different MPI rank into the temporary particle communication buffer.
The routine processes atmospheric entries in the index range [atm->np, atm->np + npart). A particle is selected for export if all of the following conditions are met:
atm->q[ctl->qnt_target_subdomain][ip] != rank,atm->q[ctl->qnt_target_subdomain][ip] >= 0,atm->q[ctl->qnt_current_subdomain][ip] >= 0.For each selected entry, the particle coordinates (time, lon, lat, p) and all quantities q[iq] are copied into particles[ip - atm->np].
After copying, the atmospheric entry is marked as inactive for local ownership by setting atm->q[ctl->qnt_current_subdomain][ip] = -1, and the cached time step cache->dt[ip] is reset to zero.
| [in] | ctl | Pointer to the control structure containing the quantity indices qnt_target_subdomain, qnt_current_subdomain, and the number of quantities nq. |
| [in,out] | cache | Pointer to the cache structure. The cached time-step values dt[ip] of exported particles are reset to zero. |
| [in,out] | atm | Pointer to the atmospheric state. Particle data are read from the domain-decomposition particle region, and exported entries are marked by setting their subdomain index to -1. |
| [out] | particles | Buffer receiving the particles selected for communication. Selected atmospheric entries are copied to index ip - atm->np. |
| [in] | npart | Number of candidate particle entries in the atmospheric communication region. |
MPI_COMM_WORLD._OPENACC, the selection and copy loop may run on the accelerator.| int dd_calc_subdomain_from_coords | ( | const ctl_t * | ctl, |
| const dd_t * | dd, | ||
| const double | lon, | ||
| const double | lat | ||
| ) |
Determine MPI subdomain from particle coordinates.
This function maps a particle position to the corresponding MPI rank in the domain-decomposed grid. The longitude and latitude are first normalized using dd_normalize_lon_lat(). The corresponding global grid cell indices are then determined using locate_reg() for longitude and locate_irr() for latitude.
The grid indices are mapped to the zonal and meridional subdomain blocks defined by ctl->dd_subdomains_zonal and ctl->dd_subdomains_meridional. The resulting subdomain indices are converted to an MPI rank assuming the layout:
rank = zonal_rank * dd_subdomains_meridional + merid_rank
This implementation works for both regular and irregular latitude grids.
| ctl | Pointer to control structure. |
| dd | Pointer to domain decomposition structure. |
| lon | Longitude of particle [deg]. |
| lat | Latitude of particle [deg]. |
| void dd_communicate_particles | ( | const ctl_t * | ctl, |
| const dd_t * | dd, | ||
| particle_t ** | particles, | ||
| int * | npart, | ||
| int * | capacity | ||
| ) |
Exchange particles between MPI ranks according to their destination rank.
This routine migrates particles to their assigned destination MPI ranks using collective communication. The destination rank of each particle is read from (*particles)[ip].q[ctl->qnt_target_subdomain].
First, the routine counts how many local particles must be sent to each MPI rank, excluding particles whose destination equals the current rank. These counts are exchanged using MPI_Alltoall to determine the number of incoming particles. Temporary send and receive buffers are then allocated, the outgoing particles are packed into the send buffer, and particle data are exchanged with MPI_Alltoallv using the MPI datatype stored in dd->MPI_Particle.
After communication, the received particles replace the leading entries of the particle buffer. For each received particle, both ctl->qnt_target_subdomain and ctl->qnt_current_subdomain are reset to the current MPI rank.
If the receive buffer exceeds the currently allocated particle capacity, the particle array is automatically resized using realloc().
| [in] | ctl | Pointer to the control structure containing the quantity indices qnt_target_subdomain and qnt_current_subdomain. |
| [in] | dd | Pointer to the domain decomposition structure containing the MPI datatype MPI_Particle used for particle communication. |
| [in,out] | particles | Pointer to the particle buffer. The pointer may be updated if the buffer is resized. On output, the leading entries contain the particles received by the current MPI rank. |
| [in,out] | npart | Number of local particles. On input, it specifies the number of valid particles stored in particles. On output, it is replaced by the number of particles received by the current MPI rank. |
| [in,out] | capacity | Allocated capacity of the particle buffer. If necessary, the routine increases this value when resizing the particle array. |
Initialize the domain decomposition infrastructure.
This routine initializes the MPI-related data structures required for domain decomposition. It first verifies that the number of MPI ranks in MPI_COMM_WORLD matches the total number of requested subdomains,
ctl->dd_subdomains_meridional * ctl->dd_subdomains_zonal.
It then constructs and commits the MPI datatype dd->MPI_Particle corresponding to the particle_t structure. This datatype is used for particle exchange between MPI ranks.
Finally, the routine performs the initial particle-to-subdomain assignment by calling dd_assign_subdomains() with init = 1.
| [in] | ctl | Pointer to the control structure containing the domain-decomposition settings, in particular the number of zonal and meridional subdomains. |
| [in,out] | dd | Pointer to the domain decomposition structure. On return, the MPI datatype MPI_Particle is initialized and committed. |
| [in,out] | atm | Pointer to the atmospheric state. Its domain-decomposition bookkeeping quantities are initialized by the call to dd_assign_subdomains(). |
dd->MPI_Particle must be released later with MPI_Type_free() when the domain decomposition is finalized.| void dd_normalize_lon_lat | ( | const dd_t * | dd, |
| double * | lon, | ||
| double * | lat | ||
| ) |
Normalize geographic coordinates to the global grid convention.
This helper routine normalizes a longitude/latitude pair so that it is consistent with the longitude convention of the global meteorological grid stored in the domain decomposition structure.
The function performs the following steps:
dd->lon_glob (either [-180,180) or [0,360)).This ensures that coordinates used for domain-decomposition checks and subdomain index calculations are always mapped consistently to the global grid representation.
| [in] | dd | Pointer to the domain decomposition structure containing the global longitude grid (dd->lon_glob). |
| [in,out] | lon | Pointer to longitude [deg]. The value is normalized in-place to match the grid convention. |
| [in,out] | lat | Pointer to latitude [deg]. Values outside [-90,90] are wrapped across the poles in-place. |
FMOD macro instead of fmod() to ensure compatibility with GPU/OpenACC builds.#pragma acc routine seq so it can be called from within OpenACC parallel loops.| void dd_particles2atm | ( | const ctl_t * | ctl, |
| cache_t * | cache, | ||
| const particle_t * | particles, | ||
| const int | npart, | ||
| atm_t * | atm | ||
| ) |
Copy received particles from the communication buffer into the atmospheric state.
This routine appends particles from the temporary domain-decomposition particle buffer to the atmospheric particle arrays. The particles are written into the index range
[atm->np, atm->np + npart)
of the atmospheric state.
For each transferred particle, the routine copies the particle coordinates (time, lon, lat, p) and all particle quantities q[iq] into the atmospheric arrays. The cached time-step value cache->dt[ip] is initialized to ctl->dt_mod for each appended particle.
After all particles have been copied, the total number of atmospheric particles atm->np is increased by npart.
| [in] | ctl | Pointer to the control structure containing the number of quantities nq and the model time step dt_mod. |
| [in,out] | cache | Pointer to the cache structure. The cached time-step values for appended particles are initialized to ctl->dt_mod. |
| [in] | particles | Buffer containing the particles to be appended to the atmospheric state. |
| [in] | npart | Number of valid particles stored in particles. |
| [in,out] | atm | Pointer to the atmospheric state. The particle data are appended to its arrays, and atm->np is incremented accordingly. |
atm->np.NP._OPENACC, particle data may be copied on the accelerator and the updated value of atm->np is synchronized to the device.Compact and stage non-local particles for domain decomposition transfer.
This routine scans all local atmospheric particles and moves those that need to be transferred to other MPI ranks to a contiguous region at the end of the particle arrays (starting at index atm->np).
A particle is selected for transfer if:
atm->q[ctl->qnt_current_subdomain][ip] != -1 (particle is active), ANDatm->q[ctl->qnt_target_subdomain][ip] != rank (destination is not current rank)For each selected particle:
q[ctl->qnt_target_subdomain][ip] = -1 and q[ctl->qnt_current_subdomain][ip] = -1cache->dt[ip] is copied to the new location and reset to zero at the original position| [in] | ctl | Pointer to the control structure containing the quantity indices qnt_target_subdomain and qnt_current_subdomain. |
| [in,out] | atm | Pointer to the atmospheric state. Selected particles are copied to indices [atm->np, atm->np+npart) and original entries are invalidated. |
| [in,out] | cache | Pointer to the cache structure. Time-step values for moved particles are preserved; original positions are reset to zero. |
| [out] | npart | Number of particles staged for transfer (i.e., the number of valid particles in the range [atm->np, atm->np + npart)). |
MPI_COMM_WORLD._OPENACC, the selection and copy operations execute on the accelerator using atomic operations for the counter.Read meteorological grid information and construct the domain-decomposed grid with halo regions.
This routine reads the global longitude and latitude grid from a NetCDF meteorological input file and initializes the domain decomposition used for parallel processing. The global grid is partitioned into zonal and meridional subdomains according to ctl->dd_subdomains_zonal and ctl->dd_subdomains_meridional.
For the current MPI rank, the routine determines the core subdomain in index space and constructs a hyperslab describing the portion of the meteorological grid that must be read from the NetCDF file. The hyperslab is extended with inner halo cells to overlap with neighboring subdomains. For ranks located at the zonal boundaries, additional periodic boundary halos are created.
The function populates the local longitude and latitude arrays in the met_t structure for the subdomain including halos, and stores the hyperslab definitions and halo offsets in the dd_t structure.
| [out] | dd | Domain decomposition structure. On return it contains global grid dimensions, hyperslab definitions for the local subdomain, and halo metadata. |
| [in] | ctl | Control structure specifying the number of zonal and meridional subdomains and the halo size (dd_subdomains_zonal, dd_subdomains_meridional, dd_halos_size). |
| [out] | met | Meteorological grid structure. On return it contains the local longitude and latitude arrays and the grid dimensions (nx, ny, np) of the halo-extended subdomain. |
| [in] | ncid | NetCDF file identifier returned by nc_open(). |
-DMPI.Definition at line 11484 of file mptrac.c.
Sort local atmospheric particles and determine export counts for domain decomposition.
This routine reorders the local atmospheric particle arrays so that particles remaining on the current MPI rank are placed first, followed by particles that must be sent to other ranks, and finally particles already marked invalid with subdomain index -1.
For each particle a sorting key is constructed:
-1 receive the largest keys and are moved to the end of the arrays.The resulting permutation is applied to all particle arrays (time, p, lon, lat, and all quantities q[iq]) using dd_sort_help().
After sorting, the routine
npart,atm->np to the number of particles kept locally.Particles with subdomain index -1 are removed from the active particle set. A warning is issued if such particles are encountered.
| [in] | ctl | Control structure containing the quantity indices qnt_current_subdomain, qnt_target_subdomain, and the number of quantities nq. |
| [in] | met0 | Meteorological grid used to compute sorting keys from particle longitude, latitude, and pressure. |
| [in,out] | atm | Atmospheric state. Particle arrays are reordered in place and atm->np is reduced to the number of particles retained locally. |
| [in,out] | dd | Domain decomposition structure providing the temporary arrays sort_key and perm used during sorting. |
| [out] | npart | Number of particles that must be sent to other MPI ranks. |
MPI_COMM_WORLD._OPENACC, key generation, counting, and selected data movement may execute on the accelerator.THRUST, sorting may use the GPU-based thrustSortWrapper; otherwise a GSL-based index sort is used.| void dd_sort_help | ( | double * | a, |
| dd_t * | dd, | ||
| const int | np | ||
| ) |
Apply the sorting permutation to a particle data array.
This helper routine reorders a one-dimensional particle data array according to the permutation indices stored in dd->perm. The reordering is performed in two passes using the temporary buffer dd->tmp:
This routine is used by dd_sort() to apply the same permutation to all atmospheric particle arrays, ensuring that time, position, pressure, and tracer quantities remain consistently aligned.
| [in,out] | a | Pointer to the array to be reordered in place. |
| [in,out] | dd | Pointer to the domain decomposition structure containing:
|
| [in] | np | Number of elements in the array. |
_OPENACC, the reordering loops may execute on the accelerator; otherwise OpenMP is used.| 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 2532 of file mptrac.c.
| void fft_help | ( | double * | fcReal, |
| double * | fcImag, | ||
| const 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 2562 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. | void get_met_filename | ( | const ctl_t * | ctl, |
| const double | t, | ||
| const int | direct, | ||
| const char * | metbase, | ||
| const double | dt_met, | ||
| char * | filename | ||
| ) |
Generates a formatted filename for meteorological data files based on the input parameters.
This function determines a rounded time interval, decodes the time components (year, month, day, hour, minute, second), and constructs a filename string for meteorological data files in various formats. The filename is adjusted based on the input control settings.
| [in] | ctl | Pointer to the control structure containing configuration settings. |
| [in] | t | The time value in seconds since a reference epoch. |
| [in] | direct | Direction to round the time value. Use -1 for rounding down and 1 for rounding up. |
| [in] | metbase | Base string for the filename, representing the dataset. |
| [in] | dt_met | Time interval for rounding in seconds. |
| [out] | filename | Output buffer to store the generated filename. |
Definition at line 2619 of file mptrac.c.
| void get_met_replace | ( | char * | orig, |
| const char * | search, | ||
| const 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 2688 of file mptrac.c.
| void get_tropo | ( | const int | met_tropo, |
| ctl_t * | ctl, | ||
| const clim_t * | clim, | ||
| met_t * | met, | ||
| const double * | lons, | ||
| const int | nx, | ||
| const double * | lats, | ||
| const 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 2712 of file mptrac.c.
| void intpol_check_lon_lat | ( | const double * | lons, |
| const int | nlon, | ||
| const double * | lats, | ||
| const int | nlat, | ||
| const double | lon, | ||
| const double | lat, | ||
| double * | lon2, | ||
| double * | lat2 | ||
| ) |
Adjusts longitude and latitude to ensure they fall within valid bounds.
This function checks and modifies the given longitude and latitude values to fit within the specified longitude and latitude arrays. The longitude is wrapped within a 360-degree range, and the latitude is clamped within the valid range defined by the latitude array.
| [in] | lons | Pointer to an array of valid longitude values. |
| [in] | nlon | Number of elements in the longitude array. |
| [in] | lats | Pointer to an array of valid latitude values. |
| [in] | nlat | Number of elements in the latitude array. |
| [in] | lon | Input longitude to be checked and adjusted. |
| [in] | lat | Input latitude to be checked and adjusted. |
| [out] | lon2 | Pointer to the adjusted longitude. |
| [out] | lat2 | Pointer to the adjusted latitude. |
Definition at line 2755 of file mptrac.c.
| void intpol_check_cartesian | ( | const double * | lons, |
| const int | nlon, | ||
| const double * | lats, | ||
| const int | nlat, | ||
| const double | lon, | ||
| const double | lat, | ||
| double * | lon2, | ||
| double * | lat2 | ||
| ) |
Clamps UTM coordinates to the valid bounds.
This function constrains the given UTM easting and northing values so they remain within the limits defined by the provided x and y coordinate arrays.
| [in] | lons | Pointer to an array of valid x/easting values. |
| [in] | nlon | Number of elements in the x/easting array. |
| [in] | lats | Pointer to an array of valid y/northing values. |
| [in] | nlat | Number of elements in the y/northing array. |
| [in] | lon | Input x/easting coordinate to be checked and adjusted. |
| [in] | lat | Input y/northing coordinate to be checked and adjusted. |
| [out] | lon2 | Pointer to the adjusted x/easting coordinate. |
| [out] | lat2 | Pointer to the adjusted y/northing coordinate. |
Definition at line 2782 of file mptrac.c.
| void intpol_met_4d_zeta | ( | const met_t * | met0, |
| float | height0[EX][EY][EP], | ||
| float | array0[EX][EY][EP], | ||
| const met_t * | met1, | ||
| float | height1[EX][EY][EP], | ||
| float | array1[EX][EY][EP], | ||
| const double | ts, | ||
| const double | height, | ||
| const double | lon, | ||
| const double | lat, | ||
| double * | var, | ||
| int * | ci, | ||
| double * | cw, | ||
| const 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 2808 of file mptrac.c.
| void intpol_met_space_3d | ( | const met_t * | met, |
| float | array[EX][EY][EP], | ||
| const double | p, | ||
| const double | lon, | ||
| const double | lat, | ||
| double * | var, | ||
| int * | ci, | ||
| double * | cw, | ||
| const 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 2985 of file mptrac.c.
| void intpol_met_space_2d | ( | const met_t * | met, |
| float | array[EX][EY], | ||
| const double | lon, | ||
| const double | lat, | ||
| double * | var, | ||
| int * | ci, | ||
| double * | cw, | ||
| const 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 3048 of file mptrac.c.
| void intpol_met_time_3d | ( | const met_t * | met0, |
| float | array0[EX][EY][EP], | ||
| const met_t * | met1, | ||
| float | array1[EX][EY][EP], | ||
| const double | ts, | ||
| const double | p, | ||
| const double | lon, | ||
| const double | lat, | ||
| double * | var, | ||
| int * | ci, | ||
| double * | cw, | ||
| const 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 3112 of file mptrac.c.
| void intpol_met_time_2d | ( | const met_t * | met0, |
| float | array0[EX][EY], | ||
| const met_t * | met1, | ||
| float | array1[EX][EY], | ||
| const double | ts, | ||
| const double | lon, | ||
| const double | lat, | ||
| double * | var, | ||
| int * | ci, | ||
| double * | cw, | ||
| const 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 3141 of file mptrac.c.
| void intpol_tropo_3d | ( | const double | time0, |
| float | array0[EX][EY], | ||
| const double | time1, | ||
| float | array1[EX][EY], | ||
| const double | lons[EX], | ||
| const double | lats[EY], | ||
| const int | nlon, | ||
| const int | nlat, | ||
| const double | time, | ||
| const double | lon, | ||
| const double | lat, | ||
| const 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 3174 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 3265 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 3298 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 3324 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 met_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 met_press_level_def in the control structure ctl. It initializes the met_np and met_p fields accordingly.
met_press_level_def are:Definition at line 3342 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 3495 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 3525 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 3559 of file mptrac.c.
| void locate_vert | ( | float | profiles[EX][EY][EP], |
| const int | np, | ||
| const int | lon_ap_ind, | ||
| const int | lat_ap_ind, | ||
| const 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 3578 of file mptrac.c.
| void module_advect | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
Advances particle positions using different advection schemes.
This function updates the positions of atmospheric particles using different advection schemes based on vertical velocity formulations (omega or zetadot). The advection is performed over a number of integration nodes, using meteorological data interpolated in time and space.
| [in] | ctl | Pointer to the control structure containing configuration settings. |
| [in] | cache | Pointer to the cache structure storing precomputed time step values. |
| [in] | met0 | Pointer to the meteorological data structure at the initial time. |
| [in] | met1 | Pointer to the meteorological data structure at the next time step. |
| [in,out] | atm | Pointer to the atmospheric data structure containing particle states. |
ctl->advect_vert_coord is 0 or 2, the function uses omega vertical velocity.ctl->advect_vert_coord is 1, the function uses zetadot vertical velocity.zeta values remain non-negative.Definition at line 3598 of file mptrac.c.
| void module_advect_init | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| 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 3762 of file mptrac.c.
| void module_bound_cond | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| const clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| 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. |
Definition at line 3789 of file mptrac.c.
| void module_chem_grid | ( | const ctl_t * | ctl, |
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| const double | t | ||
| ) |
Computes gridded chemical tracer concentrations (volume mixing ratio) from individual air parcel mass data and assigns them back to the parcels.
This function aggregates the mass of tracer particles onto a 3D chemical grid (longitude × latitude × altitude), accounting for either single or ensemble simulations depending on ctl->nens. It then interpolates meteorological temperature fields and computes volume mixing ratios, storing them in the specified tracer quantity (e.g., ctl->qnt_Cx).
If the molar mass is undefined or required quantity indices are missing, the function exits early.
Parallelization is supported via OpenMP or OpenACC.
| [in] | ctl | Pointer to the control structure containing configuration parameters, including grid dimensions, tracer indices, and simulation mode. |
| [in] | met0 | Pointer to the meteorological data at the beginning of the interpolation interval. |
| [in] | met1 | Pointer to the meteorological data at the end of the interpolation interval. |
| [in,out] | atm | Pointer to the atmospheric state, including parcel coordinates, time, mass, and output tracer fields. |
| [in] | t | Central time step used for output and interpolation. |
ctl->molmass > 0, and ctl->qnt_m and ctl->qnt_Cx to be set.ctl->nens > 0 and assigns each parcel to its ensemble member via ctl->qnt_ens.qnt_Cx) is given in ppbv.Definition at line 3885 of file mptrac.c.
| void module_chem_init | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| const clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| 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 4065 of file mptrac.c.
| void module_convection | ( | const ctl_t * | ctl, |
| cache_t * | cache, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
Performs convective mixing of atmospheric particles.
This function adjusts the pressure of atmospheric particles based on boundary layer (PBL) mixing and convective conditions driven by CAPE (Convective Available Potential Energy) and CIN (Convective Inhibition). It uses meteorological data and random numbers for vertical mixing calculations.
| [in] | ctl | Pointer to the control structure with simulation settings. |
| [in,out] | cache | Pointer to the cache structure for temporary data and random numbers. |
| [in,out] | met0 | Pointer to the meteorological data at the initial timestep. |
| [in,out] | met1 | Pointer to the meteorological data at the subsequent timestep. |
| [in,out] | atm | Pointer to the atmospheric data structure with particle properties. |
atm structure in place.Definition at line 4109 of file mptrac.c.
| void module_dd | ( | double | t, |
| const ctl_t * | ctl, | ||
| cache_t * | cache, | ||
| dd_t * | dd, | ||
| atm_t * | atm, | ||
| met_t ** | met | ||
| ) |
Perform domain decomposition and exchange particles between MPI ranks.
This routine redistributes atmospheric particles across MPI processes according to the current spatial domain decomposition so that each process owns the particles located inside its subdomain.
The routine performs the following steps:
A temporary particle buffer is maintained locally and resized dynamically when necessary to hold outgoing or incoming particles. After the communication step, the atmospheric particle arrays are updated to contain only the particles owned by the current MPI rank.
| [in] | t | Current simulation time in seconds since 2000-01-01. |
| [in] | ctl | Control structure containing simulation settings and quantity indices. |
| [in,out] | cache | Cache structure used for temporary particle data during data layout transformations. |
| [in,out] | dd | Domain decomposition structure containing grid metadata and auxiliary buffers. |
| [in,out] | atm | Atmospheric particle state. Particle ownership and particle arrays are updated in place. |
| [in] | met | Pointer to the meteorological data used for subdomain assignment and sorting. |
DD) and MPI enabled.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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| clim | Pointer to the climate data structure containing atmospheric data. |
| atm | Pointer to the atmospheric data structure containing particle information. |
Definition at line 4234 of file mptrac.c.
| void module_diff_meso | ( | const ctl_t * | ctl, |
| cache_t * | cache, | ||
| const met_t * | met0, | ||
| const met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| 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. |
TURB_MESOX and TURB_MESOZ define the subgrid-scale variability as a fraction of the grid-scale variance. Stohl et al. (2005) recommend a default value of 0.16 for both parameters, providing a standard approach for turbulence representation. However, recent findings by Bakels et al. (2024) suggest disabling this approach to improve model accuracy under certain conditions. It is advised to evaluate the applicability of these recommendations based on the specific simulation context and objectives.Definition at line 4273 of file mptrac.c.
| void module_diff_pbl | ( | const ctl_t * | ctl, |
| cache_t * | cache, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
Computes particle diffusion within the planetary boundary layer (PBL).
This function handles the effects of turbulence on particles within the PBL. It calculates turbulent velocity variances, Lagrangian timescales, and updates particle positions and perturbations based on random fluctuations and boundary layer physics. This module adapts the approach of Ryall and Maryon (1998) and Stohl et al. (2005).
| ctl | Pointer to the control structure containing model settings. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| met0 | Pointer to the meteorological data structure for the current timestep. |
| met1 | Pointer to the meteorological data structure for the next timestep. |
| atm | Pointer to the atmospheric data structure containing particle states. |
The function:
module_rng.sig_u, sig_w), their vertical derivatives, and Lagrangian timescales (tau_u, tau_w).The function uses OpenACC directives for GPU acceleration.
Definition at line 4350 of file mptrac.c.
| void module_diff_turb | ( | const ctl_t * | ctl, |
| cache_t * | cache, | ||
| const clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
Applies turbulent diffusion processes to atmospheric particles.
This function calculates and applies turbulent diffusion effects, including horizontal and vertical diffusion, as well as vertical mixing in the planetary boundary layer (PBL), to a set of atmospheric particles based on input parameters and environmental conditions.
| [in] | ctl | Pointer to the control structure containing simulation parameters. |
| [in,out] | cache | Pointer to the cache structure for temporary data and random numbers. |
| [in] | clim | Pointer to the climate structure containing climatological data. |
| [in,out] | met0 | Pointer to the meteorological data structure for the initial timestep. |
| [in,out] | met1 | Pointer to the meteorological data structure for the next timestep. |
| [in,out] | atm | Pointer to the atmospheric structure containing particle data. |
The function performs the following operations:
Turbulent diffusivity parameters are derived from control inputs and weighted based on atmospheric layer influences (PBL, troposphere, stratosphere).
TURB_DX_PBL, TURB_DX_TROP, TURB_DX_STRAT, TURB_DZ_PBL, TURB_DZ_TROP, and TURB_DZ_STRAT define horizontal and vertical diffusivities (in units of m**2 s**-1) in the PBL, troposphere, and stratosphere, respectively.TURB_PBL_SCHEME activates an optional PBL-specific closure scheme. If a PBL scheme is active, the fixed-K module is only applied outside the PBL.Definition at line 4595 of file mptrac.c.
| void module_dry_depo | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| 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. |
Definition at line 4745 of file mptrac.c.
| void module_h2o2_chem | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| const clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| 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. |
Definition at line 4808 of file mptrac.c.
| void module_isosurf_init | ( | const ctl_t * | ctl, |
| cache_t * | cache, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
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 4893 of file mptrac.c.
| void module_isosurf | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| 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. |
Definition at line 4963 of file mptrac.c.
| void module_kpp_chem | ( | ctl_t * | ctl, |
| cache_t * | cache, | ||
| clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| 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. |
| void module_meteo | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| const clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| 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. |
Definition at line 5069 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 5195 of file mptrac.c.
| void module_mixing_help | ( | const ctl_t * | ctl, |
| atm_t * | atm, | ||
| const int * | ixs, | ||
| const int * | iys, | ||
| const int * | izs, | ||
| const double * | mixparam, | ||
| const int | qnt_idx, | ||
| const int | use_ensemble | ||
| ) |
Perform subgrid-scale interparcel mixing of a given quantity.
This function computes the mixing-parameter-weighted mean of a specified quantity within each subgrid box (and optionally for each ensemble member) and applies a mixing adjustment to particle values based on the computed local mean. The weighted mean conserves the sum of the quantity when the mixing parameter varies between particles.
The mixing accounts for differences in tropospheric and stratospheric mixing via a weighted parameterization. It supports both ensemble and non-ensemble modes using the use_ensemble flag.
| [in] | ctl | Pointer to control/configuration structure. |
| [in,out] | atm | Pointer to atmospheric state (includes particles). |
| [in] | ixs | Array of x-grid indices for each particle. |
| [in] | iys | Array of y-grid indices for each particle. |
| [in] | izs | Array of z-grid indices for each particle (-1 for invalid). |
| [in] | mixparam | Mixing parameter for each particle. |
| [in] | qnt_idx | Index of the quantity in atm->q to be mixed. |
| [in] | use_ensemble | Flag indicating whether to use ensemble-aware logic (0 = no, 1 = yes). |
izs[ip] < 0 are excluded from mixing. ctl->qnt_ens to be valid if use_ensemble is true.Definition at line 5287 of file mptrac.c.
| void module_oh_chem | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| const clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| 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. |
Definition at line 5390 of file mptrac.c.
| void module_position | ( | const cache_t * | cache, |
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| const int | reflect | ||
| ) |
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:
reflect.| cache | Pointer to the cache structure for temporary data and random numbers. |
| 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. |
| reflect | Select reflection (1) or clamping (0) at vertical boundaries. |
Definition at line 5474 of file mptrac.c.
| void module_rng_init | ( | const 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 5763 of file mptrac.c.
| void module_rng | ( | const ctl_t * | ctl, |
| double * | rs, | ||
| const size_t | n, | ||
| const 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 5794 of file mptrac.c.
Apply radioactive decay to atmospheric tracer species.
This routine updates the concentrations of radioactive tracers carried by atmospheric particles by applying exponential decay over the current particle timestep. The decay constants are derived from the half-lives of the respective isotopes.
Implemented isotopes:
For each particle, the tracer mixing ratios are reduced according to \( q(t+\Delta t) = q(t) \exp(-\lambda \Delta t) \), where \(\lambda\) is the decay constant and \(\Delta t\) is the particle timestep.
Additionally, the decay of Rn-222 contributes to the production of Pb-210 via a simplified parent–daughter relationship.
The update is performed only if the corresponding tracer index in the control structure is non-negative.
| [in] | ctl | Control structure containing tracer indices. |
| [in] | cache | Cache structure providing particle timesteps. |
| [in,out] | atm | Atmospheric state containing particle tracer fields that are updated in place. |
Definition at line 5533 of file mptrac.c.
| void module_radio_depo | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| depo_t * | depo | ||
| ) |
Deposit supported radionuclides from air parcels onto the ground grid.
Applies fixed dry-deposition velocities in the surface layer and precipitation-dependent wet scavenging to Apb210, Abe7, Acs137, and Ai131. Air-parcel activities are reduced exponentially and the removed activity is accumulated in depo. When RADIO_DECAY is enabled, ground inventories are referenced to the simulation start for decay-correct output. Backward or zero-length particle timesteps are ignored.
| [in] | ctl | Control parameters, quantity indices, and deposition grid. |
| [in] | cache | Per-particle model timesteps. |
| [in] | met0 | Meteorological field at the earlier bracketing time. |
| [in] | met1 | Meteorological field at the later bracketing time. |
| [in,out] | atm | Air-parcel positions and activities [Bq]. |
| [in,out] | depo | Cumulative ground inventories [Bq per grid cell]. |
Definition at line 5591 of file mptrac.c.
| void module_sedi | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| 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. |
sedi function, which takes atmospheric pressure, temperature, particle radius, and particle density as inputs. DZ2DP function.Definition at line 5899 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 5927 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 6001 of file mptrac.c.
| void module_timesteps | ( | const ctl_t * | ctl, |
| cache_t * | cache, | ||
| const met_t * | met0, | ||
| atm_t * | atm, | ||
| const double | t | ||
| ) |
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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| met0 | Pointer to the initial meteorological data structure. |
| atm | Pointer to the atmospheric data structure containing air parcel information. |
| t | The target time for which time steps are calculated. |
t. Definition at line 6039 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 6086 of file mptrac.c.
| void module_tracer_chem | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| const clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| 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. |
Definition at line 6117 of file mptrac.c.
| void module_wet_depo | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm | ||
| ) |
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. |
| cache | Pointer to the cache structure for temporary data and random numbers. |
| 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. |
Definition at line 6195 of file mptrac.c.
| void mptrac_alloc | ( | ctl_t ** | ctl, |
| cache_t ** | cache, | ||
| clim_t ** | clim, | ||
| met_t ** | met0, | ||
| met_t ** | met1, | ||
| atm_t ** | atm, | ||
| depo_t ** | depo, | ||
| dd_t ** | dd | ||
| ) |
Allocates and initializes memory resources for MPTRAC.
This function handles memory allocation for various data structures and sets up GPU resources if available. It also creates the necessary data regions on GPUs for OpenACC-enabled execution. Each output parameter is optional; pass NULL for structures that are not needed.
| [out] | ctl | Pointer to the control structure (ctl_t). |
| [out] | cache | Pointer to the cache structure (cache_t). |
| [out] | clim | Pointer to the climatology structure (clim_t). |
| [out] | met0 | Pointer to the first meteorology structure (met_t). |
| [out] | met1 | Pointer to the second meteorology structure (met_t). |
| [out] | atm | Pointer to the atmospheric structure (atm_t). |
| [out] | depo | Pointer to the radionuclide ground-inventory structure (depo_t). |
| [out] | dd | pointer to an dd_t structure containing MPI information, including rank and neighbours. |
| Runtime | error if no GPU devices are available when OpenACC is enabled. |
Definition at line 6334 of file mptrac.c.
| void mptrac_free | ( | ctl_t * | ctl, |
| cache_t * | cache, | ||
| clim_t * | clim, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| depo_t * | depo, | ||
| dd_t * | dd | ||
| ) |
Frees memory resources allocated for MPTRAC.
This function releases the memory allocated for various data structures and deletes any associated data regions on GPUs if OpenACC is enabled. Every argument may be NULL.
| [in] | ctl | Pointer to the control structure (ctl_t) to be freed. |
| [in] | cache | Pointer to the cache structure (cache_t) to be freed. |
| [in] | clim | Pointer to the climatology structure (clim_t) to be freed. |
| [in] | met0 | Pointer to the first meteorology structure (met_t) to be freed. |
| [in] | met1 | Pointer to the second meteorology structure (met_t) to be freed. |
| [in] | atm | Pointer to the atmospheric structure (atm_t) to be freed. |
| [in] | depo | Pointer to the radionuclide ground-inventory structure (depo_t) to be freed. |
| [in] | dd | Pointer to an dd_t structure containing MPI information, including rank and neighbours. |
Definition at line 6417 of file mptrac.c.
| void mptrac_get_met | ( | ctl_t * | ctl, |
| clim_t * | clim, | ||
| const double | t, | ||
| met_t ** | met0, | ||
| met_t ** | met1, | ||
| dd_t * | dd | ||
| ) |
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. |
| dd | A pointer to an dd_t structure containing MPI information, including rank and neighbours. |
The function performs the following steps:
ctl, clim, met0, and met1 are properly initialized before calling this function.Definition at line 6478 of file mptrac.c.
| void mptrac_init | ( | ctl_t * | ctl, |
| cache_t * | cache, | ||
| clim_t * | clim, | ||
| atm_t * | atm, | ||
| depo_t * | depo, | ||
| const int | ntask | ||
| ) |
Initializes the MPTRAC model and its associated components.
This function sets up the necessary components and subsystems for the MPTRAC module, including timesteps, random number generation, and GPU memory updates.
| ctl | Pointer to the control structure containing configuration and state information. |
| cache | Pointer to the cache structure used for data storage and retrieval. |
| clim | Pointer to the climatology structure containing climate-related data. |
| atm | Pointer to the atmospheric structure containing atmospheric state data. |
| depo | Pointer to inventories allocated and zero-initialized by mptrac_alloc(). |
| ntask | Number of tasks or threads to initialize for the random number generator. |
The function performs the following operations:
module_timesteps_init function.module_rng_init function.mptrac_update_device function.Definition at line 6603 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 6628 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 6703 of file mptrac.c.
| void mptrac_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 6763 of file mptrac.c.
| int mptrac_read_met | ( | const char * | filename, |
| const ctl_t * | ctl, | ||
| const clim_t * | clim, | ||
| met_t * | met, | ||
| dd_t * | dd | ||
| ) |
Reads meteorological data from a file, supporting multiple formats and MPI broadcasting.
This function reads meteorological data from a file specified by the filename parameter. It supports both NetCDF and binary formats based on the met_type field in the ctl_t structure. The function can also handle parallel processing with MPI, broadcasting the data across ranks if required by the configuration.
| filename | A constant character pointer representing the name of the file to read the meteorological data from. |
| ctl | A pointer to a ctl_t structure, which holds control parameters including the type of meteorological data, MPI sharing flags, and configuration details. |
| clim | A pointer to a clim_t structure, which contains climatological data to be used in the process, if applicable. |
| met | A pointer to a met_t structure that will store the meteorological data read from the file. |
| dd | A pointer to an dd_t structure containing MPI information, including rank and neighbours. |
met_mpi_share flag is set in the control structure.ctl->met_type is 0, the data is read from a NetCDF file using the read_met_nc function.ctl->met_type is between 1 and 5, or equals 7 or 8, the data is read from a binary file using the read_met_bin function.ctl->met_type is 6, the data is read from grib files using the read_met_grib function.met_type is not recognized, an error message is generated.Definition at line 7802 of file mptrac.c.
| void mptrac_run_timestep | ( | ctl_t * | ctl, |
| cache_t * | cache, | ||
| clim_t * | clim, | ||
| met_t ** | met0, | ||
| met_t ** | met1, | ||
| atm_t * | atm, | ||
| depo_t * | depo, | ||
| double | t, | ||
| dd_t * | dd | ||
| ) |
Executes a single timestep of the MPTRAC model simulation.
This function performs all operations required to advance the model simulation by one timestep. It includes updating air parcel positions, applying advection, diffusion, convection, and other processes such as sedimentation, chemistry, and deposition. Each process is conditionally executed based on the control settings provided in the ctl structure.
| ctl | Pointer to the control structure containing model parameters and settings. |
| cache | Pointer to the cache structure used for intermediate calculations. |
| clim | Pointer to the climatology structure containing climatological data. |
| met0 | Pointer to the current meteorological data structure. |
| met1 | Pointer to the next meteorological data structure. |
| atm | Pointer to the atmosphere structure containing air parcel data. |
| depo | Pointer to cumulative radionuclide ground inventories. |
| t | Current simulation time in seconds. |
| dd | MPI information required for the domain decomposition. |
Definition at line 7911 of file mptrac.c.
| void mptrac_write_atm | ( | const char * | filename, |
| const ctl_t * | ctl, | ||
| const atm_t * | atm, | ||
| const double | t | ||
| ) |
Writes air parcel data to a file in various formats.
The mptrac_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 8177 of file mptrac.c.
Writes meteorological data to a file, supporting multiple formats and compression options.
This function handles writing meteorological data based on the specified control (ctl_t) and meteorological data (met_t) structures. The file format and compression type are determined by the met_type in the control structure. The function supports netCDF, binary output, and various compression methods (ZFP, ZSTD, CMS), while providing error handling for unsupported configurations.
| filename | A constant character pointer representing the name of the file to write the meteorological data to. |
| ctl | A pointer to a ctl_t structure, which holds the configuration and control parameters for the output, including the type of meteorological data and compression method. |
| met | A pointer to a met_t structure that holds the meteorological data to be written to the file. |
ctl->met_type is 3, ZFP compression is required, and the function will generate an error if compiled without ZFP support.ctl->met_type is 4, ZSTD compression is required, and the function will generate an error if compiled without ZSTD support.ctl->met_type is 5, CMS compression is required, and the function will generate an error if compiled without CMS support.ctl->met_type is 7, SZ3 compression is required, and the function will generate an error if compiled without SZ3 support.ctl->met_type is 8, LZ4 compression is required, and the function will generate an error if compiled without LZ4 support.ctl->met_type is 0, the function writes data in netCDF format via write_met_nc.ctl->met_type is between 1 and 5, or equals 7 or 8, the function writes data in binary format via write_met_bin.ctl->met_type is not recognized, an error message is generated.Definition at line 8241 of file mptrac.c.
| void mptrac_write_output | ( | const char * | dirname, |
| const ctl_t * | ctl, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| atm_t * | atm, | ||
| depo_t * | depo, | ||
| const double | t | ||
| ) |
Writes various types of output data to files in a specified directory.
The mptrac_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. |
| depo | A pointer to the cumulative radionuclide ground inventories. |
| 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 8290 of file mptrac.c.
| void mptrac_update_device | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| const clim_t * | clim, | ||
| met_t ** | met0, | ||
| met_t ** | met1, | ||
| const atm_t * | atm | ||
| ) |
Updates device memory for specified data structures.
This function updates the GPU memory with the data from the provided host data structures (ctl, cache, clim, atm) using OpenACC directives. It ensures that the host data is transferred to the device for further computation.
| [in] | ctl | Pointer to the ctl_t structure. If not NULL, the corresponding device memory for ctl is updated. |
| [in] | cache | Pointer to the cache_t structure. If not NULL, the corresponding device memory for cache is updated. |
| [in] | clim | Pointer to the clim_t structure. If not NULL, the corresponding device memory for clim is updated. |
| [in] | met0 | Pointer to the first met_t structure. If not NULL, the corresponding device memory for met0 is updated. |
| [in] | met1 | Pointer to the second met_t structure. If not NULL, the corresponding device memory for met1 is updated. |
| [in] | atm | Pointer to the atm_t structure. If not NULL, the corresponding device memory for atm is updated. |
#pragma acc update directive for device memory synchronization. Each update operation is wrapped with a timer labeled as "UPDATE_DEVICE" for performance tracking.Definition at line 8065 of file mptrac.c.
| void mptrac_update_host | ( | const ctl_t * | ctl, |
| const cache_t * | cache, | ||
| const clim_t * | clim, | ||
| met_t ** | met0, | ||
| met_t ** | met1, | ||
| const atm_t * | atm | ||
| ) |
Updates host memory for specified data structures.
This function transfers data from the device (GPU) memory back to the host memory for the provided data structures (ctl, cache, clim, atm) using OpenACC directives. It ensures that the latest data from the device is synchronized with the host.
| [in] | ctl | Pointer to the ctl_t structure. If not NULL, the corresponding host memory for ctl is updated from the device. |
| [in] | cache | Pointer to the cache_t structure. If not NULL, the corresponding host memory for cache is updated from the device. |
| [in] | clim | Pointer to the clim_t structure. If not NULL, the corresponding host memory for clim is updated from the device. |
| [in] | met0 | Pointer to the first met_t structure. If not NULL, the corresponding host memory for met0 is updated. |
| [in] | met1 | Pointer to the second met_t structure. If not NULL, the corresponding host memory for met1 is updated. |
| [in] | atm | Pointer to the atm_t structure. If not NULL, the corresponding host memory for atm is updated from the device. |
#pragma acc update directive for host memory synchronization. Each update operation is wrapped with a timer labeled as "UPDATE_HOST" for performance tracking.Definition at line 8121 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 8403 of file mptrac.c.
| double pbl_weight | ( | const ctl_t * | ctl, |
| const atm_t * | atm, | ||
| const int | ip, | ||
| const double | pbl, | ||
| const double | ps | ||
| ) |
Computes a weighting factor based on planetary boundary layer pressure.
This function calculates a weighting factor that determines the contribution of a pressure level to processes within the planetary boundary layer. The factor is based on the relative position of the pressure within a linear transition range defined by pbl, ps, and TURB_PBL_TRANS, expressed as a fraction of the PBL pressure thickness (ps - pbl).
| ctl | Pointer to the control structure containing configuration parameters. |
| atm | Pointer to the atmospheric data structure containing pressure levels. |
| ip | Index of the pressure level in the atmospheric data array. |
| pbl | Pressure at the planetary boundary layer. |
| ps | Surface pressure. |
p0).p1).Definition at line 8427 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 8449 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 8491 of file mptrac.c.
Reads atmospheric data from a CLAMS NetCDF file.
This function opens a NetCDF file, reads various atmospheric parameters, and stores them in the provided atm_t structure. It handles both zeta and pressure coordinate systems depending on the control settings.
| [in] | filename | Path to the NetCDF file containing atmospheric data. |
| [in] | ctl | Pointer to the control structure containing configuration settings. |
| [out] | atm | Pointer to the atmospheric data structure where the data will be stored. |
NPARTS).TIME_INIT) or falls back to time if unavailable.ctl->advect_vert_coord, reads ZETA and optionally PRESS, or reads PRESS_INIT with fallback to PRESS.LON) and latitude (LAT).Definition at line 8547 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 8610 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 8643 of file mptrac.c.
| void read_clim_photo_help | ( | const int | ncid, |
| const char * | varname, | ||
| const clim_photo_t * | photo, | ||
| double | var[CP][CSZA][CO3] | ||
| ) |
Reads a 3D climatological photochemistry variable from a NetCDF file.
This function reads a variable from a NetCDF file into a 3D array based on the dimensions provided by the clim_photo_t structure.
| [in] | ncid | NetCDF file ID. |
| [in] | varname | Name of the variable to read from the NetCDF file. |
| [in] | photo | Pointer to a structure defining the data dimensions (np, nsza, no3c). |
| [out] | var | 3D array to store the read data, with dimensions [CP][CSZA][CO3]. |
Definition at line 8734 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 8762 of file mptrac.c.
| void read_clim_zm | ( | const char * | filename, |
| const 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 8816 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 8915 of file mptrac.c.
Reads meteorological data from a binary file.
This function reads meteorological data from a binary file and populates the provided met_t structure with the data. It checks the binary file's format version and met_type, ensuring compatibility with the control structure (ctl_t). The function reads time, grid, surface data, and multi-level data, and supports different binary file versions.
| filename | A constant character pointer representing the name of the binary file to read the meteorological data from. |
| ctl | A pointer to a ctl_t structure that holds control parameters such as the expected met_type and other configuration options. |
| met | A pointer to a met_t structure that will store the meteorological data read from the binary file. |
FREAD macro for safe binary reading operations, which checks the integrity of the read operation.met_type and binary file version to ensure compatibility.LSM, SST, RWC, SWC, and CC).met_type in the file does not match the ctl->met_type.Definition at line 8956 of file mptrac.c.
| void read_met_bin_2d | ( | FILE * | in, |
| const met_t * | met, | ||
| float | var[EX][EY], | ||
| const 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 9116 of file mptrac.c.
| void read_met_bin_3d | ( | FILE * | in, |
| const ctl_t * | ctl, | ||
| const met_t * | met, | ||
| float | var[EX][EY][EP], | ||
| const char * | varname, | ||
| const float | bound_min, | ||
| const 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 9145 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.
| ctl | Pointer to the control structure that contains metadata about the type of data and how it is stored. |
| 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 9255 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 9375 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 9432 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 9539 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 9579 of file mptrac.c.
Reads meteorological data from a grib file and processes it.
This function reads meteorological data from a grib file specified by the filename parameter, using the ECCODES library. It reads grid, surface, and vertical level data, processes the data and calculates various derived meteorological fields such as geopotential heights, potential vorticity, cloud properties, and convective available potential energy (CAPE).
| filename | A constant character pointer representing the name of the grib files to read the meteorological data from. The sf or ml suffixes indicating surface or multi-level data should be replaced with XX. |
| ctl | A pointer to a ctl_t structure, which contains control parameters for reading and processing the meteorological data. |
| met | A pointer to a met_t structure that will store the meteorological data read and processed from the NetCDF file. |
| void read_met_grib_grid | ( | codes_handle ** | handles, |
| int | count_handles, | ||
| met_t * | met | ||
| ) |
Reads global meteorological information from a grib file.
This function reads meteorological grid information from a grib file, including time and spatial dimensions. The function computes the latitude and longitude grid based on the provided boundaries and increment values.
| handles | A pointer to an array of codes_handle pointers representing the grib messages. |
| count_handles | The total number of grib messages in the handles array. |
| met | A pointer to a structure to store meteorological data. |
The function performs the following steps:
| void read_met_grib_levels | ( | codes_handle ** | handles, |
| const int | num_messages, | ||
| const ctl_t * | ctl, | ||
| met_t * | met | ||
| ) |
Reads meteorological variables at different vertical levels from a grib 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 grib file.
| handles | A pointer to an array of codes_handle pointers representing the grib messages. |
| num_messages | The total number of grib messages in the handles array. |
| 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:
| void read_met_grib_surface | ( | codes_handle ** | handles, |
| const int | num_messages, | ||
| const ctl_t * | ctl, | ||
| met_t * | met | ||
| ) |
Reads surface meteorological data from a grib file and stores it in the meteorological data structure.
This function reads various surface meteorological variables from a grib 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.
| handles | A pointer to an array of codes_handle pointers representing the grib messages. |
| num_messages | The total number of grib messages in the handles array. |
| 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:
| void read_met_ml2pl | ( | const ctl_t * | ctl, |
| const met_t * | met, | ||
| float | var[EX][EY][EP], | ||
| const char * | varname | ||
| ) |
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 11314 of file mptrac.c.
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.
| ctl | A pointer to a control parameter structure. |
| met | A pointer to a structure containing meteorological data. |
The function performs the following steps:
Definition at line 11356 of file mptrac.c.
Reads meteorological data from a NetCDF file and processes it.
This function reads meteorological data from a NetCDF file specified by the filename parameter, using the NetCDF library. It reads grid, surface, and vertical level data, processes the data (including extrapolation, boundary conditions, and downsampling), and calculates various derived meteorological fields such as geopotential heights, potential vorticity, cloud properties, and convective available potential energy (CAPE).
| filename | A constant character pointer representing the name of the NetCDF file to read the meteorological data from. |
| ctl | A pointer to a ctl_t structure, which contains control parameters for reading and processing the meteorological data. |
| met | A pointer to a met_t structure that will store the meteorological data read and processed from the NetCDF file. |
| dd | A pointer to an dd_t structure containing MPI information, including rank and neighbours. |
nc_open and handles any errors during the file opening process.Definition at line 11441 of file mptrac.c.
| void read_met_nc_grid | ( | const char * | filename, |
| const int | ncid, | ||
| const ctl_t * | ctl, | ||
| met_t * | met, | ||
| dd_t * | dd | ||
| ) |
Reads meteorological grid data from NetCDF files with domain decomposition.
The read_met_nc_grid function reads meteorological data from NetCDF files and processes it with domain decomposition for parallel processing. It extracts time information, grid dimensions, and coordinates, and sets up hyperslabs for subdomains and halos. It also reads pressure levels and handles model level and surface data.
| filename | A string representing the filename of the NetCDF file to read. |
| ncid | A NetCDF file ID. |
| ctl | A pointer to a ctl_t structure containing control parameters. |
| met | A pointer to a met_t structure where meteorological data will be stored. |
| dd | A pointer to an dd_t structure containing MPI information, including rank and neighbours. |
The function performs the following steps:
Definition at line 9707 of file mptrac.c.
Reads and processes meteorological level data from NetCDF files with domain decomposition.
The read_met_nc_levels function reads meteorological level data from a NetCDF file and processes it for use in a domain decomposition context. It handles various meteorological parameters such as temperature, wind components, humidity, ozone, cloud data, and vertical velocity. The function also processes pressure levels and interpolates data between model and pressure levels as needed.
| ncid | An integer representing the NetCDF file ID. |
| ctl | A pointer to a ctl_t structure containing control parameters and settings. |
| met | A pointer to a met_t structure where meteorological level data will be stored. |
| dd | A pointer to an dd_t structure containing MPI information, including rank and neighbours. |
The function performs the following steps:
Definition at line 10049 of file mptrac.c.
Reads and processes surface meteorological data from NetCDF files with domain decomposition.
The read_met_nc_surface function reads surface meteorological data from a NetCDF file and processes it for use in a domain decomposition context. It handles various surface parameters such as pressure, geopotential height, temperature, wind components, and other relevant meteorological data. The function is designed to work with different meteorological data formats and configurations.
| ncid | An integer representing the NetCDF file ID. |
| ctl | A pointer to a ctl_t structure containing control parameters and settings. |
| met | A pointer to a met_t structure where surface meteorological data will be stored. |
| dd | A pointer to an dd_t structure containing MPI information, including rank and neighbours. |
The function performs the following steps:
zm surface geopotential data from m^2 s^-2 to km when MET_GP2Z = 1.Definition at line 9894 of file mptrac.c.
| int read_met_nc_2d | ( | const int | ncid, |
| const char *const | varnames[], | ||
| const size_t | nvarnames, | ||
| const ctl_t * | ctl, | ||
| const met_t * | met, | ||
| dd_t * | dd, | ||
| float | dest[EX][EY], | ||
| const float | scl, | ||
| const 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. |
| varnames | Array of variable names in order of preference. |
| nvarnames | Number of entries in varnames. |
| ctl | A pointer to a structure containing control parameters. |
| met | A pointer to a structure containing meteorological data. |
| dd | A pointer to an dd_t structure containing MPI information, including rank and neighbours. |
| 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 10267 of file mptrac.c.
| int read_met_nc_3d | ( | const int | ncid, |
| const char *const | varnames[], | ||
| const size_t | nvarnames, | ||
| const ctl_t * | ctl, | ||
| const met_t * | met, | ||
| dd_t * | dd, | ||
| float | dest[EX][EY][EP], | ||
| const 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. |
| varnames | Array of variable names in order of preference. |
| nvarnames | Number of entries in varnames. |
| ctl | A pointer to a structure containing control parameters. |
| met | A pointer to a structure containing meteorological data. |
| dd | A pointer to an dd_t structure containing MPI information, including rank and neighbours. |
| 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 10573 of file mptrac.c.
Computes the planetary boundary layer (PBL) pressure based on meteorological data.
This function determines the PBL pressure for each grid point using one of four methods: 0. Read PBL pressure from meteo data file.
| [in] | ctl | Pointer to the control structure (ctl_t), which contains parameters controlling the PBL calculation. |
| [in,out] | met | Pointer to the meteorological data structure (met_t), which contains grid and atmospheric data. The met->pbl array is updated with the calculated PBL pressure. |
Method 0 (Precomputed PBL pressure from file):
Method 1 (Precomputed PBL height from file):
Method 2 (Richardson number criterion):
Method 3 (Potential temperature difference):
Final Adjustments:
Definition at line 11660 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 11797 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 11858 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 11920 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 12040 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 12069 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 12242 of file mptrac.c.
| void read_obs | ( | const char * | filename, |
| const 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 12416 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 12460 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 12488 of file mptrac.c.
| double scan_ctl | ( | const char * | filename, |
| int | argc, | ||
| char * | argv[], | ||
| const char * | varname, | ||
| const 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 12517 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 12589 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 12622 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:
| 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 12690 of file mptrac.c.
| void timer | ( | const char * | name, |
| const char * | group, | ||
| const 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 12721 of file mptrac.c.
| double time_from_filename | ( | const char * | filename, |
| const int | offset, | ||
| const int | with_seconds | ||
| ) |
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. |
| with_seconds | An integer indicating if the filename contains seconds or not. |
The function performs the following steps:
time2jsec function.YYYY-MM-DD_HH-MM (e.g., "2023-05-27_14-45")YYYY-MM-DD_HH-MM-SS (e.g., "2023-05-27_14-45-12")Definition at line 12789 of file mptrac.c.
Computes a weighting factor based on tropopause pressure.
This function calculates a weighting factor for a given pressure value in relation to the tropopause pressure. The weighting factor is determined as follows:
| [in] | ctl | Pointer to the control data structure. |
| [in] | clim | Pointer to the climatology data structure. |
| [in] | atm | Pointer to the atmospheric data structure. |
| [in] | ip | Index of the pressure value to evaluate within the atmospheric data. |
Definition at line 12831 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 12857 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 12955 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:
NPARTS).Definition at line 13005 of file mptrac.c.
| void write_atm_clams_traj | ( | const char * | dirname, |
| const ctl_t * | ctl, | ||
| const atm_t * | atm, | ||
| const double | t | ||
| ) |
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 13061 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 13222 of file mptrac.c.
Writes a time series of the atmospheric tracer mass budget.
The output contains the number of active air parcels, their total mass, and sums of the configured cumulative mass-loss quantities. If ensembles are enabled, one row per ensemble and an additional total row are written. Domain-decomposed runs report the local budget of each subdomain. Missing mass-loss quantities are reported as NaN.
| filename | Path to the budget output file. |
| ctl | Pointer to the control structure. |
| atm | Pointer to the atmospheric data structure. |
| t | Current model time [s]. |
Definition at line 13271 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 13365 of file mptrac.c.
Writes ensemble-based Critical Success Index (CSI) and other verification statistics to an output file.
This function computes and writes various statistical verification metrics that assess the performance of ensemble forecasts compared to observations. The output includes, for each ensemble member:
On the first invocation (when t == ctl->t_start), the function loads observation data and kernel weights, allocates necessary arrays, and creates a new output file with a descriptive header. On the final call (when t == ctl->t_stop), the function closes the file and frees all persistent memory.
Output is written at time steps divisible by ctl->csi_dt_out.
| [in] | filename | Path to the output file where statistics will be written. |
| [in] | ctl | Pointer to control structure containing configuration, ensemble settings, and spatial/time grid info. |
| [in] | atm | Pointer to atmospheric data structure holding model output. |
| [in] | t | Current simulation/model time. |
| If | required quantities (mass or ensemble IDs) are undefined, if ensemble IDs are out of bounds, or if the output file cannot be created. |
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 13652 of file mptrac.c.
Convert cumulative ground inventories to Bq m^-2 and write them.
Cell centers and areas are derived from GRID_LON*, GRID_LAT*, GRID_NX, and GRID_NY. If RADIO_DECAY is enabled, inventories are decayed to t. Output format is selected by DEPO_TYPE.
| [in] | filename | Output filename. |
| [in] | ctl | Grid, format, compression, and decay settings. |
| [in] | depo | Cumulative radionuclide inventories. |
| [in] | t | Output validity time [s since 2000-01-01 00:00 UTC]. |
Definition at line 13752 of file mptrac.c.
| void write_depo_asc | ( | const char * | filename, |
| const ctl_t * | ctl, | ||
| const double * | data, | ||
| const double | t, | ||
| const double * | lon, | ||
| const double * | lat, | ||
| const double * | area | ||
| ) |
Write radioactive deposition densities as a gnuplot-compatible table.
| [in] | filename | Output filename. |
| [in] | ctl | Deposition grid dimensions. |
| [in] | data | Four consecutive fields in Bq m^-2. |
| [in] | t | Output validity time [s since 2000-01-01 00:00 UTC]. |
| [in] | lon | Grid-cell center longitudes [deg]. |
| [in] | lat | Grid-cell center latitudes [deg]. |
| [in] | area | Grid-cell areas by latitude [m^2]. |
Definition at line 13823 of file mptrac.c.
| void write_depo_nc | ( | const char * | filename, |
| const ctl_t * | ctl, | ||
| const double * | data, | ||
| const double | t, | ||
| const double * | lon, | ||
| const double * | lat, | ||
| const double * | area | ||
| ) |
Write radioactive deposition densities as a CF-style netCDF file.
| [in] | filename | Output filename. |
| [in] | ctl | Deposition grid dimensions and netCDF compression setting. |
| [in] | data | Four consecutive fields in Bq m^-2. |
| [in] | t | Output validity time [s since 2000-01-01 00:00 UTC]. |
| [in] | lon | Grid-cell center longitudes [deg]. |
| [in] | lat | Grid-cell center latitudes [deg]. |
| [in] | area | Grid-cell areas by latitude [m^2]. |
Definition at line 13863 of file mptrac.c.
| void write_grid | ( | const char * | filename, |
| const ctl_t * | ctl, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| const atm_t * | atm, | ||
| const 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 13929 of file mptrac.c.
| void write_grid_asc | ( | const char * | filename, |
| const ctl_t * | ctl, | ||
| const double * | cd, | ||
| double * | mean[NQ], | ||
| double * | sigma[NQ], | ||
| const double * | vmr_impl, | ||
| const double | t, | ||
| const double * | z, | ||
| const double * | lon, | ||
| const double * | lat, | ||
| const double * | area, | ||
| const double | dz, | ||
| const 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 14131 of file mptrac.c.
| void write_grid_nc | ( | const char * | filename, |
| const ctl_t * | ctl, | ||
| const double * | cd, | ||
| double * | mean[NQ], | ||
| double * | sigma[NQ], | ||
| const double * | vmr_impl, | ||
| const double | t, | ||
| const double * | z, | ||
| const double * | lon, | ||
| const double * | lat, | ||
| const double * | area, | ||
| const double | dz, | ||
| const 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 14235 of file mptrac.c.
Writes meteorological data in binary format to a specified file.
This function writes meteorological data from the met_t structure to a binary file. The data includes grid and surface data, as well as multi-level data such as temperature, velocity components, and atmospheric properties. The compression options for multi-level data (ZFP) are controlled via the ctl_t structure. The function supports multiple variables, such as surface pressure, temperature, wind components, and cloud properties.
| filename | A constant character pointer representing the name of the file to write the binary data to. |
| ctl | A pointer to a ctl_t structure, which holds control parameters including the type of meteorological data, compression settings, and grid dimensions. |
| met | A pointer to a met_t structure that contains the meteorological data to be written to the binary file. |
ctl->met_type) and the version of the binary format are written at the beginning of the file.write_met_bin_2d helper function.write_met_bin_3d function with optional ZFP compression settings.Definition at line 14365 of file mptrac.c.
| void write_met_bin_2d | ( | FILE * | out, |
| met_t * | met, | ||
| float | var[EX][EY], | ||
| const 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 14495 of file mptrac.c.
| void write_met_bin_3d | ( | FILE * | out, |
| const ctl_t * | ctl, | ||
| met_t * | met, | ||
| float | var[EX][EY][EP], | ||
| const char * | varname, | ||
| const int | metvar, | ||
| FILE * | level_log | ||
| ) |
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. |
| metvar | Index of the meteorological variable, used to select codec-specific compression settings. |
| level_log | Optional file stream for per-level compression diagnostics, or NULL to disable logging. |
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 14524 of file mptrac.c.
Writes meteorological data to a NetCDF file.
This function creates and writes meteorological data to a NetCDF file in the NetCDF-4 format. It defines the required dimensions, grid, surface variables, and level data within the NetCDF structure and writes the corresponding values from the met_t structure. The function uses helper functions to write 2D surface and 3D level data.
| filename | A constant character pointer representing the name of the NetCDF file to create and write the data to. |
| ctl | A pointer to a ctl_t structure that contains control parameters, such as the NetCDF level and quantization settings. |
| met | A pointer to a met_t structure that contains the meteorological data to be written to the NetCDF file. |
Definition at line 14617 of file mptrac.c.
| void write_met_nc_2d | ( | const int | ncid, |
| const char * | varname, | ||
| met_t * | met, | ||
| float | var[EX][EY], | ||
| const float | scl | ||
| ) |
Writes a 2D meteorological variable to a NetCDF file.
This function writes a 2D meteorological variable, stored in the array var, to a NetCDF file with the specified variable name. The data is scaled by a factor scl before being written. The function handles memory allocation for the data copy, scaling, and freeing the allocated memory after writing the data to the NetCDF file.
| ncid | The NetCDF file ID. This is an integer that identifies the NetCDF file where the data will be written. It is assumed that this file has already been opened for writing. |
| varname | A pointer to a string containing the name of the variable in the NetCDF file where the data will be stored. |
| met | A pointer to a structure of type met_t that contains metadata about the meteorological field, including the dimensions nx (number of points in x-direction) and ny (number of points in y-direction). |
| var | A 2D array of dimensions EX x EY containing the meteorological data to be written. The data is provided in the format var[ix][iy], where ix is the index in the x-direction and iy is the index in the y-direction. |
| scl | A scaling factor applied to each element in the var array before writing to the NetCDF file. |
Definition at line 14799 of file mptrac.c.
| void write_met_nc_3d | ( | const int | ncid, |
| const char * | varname, | ||
| met_t * | met, | ||
| float | var[EX][EY][EP], | ||
| const float | scl | ||
| ) |
Writes a 3D meteorological variable to a NetCDF file.
This function writes a 3D meteorological variable, stored in the array var, to a NetCDF file with the specified variable name. The data is scaled by a factor scl before being written. The function handles memory allocation for the data copy, scaling, and freeing the allocated memory after writing the data to the NetCDF file.
| ncid | The NetCDF file ID. This is an integer that identifies the NetCDF file where the data will be written. It is assumed that this file has already been opened for writing. |
| varname | A pointer to a string containing the name of the variable in the NetCDF file where the data will be stored. |
| met | A pointer to a structure of type met_t that contains metadata about the meteorological field, including the dimensions nx (number of points in the x-direction), ny (number of points in the y-direction), and np (number of points in the third dimension, e.g., pressure levels). |
| var | A 3D array of dimensions EX x EY x EP containing the meteorological data to be written. The data is provided in the format var[ix][iy][ip], where ix is the index in the x-direction, iy is the index in the y-direction, and ip is the index in the third dimension (e.g., vertical levels). |
| scl | A scaling factor applied to each element in the var array before writing to the NetCDF file. |
Definition at line 14829 of file mptrac.c.
| void write_prof | ( | const char * | filename, |
| const ctl_t * | ctl, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| const atm_t * | atm, | ||
| const 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 14860 of file mptrac.c.
| void write_sample | ( | const char * | filename, |
| const ctl_t * | ctl, | ||
| met_t * | met0, | ||
| met_t * | met1, | ||
| const atm_t * | atm, | ||
| const 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 15095 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 15260 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 15349 of file mptrac.c.
1.9.4