CrIS Code Collection
Macros | Functions
perturbation.c File Reference
#include "libcris.h"

Go to the source code of this file.

Macros

#define N4   54
 Number of 4.3 micron channels: More...
 
#define N15_LOW   15
 Number of 15 micron low channels: More...
 
#define N15_HIGH   1
 Number of 15 micron high channels: More...
 
#define N8   1
 Number of 8.1 micron channels: More...
 
#define N10   1
 Number of 10.4 micron channels: More...
 
#define LOOP_ALL(ntrack)
 Loop over all footprints. More...
 

Functions

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

Macro Definition Documentation

◆ N4

#define N4   54

Number of 4.3 micron channels:

Definition at line 8 of file perturbation.c.

◆ N15_LOW

#define N15_LOW   15

Number of 15 micron low channels:

Definition at line 11 of file perturbation.c.

◆ N15_HIGH

#define N15_HIGH   1

Number of 15 micron high channels:

Definition at line 14 of file perturbation.c.

◆ N8

#define N8   1

Number of 8.1 micron channels:

Definition at line 17 of file perturbation.c.

◆ N10

#define N10   1

Number of 10.4 micron channels:

Definition at line 20 of file perturbation.c.

◆ LOOP_ALL

#define LOOP_ALL (   ntrack)
Value:
for (int track = 0; track < ntrack; track++) \
for (int xtrack = 0; xtrack < L1_NXTRACK; xtrack++) \
for (int ifov = 0; ifov < L1_NFOV; ifov++)
L1_NXTRACK
#define L1_NXTRACK
Across-track size of CrIS radiance granule.
Definition: libcris.h:24
L1_NFOV
#define L1_NFOV
Number of field of views of CrIS radiance granule.
Definition: libcris.h:27

Loop over all footprints.

Definition at line 27 of file perturbation.c.

Function Documentation

◆ main()

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

Definition at line 36 of file perturbation.c.

