MPTRAC
Functions
tropo.c File Reference

Create tropopause data set from meteorological data. More...

#include "mptrac.h"

Go to the source code of this file.

Functions

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

Detailed Description

Create tropopause data set from meteorological data.

Definition in file tropo.c.

Function Documentation

◆ main()

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

Definition at line 27 of file tropo.c.

29 {
30
31 ctl_t ctl;
32
33 clim_t *clim;
34
35 met_t *met;
36
37 static double ps[EX * EY], pt[EX * EY], qt[EX * EY], o3t[EX * EY],
38 zs[EX * EY], zt[EX * EY], tt[EX * EY], lon, lons[EX], lat, lats[EY];
39
40 static int init, nx, ny, nt, ncid, varid, dims[3];
41
42 static size_t count[10], start[10];
43
44 /* Allocate... */
45 ALLOC(clim, clim_t, 1);
46 ALLOC(met, met_t, 1);
47
48 /* Check arguments... */
49 if (argc < 4)
50 ERRMSG("Give parameters: <ctl> <tropo.nc> <met0> [ <met1> ... ]");
51
52 /* Read control parameters... */
53 read_ctl(argv[1], argc, argv, &ctl);
54 double lon0 = scan_ctl(argv[1], argc, argv, "TROPO_LON0", -1, "-180", NULL);
55 double lon1 = scan_ctl(argv[1], argc, argv, "TROPO_LON1", -1, "180", NULL);
56 double dlon = scan_ctl(argv[1], argc, argv, "TROPO_DLON", -1, "-999", NULL);
57 double lat0 = scan_ctl(argv[1], argc, argv, "TROPO_LAT0", -1, "-90", NULL);
58 double lat1 = scan_ctl(argv[1], argc, argv, "TROPO_LAT1", -1, "90", NULL);
59 double dlat = scan_ctl(argv[1], argc, argv, "TROPO_DLAT", -1, "-999", NULL);
60 int h2o = (int) scan_ctl(argv[1], argc, argv, "TROPO_H2O", -1, "1", NULL);
61 int o3 = (int) scan_ctl(argv[1], argc, argv, "TROPO_O3", -1, "1", NULL);
62
63 /* Read climatological data... */
64 read_clim(&ctl, clim);
65
66 /* Loop over files... */
67 for (int i = 3; i < argc; i++) {
68
69 /* Set control parameters... */
70 ctl.met_tropo = 0;
71
72 /* Read meteorological data... */
73 if (!read_met(argv[i], &ctl, clim, met))
74 continue;
75
76 /* Set horizontal grid... */
77 if (!init) {
78 init = 1;
79
80 /* Get grid... */
81 if (dlon <= 0)
82 dlon = fabs(met->lon[1] - met->lon[0]);
83 if (dlat <= 0)
84 dlat = fabs(met->lat[1] - met->lat[0]);
85 if (lon0 < -360 && lon1 > 360) {
86 lon0 = gsl_stats_min(met->lon, 1, (size_t) met->nx);
87 lon1 = gsl_stats_max(met->lon, 1, (size_t) met->nx);
88 }
89 nx = ny = 0;
90 for (lon = lon0; lon <= lon1; lon += dlon) {
91 lons[nx] = lon;
92 if ((++nx) >= EX)
93 ERRMSG("Too many longitudes!");
94 }
95 if (lat0 < -90 && lat1 > 90) {
96 lat0 = gsl_stats_min(met->lat, 1, (size_t) met->ny);
97 lat1 = gsl_stats_max(met->lat, 1, (size_t) met->ny);
98 }
99 for (lat = lat0; lat <= lat1; lat += dlat) {
100 lats[ny] = lat;
101 if ((++ny) >= EY)
102 ERRMSG("Too many latitudes!");
103 }
104
105 /* Create netCDF file... */
106 LOG(1, "Write tropopause data file: %s", argv[2]);
107 NC(nc_create(argv[2], NC_CLOBBER, &ncid));
108
109 /* Create dimensions... */
110 NC(nc_def_dim(ncid, "time", (size_t) NC_UNLIMITED, &dims[0]));
111 NC(nc_def_dim(ncid, "lat", (size_t) ny, &dims[1]));
112 NC(nc_def_dim(ncid, "lon", (size_t) nx, &dims[2]));
113
114 /* Create variables... */
115 NC_DEF_VAR("time", NC_DOUBLE, 1, &dims[0], "time",
116 "seconds since 2000-01-01 00:00:00 UTC");
117 NC_DEF_VAR("lat", NC_DOUBLE, 1, &dims[1], "latitude", "degrees_north");
118 NC_DEF_VAR("lon", NC_DOUBLE, 1, &dims[2], "longitude", "degrees_east");
119
120 NC_DEF_VAR("clp_z", NC_FLOAT, 3, &dims[0], "cold point height", "km");
121 NC_DEF_VAR("clp_p", NC_FLOAT, 3, &dims[0], "cold point pressure",
122 "hPa");
123 NC_DEF_VAR("clp_t", NC_FLOAT, 3, &dims[0], "cold point temperature",
124 "K");
125 if (h2o)
126 NC_DEF_VAR("clp_q", NC_FLOAT, 3, &dims[0], "cold point water vapor",
127 "ppv");
128 if (o3)
129 NC_DEF_VAR("clp_o3", NC_FLOAT, 3, &dims[0], "cold point ozone",
130 "ppv");
131
132 NC_DEF_VAR("dyn_z", NC_FLOAT, 3, &dims[0],
133 "dynamical tropopause height", "km");
134 NC_DEF_VAR("dyn_p", NC_FLOAT, 3, &dims[0],
135 "dynamical tropopause pressure", "hPa");
136 NC_DEF_VAR("dyn_t", NC_FLOAT, 3, &dims[0],
137 "dynamical tropopause temperature", "K");
138 if (h2o)
139 NC_DEF_VAR("dyn_q", NC_FLOAT, 3, &dims[0],
140 "dynamical tropopause water vapor", "ppv");
141 if (o3)
142 NC_DEF_VAR("dyn_o3", NC_FLOAT, 3, &dims[0],
143 "dynamical tropopause ozone", "ppv");
144
145 NC_DEF_VAR("wmo_1st_z", NC_FLOAT, 3, &dims[0],
146 "WMO 1st tropopause height", "km");
147 NC_DEF_VAR("wmo_1st_p", NC_FLOAT, 3, &dims[0],
148 "WMO 1st tropopause pressure", "hPa");
149 NC_DEF_VAR("wmo_1st_t", NC_FLOAT, 3, &dims[0],
150 "WMO 1st tropopause temperature", "K");
151 if (h2o)
152 NC_DEF_VAR("wmo_1st_q", NC_FLOAT, 3, &dims[0],
153 "WMO 1st tropopause water vapor", "ppv");
154 if (o3)
155 NC_DEF_VAR("wmo_1st_o3", NC_FLOAT, 3, &dims[0],
156 "WMO 1st tropopause ozone", "ppv");
157
158 NC_DEF_VAR("wmo_2nd_z", NC_FLOAT, 3, &dims[0],
159 "WMO 2nd tropopause height", "km");
160 NC_DEF_VAR("wmo_2nd_p", NC_FLOAT, 3, &dims[0],
161 "WMO 2nd tropopause pressure", "hPa");
162 NC_DEF_VAR("wmo_2nd_t", NC_FLOAT, 3, &dims[0],
163 "WMO 2nd tropopause temperature", "K");
164 if (h2o)
165 NC_DEF_VAR("wmo_2nd_q", NC_FLOAT, 3, &dims[0],
166 "WMO 2nd tropopause water vapor", "ppv");
167 if (o3)
168 NC_DEF_VAR("wmo_2nd_o3", NC_FLOAT, 3, &dims[0],
169 "WMO 2nd tropopause ozone", "ppv");
170
171 NC_DEF_VAR("ps", NC_FLOAT, 3, &dims[0], "surface pressure", "hPa");
172 NC_DEF_VAR("zs", NC_FLOAT, 3, &dims[0], "surface height", "km");
173
174 /* End definition... */
175 NC(nc_enddef(ncid));
176
177 /* Write longitude and latitude... */
178 NC_PUT_DOUBLE("lat", lats, 0);
179 NC_PUT_DOUBLE("lon", lons, 0);
180 }
181
182 /* Write time... */
183 start[0] = (size_t) nt;
184 count[0] = 1;
185 start[1] = 0;
186 count[1] = (size_t) ny;
187 start[2] = 0;
188 count[2] = (size_t) nx;
189 NC_PUT_DOUBLE("time", &met->time, 1);
190
191 /* Get cold point... */
192 get_tropo(2, &ctl, clim, met, lons, nx, lats, ny, pt, zt, tt, qt, o3t, ps,
193 zs);
194 NC_PUT_DOUBLE("clp_z", zt, 1);
195 NC_PUT_DOUBLE("clp_p", pt, 1);
196 NC_PUT_DOUBLE("clp_t", tt, 1);
197 if (h2o)
198 NC_PUT_DOUBLE("clp_q", qt, 1);
199 if (o3)
200 NC_PUT_DOUBLE("clp_o3", o3t, 1);
201
202 /* Get dynamical tropopause... */
203 get_tropo(5, &ctl, clim, met, lons, nx, lats, ny, pt, zt, tt, qt, o3t, ps,
204 zs);
205 NC_PUT_DOUBLE("dyn_z", zt, 1);
206 NC_PUT_DOUBLE("dyn_p", pt, 1);
207 NC_PUT_DOUBLE("dyn_t", tt, 1);
208 if (h2o)
209 NC_PUT_DOUBLE("dyn_q", qt, 1);
210 if (o3)
211 NC_PUT_DOUBLE("dyn_o3", o3t, 1);
212
213 /* Get WMO 1st tropopause... */
214 get_tropo(3, &ctl, clim, met, lons, nx, lats, ny, pt, zt, tt, qt, o3t, ps,
215 zs);
216 NC_PUT_DOUBLE("wmo_1st_z", zt, 1);
217 NC_PUT_DOUBLE("wmo_1st_p", pt, 1);
218 NC_PUT_DOUBLE("wmo_1st_t", tt, 1);
219 if (h2o)
220 NC_PUT_DOUBLE("wmo_1st_q", qt, 1);
221 if (o3)
222 NC_PUT_DOUBLE("wmo_1st_o3", o3t, 1);
223
224 /* Get WMO 2nd tropopause... */
225 get_tropo(4, &ctl, clim, met, lons, nx, lats, ny, pt, zt, tt, qt, o3t, ps,
226 zs);
227 NC_PUT_DOUBLE("wmo_2nd_z", zt, 1);
228 NC_PUT_DOUBLE("wmo_2nd_p", pt, 1);
229 NC_PUT_DOUBLE("wmo_2nd_t", tt, 1);
230 if (h2o)
231 NC_PUT_DOUBLE("wmo_2nd_q", qt, 1);
232 if (o3)
233 NC_PUT_DOUBLE("wmo_2nd_o3", o3t, 1);
234
235 /* Write surface data... */
236 NC_PUT_DOUBLE("ps", ps, 1);
237 NC_PUT_DOUBLE("zs", zs, 1);
238
239 /* Increment time step counter... */
240 nt++;
241 }
242
243 /* Close file... */
244 NC(nc_close(ncid));
245
246 /* Free... */
247 free(clim);
248 free(met);
249
250 return EXIT_SUCCESS;
251}
read_clim
void read_clim(ctl_t *ctl, clim_t *clim)
Reads various climatological data and populates the given climatology structure.
Definition: mptrac.c:4441
get_tropo
void get_tropo(int met_tropo, ctl_t *ctl, clim_t *clim, met_t *met, double *lons, int nx, double *lats, int ny, double *pt, double *zt, double *tt, double *qt, double *o3t, double *ps, double *zs)
Calculate tropopause data.
Definition: mptrac.c:1190
read_ctl
void read_ctl(const char *filename, int argc, char *argv[], ctl_t *ctl)
Reads control parameters from a configuration file and populates the given structure.
Definition: mptrac.c:4789
read_met
int read_met(const char *filename, ctl_t *ctl, clim_t *clim, met_t *met)
Reads meteorological data from a file and populates the provided structures.
Definition: mptrac.c:5611
scan_ctl
double scan_ctl(const char *filename, int argc, char *argv[], const char *varname, int arridx, const char *defvalue, char *value)
Scans a control file or command-line arguments for a specified variable.
Definition: mptrac.c:7943
NC
#define NC(cmd)
Execute a NetCDF command and check for errors.
Definition: mptrac.h:981
ERRMSG
#define ERRMSG(...)
Print an error message with contextual information and terminate the program.
Definition: mptrac.h:1881
EY
#define EY
Maximum number of latitudes for meteo data.
Definition: mptrac.h:266
EX
#define EX
Maximum number of longitudes for meteo data.
Definition: mptrac.h:261
ALLOC
#define ALLOC(ptr, type, n)
Allocate memory for a pointer with error handling.
Definition: mptrac.h:344
LOG
#define LOG(level,...)
Print a log message with a specified logging level.
Definition: mptrac.h:1811
NC_PUT_DOUBLE
#define NC_PUT_DOUBLE(varname, ptr, hyperslab)
Write double precision data to a NetCDF variable.
Definition: mptrac.h:1077
NC_DEF_VAR
#define NC_DEF_VAR(varname, type, ndims, dims, long_name, units)
Define a NetCDF variable with attributes.
Definition: mptrac.h:1003
clim_t
Climatological data.
Definition: mptrac.h:3230
ctl_t
Control parameters.
Definition: mptrac.h:2135
ctl_t::met_tropo
int met_tropo
Tropopause definition (0=none, 1=clim, 2=cold point, 3=WMO_1st, 4=WMO_2nd, 5=dynamical).
Definition: mptrac.h:2560
met_t
Meteo data structure.
Definition: mptrac.h:3289
met_t::nx
int nx
Number of longitudes.
Definition: mptrac.h:3295
met_t::ny
int ny
Number of latitudes.
Definition: mptrac.h:3298
met_t::lon
double lon[EX]
Longitude [deg].
Definition: mptrac.h:3307
met_t::time
double time
Time [s].
Definition: mptrac.h:3292
met_t::lat
double lat[EY]
Latitude [deg].
Definition: mptrac.h:3310
Here is the call graph for this function:
Generated on Mon Sep 16 2024 17:36:26 for MPTRAC by doxygen 1.9.4