mexcdf.r2859.zip_matlab 2009b_matlab isvar_matlab r2009 mexcdf_m

/********************************************************************** * * netcdf3.c * * This file contains code to handle the mexfile interface to NetCDF-3 * API calls. * *********************************************************************/ /* * $Id: netcdf3.c 2828 2010-01-15 14:48:55Z johnevans007 $ * */ # include <ctype.h> # include <errno.h> # include <stdio.h> # include <stdlib.h> # include <string.h> # include "netcdf.h" # include "mex.h" # include "mexnc.h" # include "netcdf3.h" /* * Let each function use this over and over again. * */ static char error_message[1000]; /* * DETERMINE_VARM_OUTPUT_SIZE * * The get_varm and put_varm routines are unusual in that the size of * the resulting output array do not necessarily match what that of * the nc_count_coord array. * * The size of each edge depends upon computing the ratio of the two * smallest edges of the imap vector. It has to be that way, if you think * about the definition of the imap vector. For example, suppose that the * netcdf variable is 8x6x4 in size and that the imap vector is [24 4 1]. * This trivially maps the output into an 8x6x4 output array, but it's a * good example, nevertheless. * * There are 192 elements to be retrieved, but the first imap element of * 24 means that there is a distance of * So the output size is [8 6 4]. * * As a second example, suppose we have the same 8x6x4 array, but the * imap vector is [1 8 48]. * * Ex. #2 NetCDF size = [5 4 3 2]. * Want the trivial matlab size of [5 4 3 2]. * Think in terms of the transpose, [2 3 4 5]. * Imap vector is [1 5 20 60] * * Not sure why this works, exactly. It makes sense for the trivial * mapping, but hard to figure otherwise. * * */ void determine_varm_output_size ( int ndims, int num_requested_elements, size_t *nc_count_coord, ptrdiff_t *nc_stride_coord, ptrdiff_t *nc_imap_coord, int *result_size ) { /* * Loop index * */ int j, k; char error_message[1000]; /* * If an element in this array is flagged, it means that we * have already figured its contribution. * */ int still_unused[MAX_NC_DIMS]; /* * Keeps track of the largest remaining imap coordinate for * each dimension. * */ ptrdiff_t max_imap_element_size; ptrdiff_t max_imap_element_index; /* * Initialize the flag array. * */ for ( j = 0; j < MAX_NC_DIMS; ++j ) { still_unused[j] = 1; } /* * For each dimension, figure the contribution. * */ for ( j = 0; j < ndims; ++j ) { /* * Find the largest remaining imap element. * */ max_imap_element_size = -1; for ( k = 0; k < ndims; ++k ) { if ( (nc_imap_coord[k] > max_imap_element_size) && ( still_unused[k] ) ) { max_imap_element_index = k; max_imap_element_size = nc_imap_coord[k]; } } /* * Figure the contribution for this dimension. * Remember to reverse the order, otherwise the row-major order * and column-major order issue will get us.:W * */ result_size[j] = num_requested_elements / nc_imap_coord [ max_imap_element_index ]; /* * Reduce the dimensionality. * */ num_requested_elements = nc_imap_coord [ max_imap_element_index ]; /* * We are done with this imap coordinate. Mark it as used. * */ still_unused[max_imap_element_index] = 0; } /* * Check that no elements in result_size are non-positive. * That can crash matlab hard. * */ for ( j = 0; j < ndims; ++j ) { if ( result_size[j] < 1 ) { sprintf ( error_message, "Requested data extent is invalid.\n\n" ); mexPrintf ( error_message ); mexPrintf ( "count " ); for ( j = 0; j < ndims; ++j ) { printf ( "[%d]", nc_count_coord[j] ); } mexPrintf ( "\n\n" ); mexPrintf ( "stride " ); for ( j = 0; j < ndims; ++j ) { printf ( "[%d]", nc_stride_coord[j] ); } mexPrintf ( "\n\n" ); mexPrintf ( "imap " ); for ( j = 0; j < ndims; ++j ) { printf ( "[%d]", nc_imap_coord[j] ); } mexPrintf ( "\n\n" ); mexPrintf ( "result_size " ); for ( j = 0; j < ndims; ++j ) { printf ( "[%d]", result_size[j] ); } mexPrintf ( "\n\n" ); mexErrMsgTxt ( "\n\n" ); } } } /*********************************************************************** * * HANDLE_NC_ABORT * * code for handling the nc_abort routine. * * PARAMETERS: * nlhs: * number of output arguments defined in matlab * plhs: * Array of output matlab arrays. * plhs[0] = return status of the nc_abort function call * nrhs: * number of input arguments defined in matlab * prhs: * Array of input matlab arrays. * prhs[0] = 'ABORT' * prhs[1] = file ID (ncid) of the netcdf file * nc_op: * Pointer to the current mexnc operation structure. * **********************************************************************/ void handle_nc_abort ( int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[], op *nc_op ) { /* * Pointer shortcut to matrix data. * */ double *pr; /* * file and variable IDs used in nc_copy_att * */ int ncid; /* * Return status from netcdf operation. * */ int status; /* * Make sure that the inputs are the right type. * */ check_numeric_argument_type ( prhs, nc_op->opname, 1 ); pr = mxGetData ( prhs[1] ); ncid = (int)(pr[0]); status = nc_abort(ncid); plhs[0] = mexncCreateDoubleScalar(status); return; } /*********************************************************************** * * HANDLE_NC_CLOSE * * code for handling the nc_close routine. * * PARAMETERS: * nlhs: * number of output arguments defined in matlab * plhs: * Array of output matlab arrays. * plhs[0] = the return status of the nc_close function call * nrhs: * number of input arguments defined in matlab * prhs: * Array of input matlab arrays. * prhs[0] = 'INQ_CLOSE' * prhs[1] = file ID (ncid) of the netcdf file * nc_op: * Pointer to the current mexnc operation structure. * **********************************************************************/ void handle_nc_close ( int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[], op *nc_op ) { /* * Pointer shortcut to matrix data. * */ double *pr; /* * file and variable IDs used in nc_copy_att * */ int ncid; /* * Return status from netcdf operation. * */ int status; /* * Make sure that the inputs are the right type. * */ check_numeric_argument_type ( prhs, nc_op->opname, 1 ); /* * Unpack the parameters. * */ pr = mxGetData ( prhs[1] ); ncid = (int)(pr[0]); status = nc_close(ncid); plhs[0] = mexncCreateDoubleScalar ( status ); return; } /*********************************************************************** * * HANDLE_COPY_ATT: * * code for handling the nc_copy_att routine. * * PARAMETERS: * nlhs: * number of output arguments defined in matlab * plhs: * Array of output matlab arrays. The return status of the * nc_copy_att function call is placed into plhs[0]. * nrhs: * numbe