33 {
34
35 static atm_t atm_i, atm_apr;
37 static obs_t obs_i, obs_meas;
39
40 FILE *dirlist;
41
42
43 if (argc < 3)
44 ERRMSG(
"Give parameters: <ctl> <dirlist>");
45
46
48
49
52
53
55
56
57 if (!(dirlist = fopen(argv[2], "r")))
58 ERRMSG(
"Cannot open directory list!");
59
60
61 while (fscanf(dirlist,
"%4999s", ret.
dir) != EOF) {
62
63
64 LOG(1,
"\nRetrieve in directory %s...\n", ret.
dir);
65
66
68
69
71
72
73 double chisq;
75 &chisq);
76
77
79 }
80
81
82 LOG(1,
"\nRetrieval done...");
83
84
86
87
88 free(tbl);
89
90 return EXIT_SUCCESS;
91}
void read_ctl(int argc, char *argv[], ctl_t *ctl)
Read model control parameters from command-line and configuration input.
void optimal_estimation(ret_t *ret, ctl_t *ctl, tbl_t *tbl, obs_t *obs_meas, obs_t *obs_i, atm_t *atm_apr, atm_t *atm_i, double *chisq)
Perform optimal estimation retrieval using Levenberg–Marquardt minimization.
void read_atm(const char *dirname, const char *filename, const ctl_t *ctl, atm_t *atm)
Read atmospheric input data from a file.
void read_obs(const char *dirname, const char *filename, const ctl_t *ctl, obs_t *obs)
Read observation data from an input file.
void read_ret(int argc, char *argv[], const ctl_t *ctl, ret_t *ret)
Read retrieval configuration and error parameters.
tbl_t * read_tbl(const ctl_t *ctl)
Read all emissivity lookup tables for all gases and frequencies.
#define ERRMSG(...)
Print an error message with contextual information and terminate the program.
#define TIMER(name, mode)
Start or stop a named timer.
#define LOG(level,...)
Print a log message with a specified logging level.
Atmospheric profile data.
Observation geometry and radiance data.
Retrieval control parameters.
char dir[LEN]
Working directory.
Emissivity look-up tables.
1.9.4