1 | /********************************************************************* |
---|
2 | * Copyright 1993, UCAR/Unidata |
---|
3 | * See netcdf/COPYRIGHT file for copying and redistribution conditions. |
---|
4 | * $Header: /upc/share/CVS/netcdf-3/libncdap3/daputil.c,v 1.47 2010/05/21 23:24:15 dmh Exp $ |
---|
5 | *********************************************************************/ |
---|
6 | |
---|
7 | #include "config.h" |
---|
8 | |
---|
9 | #include <sys/time.h> |
---|
10 | |
---|
11 | #include "oc.h" |
---|
12 | extern int oc_dumpnode(OClink, OCobject); |
---|
13 | |
---|
14 | #include "ncdap3.h" |
---|
15 | #include "dapalign.h" |
---|
16 | #include "dapodom.h" |
---|
17 | |
---|
18 | #define LBRACKET '[' |
---|
19 | #define RBRACKET ']' |
---|
20 | |
---|
21 | |
---|
22 | /**************************************************/ |
---|
23 | /** |
---|
24 | * Provide a hidden interface to allow utilities |
---|
25 | * to check if a given path name is really an ncdap3 url. |
---|
26 | * If no, return null, else return basename of the url |
---|
27 | * minus any extension. |
---|
28 | */ |
---|
29 | |
---|
30 | int |
---|
31 | nc__testurl(const char* path, char** basenamep) |
---|
32 | { |
---|
33 | NC_URI* uri; |
---|
34 | int ok = nc_uriparse(path,&uri); |
---|
35 | if(ok) { |
---|
36 | char* slash = strrchr(uri->file, '/'); |
---|
37 | char* dot; |
---|
38 | if(slash == NULL) slash = (char*)path; else slash++; |
---|
39 | slash = nulldup(slash); |
---|
40 | dot = strrchr(slash, '.'); |
---|
41 | if(dot != NULL && dot != slash) *dot = '\0'; |
---|
42 | if(basenamep) *basenamep=slash ; else free(slash); |
---|
43 | nc_urifree(uri); |
---|
44 | } |
---|
45 | return ok; |
---|
46 | } |
---|
47 | |
---|
48 | /**************************************************/ |
---|
49 | |
---|
50 | /* |
---|
51 | Given a legal dap name with arbitrary characters, |
---|
52 | convert to equivalent legal cdf name |
---|
53 | With the new name policy for netcdf, this procedure |
---|
54 | does nothing. |
---|
55 | */ |
---|
56 | |
---|
57 | char* |
---|
58 | cdflegalname3(char* dapname) |
---|
59 | { |
---|
60 | return nulldup(dapname); |
---|
61 | } |
---|
62 | |
---|
63 | /* Define the type conversion of the DAP variables |
---|
64 | to the external netCDF variable type. |
---|
65 | The proper way is to, for example, convert unsigned short |
---|
66 | to an int to maintain the values. |
---|
67 | Unfortuneately, libnc-dap does not do this: |
---|
68 | it translates the types directly. For example |
---|
69 | libnc-dap upgrades the DAP byte type, which is unsigned char, |
---|
70 | to NC_BYTE, which signed char. |
---|
71 | Oh well. |
---|
72 | For netcdf-4, we can do proper type conversion. |
---|
73 | */ |
---|
74 | nc_type |
---|
75 | nctypeconvert(NCDAPCOMMON* drno, nc_type nctype) |
---|
76 | { |
---|
77 | nc_type upgrade = NC_NAT; |
---|
78 | if(drno->controls.flags & NCF_NC3) { |
---|
79 | /* libnc-dap mimic invariant is to maintain type size */ |
---|
80 | switch (nctype) { |
---|
81 | case NC_CHAR: upgrade = NC_CHAR; break; |
---|
82 | case NC_BYTE: upgrade = NC_BYTE; break; |
---|
83 | case NC_UBYTE: upgrade = NC_BYTE; break; |
---|
84 | case NC_SHORT: upgrade = NC_SHORT; break; |
---|
85 | case NC_USHORT: upgrade = NC_SHORT; break; |
---|
86 | case NC_INT: upgrade = NC_INT; break; |
---|
87 | case NC_UINT: upgrade = NC_INT; break; |
---|
88 | case NC_INT64: upgrade = NC_INT64; break; |
---|
89 | case NC_UINT64: upgrade = NC_UINT64; break; |
---|
90 | case NC_FLOAT: upgrade = NC_FLOAT; break; |
---|
91 | case NC_DOUBLE: upgrade = NC_DOUBLE; break; |
---|
92 | case NC_URL: |
---|
93 | case NC_STRING: upgrade = NC_CHAR; break; |
---|
94 | default: break; |
---|
95 | } |
---|
96 | } else if(drno->controls.flags & NCF_NC4) { |
---|
97 | /* netcdf-4 conversion is more correct */ |
---|
98 | switch (nctype) { |
---|
99 | case NC_CHAR: upgrade = NC_CHAR; break; |
---|
100 | case NC_BYTE: upgrade = NC_BYTE; break; |
---|
101 | case NC_UBYTE: upgrade = NC_UBYTE; break; |
---|
102 | case NC_SHORT: upgrade = NC_SHORT; break; |
---|
103 | case NC_USHORT: upgrade = NC_USHORT; break; |
---|
104 | case NC_INT: upgrade = NC_INT; break; |
---|
105 | case NC_UINT: upgrade = NC_UINT; break; |
---|
106 | case NC_INT64: upgrade = NC_INT64; break; |
---|
107 | case NC_UINT64: upgrade = NC_UINT64; break; |
---|
108 | case NC_FLOAT: upgrade = NC_FLOAT; break; |
---|
109 | case NC_DOUBLE: upgrade = NC_DOUBLE; break; |
---|
110 | case NC_URL: |
---|
111 | case NC_STRING: upgrade = NC_STRING; break; |
---|
112 | default: break; |
---|
113 | } |
---|
114 | } |
---|
115 | return upgrade; |
---|
116 | } |
---|
117 | |
---|
118 | nc_type |
---|
119 | octypetonc(OCtype etype) |
---|
120 | { |
---|
121 | switch (etype) { |
---|
122 | case OC_Char: return NC_CHAR; |
---|
123 | case OC_Byte: return NC_UBYTE; |
---|
124 | case OC_UByte: return NC_UBYTE; |
---|
125 | case OC_Int16: return NC_SHORT; |
---|
126 | case OC_UInt16: return NC_USHORT; |
---|
127 | case OC_Int32: return NC_INT; |
---|
128 | case OC_UInt32: return NC_UINT; |
---|
129 | case OC_Int64: return NC_INT64; |
---|
130 | case OC_UInt64: return NC_UINT64; |
---|
131 | case OC_Float32: return NC_FLOAT; |
---|
132 | case OC_Float64: return NC_DOUBLE; |
---|
133 | case OC_String: return NC_STRING; |
---|
134 | case OC_URL: return NC_STRING; |
---|
135 | case OC_Dataset: return NC_Dataset; |
---|
136 | case OC_Sequence: return NC_Sequence; |
---|
137 | case OC_Structure: return NC_Structure; |
---|
138 | case OC_Grid: return NC_Grid; |
---|
139 | case OC_Dimension: return NC_Dimension; |
---|
140 | case OC_Primitive: return NC_Primitive; |
---|
141 | default: break; |
---|
142 | } |
---|
143 | return NC_NAT; |
---|
144 | } |
---|
145 | |
---|
146 | OCtype |
---|
147 | nctypetodap(nc_type nctype) |
---|
148 | { |
---|
149 | switch (nctype) { |
---|
150 | case NC_CHAR: return OC_Char; |
---|
151 | case NC_BYTE: return OC_Byte; |
---|
152 | case NC_UBYTE: return OC_UByte; |
---|
153 | case NC_SHORT: return OC_Int16; |
---|
154 | case NC_USHORT: return OC_UInt16; |
---|
155 | case NC_INT: return OC_Int32; |
---|
156 | case NC_UINT: return OC_UInt32; |
---|
157 | case NC_INT64: return OC_Int64; |
---|
158 | case NC_UINT64: return OC_UInt64; |
---|
159 | case NC_FLOAT: return OC_Float32; |
---|
160 | case NC_DOUBLE: return OC_Float64; |
---|
161 | case NC_STRING: return OC_String; |
---|
162 | default : break; |
---|
163 | } |
---|
164 | return OC_NAT; |
---|
165 | } |
---|
166 | |
---|
167 | size_t |
---|
168 | nctypesizeof(nc_type nctype) |
---|
169 | { |
---|
170 | switch (nctype) { |
---|
171 | case NC_CHAR: return sizeof(char); |
---|
172 | case NC_BYTE: return sizeof(signed char); |
---|
173 | case NC_UBYTE: return sizeof(unsigned char); |
---|
174 | case NC_SHORT: return sizeof(short); |
---|
175 | case NC_USHORT: return sizeof(unsigned short); |
---|
176 | case NC_INT: return sizeof(int); |
---|
177 | case NC_UINT: return sizeof(unsigned int); |
---|
178 | case NC_INT64: return sizeof(long long); |
---|
179 | case NC_UINT64: return sizeof(unsigned long long); |
---|
180 | case NC_FLOAT: return sizeof(float); |
---|
181 | case NC_DOUBLE: return sizeof(double); |
---|
182 | case NC_STRING: return sizeof(char*); |
---|
183 | default: PANIC("nctypesizeof"); |
---|
184 | } |
---|
185 | return 0; |
---|
186 | } |
---|
187 | |
---|
188 | char* |
---|
189 | nctypetostring(nc_type nctype) |
---|
190 | { |
---|
191 | switch (nctype) { |
---|
192 | case NC_NAT: return "NC_NAT"; |
---|
193 | case NC_BYTE: return "NC_BYTE"; |
---|
194 | case NC_CHAR: return "NC_CHAR"; |
---|
195 | case NC_SHORT: return "NC_SHORT"; |
---|
196 | case NC_INT: return "NC_INT"; |
---|
197 | case NC_FLOAT: return "NC_FLOAT"; |
---|
198 | case NC_DOUBLE: return "NC_DOUBLE"; |
---|
199 | case NC_UBYTE: return "NC_UBYTE"; |
---|
200 | case NC_USHORT: return "NC_USHORT"; |
---|
201 | case NC_UINT: return "NC_UINT"; |
---|
202 | case NC_INT64: return "NC_INT64"; |
---|
203 | case NC_UINT64: return "NC_UINT64"; |
---|
204 | case NC_STRING: return "NC_STRING"; |
---|
205 | case NC_VLEN: return "NC_VLEN"; |
---|
206 | case NC_OPAQUE: return "NC_OPAQUE"; |
---|
207 | case NC_ENUM: return "NC_ENUM"; |
---|
208 | case NC_COMPOUND: return "NC_COMPOUND"; |
---|
209 | case NC_URL: return "NC_URL"; |
---|
210 | case NC_SET: return "NC_SET"; |
---|
211 | case NC_Dataset: return "NC_Dataset"; |
---|
212 | case NC_Sequence: return "NC_Sequence"; |
---|
213 | case NC_Structure: return "NC_Structure"; |
---|
214 | case NC_Grid: return "NC_Grid"; |
---|
215 | case NC_Dimension: return "NC_Dimension"; |
---|
216 | case NC_Primitive: return "NC_Primitive"; |
---|
217 | default: break; |
---|
218 | } |
---|
219 | return NULL; |
---|
220 | } |
---|
221 | |
---|
222 | |
---|
223 | /* Pad a buffer */ |
---|
224 | int |
---|
225 | alignbuffer3(NCbytes* buf, int alignment) |
---|
226 | { |
---|
227 | int pad; |
---|
228 | unsigned long len; |
---|
229 | if(buf == NULL) return 0; |
---|
230 | len = ncbyteslength(buf); |
---|
231 | pad = nccpadding(len,alignment); |
---|
232 | |
---|
233 | #ifdef TEST |
---|
234 | for(;pad > 0;pad--) |
---|
235 | ncbytesappend(buf,0x3a); /* 0x3a was chosen at random */ |
---|
236 | #else |
---|
237 | ncbytessetlength(buf,len+pad); |
---|
238 | #endif |
---|
239 | return 1; |
---|
240 | } |
---|
241 | |
---|
242 | size_t |
---|
243 | dimproduct3(NClist* dimensions) |
---|
244 | { |
---|
245 | size_t size = 1; |
---|
246 | unsigned int i; |
---|
247 | if(dimensions == NULL) return size; |
---|
248 | for(i=0;i<nclistlength(dimensions);i++) { |
---|
249 | CDFnode* dim = (CDFnode*)nclistget(dimensions,i); |
---|
250 | size *= dim->dim.declsize; |
---|
251 | } |
---|
252 | return size; |
---|
253 | } |
---|
254 | |
---|
255 | |
---|
256 | /* Return value of param or NULL if not found */ |
---|
257 | const char* |
---|
258 | paramvalue34(NCDAPCOMMON* nccomm, const char* key) |
---|
259 | { |
---|
260 | const char* value; |
---|
261 | |
---|
262 | if(nccomm == NULL || key == NULL) return 0; |
---|
263 | if(!nc_urilookup(nccomm->oc.url,key,&value)) |
---|
264 | return NULL; |
---|
265 | return value; |
---|
266 | } |
---|
267 | |
---|
268 | static const char* checkseps = "+,:;"; |
---|
269 | |
---|
270 | /* Search for substring in value of param. If substring == NULL; then just |
---|
271 | check if param is defined. |
---|
272 | */ |
---|
273 | int |
---|
274 | paramcheck34(NCDAPCOMMON* nccomm, const char* key, const char* subkey) |
---|
275 | { |
---|
276 | const char* value; |
---|
277 | char* p; |
---|
278 | |
---|
279 | if(nccomm == NULL || key == NULL) return 0; |
---|
280 | if(!nc_urilookup(nccomm->oc.url,key,&value)) |
---|
281 | return 0; |
---|
282 | if(subkey == NULL) return 1; |
---|
283 | p = strstr(value,subkey); |
---|
284 | if(p == NULL) return 0; |
---|
285 | p += strlen(subkey); |
---|
286 | if(*p != '\0' && strchr(checkseps,*p) == NULL) return 0; |
---|
287 | return 1; |
---|
288 | } |
---|
289 | |
---|
290 | |
---|
291 | /* This is NOT UNION */ |
---|
292 | int |
---|
293 | nclistconcat(NClist* l1, NClist* l2) |
---|
294 | { |
---|
295 | unsigned int i; |
---|
296 | for(i=0;i<nclistlength(l2);i++) nclistpush(l1,nclistget(l2,i)); |
---|
297 | return 1; |
---|
298 | } |
---|
299 | |
---|
300 | int |
---|
301 | nclistminus(NClist* l1, NClist* l2) |
---|
302 | { |
---|
303 | unsigned int i,len,found; |
---|
304 | len = nclistlength(l2); |
---|
305 | found = 0; |
---|
306 | for(i=0;i<len;i++) { |
---|
307 | if(nclistdeleteall(l1,nclistget(l2,i))) found = 1; |
---|
308 | } |
---|
309 | return found; |
---|
310 | } |
---|
311 | |
---|
312 | int |
---|
313 | nclistdeleteall(NClist* l, ncelem elem) |
---|
314 | { |
---|
315 | int i; /* do not make unsigned */ |
---|
316 | unsigned int len,found; |
---|
317 | found = 0; |
---|
318 | len = nclistlength(l); |
---|
319 | for(i=len-1;i>=0;i--) { |
---|
320 | ncelem test = nclistget(l,i); |
---|
321 | if(test==elem) { |
---|
322 | nclistremove(l,i); |
---|
323 | found=1; |
---|
324 | } |
---|
325 | } |
---|
326 | return found; |
---|
327 | } |
---|
328 | |
---|
329 | /* Collect the set of container nodes ending in "container"*/ |
---|
330 | void |
---|
331 | collectnodepath3(CDFnode* node, NClist* path, int withdataset) |
---|
332 | { |
---|
333 | if(node == NULL) return; |
---|
334 | nclistpush(path,(ncelem)node); |
---|
335 | while(node->container != NULL) { |
---|
336 | node = node->container; |
---|
337 | if(!withdataset && node->nctype == NC_Dataset) break; |
---|
338 | nclistinsert(path,0,(ncelem)node); |
---|
339 | } |
---|
340 | } |
---|
341 | |
---|
342 | /* Like collectnodepath3, but in ocspace */ |
---|
343 | void |
---|
344 | collectocpath(OCconnection conn, OCobject node, NClist* path) |
---|
345 | { |
---|
346 | OCobject container; |
---|
347 | OCtype octype; |
---|
348 | if(node == OCNULL) return; |
---|
349 | oc_inq_class(conn,node,&octype); |
---|
350 | if(octype == OC_Dataset) return; |
---|
351 | oc_inq_container(conn,node,&container); |
---|
352 | if(container != OCNULL) |
---|
353 | collectocpath(conn,container,path); |
---|
354 | nclistpush(path,(ncelem)node); |
---|
355 | } |
---|
356 | |
---|
357 | char* |
---|
358 | makeocpathstring3(OCconnection conn, OCobject node, const char* sep) |
---|
359 | { |
---|
360 | int slen,i,len,first,seplen; |
---|
361 | char* pathname; |
---|
362 | OCtype octype; |
---|
363 | NClist* ocpath = nclistnew(); |
---|
364 | |
---|
365 | collectocpath(conn,node,ocpath); |
---|
366 | len = nclistlength(ocpath); |
---|
367 | assert(len > 0); /* dataset at least */ |
---|
368 | |
---|
369 | oc_inq_type(conn,node,&octype); |
---|
370 | if(octype == OC_Dataset) |
---|
371 | {pathname = nulldup(""); goto done;} /* Dataset */ |
---|
372 | |
---|
373 | seplen = strlen(sep); |
---|
374 | for(slen=0,i=0;i<len;i++) { |
---|
375 | OCobject node = (OCobject)nclistget(ocpath,i); |
---|
376 | char* name; |
---|
377 | oc_inq_type(conn,node,&octype); |
---|
378 | if(octype == OC_Dataset) continue; |
---|
379 | oc_inq_name(conn,node,&name); |
---|
380 | slen += (name == NULL? 0 : strlen(name)); |
---|
381 | slen += seplen; |
---|
382 | nullfree(name); |
---|
383 | } |
---|
384 | slen += 1; /* for null terminator*/ |
---|
385 | pathname = (char*)malloc(slen); |
---|
386 | MEMCHECK(pathname,NULL); |
---|
387 | pathname[0] = '\0'; |
---|
388 | for(first=1,i=0;i<len;i++) { |
---|
389 | OCobject node = (OCobject)nclistget(ocpath,i); |
---|
390 | char* name; |
---|
391 | oc_inq_type(conn,node,&octype); |
---|
392 | if(octype == OC_Dataset) continue; |
---|
393 | oc_inq_name(conn,node,&name); |
---|
394 | if(!first) strcat(pathname,sep); |
---|
395 | if(name != NULL) strcat(pathname,name); |
---|
396 | nullfree(name); |
---|
397 | first = 0; |
---|
398 | } |
---|
399 | done: |
---|
400 | nclistfree(ocpath); |
---|
401 | return pathname; |
---|
402 | } |
---|
403 | |
---|
404 | char* |
---|
405 | makepathstring3(NClist* path, const char* separator, int flags) |
---|
406 | { |
---|
407 | int slen,i,len,first,seplen; |
---|
408 | char* pathname; |
---|
409 | |
---|
410 | len = nclistlength(path); |
---|
411 | ASSERT(len > 0); /* dataset at least */ |
---|
412 | seplen = strlen(separator); |
---|
413 | ASSERT(seplen > 0); |
---|
414 | for(slen=0,i=0;i<len;i++) { |
---|
415 | CDFnode* node = (CDFnode*)nclistget(path,i); |
---|
416 | if(node->nctype == NC_Dataset) continue; |
---|
417 | slen += strlen(node->ncbasename); |
---|
418 | slen += seplen; |
---|
419 | } |
---|
420 | slen += 1; /* for null terminator*/ |
---|
421 | pathname = (char*)malloc(slen); |
---|
422 | MEMCHECK(pathname,NULL); |
---|
423 | pathname[0] = '\0'; |
---|
424 | for(first=1,i=0;i<len;i++) { |
---|
425 | CDFnode* node = (CDFnode*)nclistget(path,i); |
---|
426 | char* name; |
---|
427 | if(!node->elided || (flags & PATHELIDE)==0) { |
---|
428 | if(node->nctype != NC_Dataset) { |
---|
429 | name = node->ncbasename; |
---|
430 | if(!first) strcat(pathname,separator); |
---|
431 | strcat(pathname,name); |
---|
432 | first = 0; |
---|
433 | } |
---|
434 | } |
---|
435 | } |
---|
436 | return pathname; |
---|
437 | } |
---|
438 | |
---|
439 | |
---|
440 | /* convert path to string using the ncname field */ |
---|
441 | char* |
---|
442 | makecdfpathstring3(CDFnode* var, const char* separator) |
---|
443 | { |
---|
444 | char* spath; |
---|
445 | NClist* path = nclistnew(); |
---|
446 | collectnodepath3(var,path,WITHDATASET); /* <= note */ |
---|
447 | spath = makepathstring3(path,separator,PATHNC); |
---|
448 | nclistfree(path); |
---|
449 | return spath; |
---|
450 | } |
---|
451 | |
---|
452 | /* Collect the set names of container nodes ending in "container"*/ |
---|
453 | void |
---|
454 | clonenodenamepath3(CDFnode* node, NClist* path, int withdataset) |
---|
455 | { |
---|
456 | if(node == NULL) return; |
---|
457 | /* stop at the dataset container as well*/ |
---|
458 | if(node->nctype != NC_Dataset) |
---|
459 | clonenodenamepath3(node->container,path,withdataset); |
---|
460 | if(node->nctype != NC_Dataset || withdataset) |
---|
461 | nclistpush(path,(ncelem)nulldup(node->ncbasename)); |
---|
462 | } |
---|
463 | |
---|
464 | char* |
---|
465 | simplepathstring3(NClist* names, char* separator) |
---|
466 | { |
---|
467 | int i; |
---|
468 | size_t len; |
---|
469 | char* result; |
---|
470 | if(nclistlength(names) == 0) return nulldup(""); |
---|
471 | for(len=0,i=0;i<nclistlength(names);i++) { |
---|
472 | char* name = (char*)nclistget(names,i); |
---|
473 | len += strlen(name); |
---|
474 | len += strlen(separator); |
---|
475 | } |
---|
476 | len++; /* null terminator */ |
---|
477 | result = (char*)malloc(len); |
---|
478 | result[0] = '\0'; |
---|
479 | for(i=0;i<nclistlength(names);i++) { |
---|
480 | char* segment = (char*)nclistget(names,i); |
---|
481 | if(i > 0) strcat(result,separator); |
---|
482 | strcat(result,segment); |
---|
483 | } |
---|
484 | return result; |
---|
485 | } |
---|
486 | |
---|
487 | /* Define a number of location tests */ |
---|
488 | |
---|
489 | /* Is node contained (transitively) in a sequence ? */ |
---|
490 | BOOL |
---|
491 | dapinsequence(CDFnode* node) |
---|
492 | { |
---|
493 | if(node == NULL || node->container == NULL) return TRUE; |
---|
494 | for(node=node->container;node->nctype != NC_Dataset;node=node->container) { |
---|
495 | if(node->nctype == NC_Sequence) return TRUE; |
---|
496 | } |
---|
497 | return FALSE; |
---|
498 | } |
---|
499 | |
---|
500 | /* Is node a map field of a grid? */ |
---|
501 | BOOL |
---|
502 | dapgridmap(CDFnode* node) |
---|
503 | { |
---|
504 | if(node != NULL && node->container != NULL |
---|
505 | && node->container->nctype == NC_Grid) { |
---|
506 | CDFnode* array = (CDFnode*)nclistget(node->container->subnodes,0); |
---|
507 | return (node != array); |
---|
508 | } |
---|
509 | return FALSE; |
---|
510 | } |
---|
511 | |
---|
512 | /* Is node an array field of a grid? */ |
---|
513 | BOOL |
---|
514 | dapgridarray(CDFnode* node) |
---|
515 | { |
---|
516 | if(node != NULL && node->container != NULL |
---|
517 | && node->container->nctype == NC_Grid) { |
---|
518 | CDFnode* array = (CDFnode*)nclistget(node->container->subnodes,0); |
---|
519 | return (node == array); |
---|
520 | } |
---|
521 | return FALSE; |
---|
522 | } |
---|
523 | |
---|
524 | BOOL |
---|
525 | dapgridelement(CDFnode* node) |
---|
526 | { |
---|
527 | return dapgridarray(node) |
---|
528 | || dapgridmap(node); |
---|
529 | } |
---|
530 | |
---|
531 | /* Is node a top-level grid node? */ |
---|
532 | BOOL |
---|
533 | daptopgrid(CDFnode* grid) |
---|
534 | { |
---|
535 | if(grid == NULL || grid->nctype != NC_Grid) return FALSE; |
---|
536 | return daptoplevel(grid); |
---|
537 | } |
---|
538 | |
---|
539 | /* Is node a top-level sequence node? */ |
---|
540 | BOOL |
---|
541 | daptopseq(CDFnode* seq) |
---|
542 | { |
---|
543 | if(seq == NULL || seq->nctype != NC_Sequence) return FALSE; |
---|
544 | return daptoplevel(seq); |
---|
545 | } |
---|
546 | |
---|
547 | /* Is node a top-level node? */ |
---|
548 | BOOL |
---|
549 | daptoplevel(CDFnode* node) |
---|
550 | { |
---|
551 | if(node->container == NULL |
---|
552 | || node->container->nctype != NC_Dataset) return FALSE; |
---|
553 | return TRUE; |
---|
554 | } |
---|
555 | |
---|
556 | unsigned int |
---|
557 | modeldecode(int translation, const char* smodel, |
---|
558 | const struct NCTMODEL* models, |
---|
559 | unsigned int dfalt) |
---|
560 | { |
---|
561 | for(;models->translation;models++) { |
---|
562 | if(translation != models->translation) continue; |
---|
563 | if(smodel == models->model |
---|
564 | || (models->model != NULL && strcasecmp(smodel,models->model)==0)) { |
---|
565 | /* We have a match */ |
---|
566 | return models->flags; |
---|
567 | } |
---|
568 | } |
---|
569 | return dfalt; |
---|
570 | } |
---|
571 | |
---|
572 | unsigned long |
---|
573 | getlimitnumber(const char* limit) |
---|
574 | { |
---|
575 | size_t slen; |
---|
576 | unsigned long multiplier = 1; |
---|
577 | unsigned long lu; |
---|
578 | |
---|
579 | if(limit == NULL) return 0; |
---|
580 | slen = strlen(limit); |
---|
581 | if(slen == 0) return 0; |
---|
582 | switch (limit[slen-1]) { |
---|
583 | case 'G': case 'g': multiplier = GIGBYTE; break; |
---|
584 | case 'M': case 'm': multiplier = MEGBYTE; break; |
---|
585 | case 'K': case 'k': multiplier = KILBYTE; break; |
---|
586 | default: break; |
---|
587 | } |
---|
588 | sscanf(limit,"%lu",&lu); |
---|
589 | return (lu*multiplier); |
---|
590 | } |
---|
591 | |
---|
592 | void |
---|
593 | dapexpandescapes(char *termstring) |
---|
594 | { |
---|
595 | char *s, *t, *endp; |
---|
596 | |
---|
597 | /* expand "\" escapes, e.g. "\t" to tab character; |
---|
598 | will only shorten string length, never increase it |
---|
599 | */ |
---|
600 | s = termstring; |
---|
601 | t = termstring; |
---|
602 | while(*t) { |
---|
603 | if (*t == '\\') { |
---|
604 | t++; |
---|
605 | switch (*t) { |
---|
606 | case 'a': |
---|
607 | *s++ = '\007'; t++; /* will use '\a' when STDC */ |
---|
608 | break; |
---|
609 | case 'b': |
---|
610 | *s++ = '\b'; t++; |
---|
611 | break; |
---|
612 | case 'f': |
---|
613 | *s++ = '\f'; t++; |
---|
614 | break; |
---|
615 | case 'n': |
---|
616 | *s++ = '\n'; t++; |
---|
617 | break; |
---|
618 | case 'r': |
---|
619 | *s++ = '\r'; t++; |
---|
620 | break; |
---|
621 | case 't': |
---|
622 | *s++ = '\t'; t++; |
---|
623 | break; |
---|
624 | case 'v': |
---|
625 | *s++ = '\v'; t++; |
---|
626 | break; |
---|
627 | case '\\': |
---|
628 | *s++ = '\\'; t++; |
---|
629 | break; |
---|
630 | case '?': |
---|
631 | *s++ = '\177'; t++; |
---|
632 | break; |
---|
633 | case 'x': |
---|
634 | t++; /* now t points to one or more hex digits */ |
---|
635 | *s++ = (char) strtol(t, &endp, 16); |
---|
636 | t = endp; |
---|
637 | break; |
---|
638 | case '0': |
---|
639 | case '1': |
---|
640 | case '2': |
---|
641 | case '3': |
---|
642 | case '4': |
---|
643 | case '5': |
---|
644 | case '6': |
---|
645 | case '7': { |
---|
646 | /* t should now point to 3 octal digits */ |
---|
647 | int c; |
---|
648 | c = t[0]; |
---|
649 | if(c == 0 || c < '0' || c > '7') goto normal; |
---|
650 | c = t[1]; |
---|
651 | if(c == 0 || c < '0' || c > '7') goto normal; |
---|
652 | c = t[2]; |
---|
653 | if(c == 0 || c < '0' || c > '7') goto normal; |
---|
654 | c = ((t[0]-'0')<<6)+((t[1]-'0')<<3)+(t[2]-'0'); |
---|
655 | *s++ = (char)c; |
---|
656 | t += 3; |
---|
657 | } break; |
---|
658 | default: |
---|
659 | if(*t == 0) |
---|
660 | *s++ = '\\'; |
---|
661 | else |
---|
662 | *s++ = *t++; |
---|
663 | break; |
---|
664 | } |
---|
665 | } else { |
---|
666 | normal: *s++ = *t++; |
---|
667 | } |
---|
668 | } |
---|
669 | *s = '\0'; |
---|
670 | return; |
---|
671 | } |
---|
672 | |
---|
673 | #ifdef HAVE_GETTIMEOFDAY |
---|
674 | static struct timeval time0; |
---|
675 | static struct timeval time1; |
---|
676 | |
---|
677 | static double |
---|
678 | deltatime() |
---|
679 | { |
---|
680 | double t0, t1; |
---|
681 | t0 = ((double)time0.tv_sec); |
---|
682 | t0 += ((double)time0.tv_usec) / 1000000.0; |
---|
683 | t1 = ((double)time1.tv_sec); |
---|
684 | t1 += ((double)time1.tv_usec) / 1000000.0; |
---|
685 | return (t1 - t0); |
---|
686 | } |
---|
687 | #endif |
---|
688 | |
---|
689 | /* Provide a wrapper for oc_fetch so we can log what it does */ |
---|
690 | OCerror |
---|
691 | dap_fetch(NCDAPCOMMON* nccomm, OCconnection conn, const char* ce, |
---|
692 | OCdxd dxd, OCobject* rootp) |
---|
693 | { |
---|
694 | OCerror ocstat; |
---|
695 | char* ext; |
---|
696 | OCflags flags = 0; |
---|
697 | |
---|
698 | if(dxd == OCDDS) ext = ".dds"; |
---|
699 | else if(dxd == OCDAS) ext = ".das"; |
---|
700 | else ext = ".dods"; |
---|
701 | |
---|
702 | if(ce != NULL && strlen(ce) == 0) |
---|
703 | ce = NULL; |
---|
704 | |
---|
705 | if(FLAGSET(nccomm->controls,NCF_UNCONSTRAINABLE)) { |
---|
706 | ce = NULL; |
---|
707 | } |
---|
708 | |
---|
709 | if(FLAGSET(nccomm->controls,NCF_ONDISK)) { |
---|
710 | flags |= OCONDISK; |
---|
711 | } |
---|
712 | |
---|
713 | if(SHOWFETCH) { |
---|
714 | /* Build uri string minus the constraint */ |
---|
715 | char* baseurl = nc_uribuild(nccomm->oc.url,NULL,ext,0); |
---|
716 | if(ce == NULL) |
---|
717 | LOG1(NCLOGNOTE,"fetch: %s\n",baseurl); |
---|
718 | else |
---|
719 | LOG2(NCLOGNOTE,"fetch: %s?%s\n",baseurl,ce); |
---|
720 | nullfree(baseurl); |
---|
721 | #ifdef HAVE_GETTIMEOFDAY |
---|
722 | gettimeofday(&time0,NULL); |
---|
723 | #endif |
---|
724 | } |
---|
725 | ocstat = oc_fetchf(conn,ce,dxd,flags,rootp); |
---|
726 | if(FLAGSET(nccomm->controls,NCF_SHOWFETCH)) { |
---|
727 | #ifdef HAVE_GETTIMEOFDAY |
---|
728 | double secs; |
---|
729 | gettimeofday(&time1,NULL); |
---|
730 | secs = deltatime(); |
---|
731 | nclog(NCLOGNOTE,"fetch complete: %0.3f secs",secs); |
---|
732 | #else |
---|
733 | nclog(NCLOGNOTE,"fetch complete."); |
---|
734 | #endif |
---|
735 | } |
---|
736 | #ifdef DEBUG2 |
---|
737 | fprintf(stderr,"fetch: dds:\n"); |
---|
738 | oc_dumpnode(conn,*rootp); |
---|
739 | #endif |
---|
740 | return ocstat; |
---|
741 | } |
---|
742 | |
---|
743 | /* Check a name to see if it contains illegal dap characters |
---|
744 | */ |
---|
745 | |
---|
746 | static char* badchars = "./"; |
---|
747 | |
---|
748 | |
---|
749 | int |
---|
750 | dap_badname(char* name) |
---|
751 | { |
---|
752 | char* p; |
---|
753 | if(name == NULL) return 0; |
---|
754 | for(p=badchars;*p;p++) { |
---|
755 | if(strchr(name,*p) != NULL) return 1; |
---|
756 | } |
---|
757 | return 0; |
---|
758 | } |
---|
759 | |
---|
760 | /* Check a name to see if it contains illegal dap characters |
---|
761 | and repair them |
---|
762 | */ |
---|
763 | |
---|
764 | char* |
---|
765 | dap_repairname(char* name) |
---|
766 | { |
---|
767 | char* newname; |
---|
768 | char *p, *q; int c; |
---|
769 | |
---|
770 | if(name == NULL) return NULL; |
---|
771 | /* assume that dap_badname was called on this name and returned 1 */ |
---|
772 | newname = (char*)malloc(1+(3*strlen(name))); /* max needed */ |
---|
773 | newname[0] = '\0'; /* so we can use strcat */ |
---|
774 | for(p=name,q=newname;(c=*p);p++) { |
---|
775 | if(strchr(badchars,c) != NULL) { |
---|
776 | char newchar[8]; |
---|
777 | snprintf(newchar,sizeof(newchar),"%%%hhx",c); |
---|
778 | strcat(newname,newchar); |
---|
779 | q += 3; /*strlen(newchar)*/ |
---|
780 | } else |
---|
781 | *q++ = c; |
---|
782 | *q = '\0'; /* so we can always do strcat */ |
---|
783 | } |
---|
784 | *q = '\0'; /* ensure trailing null */ |
---|
785 | return newname; |
---|
786 | } |
---|