Due to the lapse in federal government funding, NASA is not updating this website. We sincerely regret this inconvenience.
A sphere with three ocean waves in differing shades of blue. Next to the sphere, there is the bolded text "Ocean Color". Under that, there are two acronyms, separated by a vertical pipe: OB.DAAC (Ocean Biology Distributed Active Archive Center) and OBPG (Ocean Biology Processing Group). Lorem ipsum dolor sit amet
NASA Logo
Ocean Color Science Software

ocssw V2022
rdatreminfo.c
Go to the documentation of this file.
1 /* ========================================================== */
2 /* rdatreminfo() - reads parameters from ATREM sensor files */
3 /* */
4 /* Inputs: */
5 /* sensorID - sensor ID as defined in l12_parms.h */
6 /* pname - parameter name from sensor table */
7 /* */
8 /* Outputs: */
9 /* pval - pointer to scalar or array containing param */
10 /* */
11 /* Written By: */
12 /* R. Healy, March 2015 */
13 /* Based on rdsensorinfo.c */
14 /* */
15 /* ========================================================== */
16 
17 #include <stdio.h>
18 #include <string.h>
19 #include <stdlib.h>
20 #include <ctype.h>
21 #include "sensorInfo.h"
22 #include "genutils.h"
23 
24 // this is really defined in l2gen/l12_params.h
25 // but I need it here also.
26 #define NEWSENSINFO 32 /* use test sensor info file */
27 void parseline(char *line, char *name, char *value);
28 
38 int32_t rdatreminfo(int32_t sensorID, int32_t evalmask, const char *pname, void **pval) {
39  FILE *fp;
40  char *filedir;
41  char filename[FILENAME_MAX];
42  char line[80];
43  char name[80];
44  char value[80];
45  char param[80];
46  char *p;
47  int i;
48  int status;
49 
50  static int32_t sensorID_s = -999;
51  static int32_t nbands = 0;
52 
53  static float *fwhm, *fwave;
54  static int32_t h2o, co2, co, ch4, no2, n2o, o3, o2;
55  static float vrto3, sno2;
56  static float window1, window2, window3, window4, wp94c, w1p14c;
57  static int32_t nb1, nb2, nb3, nb4, nbp94, nb1p14, full_calc, dogeom;
58 
59  want_verbose = 0;
60  if (sensorID != sensorID_s) {
61 
62  if (want_verbose)
63  printf("\nLoading atrem (atmospheric removal - water vapor,etc.) information for %s\n",
64  sensorId2SensorName(sensorID));
65 
66  /* */
67  /* Locate atmocor datafile using environment variable */
68  /* */
69  if ((filedir = getenv("OCDATAROOT")) == NULL) {
70  printf("-E- %s: OCDATAROOT env variable undefined.\n", __FILE__);
71  return (-1);
72  }
73 
74  strcpy(filename, filedir);
75  if ((evalmask & NEWSENSINFO) != 0)
76  strcat(filename, "/eval/");
77  else
78  strcat(filename, "/");
79  strcat(filename, sensorId2SensorDir(sensorID));
80 
81  strcat(filename, "/msl12_atrem_info.dat");
82 
83  if (want_verbose)
84  printf("Opening atrem information file %s\n Looking for %s\n", filename, pname);
85 
86  if ((fp = fopen(filename, "r")) == NULL) {
87  fprintf(stderr, "-E- %s line %d: unable to open %s for reading\n",
88  __FILE__, __LINE__, filename);
89  return (-1);
90  }
91 
92  // read file to determine number of bands
93 
94  while (fgets(line, 80, fp)) {
95 
96  memset(name, '\0', sizeof (name));
97  memset(value, '\0', sizeof (value));
98 
99  // skip comment lines, empty lines, and lines without a name = value pair.
100 
101  if (line[0] == '#' || line[0] == '\n')
102  continue;
103  if (!(p = strchr(line, '=')))
104  continue;
105 
106  // parse parameter name and value, looking for Nbands
107 
108  parseline(line, name, value);
109 
110  if (strcmp(name, "Nbands") == 0) {
111  nbands = (int32_t) atoi(value);
112  break;
113  }
114  }
115 
116  if (nbands <= 0) {
117  printf("-E- %s: error find Nbands in %s.\n", __FILE__, filename);
118  return (-1);
119  }
120 
121  // allocate space for static data
122 
123  status = 0;
124  if ((fwhm = (float*) calloc(nbands, sizeof (float))) == NULL)
125  status = 1;
126  if ((fwave = (float*) calloc(nbands, sizeof (float))) == NULL)
127  status = 1;
128  if (status == 1) {
129  printf(
130  "-E- %s: error allocating space for %d bands in sensor info.\n",
131  __FILE__, nbands);
132  return (-1);
133  }
134 
135 
136  // Loop through each line of file looking for identifiers
137 
138  while (fgets(line, 80, fp)) {
139 
140  memset(name, '\0', sizeof (name));
141  memset(value, '\0', sizeof (value));
142 
143  // skip comment lines, empty lines, and lines without a name = value pair.
144 
145  if (line[0] == '#' || line[0] == '\n')
146  continue;
147  if (!(p = strchr(line, '=')))
148  continue;
149 
150  // parse parameter name and value
151 
152  parseline(line, name, value);
153 
154  // Copy value to appropriate variable
155 
156  if (strcmp(name, "Nbands") == 0)
157  nbands = (int32_t) atoi(value);
158  else if (strcmp(name, "h2o") == 0)
159  h2o = (int32_t) atoi(value);
160  else if (strcmp(name, "o2") == 0)
161  o2 = (int32_t) atoi(value);
162  else if (strcmp(name, "o3") == 0)
163  o3 = (int32_t) atoi(value);
164  else if (strcmp(name, "ch4") == 0)
165  ch4 = (int32_t) atoi(value);
166  else if (strcmp(name, "co") == 0)
167  co = (int32_t) atoi(value);
168  else if (strcmp(name, "n2o") == 0)
169  n2o = (int32_t) atoi(value);
170  else if (strcmp(name, "no2") == 0)
171  no2 = (int32_t) atoi(value);
172  else if (strcmp(name, "co2") == 0)
173  co2 = (int32_t) atoi(value);
174  else if (strcmp(name, "vrto3") == 0)
175  vrto3 = (float) atof(value);
176  else if (strcmp(name, "sno2") == 0)
177  sno2 = (float) atof(value);
178  else if (strcmp(name, "window1") == 0)
179  window1 = (float) atof(value);
180  else if (strcmp(name, "window2") == 0)
181  window2 = (float) atof(value);
182  else if (strcmp(name, "window3") == 0)
183  window3 = (float) atof(value);
184  else if (strcmp(name, "window4") == 0)
185  window4 = (float) atof(value);
186  else if (strcmp(name, "w1p14c") == 0)
187  w1p14c = (float) atof(value);
188  else if (strcmp(name, "wp94c") == 0)
189  wp94c = (float) atof(value);
190  else if (strcmp(name, "nb1") == 0)
191  nb1 = (int32_t) atoi(value);
192  else if (strcmp(name, "nb2") == 0)
193  nb2 = (int32_t) atoi(value);
194  else if (strcmp(name, "nb3") == 0)
195  nb3 = (int32_t) atoi(value);
196  else if (strcmp(name, "nb4") == 0)
197  nb4 = (int32_t) atoi(value);
198  else if (strcmp(name, "nbp94") == 0)
199  nbp94 = (int32_t) atoi(value);
200  else if (strcmp(name, "nb1p14") == 0)
201  nb1p14 = (int32_t) atoi(value);
202  else if (strcmp(name, "full_calc") == 0)
203  full_calc = (int32_t) atoi(value);
204  else if (strcmp(name, "dogeom") == 0)
205  dogeom = (int32_t) atoi(value);
206  else {
207  for (i = 0; i < nbands; i++) {
208  sprintf(param, "fwhm(%d)", i + 1);
209  if (strcmp(name, param) == 0) {
210  fwhm[i] = (float) atof(value);
211  break;
212  }
213  sprintf(param, "Lambda(%d)", i + 1);
214  if (strcmp(name, param) == 0) {
215  fwave[i] = (float) atof(value);
216  break;
217  }
218 
219  }
220  }
221 
222  }
223 
224 
225  /* */
226  /* Write-out what was read, for informational purposes */
227  /* */
228  if (want_verbose) {
229 
230  printf("window1=%f\n", window1);
231  printf(" %3s %8s\n", "Bnd",
232  "Fwhm");
233  for (i = 0; i < nbands; i++)
234  printf(
235  " %3d %8.3f\n",
236  i, fwhm[i]);
237  printf("\n");
238 
239  printf("\n");
240  } // want_verbose
241  }
242 
243  /* */
244  /* Return a pointer to the requested parameter */
245  /* */
246  if (pval != NULL) {
247  if (strcmp(pname, "Nbands") == 0)
248  *pval = (void *) &nbands;
249  else if (strcmp(pname, "fwhm") == 0)
250  *pval = (void *) fwhm;
251  else if (strcmp(pname, "Lambda") == 0)
252  *pval = (void *) fwave;
253  else if (strcmp(pname, "o2") == 0)
254  *pval = (void *) &o2;
255  else if (strcmp(pname, "o3") == 0)
256  *pval = (void *) &o3;
257  else if (strcmp(pname, "n2o") == 0)
258  *pval = (void *) &n2o;
259  else if (strcmp(pname, "ch4") == 0)
260  *pval = (void *) &ch4;
261  else if (strcmp(pname, "co") == 0)
262  *pval = (void *) &co;
263  else if (strcmp(pname, "no2") == 0)
264  *pval = (void *) &no2;
265  else if (strcmp(pname, "h2o") == 0)
266  *pval = (void *) &h2o;
267  else if (strcmp(pname, "co2") == 0)
268  *pval = (void *) &co2;
269  else if (strcmp(pname, "vrto3") == 0)
270  *pval = (void *) &vrto3;
271  else if (strcmp(pname, "sno2") == 0)
272  *pval = (void *) &sno2;
273  else if (strcmp(pname, "window1") == 0)
274  *pval = (void *) &window1;
275  else if (strcmp(pname, "window2") == 0)
276  *pval = (void *) &window2;
277  else if (strcmp(pname, "window3") == 0)
278  *pval = (void *) &window3;
279  else if (strcmp(pname, "window4") == 0)
280  *pval = (void *) &window4;
281  else if (strcmp(pname, "wp94c") == 0)
282  *pval = (void *) &wp94c;
283  else if (strcmp(pname, "w1p14c") == 0)
284  *pval = (void *) &w1p14c;
285  else if (strcmp(pname, "nb1") == 0)
286  *pval = (void *) &nb1;
287  else if (strcmp(pname, "nb2") == 0)
288  *pval = (void *) &nb2;
289  else if (strcmp(pname, "nb3") == 0)
290  *pval = (void *) &nb3;
291  else if (strcmp(pname, "nb4") == 0)
292  *pval = (void *) &nb4;
293  else if (strcmp(pname, "nbp94") == 0)
294  *pval = (void *) &nbp94;
295  else if (strcmp(pname, "nb1p14") == 0)
296  *pval = (void *) &nb1p14;
297  else if (strcmp(pname, "full_calc") == 0)
298  *pval = (void *) &full_calc;
299  else if (strcmp(pname, "dogeom") == 0)
300  *pval = (void *) &dogeom;
301 
302  else
303  return (-1);
304  }
305 
306  sensorID_s = sensorID;
307 
308  if (pname != NULL)
309  return (nbands);
310 
311  return (0);
312 }
313 
value
int32 value
Definition: Granule.c:1235
sensorId2SensorDir
const char * sensorId2SensorDir(int sensorId)
Definition: sensorInfo.c:315
co
int32_t co
Definition: atrem_corl1.h:118
full_calc
int32_t full_calc
Definition: atrem_corl1.h:148
status
int status
Definition: l1_czcs_hdf.c:32
o3
int32_t o3
Definition: atrem_corl1.h:118
co2
int32_t co2
Definition: atrem_corl1.h:118
NULL
#define NULL
Definition: decode_rs.h:63
sno2
float sno2
Definition: atrem_corl1.h:133
rdatreminfo
int32_t rdatreminfo(int32_t sensorID, int32_t evalmask, const char *pname, void **pval)
Definition: rdatreminfo.c:38
nb4
int32_t nb4
Definition: atrem_cor.h:98
nbp94
int32_t nbp94
Definition: atrem_cor.h:98
fwhm
float fwhm[NBANDS]
Definition: atrem_corl1.h:144
vrto3
float vrto3
Definition: atrem_corl1.h:133
parseline
void parseline(char *line, char *name, char *value)
Definition: rdsensorinfo.c:32
h2o
int32_t h2o
Definition: atrem_corl1.h:118
no2
int32_t no2
Definition: atrem_corl1.h:118
nb1
int32_t nb1
Definition: atrem_cor.h:98
nb2
int32_t nb2
Definition: atrem_cor.h:98
w1p14c
float w1p14c
Definition: atrem_corl1.h:153
want_verbose
int want_verbose
Definition: genutils_globals.c:3
filename
char filename[FILENAME_MAX]
Definition: atrem_corl1.h:122
o2
int32_t o2
Definition: atrem_corl1.h:118
fix_mac_rpath.line
line
Definition: fix_mac_rpath.py:56
wp94c
float wp94c
Definition: atrem_corl1.h:153
NEWSENSINFO
#define NEWSENSINFO
Definition: rdatreminfo.c:26
nbands
int32_t nbands
Definition: atrem_corl1v3.h:190
sensorId2SensorName
const char * sensorId2SensorName(int sensorId)
Definition: sensorInfo.c:273
ch4
int32_t ch4
Definition: atrem_corl1.h:118
i
int i
Definition: decode_rs.h:71
strcpy
How many dimensions is the output array Default is Not sure if anything above will work correctly strcpy(l2prod->title, "no title yet")
sensorID
int32_t sensorID[MAXNFILES]
Definition: l2bin.cpp:91
nb1p14
int32_t nb1p14
Definition: atrem_cor.h:98
nb3
int32_t nb3
Definition: atrem_cor.h:98
p
float p[MODELMAX]
Definition: atrem_corl1.h:131
n2o
int32_t n2o
Definition: atrem_corl1.h:118


ABOUT

Our Team

Contact Us

DATA

Get Started

Find Data

Use Data

TOOLS

SeaBASS

SeaDAS

Giovanni Data Visualizer

Worldview Mapping Interface

Mission Quality Monitor

HELP

How to Use Earthdata Search

Level 3 & 4 Browser Help

File Search Help

Search & Download Help

Earthdata Forum

Join Mailing List

How to Cite Data

CONNECT WITH US
COMPLIANCE

Web Privacy Policy

Data & Information Policy

Communications Policy

Accessibility

Freedom of Information Act

USA.gov

Responsible NASA Official: Sean Bailey
Curator: OceanColor Webmaster