int main( )
{
/************************* Variable declaration **************************/
intn status_n; /* returned status for functions returning an intn */
int32 status_32, /* returned status for functions returning an int32 */
file_id, vdata_id, vdata_ref,
index = 0; /* index of the vdata in the file - manually kept */
int8 found_fields; /* TRUE if the specified fields exist in the vdata */
/********************** End of variable declaration **********************/
/*
* Open the HDF file for reading.
*/
file_id = Hopen (FILE_NAME, DFACC_READ, 0);
/*
* Initialize the VS interface.
*/
status_n = Vstart (file_id);
/*
* Set the reference number to -1 to start the search from
* the beginning of file.
*/
vdata_ref = -1;
/*
* Assume that the specified fields are not found in the current vdata.
*/
found_fields = FALSE;
/*
* Use VSgetid to obtain each vdata by its reference number then
* attach to the vdata and search for the fields. The loop
* terminates when the last vdata is reached or when a vdata which
* contains the fields listed in SEARCHED_FIELDS is found.
*/
while ((vdata_ref = VSgetid (file_id, vdata_ref)) != FAIL)
{
vdata_id = VSattach (file_id, vdata_ref, "r");
if ((status_n = VSfexist (vdata_id, SEARCHED_FIELDS)) != FAIL)
{
found_fields = TRUE;
break;
}
/*
* Detach from the current vdata before continuing searching.
*/
status_32 = VSdetach (vdata_id);
index++; /* advance the index by 1 for the next vdata */
}
/*
* Print the index of the vdata containing the fields or a "not found"
* message if no such vdata is found. Also detach from the vdata found.
*/
if (!found_fields)
printf ("Fields Position and Temperature were not found.\n");
else
{
printf
("Fields Position and Temperature found in the vdata at position %d\n",
index);
status_32 = VSdetach (vdata_id);
}
/*
* Terminate access to the VS interface and close the HDF file.
*/
status_n = Vend (file_id);
status_32 = Hclose (file_id);
return 0;
}
int main(void)
{
int32 sd1_id; /* SD interface identifier */
int32 sd2_id; /* SD interface identifier */
int32 sds1_id; /* SDS identifier */
int32 sds2_id; /* SDS identifier */
int32 dim_sizes[2]; /* sizes of the SDS dimensions */
int32 start[2]; /* start location to write */
int32 edges[2]; /* number of elements to write */
int32 n_values;
int32 buf1a[Y_LENGTH][X_LENGTH] = {{1,1},{1,1},{5,6}};
int32 buf1b[Y_LENGTH][X_LENGTH] = {{1,2},{3,4},{5,6}};
/* percent test: compare divide by zero, both zero */
int32 buf2a[Y_LENGTH][X_LENGTH] = {{100,100},{100,0},{0,100}};
int32 buf2b[Y_LENGTH][X_LENGTH] = {{120,80},{0,100},{0,50}};
/* global attributes */
char8 bufga1[] = "Storm_track_data1";
char8 bufga2[] = "Storm_track_data2";
float32 bufa1[2] = {1., 1.};
float32 bufa2[2] = {1., 2.};
/*vdata*/
int32 file1_id;
int32 file2_id;
int32 vdata1_id;
int32 vdata2_id;
char8 vdata1_buf1 [N_RECORDS_1] = {'V', 'D', 'A', 'T', 'A'};
char8 vdata1_buf2 [N_RECORDS_1] = {'X', 'D', 'A', 'T', 'A'};
int32 vdata2_buf1 [N_RECORDS_2][ORDER_2] = {{1, 2, 3, 4}, {5, 6, 7, 8}};
int32 vdata2_buf2 [N_RECORDS_2][ORDER_2] = {{1, 1, 1, 1}, {5, 6, 7, 8}};
float32 vdata3_buf1[N_RECORDS_2][N_VALS_PER_REC]={{1,2,3,4,5,6},{7,8,9,10,11,12}};
float32 vdata3_buf2[N_RECORDS_2][N_VALS_PER_REC]={{1,1,1,1,1,1},{7,8,9,10,11,12}};
/* Define the location and size of the data to be written to the data set*/
start[0] = 0;
start[1] = 0;
edges[0] = Y_LENGTH;
edges[1] = X_LENGTH;
/* Define the dimensions of the array to be created */
dim_sizes[0] = Y_LENGTH;
dim_sizes[1] = X_LENGTH;
/*-------------------------------------------------------------------------
* SD data
*-------------------------------------------------------------------------
*/
/* Create the files and initialize the SD interface */
if ((sd1_id = SDstart (FILE1_NAME, DFACC_CREATE))==FAIL)
goto error;
if ((sd2_id = SDstart (FILE2_NAME, DFACC_CREATE))==FAIL)
goto error;
/* Set a global attribute */
n_values = sizeof(bufga1);
if ( SDsetattr (sd1_id, FILE_ATTR_NAME, DFNT_CHAR8, n_values, (VOIDP)bufga1)==FAIL)
goto error;
if ( SDsetattr (sd2_id, FILE_ATTR_NAME, DFNT_CHAR8, n_values, (VOIDP)bufga2)==FAIL)
goto error;
/* Create the data sets */
if ((sds1_id = SDcreate (sd1_id, "dset1", DFNT_INT32, RANK, dim_sizes))==FAIL)
goto error;
if ((sds2_id = SDcreate (sd2_id, "dset1", DFNT_INT32, RANK, dim_sizes))==FAIL)
goto error;
/* Assign attribute */
n_values = 2;
if ( SDsetattr (sds1_id, SDS_ATTR_NAME, DFNT_FLOAT32, n_values, (VOIDP)bufa1)==FAIL)
goto error;
if ( SDsetattr (sds2_id, SDS_ATTR_NAME, DFNT_FLOAT32, n_values, (VOIDP)bufa2)==FAIL)
goto error;
/* Write the stored data to the data set */
if ( SDwritedata (sds1_id, start, NULL, edges, (VOIDP)buf1a)==FAIL)
goto error;
if ( SDwritedata (sds2_id, start, NULL, edges, (VOIDP)buf1b)==FAIL)
goto error;
/* Terminate access to the data set */
if ( SDendaccess (sds1_id)==FAIL)
goto error;
if ( SDendaccess (sds2_id)==FAIL)
goto error;
/* Create another data set */
if (( sds1_id = SDcreate (sd1_id, "dset2", DFNT_INT32, RANK, dim_sizes))==FAIL)
goto error;
if (( sds2_id = SDcreate (sd2_id, "dset2", DFNT_INT32, RANK, dim_sizes))==FAIL)
goto error;
if ( SDwritedata (sds1_id, start, NULL, edges, (VOIDP)buf1a)==FAIL)
goto error;
if ( SDwritedata (sds2_id, start, NULL, edges, (VOIDP)buf1b)==FAIL)
goto error;
if ( SDendaccess (sds1_id)==FAIL)
goto error;
if ( SDendaccess (sds2_id)==FAIL)
goto error;
/* data sets for -p test */
if (( sds1_id = SDcreate (sd1_id, "dset3", DFNT_INT32, RANK, dim_sizes))==FAIL)
goto error;
if (( sds2_id = SDcreate (sd2_id, "dset3", DFNT_INT32, RANK, dim_sizes))==FAIL)
goto error;
if ( SDwritedata (sds1_id, start, NULL, edges, (VOIDP)buf2a)==FAIL)
goto error;
if ( SDwritedata (sds2_id, start, NULL, edges, (VOIDP)buf2b)==FAIL)
goto error;
if ( SDendaccess (sds1_id)==FAIL)
goto error;
if ( SDendaccess (sds2_id)==FAIL)
goto error;
/*-------------------------------------------------------------------------
* end SD
*-------------------------------------------------------------------------
*/
/* Terminate access to the SD interface and close the file */
if ( SDend (sd1_id)==FAIL)
goto error;
if ( SDend (sd2_id)==FAIL)
goto error;
/*-------------------------------------------------------------------------
* VD data
*-------------------------------------------------------------------------
*/
/* Open the HDF file for writing */
if ((file1_id = Hopen (FILE1_NAME, DFACC_WRITE, 0))==FAIL)
goto error;
if ((file2_id = Hopen (FILE2_NAME, DFACC_WRITE, 0))==FAIL)
goto error;
/* Initialize the VS interface */
if ( Vstart (file1_id)==FAIL)
goto error;
if ( Vstart (file2_id)==FAIL)
goto error;
/*-------------------------------------------------------------------------
* VD data one field
*-------------------------------------------------------------------------
*/
/* Create the first vdata and populate it with data from vdata1_buf */
if (VHstoredata (file1_id, FIELD1_NAME, (uint8 *)vdata1_buf1,
N_RECORDS_1, DFNT_CHAR8, VDATA1_NAME, CLASS1_NAME)==FAIL)
goto error;
if (VHstoredata (file2_id, FIELD1_NAME, (uint8 *)vdata1_buf2,
N_RECORDS_1, DFNT_CHAR8, VDATA1_NAME, CLASS1_NAME)==FAIL)
goto error;
/*-------------------------------------------------------------------------
* VD data one field, order 4
*-------------------------------------------------------------------------
*/
/* Create the second vdata and populate it with data from vdata2_buf */
if ( VHstoredatam (file1_id, FIELD2_NAME, (uint8 *)vdata2_buf1,
N_RECORDS_2, DFNT_INT32, VDATA2_NAME, CLASS2_NAME, ORDER_2)==FAIL)
goto error;
if ( VHstoredatam (file2_id, FIELD2_NAME, (uint8 *)vdata2_buf2,
N_RECORDS_2, DFNT_INT32, VDATA2_NAME, CLASS2_NAME, ORDER_2)==FAIL)
goto error;
/*-------------------------------------------------------------------------
* VD data several fields
*-------------------------------------------------------------------------
*/
/* Create a new vdata */
if ((vdata1_id = VSattach (file1_id, -1, "w"))==FAIL)
goto error;
if ((vdata2_id = VSattach (file2_id, -1, "w"))==FAIL)
goto error;
/* Set name and class name of the vdata */
if ( VSsetname (vdata1_id, VDATA3_NAME)==FAIL)
goto error;
if ( VSsetclass (vdata1_id, CLASS3_NAME)==FAIL)
goto error;
if ( VSsetname (vdata2_id, VDATA3_NAME)==FAIL)
goto error;
if ( VSsetclass (vdata2_id, CLASS3_NAME)==FAIL)
goto error;
/* Define fields */
if ( VSfdefine (vdata1_id, FIELD3_NAME1, DFNT_FLOAT32, ORDER3_1 )==FAIL)
goto error;
if ( VSfdefine (vdata1_id, FIELD3_NAME2, DFNT_FLOAT32, ORDER3_2 )==FAIL)
goto error;
if ( VSfdefine (vdata1_id, FIELD3_NAME3, DFNT_FLOAT32, ORDER3_3 )==FAIL)
goto error;
if ( VSsetfields (vdata1_id, FIELDNAME3_LIST )==FAIL)
goto error;
if ( VSfdefine (vdata2_id, FIELD3_NAME1, DFNT_FLOAT32, ORDER3_1 )==FAIL)
goto error;
if ( VSfdefine (vdata2_id, FIELD3_NAME2, DFNT_FLOAT32, ORDER3_2 )==FAIL)
goto error;
if ( VSfdefine (vdata2_id, FIELD3_NAME3, DFNT_FLOAT32, ORDER3_3 )==FAIL)
goto error;
if ( VSsetfields (vdata2_id, FIELDNAME3_LIST)==FAIL)
goto error;
/* Write the data with full interlacing mode */
if ( VSwrite (vdata1_id, (uint8 *)vdata3_buf1, N_RECORDS_2, FULL_INTERLACE)==FAIL)
goto error;
if ( VSwrite (vdata2_id, (uint8 *)vdata3_buf2, N_RECORDS_2, FULL_INTERLACE)==FAIL)
goto error;
if ( VSdetach (vdata1_id)==FAIL)
goto error;
if ( VSdetach (vdata2_id)==FAIL)
goto error;
/*-------------------------------------------------------------------------
* end VD data
*-------------------------------------------------------------------------
*/
/* Terminate access to the VS interface and close the HDF file */
if ( Vend (file1_id)==FAIL)
goto error;
if ( Vend (file2_id)==FAIL)
goto error;
if ( Hclose (file1_id)==FAIL)
goto error;
if ( Hclose (file2_id)==FAIL)
goto error;
/*-------------------------------------------------------------------------
* write 2 big files for hyperslab reading
*-------------------------------------------------------------------------
*/
if (do_big_file()==FAIL)
goto error;
/*-------------------------------------------------------------------------
* groups
*-------------------------------------------------------------------------
*/
if (do_groups()==FAIL)
goto error;
return 0;
error:
printf("Error...Exiting...\n");
return 1;
}
void main(void)
{
char vdata_name[MAX_NC_NAME], vdata_class[MAX_NC_NAME];
char fields[60];
int32 file_id, vdata_id, istat;
int32 n_records, interlace, vdata_size, vdata_ref;
int bufsz = (2 * sizeof(float32) + sizeof(char) \
+ sizeof(int16)) * NRECORDS;
int i;
uint8 databuf[((2 * sizeof(float32)) + sizeof(char) \
+ sizeof(int16)) * NRECORDS];
VOIDP fldbufptrs[4];
float32 itemp[NRECORDS], ispeed[NRECORDS];
int16 iheight[NRECORDS];
char idents[NRECORDS];
/* Open the HDF file. */
file_id = Hopen("VD_Ex4.hdf", DFACC_READ, 0);
/* Initialize the Vset interface. */
istat = Vstart(file_id);
/*
* Get the reference number for the first Vdata in
* the file.
*/
vdata_ref = -1;
vdata_ref = VSgetid(file_id, vdata_ref);
/* Attach to the first Vdata in read mode. */
vdata_id = VSattach(file_id, vdata_ref, "r");
for (i=0; i<60; i++)
fields[i] = '\0';
/* Get the list of field names. */
istat =VSinquire(vdata_id, &n_records, &interlace,
fields, &vdata_size, vdata_name);
printf("files: %s, n_records: %d, vdata_size: %d\n",
fields, n_records, vdata_size);
/* Get the class. */
istat = VSgetclass(vdata_id, vdata_class);
/* Determine the fields that will be read. */
istat = VSsetfields(vdata_id, fields);
/* Print the Vdata information. */
printf("Current Vdata name: %s \nCurrent Vdata class: %s.\n",
vdata_name, vdata_class);
/* Read the data. */
istat = VSread(vdata_id, (VOIDP)databuf, n_records,
FULL_INTERLACE);
/* set fldbufptrs and unpack field values */
fldbufptrs[0] = &idents[0];
fldbufptrs[1] = &ispeed[0];
fldbufptrs[2] = &iheight[0];
fldbufptrs[3] = &itemp[0];
istat = VSfpack(vdata_id, _HDF_VSUNPACK, fields, databuf,
bufsz, n_records, "Ident,Speed,Height,Temp",
fldbufptrs);
printf(" Temp Height Speed Ident\n");
for (i=0; i < n_records; i++) {
printf(" %6.2f %6d %6.2f %c\n",
itemp[i],iheight[i],ispeed[i],idents[i]);
}
/* Detach from the Vdata, close the interface and the file. */
istat = VSdetach(vdata_id);
istat = Vend(file_id);
istat = Hclose(file_id);
}
uint32 diff_vs( int32 file1_id,
int32 file2_id,
int32 ref1,
int32 ref2,
diff_opt_t * opt)
{
int32 vdata1_id, /* vdata identifier */
n_records1, /* number of records */
vdata1_size,
interlace1_mode,
vdata2_id=-1, /* vdata identifier */
n_records2, /* number of records */
vdata2_size,
interlace2_mode;
char vdata1_name [VSNAMELENMAX];
char vdata1_class[VSNAMELENMAX];
char fieldname1_list[VSFIELDMAX*FIELDNAMELENMAX];
char vdata2_name [VSNAMELENMAX];
char vdata2_class[VSNAMELENMAX];
char fieldname2_list[VSFIELDMAX*FIELDNAMELENMAX];
uint32 nfound=0;
/*-------------------------------------------------------------------------
* object 1
*-------------------------------------------------------------------------
*/
if (Vstart (file1_id)==FAIL) {
printf("Error: Could not start VS interface in VS ref %d\n", ref1);
goto out;
}
if ((vdata1_id = VSattach (file1_id, ref1, "r")) == FAIL ){
printf( "Failed to attach VS ref %d\n", ref1);
goto out;
}
if (VSgetname (vdata1_id, vdata1_name) == FAIL ){
printf( "Failed to name for VS ref %d\n", ref1);
goto out;
}
if (VSgetclass (vdata1_id, vdata1_class) == FAIL ){
printf( "Failed to name for VS ref %d\n", ref1);
goto out;
}
if (VSinquire(vdata1_id, &n_records1, &interlace1_mode, fieldname1_list,
&vdata1_size, vdata1_name) == FAIL) {
printf( "Failed to get info for VS ref %d\n", ref1);
goto out;
}
if (VFnfields(vdata1_id)== FAIL ){
printf( "Failed getting fields forVS ref %d\n", ref1);
goto out;
}
/*-------------------------------------------------------------------------
* object 2
*-------------------------------------------------------------------------
*/
if (Vstart (file2_id)==FAIL) {
printf("Error: Could not start VS interface in VS ref %d\n", ref1);
goto out;
}
if ((vdata2_id = VSattach (file2_id, ref2, "r")) == FAIL ){
printf("Failed to attach VS ref %d\n", ref2);
goto out;
}
if (VSgetname (vdata2_id, vdata2_name) == FAIL ){
printf("Failed to name for VS ref %d\n", ref2);
goto out;
}
if (VSgetclass (vdata2_id, vdata2_class) == FAIL ){
printf("Failed to name for VS ref %d\n", ref2);
goto out;
}
if (VSinquire(vdata2_id, &n_records2, &interlace2_mode, fieldname2_list,
&vdata2_size, vdata2_name) == FAIL) {
printf("Failed to get info for VS ref %d\n", ref2);
goto out;
}
if (VFnfields(vdata2_id)== FAIL ){
printf("Failed getting fields forVS ref %d\n", ref2);
goto out;
}
/*-------------------------------------------------------------------------
* check for input VSs
*-------------------------------------------------------------------------
*/
if (opt->nuvars > 0) /* if specified vdata is selected */
{
int imatch = 0, j;
for (j = 0; j < opt->nuvars; j++)
{
if (strcmp(vdata1_name, opt->uvars[j]) == 0)
{
imatch = 1;
break;
}
}
if (imatch == 0)
{
goto do_nothing;
}
}
/*-------------------------------------------------------------------------
* Comparing
*-------------------------------------------------------------------------
*/
if (opt->verbose)
printf("Comparing <%s>\n",vdata1_name);
nfound=vdata_cmp(vdata1_id,vdata2_id,vdata1_name,vdata1_class,opt);
do_nothing:
/* terminate access to the VSs */
if (VSdetach (vdata1_id)==FAIL) {
printf("Failed to dettach VS ref %d\n", ref1);
goto out;
}
if (vdata2_id!=-1)
{
if (VSdetach (vdata2_id)==FAIL) {
printf("Failed to dettach VS ref %d\n", ref2);
goto out;
}
}
return nfound;
out:
opt->err_stat = 1;
return 0;
}
int main( )
{
/************************* Variable declaration **************************/
intn status_n; /* returned status for functions returning an intn */
int32 status_32, /* returned status for functions returning an int32 */
file_id, vdata_ref, vdata_id,
field_index, /* index of a field within the vdata */
n_vdattrs, /* number of vdata attributes */
n_fldattrs, /* number of field attributes */
vdata_type, /* to hold the type of vdata's attribute */
vdata_n_values,/* to hold the number of vdata's attribute values */
vdata_size, /* to hold the size of vdata's attribute values */
field_type, /* to hold the type of field's attribute */
field_n_values,/* to hold the number of field's attribute values */
field_size; /* to hold the size of field's attribute values */
char vd_attr[VATTR_N_VALUES] = {'A', 'B', 'C'};/* vdata attribute values*/
int32 fld_attr[FATTR_N_VALUES] = {2, 4, 6, 8}; /* field attribute values*/
char vattr_buf[VATTR_N_VALUES]; /* to hold vdata attribute's values */
int32 fattr_buf[FATTR_N_VALUES]; /* to hold field attribute's values */
char vattr_name[30], /* name of vdata attribute */
fattr_name[30]; /* name of field attribute */
/********************** End of variable declaration **********************/
/*
* Open the HDF file for writing.
*/
file_id = Hopen (FILE_NAME, DFACC_WRITE, 0);
/*
* Initialize the VS interface.
*/
status_n = Vstart (file_id);
/*
* Get the reference number of the vdata named VDATA_NAME.
*/
vdata_ref = VSfind (file_id, VDATA_NAME);
/*
* Attach to the vdata for writing.
*/
vdata_id = VSattach (file_id, vdata_ref, "w");
/*
* Attach an attribute to the vdata, i.e., indicated by the second parameter.
*/
status_n = VSsetattr (vdata_id, _HDF_VDATA, VATTR_NAME, DFNT_CHAR,
VATTR_N_VALUES, vd_attr);
/*
* Get the index of the field FIELD_NAME within the vdata.
*/
status_n = VSfindex (vdata_id, FIELD_NAME, &field_index);
/*
* Attach an attribute to the field field_index.
*/
status_n = VSsetattr (vdata_id, field_index, FATTR_NAME, DFNT_INT32,
FATTR_N_VALUES, fld_attr);
/*
* Get the number of attributes attached to the vdata's first
* field - should be 0.
*/
n_fldattrs = VSfnattrs (vdata_id, 0);
printf ( "Number of attributes of the first field of the vdata: %d\n",
n_fldattrs);
/*
* Get the number of attributes attached to the field specified by
* field_index - should be 1.
*/
n_fldattrs = VSfnattrs (vdata_id, field_index);
printf ( "Number of attributes of field %s: %d\n", FIELD_NAME, n_fldattrs);
/*
* Get the total number of the field's and vdata's attributes - should be 2.
*/
n_vdattrs = VSnattrs (vdata_id);
printf ( "Number of attributes of the vdata and its fields: %d\n",
n_vdattrs);
/*
* Get information about the vdata's first attribute, indicated
* by the third parameter which is the index of the attribute.
*/
status_n = VSattrinfo (vdata_id, _HDF_VDATA, 0, vattr_name,
&vdata_type, &vdata_n_values, &vdata_size);
/*
* Get information about the first attribute of the field specified by
* field_index.
*/
status_n = VSattrinfo (vdata_id, field_index, 0, fattr_name, &field_type,
&field_n_values, &field_size);
/*
* Get the vdata's first attribute.
*/
status_n = VSgetattr (vdata_id, _HDF_VDATA, 0, vattr_buf);
printf("Values of the vdata attribute = %c %c %c\n", vattr_buf[0],
vattr_buf[1], vattr_buf[2]);
/*
* Get the first attribute of the field specified by field_index.
*/
status_n = VSgetattr (vdata_id, field_index, 0, fattr_buf);
printf("Values of the field attribute = %d %d %d %d\n", fattr_buf[0],
fattr_buf[1], fattr_buf[2], fattr_buf[3]);
/*
* Terminate access to the vdata and to the VS interface, then close
* the HDF file.
*/
status_32 = VSdetach (vdata_id);
status_n = Vend (file_id);
status_32 = Hclose (file_id);
return 0;
}
//// [[Rcpp::export]]
int binlist(int x) {
int numbins = 0, *binnums = NULL;
int i;
int file_id,vdata_ref, vdata_id,numrecs,pvd_id;
char PARAM[] = "chlor_a";
initbin();
/* Open the HDF file. */
file_id = Hopen("S1998001.L3b_DAY_CHL.main", DFACC_READ, 0);
if(file_id == FAIL){
fprintf(stderr,"-E- %s line %d: Hopen(%s,DFACC_READ,0) failed.\n",
__FILE__,__LINE__,"S1998001.L3b_DAY_CHL.main");
return(EXIT_FAILURE);
}
/* Initialize the Vdata interface. */
if(Vstart(file_id) == FAIL){
fprintf(stderr,"-E- %s line %d: Vstart(%d) failed.\n",
__FILE__,__LINE__,file_id);
return(EXIT_FAILURE);
}
/* Open the "BinList" Vdata. */
vdata_ref = VSfind(file_id,"BinList");
if(vdata_ref == FAIL){
fprintf(stderr,"-E- %s line %d: VSfind(%d,\"BinList\") failed.\n",
__FILE__,__LINE__,file_id);
return(EXIT_FAILURE);
}
vdata_id = VSattach(file_id, vdata_ref, "r");
if(vdata_id == FAIL){
fprintf(stderr,"-E- %s line %d: VSattach(%d,%d,\"r\") failed.\n",
__FILE__,__LINE__,file_id,vdata_ref);
return(EXIT_FAILURE);
}
/* Find out how many bins are stored in this file. */
numrecs = VSelts(vdata_id);
if(numrecs == FAIL){
fprintf(stderr,"-E- %s line %d: VSelts(%d) failed.\n",
__FILE__,__LINE__,vdata_id);
return(EXIT_FAILURE);
}
// this doesn't seem to be necessary?
/* Set up to read the fields in the BinList Vdata records. */
// if(VSsetfields(vdata_id,BLIST_FIELDS) == FAIL){
// fprintf(stderr,"-E- %s line %d: VSsetfields(%d,%s) failed.\n",
// __FILE__,__LINE__,vdata_id,BLIST_FIELDS);
// return(EXIT_FAILURE);
// }
/* Open the parameter-specific Vdata. */
vdata_ref = VSfind(file_id,PARAM);
if(vdata_ref == 0){
fprintf(stderr,"-E- %s line %d: VSfind(%d,\"%s\") failed.\n",
__FILE__,__LINE__,file_id,PARAM);
return(EXIT_FAILURE);
}
pvd_id = VSattach(file_id, vdata_ref, "r");
if(pvd_id == FAIL){
fprintf(stderr,"-E- %s line %d: VSattach(%d,%d,\"r\") failed.\n",
__FILE__,__LINE__,file_id,vdata_ref);
return(EXIT_FAILURE);
}
/* Set up to read the fields in the parameter-specific Vdata records. */
{
int len;
len = 2*strlen(PARAM) + strlen("_sum,") + strlen("_sum_sq") + 1;
param_fields = (char *)malloc(len);
if(param_fields == NULL){
fprintf(stderr,"-E- %s line %d: Memory allocation failed.\n",
__FILE__,__LINE__);
return(EXIT_FAILURE);
}
strcpy(param_fields,PARAM);
strcat(param_fields,"_sum,");
strcat(param_fields,PARAM);
strcat(param_fields,"_sum_sq");
}
if(VSsetfields(pvd_id,param_fields) == FAIL){
fprintf(stderr,"-E- %s line %d: VSsetfields(%d,%s) failed.\n",
__FILE__,__LINE__,pvd_id,param_fields);
return(EXIT_FAILURE);
}
printf(param_fields);
printf("\n");
/* Output a header record to identify the fields written out below. */
printf("%80s%15.15s %15.15s\n"," ",PARAM,PARAM);
printf(" bin centerlat centerlon");
printf(" north south west east");
printf(" n N sum_obs sum_squared_obs weight");
printf(" time_trend_bits l2_flag_bits sel\n");
printf("------- --------- ----------");
printf(" --------- --------- ---------- ----------");
printf(" ---- --- --------------- --------------- ---------------");
printf(" ---------------- -------------------------------- ---\n");
// for(i = 0; i < numbins; i++){
for(i = 0; i < 10; i++){
int recno;
// recno = binsearch(binnums[i],vdata_id,numrecs);
recno = i;
if(recno >= 0){
double n,s,w,e,clat,clon;
/*
Read the sum and sum-of-squares for the
the specified parameter for this bin.
*/
if(VSseek(pvd_id,recno) == FAIL){
fprintf(stderr,"-E- %s line %d: VSseek(%d,%d) failed.\n",
__FILE__,__LINE__,pvd_id,recno);
return(EXIT_FAILURE);
}
if(VSread(pvd_id,paramrec,1,FULL_INTERLACE) != 1){
fprintf(stderr,"-E- %s line %d: ",__FILE__,__LINE__);
fprintf(stderr,"VSread(%d,paramrec,1,FULL_INTERLACE) failed.\n",
pvd_id);
return(EXIT_FAILURE);
}
/*
VSfpack() sets the global sum and sum_sq variables
via the paramptrs pointer array.
*/
if(
VSfpack(
pvd_id,_HDF_VSUNPACK,param_fields,paramrec,PREC_SIZE,1,NULL,paramptrs
)
== FAIL){
fprintf(stderr,"-E- %s line %d: ",__FILE__,__LINE__);
fprintf(stderr,"VSfpack(%d, ...) failed.\n", pvd_id);
return(EXIT_FAILURE);
}
/* Get the geographical coordinates associated with this bin. */
// bin2latlon(binnums[i],&clat,&clon);
// bin2bounds(binnums[i],&n,&s,&w,&e);
/* Output the results. */
// printf("%7d %9.5f %10.5f %9.5f %9.5f %10.5f %10.5f ",
// binnums[i],clat,clon,n,s,w,e);
// printf("%4d %3d ",nobs,nscenes);
// printf("% .8e % .8e % .8e ",sum,sum_sq,weights);
printf("% .8e % .8e % .8e ",sum,sum_sq,0.0);
// printf("%.16s %.32s ",bitstr16(time_rec),bitstr32(flags_set));
// printf("%3d",sel_cat);
printf("\n");
}
}
if(VSdetach(pvd_id) == FAIL){
fprintf(stderr,"-E- %s line %d: VSdetach(%d) failed.\n",
__FILE__,__LINE__,pvd_id);
return(EXIT_FAILURE);
}
if(VSdetach(vdata_id) == FAIL){
fprintf(stderr,"-E- %s line %d: VSdetach(%d) failed.\n",
__FILE__,__LINE__,vdata_id);
return(EXIT_FAILURE);
}
if(Vend(file_id) == FAIL){
fprintf(stderr,"-E- %s line %d: Vend(%d) failed.\n",
__FILE__,__LINE__,file_id);
return(EXIT_FAILURE);
}
if(Hclose(file_id) == FAIL){
fprintf(stderr,"-E- %s line %d: Hclose(%d) failed.\n",
__FILE__,__LINE__,file_id);
return(EXIT_FAILURE);
}
free(param_fields);
free(binnums);
return 2;
}
