MPTRAC
Macros | Functions
met_zm.c File Reference

Extract zonal mean from meteorological data. More...

#include "mptrac.h"

Go to the source code of this file.

Macros

#define NZ   1000
 Maximum number of altitudes. More...
 
#define NY   EY
 Maximum number of latitudes. More...
 

Functions

void usage (void)
 Print command-line help. More...
 
int main (int argc, char *argv[])
 

Detailed Description

Extract zonal mean from meteorological data.

Definition in file met_zm.c.

Macro Definition Documentation

◆ NZ

#define NZ   1000

Maximum number of altitudes.

Definition at line 32 of file met_zm.c.

◆ NY

#define NY   EY

Maximum number of latitudes.

Definition at line 35 of file met_zm.c.

Function Documentation

◆ usage()

void usage ( void  )

Print command-line help.

Definition at line 286 of file met_zm.c.

287 {
288
289 printf("\nMPTRAC met_zm tool.\n\n");
290 printf("Extract zonal means from meteorological data.\n");
291 printf("\n");
292 printf("Usage:\n");
293 printf(" met_zm <ctl> <zm.tab> <met0> [<met1> ...]\n");
294 printf("\n");
295 printf("Arguments:\n");
296 printf(" <ctl> Control file.\n");
297 printf(" <zm.tab> Output table.\n");
298 printf(" <met*> Meteorological input files.\n");
299 printf("\nFurther information:\n");
300 printf(" Manual: https://slcs-jsc.github.io/mptrac/\n");
301}

◆ main()

int main ( int  argc,
char *  argv[] 
)

Definition at line 49 of file met_zm.c.

