52 static double timem[
NZ], z, lon, lonm[
NZ], lat, latm[
NZ], t, tm[
NZ], u,
53 um[
NZ], v, vm[
NZ], w, wm[
NZ], h2o, h2om[
NZ], h2ot, h2otm[
NZ], o3, o3m[
NZ],
54 lwc, lwcm[
NZ], rwc, rwcm[
NZ], iwc, iwcm[
NZ], swc, swcm[
NZ], cc, ccm[
NZ],
55 ps, psm[
NZ], ts, tsm[
NZ], zs, zsm[
NZ], us, usm[
NZ], vs, vsm[
NZ],
56 ess, essm[
NZ], nss, nssm[
NZ], shf, shfm[
NZ], lsm, lsmm[
NZ],
57 sst, sstm[
NZ], pbl, pblm[
NZ], pt, ptm[
NZ], pct, pctm[
NZ], pcb,
58 pcbm[
NZ], cl, clm[
NZ], plcl, plclm[
NZ], plfc, plfcm[
NZ], pel, pelm[
NZ],
59 cape, capem[
NZ], cin, cinm[
NZ], o3c, o3cm[
NZ], tt, ttm[
NZ], zm[
NZ], zt,
60 ztm[
NZ], pv, pvm[
NZ], plev[
NZ], rhm[
NZ], rhicem[
NZ], tdewm[
NZ], ticem[
NZ],
61 tnatm[
NZ], hno3m[
NZ], ohm[
NZ], h2o2m[
NZ], ho2m[
NZ], o1dm[
NZ];
63 static int iz, np[
NZ], npc[
NZ], npt[
NZ], nz;
72 ERRMSG(
"Give parameters: <ctl> <prof.tab> <met0> [ <met1> ... ]");
76 double z0 =
scan_ctl(argv[1], argc, argv,
"PROF_Z0", -1,
"-999", NULL);
77 double z1 =
scan_ctl(argv[1], argc, argv,
"PROF_Z1", -1,
"-999", NULL);
78 double dz =
scan_ctl(argv[1], argc, argv,
"PROF_DZ", -1,
"-999", NULL);
79 double lon0 =
scan_ctl(argv[1], argc, argv,
"PROF_LON0", -1,
"0", NULL);
80 double lon1 =
scan_ctl(argv[1], argc, argv,
"PROF_LON1", -1,
"0", NULL);
81 double dlon =
scan_ctl(argv[1], argc, argv,
"PROF_DLON", -1,
"-999", NULL);
82 double lat0 =
scan_ctl(argv[1], argc, argv,
"PROF_LAT0", -1,
"0", NULL);
83 double lat1 =
scan_ctl(argv[1], argc, argv,
"PROF_LAT1", -1,
"0", NULL);
84 double dlat =
scan_ctl(argv[1], argc, argv,
"PROF_DLAT", -1,
"-999", NULL);
90 for (
int i = 3; i < argc; i++) {
100 z1 =
Z(met->
p[met->
np - 1]);
103 for (iz = 0; iz < met->
np; iz++)
104 if (
Z(met->
p[iz]) >= z0 &&
Z(met->
p[iz]) <= z1) {
105 plev[nz] = met->
p[iz];
107 ERRMSG(
"Too many pressure levels!");
110 for (z = z0; z <= z1; z += dz) {
113 ERRMSG(
"Too many pressure levels!");
118 dlon = fabs(met->
lon[1] - met->
lon[0]);
120 dlat = fabs(met->
lat[1] - met->
lat[0]);
123 for (iz = 0; iz < nz; iz++)
124 for (lon = lon0; lon <= lon1; lon += dlon)
125 for (lat = lat0; lat <= lat1; lat += dlat) {
132 if (isfinite(t) && isfinite(u) && isfinite(v) && isfinite(w)) {
133 timem[iz] += met->
time;
163 if (isfinite(plfc) && isfinite(pel) && cape >= ctl.
conv_cape
180 rhm[iz] +=
RH(plev[iz], t, h2o);
181 rhicem[iz] +=
RHICE(plev[iz], t, h2o);
182 tdewm[iz] +=
TDEW(plev[iz], h2o);
183 ticem[iz] +=
TICE(plev[iz], h2o);
188 ohm[iz] +=
clim_oh(&ctl, clim, met->
time, lon, lat, plev[iz]);
198 LOG(1,
"Write meteorological data file: %s", argv[2]);
199 if (!(out = fopen(argv[2],
"w")))
200 ERRMSG(
"Cannot create file!");
207 for (iz = 0; iz < nz; iz++)
209 "%.2f %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g"
210 " %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g"
211 " %g %g %g %g %g %g %g %g %g %g %g %g %g %d %d %d\n",
212 timem[iz] / np[iz],
Z(plev[iz]), lonm[iz] / np[iz],
213 latm[iz] / np[iz], plev[iz], tm[iz] / np[iz], um[iz] / np[iz],
214 vm[iz] / np[iz], wm[iz] / np[iz], h2om[iz] / np[iz],
215 o3m[iz] / np[iz], zm[iz] / np[iz], pvm[iz] / np[iz],
216 psm[iz] / np[iz], tsm[iz] / np[iz], zsm[iz] / np[iz],
217 usm[iz] / np[iz], vsm[iz] / np[iz], essm[iz] / np[iz],
218 nssm[iz] / np[iz], shfm[iz] / np[iz], lsmm[iz] / np[iz],
219 sstm[iz] / np[iz], ptm[iz] / npt[iz], ztm[iz] / npt[iz],
220 ttm[iz] / npt[iz], h2otm[iz] / npt[iz],
221 lwcm[iz] / np[iz], rwcm[iz] / np[iz], iwcm[iz] / np[iz],
222 swcm[iz] / np[iz], ccm[iz] / np[iz], clm[iz] / np[iz],
223 pctm[iz] / np[iz], pcbm[iz] / np[iz], plclm[iz] / npc[iz],
224 plfcm[iz] / npc[iz], pelm[iz] / npc[iz], capem[iz] / npc[iz],
225 cinm[iz] / npc[iz], rhm[iz] / np[iz], rhicem[iz] / np[iz],
226 tdewm[iz] / np[iz], ticem[iz] / np[iz], tnatm[iz] / np[iz],
227 hno3m[iz] / np[iz], ohm[iz] / np[iz], h2o2m[iz] / np[iz],
228 ho2m[iz] / np[iz], o1dm[iz] / np[iz], pblm[iz] / np[iz],
229 o3cm[iz] / np[iz], np[iz], npt[iz], npc[iz]);
int main(int argc, char *argv[])
#define NZ
Maximum number of altitudes.
double clim_zm(const clim_zm_t *zm, const double t, const double lat, const double p)
Interpolates monthly mean zonal mean climatological variables.
double nat_temperature(const double p, const double h2o, const double hno3)
Calculates the nitric acid trihydrate (NAT) temperature.
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.
void mptrac_read_clim(const ctl_t *ctl, clim_t *clim)
Reads various climatological data and populates the given climatology structure.
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 co...
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.
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.
MPTRAC library declarations.
#define INTPOL_SPACE_ALL(p, lon, lat)
Interpolate multiple meteorological variables in space.
#define INTPOL_INIT
Initialize arrays for interpolation.
#define ERRMSG(...)
Print an error message with contextual information and terminate the program.
#define Z(p)
Convert pressure to altitude.
#define P(z)
Compute pressure at given altitude.
#define MET_HEADER
Write header for meteorological data file.
#define TICE(p, h2o)
Calculate frost point temperature (WMO, 2018).
#define RHICE(p, t, h2o)
Compute relative humidity over ice.
#define ALLOC(ptr, type, n)
Allocate memory for a pointer with error handling.
#define RH(p, t, h2o)
Compute relative humidity over water.
#define LOG(level,...)
Print a log message with a specified logging level.
#define TDEW(p, h2o)
Calculate dew point temperature.
clim_zm_t ho2
HO2 zonal means.
clim_zm_t hno3
HNO3 zonal means.
clim_zm_t o1d
O(1D) zonal means.
clim_zm_t h2o2
H2O2 zonal means.
double conv_cape
CAPE threshold for convection module [J/kg].
double conv_cin
CIN threshold for convection module [J/kg].
Domain decomposition data structure.
int np
Number of pressure levels.
double lon[EX]
Longitudes [deg].
double lat[EY]
Latitudes [deg].
double p[EP]
Pressure levels [hPa].