38 {
39
40 static cris_l1_t l1;
41
42 static pert_t *pert_4mu, *pert_15mu_low, *pert_15mu_high;
43
44 const double var_dh2 = 100. * 100.;
45
46 static double help[PERT_NTRACK * PERT_NXTRACK * PERT_NFOV], rad, nu,
47 x[L1_NXTRACK], y[L1_NXTRACK], x0[3], x1[3];
48
49 const int list_4mu[N4]
50 = { 270, 271, 272, 273, 274, 275, 276, 277, 278, 279, 280, 281, 282, 283,
51 284, 285, 286, 288, 289, 291, 292, 294, 295, 296, 297, 298, 299, 300,
52 302, 303, 304, 305, 306, 307, 318, 319, 321, 322, 323, 324, 325, 326,
53 328, 330, 332, 333, 334, 335, 336, 337, 338, 339, 340, 341
54 };
55
56 const int list_15mu_low[N15_LOW]
57 = { 5, 8, 10, 13, 15, 17, 23, 34, 36, 39, 41, 44, 46, 49, 51 };
58
59 const int list_15mu_high[N15_HIGH] = { 31 };
60
61 const int list_8mu[N8] = { 36 };
62
63 const int list_10mu[N10] = { 502 };
64
65 const int dtrack = 3, dxtrack = 3;
66
67 static int dimid[3], n, ncid, track0,
68 time_varid, lon_varid, lat_varid, bt_4mu_varid, init,
69 bt_4mu_pt_varid, bt_4mu_var_varid, bt_8mu_varid, bt_10mu_varid,
70 bt_15mu_low_varid, bt_15mu_low_pt_varid, bt_15mu_low_var_varid,
71 bt_15mu_high_varid, bt_15mu_high_pt_varid, bt_15mu_high_var_varid;
72
73 /* Check arguments... */
74 if (argc < 4)
75 ERRMSG("Give parameters: <ctl> <pert.nc> <l1b_file1> [<l1b_file2> ...]");
76
77 /* Get control parameters... */
78 const int apo = (int) scan_ctl(argc, argv, "APO", -1, "0", NULL);
79
80 /* Allocate... */
81 ALLOC(pert_4mu, pert_t, 1);
82 ALLOC(pert_15mu_low, pert_t, 1);
83 ALLOC(pert_15mu_high, pert_t, 1);
84
85 /* ------------------------------------------------------------
86 Read CrIS Level-1B files...
87 ------------------------------------------------------------ */
88
89 /* Loop over files... */
90 for (int iarg = 3; iarg < argc; iarg++) {
91
92 /* Read CrIS data... */
93 if (!read_cris_l1(argv[iarg], &l1, apo))
94 continue;
95
96 /* Write channel information... */
97 if (!init) {
98 init = 1;
99
100 /* 4.3 micron channels... */
101 LOG(2, "4.3 micron channels:");
102 n = 0;
103 nu = 0;
104 for (int i = 0; i < N4; i++) {
105 nu += l1.nu_sw[list_4mu[i]];
106 n++;
107 LOG(2, " count= %2d | channel index= SW_%03d | nu= %.4f cm^-1", n,
108 list_4mu[i], l1.nu_sw[list_4mu[i]]);
109 }
110 LOG(2, " nu_mean= %.4f cm^-1", nu / n);
111
112 /* 15 micron low channels... */
113 LOG(2, "15 micron low channels:");
114 n = 0;
115 nu = 0;
116 for (int i = 0; i < N15_LOW; i++) {
117 nu += l1.nu_lw[list_15mu_low[i]];
118 n++;
119 LOG(2, " count= %2d | channel index= LW_%03d | nu= %.4f cm^-1", n,
120 list_15mu_low[i], l1.nu_lw[list_15mu_low[i]]);
121 }
122 LOG(2, " nu_mean= %.4f cm^-1", nu / n);
123
124 /* 15 micron high channels... */
125 LOG(2, "15 micron high channels:");
126 n = 0;
127 nu = 0;
128 for (int i = 0; i < N15_HIGH; i++) {
129 nu += l1.nu_lw[list_15mu_high[i]];
130 n++;
131 LOG(2, " count= %2d | channel index= LW_%03d | nu= %.4f cm^-1", n,
132 list_15mu_high[i], l1.nu_lw[list_15mu_high[i]]);
133 }
134 LOG(2, " nu_mean= %.4f cm^-1", nu / n);
135
136 /* 8.1 micron channels... */
137 LOG(2, "8.1 micron channels:");
138 n = 0;
139 nu = 0;
140 for (int i = 0; i < N8; i++) {
141 nu += l1.nu_mw[list_8mu[i]];
142 n++;
143 LOG(2, " count= %2d | channel index= MW_%03d | nu= %.4f cm^-1", n,
144 list_8mu[i], l1.nu_mw[list_8mu[i]]);
145 }
146 LOG(2, " nu_mean= %.4f cm^-1", nu / n);
147
148 /* 10.4 micron channels... */
149 LOG(2, "10.4 micron channels:");
150 n = 0;
151 nu = 0;
152 for (int i = 0; i < N10; i++) {
153 nu += l1.nu_lw[list_10mu[i]];
154 n++;
155 LOG(2, " count= %2d | channel index= LW_%03d | nu= %.4f cm^-1", n,
156 list_10mu[i], l1.nu_lw[list_10mu[i]]);
157 }
158 LOG(2, " nu_mean= %.4f cm^-1", nu / n);
159 }
160
161 /* Save geolocation... */
162 pert_4mu->ntrack += L1_NTRACK;
163 if (pert_4mu->ntrack > PERT_NTRACK)
164 ERRMSG("Too many granules!");
165 pert_4mu->nxtrack = L1_NXTRACK;
166 if (pert_4mu->nxtrack > PERT_NXTRACK)
167 ERRMSG("Too many tracks!");
168 pert_4mu->nfov = L1_NFOV;
169 if (pert_4mu->nfov > PERT_NFOV)
170 ERRMSG("Too many field of views!");
171 LOOP_ALL(L1_NTRACK) {
172 pert_4mu->time[track0 + track][xtrack][ifov]
173 = l1.time[track][xtrack] - 220838400.;
174 pert_4mu->lon[track0 + track][xtrack][ifov]
175 = l1.lon[track][xtrack][ifov];
176 pert_4mu->lat[track0 + track][xtrack][ifov]
177 = l1.lat[track][xtrack][ifov];
178 }
179
180 /* Get 8.1 micron brightness temperature... */
181 LOOP_ALL(L1_NTRACK)
182 pert_4mu->dc[track0 + track][xtrack][ifov]
183 = BRIGHT(l1.rad_mw[track][xtrack][ifov][list_8mu[0]] * 1e-3,
184 l1.nu_mw[list_8mu[0]]);
185
186 /* Get 10.4 micron brightness temperature... */
187 LOOP_ALL(L1_NTRACK)
188 pert_15mu_high->dc[track0 + track][xtrack][ifov]
189 = BRIGHT(l1.rad_lw[track][xtrack][ifov][list_10mu[0]] * 1e-3,
190 l1.nu_lw[list_10mu[0]]);
191
192 /* Get 4.3 micron brightness temperature... */
193 LOOP_ALL(L1_NTRACK) {
194 n = 0;
195 nu = rad = 0;
196 for (int i = 0; i < N4; i++)
197 if (gsl_finite(l1.rad_sw[track][xtrack][ifov][list_4mu[i]])) {
198 rad += l1.rad_sw[track][xtrack][ifov][list_4mu[i]] * 1e-3;
199 nu += l1.nu_sw[list_4mu[i]];
200 n++;
201 }
202 if (n > 0.9 * N4)
203 pert_4mu->bt[track0 + track][xtrack][ifov] = BRIGHT(rad / n, nu / n);
204 else
205 pert_4mu->bt[track0 + track][xtrack][ifov] = GSL_NAN;
206 }
207
208 /* Get 15 micron low brightness temperature... */
209 LOOP_ALL(L1_NTRACK) {
210 n = 0;
211 nu = rad = 0;
212 for (int i = 0; i < N15_LOW; i++)
213 if (gsl_finite(l1.rad_lw[track][xtrack][ifov][list_15mu_low[i]])) {
214 rad += l1.rad_lw[track][xtrack][ifov][list_15mu_low[i]] * 1e-3;
215 nu += l1.nu_lw[list_15mu_low[i]];
216 n++;
217 }
218 if (n > 0.9 * N15_LOW)
219 pert_15mu_low->bt[track0 + track][xtrack][ifov] =
220 BRIGHT(rad / n, nu / n);
221 else
222 pert_15mu_low->bt[track0 + track][xtrack][ifov] = GSL_NAN;
223 }
224
225 /* Get 15 micron high brightness temperature... */
226 LOOP_ALL(L1_NTRACK) {
227 n = 0;
228 nu = rad = 0;
229 for (int i = 0; i < N15_HIGH; i++)
230 if (gsl_finite(l1.rad_lw[track][xtrack][ifov][list_15mu_high[i]])) {
231 rad += l1.rad_lw[track][xtrack][ifov][list_15mu_high[i]] * 1e-3;
232 nu += l1.nu_lw[list_15mu_high[i]];
233 n++;
234 }
235 if (n > 0.9 * N15_HIGH)
236 pert_15mu_high->bt[track0 + track][xtrack][ifov] =
237 BRIGHT(rad / n, nu / n);
238 else
239 pert_15mu_high->bt[track0 + track][xtrack][ifov] = GSL_NAN;
240 }
241
242 /* Increment track counter... */
243 track0 += L1_NTRACK;
244 }
245
246 /* ------------------------------------------------------------
247 Calculate perturbations...
248 ------------------------------------------------------------ */
249
250 /* Write info... */
251 LOG(1, "Calculate perturbations...");
252
253 /* Loop over scans and field of views... */
254 for (int track = 0; track < pert_4mu->ntrack; track++)
255 for (int ifov = 0; ifov < L1_NFOV; ifov++) {
256
257 /* Skip scan edges... */
258 pert_4mu->dc[track][0][ifov] = GSL_NAN;
259 pert_4mu->dc[track][L1_NXTRACK - 1][ifov] = GSL_NAN;
260 pert_15mu_high->dc[track][0][ifov] = GSL_NAN;
261 pert_15mu_high->dc[track][L1_NXTRACK - 1][ifov] = GSL_NAN;
262 pert_4mu->bt[track][0][ifov] = GSL_NAN;
263 pert_4mu->bt[track][L1_NXTRACK - 1][ifov] = GSL_NAN;
264 pert_15mu_low->bt[track][0][ifov] = GSL_NAN;
265 pert_15mu_low->bt[track][L1_NXTRACK - 1][ifov] = GSL_NAN;
266 pert_15mu_high->bt[track][0][ifov] = GSL_NAN;
267 pert_15mu_high->bt[track][L1_NXTRACK - 1][ifov] = GSL_NAN;
268
269 /* Get 4 micron perturbations... */
270 for (int xtrack = 0; xtrack < L1_NXTRACK; xtrack++) {
271 x[xtrack] = (double) xtrack;
272 y[xtrack] = pert_4mu->bt[track][xtrack][ifov];
273 }
274 background_poly_help(x, y, L1_NXTRACK, 5);
275 for (int xtrack = 0; xtrack < L1_NXTRACK; xtrack++)
276 pert_4mu->pt[track][xtrack][ifov] =
277 pert_4mu->bt[track][xtrack][ifov] - y[xtrack];
278
279 /* Get 15 micron low perturbations... */
280 for (int xtrack = 0; xtrack < L1_NXTRACK; xtrack++) {
281 x[xtrack] = (double) xtrack;
282 y[xtrack] = pert_15mu_low->bt[track][xtrack][ifov];
283 }
284 background_poly_help(x, y, L1_NXTRACK, 7);
285 for (int xtrack = 0; xtrack < L1_NXTRACK; xtrack++)
286 pert_15mu_low->pt[track][xtrack][ifov] =
287 pert_15mu_low->bt[track][xtrack][ifov] - y[xtrack];
288
289 /* Get 15 micron high perturbations... */
290 for (int xtrack = 0; xtrack < L1_NXTRACK; xtrack++) {
291 x[xtrack] = (double) xtrack;
292 y[xtrack] = pert_15mu_high->bt[track][xtrack][ifov];
293 }
294 background_poly_help(x, y, L1_NXTRACK, 5);
295 for (int xtrack = 0; xtrack < L1_NXTRACK; xtrack++)
296 pert_15mu_high->pt[track][xtrack][ifov] =
297 pert_15mu_high->bt[track][xtrack][ifov] - y[xtrack];
298 }
299
300 /* ------------------------------------------------------------
301 Calculate variances...
302 ------------------------------------------------------------ */
303
304 /* Write info... */
305 LOG(1, "Calculate variances...");
306
307 /* Loop over footprints... */
308 LOOP_ALL(pert_4mu->ntrack) {
309
310 /* Get geolocation... */
311 geo2cart(0.0, pert_4mu->lon[track][xtrack][ifov],
312 pert_4mu->lat[track][xtrack][ifov], x0);
313
314 /* Init... */
315 int n_4mu = 0, n_15mu_low = 0, n_15mu_high = 0;
316 double mean_4mu = 0, mean_15mu_low = 0, mean_15mu_high = 0;
317 double var_4mu = 0, var_15mu_low = 0, var_15mu_high = 0;
318
319 /* Loop over neighbouring footprints... */
320 for (int track2 = track - dtrack; track2 <= track + dtrack; track2++)
321 for (int xtrack2 = xtrack - dxtrack; xtrack2 <= xtrack + dxtrack;
322 xtrack2++)
323 for (int ifov2 = 0; ifov2 < L1_NFOV; ifov2++)
324 if (track2 >= 0 && track2 < pert_4mu->ntrack
325 && xtrack2 >= 0 && xtrack2 < L1_NXTRACK) {
326
327 /* Check horizontal distance... */
328 geo2cart(0.0, pert_4mu->lon[track2][xtrack2][ifov2],
329 pert_4mu->lat[track2][xtrack2][ifov2], x1);
330 if (DIST2(x0, x1) > var_dh2)
331 continue;
332
333 /* Calculate variance... */
334 if (gsl_finite(pert_4mu->pt[track2][xtrack2][ifov2])) {
335 mean_4mu += pert_4mu->pt[track2][xtrack2][ifov2];
336 var_4mu += gsl_pow_2(pert_4mu->pt[track2][xtrack2][ifov2]);
337 n_4mu++;
338 }
339
340 if (gsl_finite(pert_15mu_low->pt[track2][xtrack2][ifov2])) {
341 mean_15mu_low += pert_15mu_low->pt[track2][xtrack2][ifov2];
342 var_15mu_low +=
343 gsl_pow_2(pert_15mu_low->pt[track2][xtrack2][ifov2]);
344 n_15mu_low++;
345 }
346
347 if (gsl_finite(pert_15mu_high->pt[track2][xtrack2][ifov2])) {
348 mean_15mu_high += pert_15mu_high->pt[track2][xtrack2][ifov2];
349 var_15mu_high +=
350 gsl_pow_2(pert_15mu_high->pt[track2][xtrack2][ifov2]);
351 n_15mu_high++;
352 }
353 }
354
355 /* Save variance data... */
356 if (n_4mu > 0 && xtrack > 0 && xtrack < L1_NXTRACK - 1)
357 pert_4mu->var[track][xtrack][ifov] =
358 var_4mu / n_4mu - gsl_pow_2(mean_4mu / n_4mu);
359 else
360 pert_4mu->var[track][xtrack][ifov] = GSL_NAN;
361
362 if (n_15mu_low > 0 && xtrack > 0 && xtrack < L1_NXTRACK - 1)
363 pert_15mu_low->var[track][xtrack][ifov] =
364 var_15mu_low / n_15mu_low - gsl_pow_2(mean_15mu_low / n_15mu_low);
365 else
366 pert_15mu_low->var[track][xtrack][ifov] = GSL_NAN;
367
368 if (n_15mu_high > 0 && xtrack > 0 && xtrack < L1_NXTRACK - 1)
369 pert_15mu_high->var[track][xtrack][ifov] =
370 var_15mu_high / n_15mu_high - gsl_pow_2(mean_15mu_high / n_15mu_high);
371 else
372 pert_15mu_high->var[track][xtrack][ifov] = GSL_NAN;
373 }
374
375 /* ------------------------------------------------------------
376 Write netCDF file...
377 ------------------------------------------------------------ */
378
379 /* Write info... */
380 LOG(1, "Write perturbation data file: %s", argv[2]);
381
382 /* Create netCDF file... */
383 NC(nc_create(argv[2], NC_CLOBBER, &ncid));
384
385 /* Set dimensions... */
386 NC(nc_def_dim(ncid, "NTRACK", (size_t) pert_4mu->ntrack, &dimid[0]));
387 NC(nc_def_dim(ncid, "NXTRACK", L1_NXTRACK, &dimid[1]));
388 NC(nc_def_dim(ncid, "NFOV", L1_NFOV, &dimid[2]));
389
390 /* Add variables... */
391 NC(nc_def_var(ncid, "time", NC_DOUBLE, 3, dimid, &time_varid));
392 add_att(ncid, time_varid, "s", "time (seconds since 2000-01-01T00:00Z)");
393 NC(nc_def_var(ncid, "lon", NC_DOUBLE, 3, dimid, &lon_varid));
394 add_att(ncid, lon_varid, "deg", "footprint longitude");
395 NC(nc_def_var(ncid, "lat", NC_DOUBLE, 3, dimid, &lat_varid));
396 add_att(ncid, lat_varid, "deg", "footprint latitude");
397
398 NC(nc_def_var(ncid, "bt_8mu", NC_FLOAT, 3, dimid, &bt_8mu_varid));
399 add_att(ncid, bt_8mu_varid, "K", "brightness temperature at 8.1 micron");
400
401 NC(nc_def_var(ncid, "bt_10mu", NC_FLOAT, 3, dimid, &bt_10mu_varid));
402 add_att(ncid, bt_10mu_varid, "K", "brightness temperature at 10.4 micron");
403
404 NC(nc_def_var(ncid, "bt_4mu", NC_FLOAT, 3, dimid, &bt_4mu_varid));
405 add_att(ncid, bt_4mu_varid, "K", "brightness temperature" " at 4.3 micron");
406 NC(nc_def_var(ncid, "bt_4mu_pt", NC_FLOAT, 3, dimid, &bt_4mu_pt_varid));
407 add_att(ncid, bt_4mu_pt_varid, "K", "brightness temperature perturbation"
408 " at 4.3 micron");
409 NC(nc_def_var(ncid, "bt_4mu_var", NC_FLOAT, 3, dimid, &bt_4mu_var_varid));
410 add_att(ncid, bt_4mu_var_varid, "K^2", "brightness temperature variance"
411 " at 4.3 micron");
412
413 NC(nc_def_var(ncid, "bt_15mu_low", NC_FLOAT, 3, dimid, &bt_15mu_low_varid));
414 add_att(ncid, bt_15mu_low_varid, "K", "brightness temperature"
415 " at 15 micron (low altitudes)");
416 NC(nc_def_var(ncid, "bt_15mu_low_pt", NC_FLOAT, 3, dimid,
417 &bt_15mu_low_pt_varid));
418 add_att(ncid, bt_15mu_low_pt_varid, "K",
419 "brightness temperature perturbation"
420 " at 15 micron (low altitudes)");
421 NC(nc_def_var
422 (ncid, "bt_15mu_low_var", NC_FLOAT, 3, dimid, &bt_15mu_low_var_varid));
423 add_att(ncid, bt_15mu_low_var_varid, "K^2",
424 "brightness temperature variance" " at 15 micron (low altitudes)");
425
426 NC(nc_def_var(ncid, "bt_15mu_high", NC_FLOAT, 3, dimid,
427 &bt_15mu_high_varid));
428 add_att(ncid, bt_15mu_high_varid, "K", "brightness temperature"
429 " at 15 micron (high altitudes)");
430 NC(nc_def_var(ncid, "bt_15mu_high_pt", NC_FLOAT, 3, dimid,
431 &bt_15mu_high_pt_varid));
432 add_att(ncid, bt_15mu_high_pt_varid, "K",
433 "brightness temperature perturbation"
434 " at 15 micron (high altitudes)");
435 NC(nc_def_var
436 (ncid, "bt_15mu_high_var", NC_FLOAT, 3, dimid, &bt_15mu_high_var_varid));
437 add_att(ncid, bt_15mu_high_var_varid, "K^2",
438 "brightness temperature variance" " at 15 micron (high altitudes)");
439
440 /* Leave define mode... */
441 NC(nc_enddef(ncid));
442
443 /* Write data... */
444 n = 0;
445 LOOP_ALL(pert_4mu->ntrack)
446 help[n++] = pert_4mu->time[track][xtrack][ifov];
447 NC(nc_put_var_double(ncid, time_varid, help));
448
449 n = 0;
450 LOOP_ALL(pert_4mu->ntrack)
451 help[n++] = pert_4mu->lon[track][xtrack][ifov];
452 NC(nc_put_var_double(ncid, lon_varid, help));
453
454 n = 0;
455 LOOP_ALL(pert_4mu->ntrack)
456 help[n++] = pert_4mu->lat[track][xtrack][ifov];
457 NC(nc_put_var_double(ncid, lat_varid, help));
458
459 n = 0;
460 LOOP_ALL(pert_4mu->ntrack)
461 help[n++] = pert_4mu->dc[track][xtrack][ifov];
462 NC(nc_put_var_double(ncid, bt_8mu_varid, help));
463
464 n = 0;
465 LOOP_ALL(pert_4mu->ntrack)
466 help[n++] = pert_15mu_high->dc[track][xtrack][ifov];
467 NC(nc_put_var_double(ncid, bt_10mu_varid, help));
468
469 n = 0;
470 LOOP_ALL(pert_4mu->ntrack)
471 help[n++] = pert_4mu->bt[track][xtrack][ifov];
472 NC(nc_put_var_double(ncid, bt_4mu_varid, help));
473
474 n = 0;
475 LOOP_ALL(pert_4mu->ntrack)
476 help[n++] = pert_4mu->pt[track][xtrack][ifov];
477 NC(nc_put_var_double(ncid, bt_4mu_pt_varid, help));
478
479 n = 0;
480 LOOP_ALL(pert_4mu->ntrack)
481 help[n++] = pert_4mu->var[track][xtrack][ifov];
482 NC(nc_put_var_double(ncid, bt_4mu_var_varid, help));
483
484 n = 0;
485 LOOP_ALL(pert_4mu->ntrack)
486 help[n++] = pert_15mu_low->bt[track][xtrack][ifov];
487 NC(nc_put_var_double(ncid, bt_15mu_low_varid, help));
488
489 n = 0;
490 LOOP_ALL(pert_4mu->ntrack)
491 help[n++] = pert_15mu_low->pt[track][xtrack][ifov];
492 NC(nc_put_var_double(ncid, bt_15mu_low_pt_varid, help));
493
494 n = 0;
495 LOOP_ALL(pert_4mu->ntrack)
496 help[n++] = pert_15mu_low->var[track][xtrack][ifov];
497 NC(nc_put_var_double(ncid, bt_15mu_low_var_varid, help));
498
499 n = 0;
500 LOOP_ALL(pert_4mu->ntrack)
501 help[n++] = pert_15mu_high->bt[track][xtrack][ifov];
502 NC(nc_put_var_double(ncid, bt_15mu_high_varid, help));
503
504 n = 0;
505 LOOP_ALL(pert_4mu->ntrack)
506 help[n++] = pert_15mu_high->pt[track][xtrack][ifov];
507 NC(nc_put_var_double(ncid, bt_15mu_high_pt_varid, help));
508
509 n = 0;
510 LOOP_ALL(pert_4mu->ntrack)
511 help[n++] = pert_15mu_high->var[track][xtrack][ifov];
512 NC(nc_put_var_double(ncid, bt_15mu_high_var_varid, help));
513
514 /* Close file... */
515 NC(nc_close(ncid));
516
517 /* Free... */
518 free(pert_4mu);
519 free(pert_15mu_low);
520 free(pert_15mu_high);
521
522 return EXIT_SUCCESS;
523}
scan_ctl
double scan_ctl(int argc, char *argv[], const char *varname, int arridx, const char *defvalue, char *value)
Search control parameter file for variable entry.
Definition: jurassic.c:5114
geo2cart
void geo2cart(const double z, const double lon, const double lat, double *x)
Convert geolocation to Cartesian coordinates.
Definition: jurassic.c:3500
BRIGHT
#define BRIGHT(rad, nu)
Compute brightness temperature.
Definition: jurassic.h:126
ERRMSG
#define ERRMSG(...)
Print error message and quit program.
Definition: jurassic.h:237
ALLOC
#define ALLOC(ptr, type, n)
Allocate memory.
Definition: jurassic.h:121
LOG
#define LOG(level,...)
Print log message.
Definition: jurassic.h:221
DIST2
#define DIST2(a, b)
Compute squared distance between two vectors.
Definition: jurassic.h:137
background_poly_help
void background_poly_help(const double *xx, double *yy, const int n, const int dim)
Get background based on polynomial fits.
Definition: libcris.c:47
read_cris_l1
int read_cris_l1(char *filename, cris_l1_t *l1, int apo)
Read CrIS Level-1 data.
Definition: libcris.c:1100
add_att
void add_att(const int ncid, const int varid, const char *unit, const char *long_name)
Add variable attributes to netCDF file.
Definition: libcris.c:5
NC
#define NC(cmd)
Execute netCDF library command and check result.
Definition: libcris.h:61
PERT_NXTRACK
#define PERT_NXTRACK
Across-track size of perturbation data.
Definition: libcris.h:42
L1_NTRACK
#define L1_NTRACK
Along-track size of CrIS radiance granule.
Definition: libcris.h:21
PERT_NTRACK
#define PERT_NTRACK
Along-track size of perturbation data.
Definition: libcris.h:39
PERT_NFOV
#define PERT_NFOV
Number of field of views of perturbation data.
Definition: libcris.h:45
N15_LOW
#define N15_LOW
Number of 15 micron low channels:
Definition: perturbation.c:11
N10
#define N10
Number of 10.4 micron channels:
Definition: perturbation.c:20
N4
#define N4
Number of 4.3 micron channels:
Definition: perturbation.c:8
N8
#define N8
Number of 8.1 micron channels:
Definition: perturbation.c:17
N15_HIGH
#define N15_HIGH
Number of 15 micron high channels:
Definition: perturbation.c:14
LOOP_ALL
#define LOOP_ALL(ntrack)
Loop over all footprints.
Definition: perturbation.c:27
cris_l1_t
CrIS Level-1 data.
Definition: libcris.h:72
cris_l1_t::nu_sw
double nu_sw[L1_NCHAN_SW]
Shortwave channel frequencies [cm^-1].
Definition: libcris.h:117
cris_l1_t::rad_mw
float rad_mw[L1_NTRACK][L1_NXTRACK][L1_NFOV][L1_NCHAN_MW]
Midwave radiance [W/(m^2 sr cm^-1)].
Definition: libcris.h:108
cris_l1_t::nu_mw
double nu_mw[L1_NCHAN_MW]
Midwave channel frequencies [cm^-1].
Definition: libcris.h:105
cris_l1_t::rad_sw
float rad_sw[L1_NTRACK][L1_NXTRACK][L1_NFOV][L1_NCHAN_SW]
Shortwave radiance [W/(m^2 sr cm^-1)].
Definition: libcris.h:120
cris_l1_t::lon
double lon[L1_NTRACK][L1_NXTRACK][L1_NFOV]
Footprint longitude [deg].
Definition: libcris.h:78
cris_l1_t::rad_lw
float rad_lw[L1_NTRACK][L1_NXTRACK][L1_NFOV][L1_NCHAN_LW]
Longwave radiance [W/(m^2 sr cm^-1)].
Definition: libcris.h:96
cris_l1_t::time
double time[L1_NTRACK][L1_NXTRACK]
Time (seconds since 2000-01-01T00:00Z).
Definition: libcris.h:75
cris_l1_t::nu_lw
double nu_lw[L1_NCHAN_LW]
Longwave channel frequencies [cm^-1].
Definition: libcris.h:93
cris_l1_t::lat
double lat[L1_NTRACK][L1_NXTRACK][L1_NFOV]
Footprint latitude [deg].
Definition: libcris.h:81
pert_t
Perturbation data.
Definition: libcris.h:131
pert_t::dc
double dc[PERT_NTRACK][PERT_NXTRACK][PERT_NFOV]
Brightness temperature (cloud channel) [K].
Definition: libcris.h:152
pert_t::lat
double lat[PERT_NTRACK][PERT_NXTRACK][PERT_NFOV]
Latitude [deg].
Definition: libcris.h:149
pert_t::nfov
int nfov
Number of field of views.
Definition: libcris.h:140
pert_t::pt
double pt[PERT_NTRACK][PERT_NXTRACK][PERT_NFOV]
Brightness temperature perturbation (4 or 15 micron) [K].
Definition: libcris.h:158
pert_t::ntrack
int ntrack
Number of along-track values.
Definition: libcris.h:134
pert_t::var
double var[PERT_NTRACK][PERT_NXTRACK][PERT_NFOV]
Brightness temperature variance (4 or 15 micron) [K].
Definition: libcris.h:161
pert_t::time
double time[PERT_NTRACK][PERT_NXTRACK][PERT_NFOV]
Time (seconds since 2000-01-01T00:00Z).
Definition: libcris.h:143
pert_t::bt
double bt[PERT_NTRACK][PERT_NXTRACK][PERT_NFOV]
Brightness temperature (4.3 or 15 micron) [K].
Definition: libcris.h:155
pert_t::nxtrack
int nxtrack
Number of across-track values.
Definition: libcris.h:137
pert_t::lon
double lon[PERT_NTRACK][PERT_NXTRACK][PERT_NFOV]
Longitude [deg].
Definition: libcris.h:146
Here is the call graph for this function:
Generated on Wed Feb 5 2025 20:00:07 for CrIS Code Collection by doxygen 1.9.4