51 {
52
53 ctl_t ctl;
54
55 clim_t *clim;
56
57 met_t *met;
58
59 dd_t *dd;
60
61 FILE *out;
62
63 static double timem[NZ][NY], psm[NZ][NY], tsm[NZ][NY], zsm[NZ][NY],
64 usm[NZ][NY], vsm[NZ][NY], essm[NZ][NY], nssm[NZ][NY], shfm[NZ][NY],
65 lsmm[NZ][NY], sstm[NZ][NY], pblm[NZ][NY],
66 ptm[NZ][NY], pctm[NZ][NY], pcbm[NZ][NY], clm[NZ][NY], plclm[NZ][NY],
67 plfcm[NZ][NY], pelm[NZ][NY], capem[NZ][NY], cinm[NZ][NY], o3cm[NZ][NY],
68 ttm[NZ][NY], ztm[NZ][NY], tm[NZ][NY], um[NZ][NY], vm[NZ][NY], wm[NZ][NY],
69 h2om[NZ][NY], h2otm[NZ][NY], pvm[NZ][NY], o3m[NZ][NY], lwcm[NZ][NY],
70 rwcm[NZ][NY], iwcm[NZ][NY], swcm[NZ][NY], ccm[NZ][NY], zm[NZ][NY],
71 rhm[NZ][NY], rhicem[NZ][NY], tdewm[NZ][NY], ticem[NZ][NY], tnatm[NZ][NY],
72 hno3m[NZ][NY], ohm[NZ][NY], h2o2m[NZ][NY], ho2m[NZ][NY], o1dm[NZ][NY], z,
73 zt, tt, plev[NZ], ps, ts, zs, us, vs, ess, nss, shf, lsm, sst, pbl,
74 pt, pct, pcb, plcl, plfc, pel, cape, cin, o3c, cl, t, u, v, w, pv,
75 h2o, h2ot, o3, lwc, rwc, iwc, swc, cc, lat, lats[NY], lonm[NZ][NY], cw[3];
76
77 static int np[NZ][NY], npc[NZ][NY], npt[NZ][NY], ny, nz, ci[3];
78
79 /* Allocate... */
80 ALLOC(clim, clim_t, 1);
81 ALLOC(met, met_t, 1);
82 ALLOC(dd, dd_t, 1);
83
84 /* Print usage information... */
85 USAGE;
86
87 /* Check arguments... */
88 if (argc < 4)
89 ERRMSG("Missing or invalid command-line arguments.\n\n"
90 "Usage: met_zm <ctl> <zm.tab> <met0> [<met1> ...]\n\n"
91 "Use -h for full help.");
92
93 /* Read control parameters... */
94 mptrac_read_ctl(argv[1], argc, argv, &ctl);
95 double z0 = scan_ctl(argv[1], argc, argv, "ZM_Z0", -1, "-999", NULL);
96 double z1 = scan_ctl(argv[1], argc, argv, "ZM_Z1", -1, "-999", NULL);
97 double dz = scan_ctl(argv[1], argc, argv, "ZM_DZ", -1, "-999", NULL);
98 double lon0 = scan_ctl(argv[1], argc, argv, "ZM_LON0", -1, "-360", NULL);
99 double lon1 = scan_ctl(argv[1], argc, argv, "ZM_LON1", -1, "360", NULL);
100 double lat0 = scan_ctl(argv[1], argc, argv, "ZM_LAT0", -1, "-90", NULL);
101 double lat1 = scan_ctl(argv[1], argc, argv, "ZM_LAT1", -1, "90", NULL);
102 double dlat = scan_ctl(argv[1], argc, argv, "ZM_DLAT", -1, "-999", NULL);
103
104 /* Read climatological data... */
105 mptrac_read_clim(&ctl, clim);
106
107 /* Loop over files... */
108 for (int i = 3; i < argc; i++) {
109
110 /* Read meteorological data... */
111 if (!mptrac_read_met(argv[i], &ctl, clim, met, dd))
112 continue;
113
114 /* Set vertical grid... */
115 if (z0 < 0)
116 z0 = Z(met->p[0]);
117 if (z1 < 0)
118 z1 = Z(met->p[met->np - 1]);
119 nz = 0;
120 if (dz < 0) {
121 for (int iz = 0; iz < met->np; iz++)
122 if (Z(met->p[iz]) >= z0 && Z(met->p[iz]) <= z1) {
123 plev[nz] = met->p[iz];
124 if ((++nz) >= NZ)
125 ERRMSG("Too many pressure levels!");
126 }
127 } else
128 for (z = z0; z <= z1; z += dz) {
129 plev[nz] = P(z);
130 if ((++nz) >= NZ)
131 ERRMSG("Too many pressure levels!");
132 }
133
134 /* Set horizontal grid... */
135 if (dlat <= 0)
136 dlat = fabs(met->lat[1] - met->lat[0]);
137 ny = 0;
138 if (lat0 < -90 && lat1 > 90) {
139 lat0 = gsl_stats_min(met->lat, 1, (size_t) met->ny);
140 lat1 = gsl_stats_max(met->lat, 1, (size_t) met->ny);
141 }
142 for (lat = lat0; lat <= lat1 + 0.001; lat += dlat) {
143 lats[ny] = round(lat * 1e3) / 1e3;
144 if ((++ny) >= NY)
145 ERRMSG("Too many latitudes!");
146 }
147
148 /* Average... */
149 for (int ix = 0; ix < met->nx; ix++)
150 if (met->lon[ix] >= lon0 && met->lon[ix] <= lon1)
151 for (int iy = 0; iy < ny; iy++)
152 for (int iz = 0; iz < nz; iz++) {
153
154 /* Interpolate meteo data... */
155 INTPOL_SPACE_ALL(plev[iz], met->lon[ix], lats[iy]);
156
157 /* Averaging... */
158 timem[iz][iy] += met->time;
159 lonm[iz][iy] += met->lon[ix];
160 zm[iz][iy] += z;
161 tm[iz][iy] += t;
162 um[iz][iy] += u;
163 vm[iz][iy] += v;
164 wm[iz][iy] += w;
165 pvm[iz][iy] += pv;
166 h2om[iz][iy] += h2o;
167 o3m[iz][iy] += o3;
168 lwcm[iz][iy] += lwc;
169 rwcm[iz][iy] += rwc;
170 iwcm[iz][iy] += iwc;
171 swcm[iz][iy] += swc;
172 ccm[iz][iy] += cc;
173 psm[iz][iy] += ps;
174 tsm[iz][iy] += ts;
175 zsm[iz][iy] += zs;
176 usm[iz][iy] += us;
177 vsm[iz][iy] += vs;
178 essm[iz][iy] += ess;
179 nssm[iz][iy] += nss;
180 shfm[iz][iy] += shf;
181 lsmm[iz][iy] += lsm;
182 sstm[iz][iy] += sst;
183 pblm[iz][iy] += pbl;
184 pctm[iz][iy] += pct;
185 pcbm[iz][iy] += pcb;
186 clm[iz][iy] += cl;
187 if (isfinite(plfc) && isfinite(pel) && cape >= ctl.conv_cape
188 && (ctl.conv_cin <= 0 || cin < ctl.conv_cin)) {
189 plclm[iz][iy] += plcl;
190 plfcm[iz][iy] += plfc;
191 pelm[iz][iy] += pel;
192 capem[iz][iy] += cape;
193 cinm[iz][iy] += cin;
194 npc[iz][iy]++;
195 }
196 if (isfinite(pt)) {
197 ptm[iz][iy] += pt;
198 ztm[iz][iy] += zt;
199 ttm[iz][iy] += tt;
200 h2otm[iz][iy] += h2ot;
201 npt[iz][iy]++;
202 }
203 o3cm[iz][iy] += o3c;
204 rhm[iz][iy] += RH(plev[iz], t, h2o);
205 rhicem[iz][iy] += RHICE(plev[iz], t, h2o);
206 tdewm[iz][iy] += TDEW(plev[iz], h2o);
207 ticem[iz][iy] += TICE(plev[iz], h2o);
208 hno3m[iz][iy] +=
209 clim_zm(&clim->hno3, met->time, lats[iy], plev[iz]);
210 tnatm[iz][iy] +=
211 nat_temperature(plev[iz], h2o,
212 clim_zm(&clim->hno3, met->time, lats[iy],
213 plev[iz]));
214 ohm[iz][iy] +=
215 clim_oh(&ctl, clim, met->time, met->lon[ix], lats[iy],
216 plev[iz]);
217 h2o2m[iz][iy]
218 += clim_zm(&clim->h2o2, met->time, lats[iy], plev[iz]);
219 ho2m[iz][iy]
220 += clim_zm(&clim->ho2, met->time, lats[iy], plev[iz]);
221 o1dm[iz][iy]
222 += clim_zm(&clim->o1d, met->time, lats[iy], plev[iz]);
223 np[iz][iy]++;
224 }
225 }
226
227 /* Create output file... */
228 LOG(1, "Write meteorological data file: %s", argv[2]);
229 if (!(out = fopen(argv[2], "w")))
230 ERRMSG("Cannot create file!");
231
232 /* Write header... */
233 MET_HEADER;
234
235 /* Write data... */
236 for (int iz = 0; iz < nz; iz++) {
237 fprintf(out, "\n");
238 for (int iy = 0; iy < ny; iy++)
239 fprintf(out,
240 "%.2f %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g"
241 " %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g"
242 " %g %g %g %g %g %g %g %g %g %g %g %g %g %d %d %d\n",
243 timem[iz][iy] / np[iz][iy], Z(plev[iz]),
244 lonm[iz][iy] / np[iz][iy], lats[iy],
245 plev[iz], tm[iz][iy] / np[iz][iy], um[iz][iy] / np[iz][iy],
246 vm[iz][iy] / np[iz][iy], wm[iz][iy] / np[iz][iy],
247 h2om[iz][iy] / np[iz][iy], o3m[iz][iy] / np[iz][iy],
248 zm[iz][iy] / np[iz][iy], pvm[iz][iy] / np[iz][iy],
249 psm[iz][iy] / np[iz][iy], tsm[iz][iy] / np[iz][iy],
250 zsm[iz][iy] / np[iz][iy], usm[iz][iy] / np[iz][iy],
251 vsm[iz][iy] / np[iz][iy], essm[iz][iy] / np[iz][iy],
252 nssm[iz][iy] / np[iz][iy], shfm[iz][iy] / np[iz][iy],
253 lsmm[iz][iy] / np[iz][iy],
254 sstm[iz][iy] / np[iz][iy], ptm[iz][iy] / npt[iz][iy],
255 ztm[iz][iy] / npt[iz][iy], ttm[iz][iy] / npt[iz][iy],
256 h2otm[iz][iy] / npt[iz][iy], lwcm[iz][iy] / np[iz][iy],
257 rwcm[iz][iy] / np[iz][iy], iwcm[iz][iy] / np[iz][iy],
258 swcm[iz][iy] / np[iz][iy], ccm[iz][iy] / np[iz][iy],
259 clm[iz][iy] / np[iz][iy], pctm[iz][iy] / np[iz][iy],
260 pcbm[iz][iy] / np[iz][iy], plclm[iz][iy] / npc[iz][iy],
261 plfcm[iz][iy] / npc[iz][iy], pelm[iz][iy] / npc[iz][iy],
262 capem[iz][iy] / npc[iz][iy], cinm[iz][iy] / npc[iz][iy],
263 rhm[iz][iy] / np[iz][iy], rhicem[iz][iy] / np[iz][iy],
264 tdewm[iz][iy] / np[iz][iy], ticem[iz][iy] / np[iz][iy],
265 tnatm[iz][iy] / np[iz][iy], hno3m[iz][iy] / np[iz][iy],
266 ohm[iz][iy] / np[iz][iy], h2o2m[iz][iy] / np[iz][iy],
267 ho2m[iz][iy] / np[iz][iy], o1dm[iz][iy] / np[iz][iy],
268 pblm[iz][iy] / np[iz][iy], o3cm[iz][iy] / np[iz][iy],
269 np[iz][iy], npt[iz][iy], npc[iz][iy]);
270 }
271
272 /* Close file... */
273 fclose(out);
274
275 /* Free... */
276 free(clim);
277 free(met);
278 free(dd);
279
280 return EXIT_SUCCESS;
281}
NY
#define NY
Maximum number of latitudes.
Definition: met_zm.c:35
NZ
#define NZ
Maximum number of altitudes.
Definition: met_zm.c:32
clim_zm
double clim_zm(const clim_zm_t *zm, const double t, const double lat, const double p)
Interpolates monthly mean zonal mean climatological variables.
Definition: mptrac.c:405
nat_temperature
double nat_temperature(const double p, const double h2o, const double hno3)
Calculates the nitric acid trihydrate (NAT) temperature.
Definition: mptrac.c:6721
scan_ctl
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.
Definition: mptrac.c:10745
mptrac_read_clim
void mptrac_read_clim(const ctl_t *ctl, clim_t *clim)
Reads various climatological data and populates the given climatology structure.
Definition: mptrac.c:5188
clim_oh
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...
Definition: mptrac.c:89
mptrac_read_met
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.
Definition: mptrac.c:6151
mptrac_read_ctl
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.
Definition: mptrac.c:5248
INTPOL_SPACE_ALL
#define INTPOL_SPACE_ALL(p, lon, lat)
Interpolate multiple meteorological variables in space.
Definition: mptrac.h:920
ERRMSG
#define ERRMSG(...)
Print an error message with contextual information and terminate the program.
Definition: mptrac.h:2102
USAGE
#define USAGE
Print usage information on -h or --help.
Definition: mptrac.h:1909
Z
#define Z(p)
Convert pressure to altitude.
Definition: mptrac.h:1939
P
#define P(z)
Compute pressure at given altitude.
Definition: mptrac.h:1480
MET_HEADER
#define MET_HEADER
Write header for meteorological data file.
Definition: mptrac.h:1088
TICE
#define TICE(p, h2o)
Calculate frost point temperature (WMO, 2018).
Definition: mptrac.h:1796
RHICE
#define RHICE(p, t, h2o)
Compute relative humidity over ice.
Definition: mptrac.h:1632
ALLOC
#define ALLOC(ptr, type, n)
Allocate memory for a pointer with error handling.
Definition: mptrac.h:453
RH
#define RH(p, t, h2o)
Compute relative humidity over water.
Definition: mptrac.h:1602
LOG
#define LOG(level,...)
Print a log message with a specified logging level.
Definition: mptrac.h:2032
TDEW
#define TDEW(p, h2o)
Calculate dew point temperature.
Definition: mptrac.h:1771
clim_t
Climatological data.
Definition: mptrac.h:3436
clim_t::ho2
clim_zm_t ho2
HO2 zonal means.
Definition: mptrac.h:3466
clim_t::hno3
clim_zm_t hno3
HNO3 zonal means.
Definition: mptrac.h:3457
clim_t::o1d
clim_zm_t o1d
O(1D) zonal means.
Definition: mptrac.h:3469
clim_t::h2o2
clim_zm_t h2o2
H2O2 zonal means.
Definition: mptrac.h:3463
ctl_t
Control parameters.
Definition: mptrac.h:2190
ctl_t::conv_cape
double conv_cape
CAPE threshold for convection module [J/kg].
Definition: mptrac.h:2765
ctl_t::conv_cin
double conv_cin
CIN threshold for convection module [J/kg].
Definition: mptrac.h:2768
dd_t
Domain decomposition data structure.
Definition: mptrac.h:3669
met_t
Meteo data structure.
Definition: mptrac.h:3495
met_t::nx
int nx
Number of longitudes.
Definition: mptrac.h:3501
met_t::ny
int ny
Number of latitudes.
Definition: mptrac.h:3504
met_t::np
int np
Number of pressure levels.
Definition: mptrac.h:3507
met_t::lon
double lon[EX]
Longitudes [deg].
Definition: mptrac.h:3513
met_t::time
double time
Time [s].
Definition: mptrac.h:3498
met_t::lat
double lat[EY]
Latitudes [deg].
Definition: mptrac.h:3516
met_t::p
double p[EP]
Pressure levels [hPa].
Definition: mptrac.h:3519
Here is the call graph for this function:
Generated by doxygen 1.9.4