The NetCDF Fortran 90 Interface Guide ************************************* This document describes the Fortran 90 interface to the netCDF library. It applies to netCDF version 4.0.1-beta1. This document was last updated in 2 October 2008. For a complete description of the netCDF format and utilities see *note The NetCDF Users Guide: (netcdf)Top. 1 Use of the NetCDF Library *************************** You can use the netCDF library without knowing about all of the netCDF interface. If you are creating a netCDF dataset, only a handful of routines are required to define the necessary dimensions, variables, and attributes, and to write the data to the netCDF dataset. (Even less are needed if you use the ncgen utility to create the dataset before running a program using netCDF library calls to write data. *Note ncgen: (netcdf)ncgen.) Similarly, if you are writing software to access data stored in a particular netCDF object, only a small subset of the netCDF library is required to open the netCDF dataset and access the data. Authors of generic applications that access arbitrary netCDF datasets need to be familiar with more of the netCDF library. In this chapter we provide templates of common sequences of netCDF calls needed for common uses. For clarity we present only the names of routines; omit declarations and error checking; omit the type-specific suffixes of routine names for variables and attributes; indent statements that are typically invoked multiple times; and use ... to represent arbitrary sequences of other statements. Full parameter lists are described in later chapters. 1.1 Creating a NetCDF Dataset ============================= Here is a typical sequence of netCDF calls used to create a new netCDF dataset: NF90_CREATE ! create netCDF dataset: enter define mode ... NF90_DEF_DIM ! define dimensions: from name and length ... NF90_DEF_VAR ! define variables: from name, type, dims ... NF90_PUT_ATT ! assign attribute values ... NF90_ENDDEF ! end definitions: leave define mode ... NF90_PUT_VAR ! provide values for variable ... NF90_CLOSE ! close: save new netCDF dataset Only one call is needed to create a netCDF dataset, at which point you will be in the first of two netCDF modes. When accessing an open netCDF dataset, it is either in define mode or data mode. In define mode, you can create dimensions, variables, and new attributes, but you cannot read or write variable data. In data mode, you can access data and change existing attributes, but you are not permitted to create new dimensions, variables, or attributes. One call to NF90_DEF_DIM is needed for each dimension created. Similarly, one call to NF90_DEF_VAR is needed for each variable creation, and one call to a member of the NF90_PUT_ATT family is needed for each attribute defined and assigned a value. To leave define mode and enter data mode, call NF90_ENDDEF. Once in data mode, you can add new data to variables, change old values, and change values of existing attributes (so long as the attribute changes do not require more storage space). Data of all types is written to a netCDF variable using the NF90_PUT_VAR subroutine. Single values, arrays, or array sections may be supplied to NF90_PUT_VAR; optional arguments allow the writing of subsampled or mapped portions of the variable. (Subsampled and mapped access are general forms of data access that are explained later.) Finally, you should explicitly close all netCDF datasets that have been opened for writing by calling NF90_CLOSE. By default, access to the file system is buffered by the netCDF library. If a program terminates abnormally with netCDF datasets open for writing, your most recent modifications may be lost. This default buffering of data is disabled by setting the NF90_SHARE flag when opening the dataset. But even if this flag is set, changes to attribute values or changes made in define mode are not written out until NF90_SYNC or NF90_CLOSE is called. 1.2 Reading a NetCDF Dataset with Known Names ============================================= Here we consider the case where you know the names of not only the netCDF datasets, but also the names of their dimensions, variables, and attributes. (Otherwise you would have to do "inquire" calls.) The order of typical C calls to read data from those variables in a netCDF dataset is: NF90_OPEN ! open existing netCDF dataset ... NF90_INQ_DIMID ! get dimension IDs ... NF90_INQ_VARID ! get variable IDs ... NF90_GET_ATT ! get attribute values ... NF90_GET_VAR ! get values of variables ... NF90_CLOSE ! close netCDF dataset First, a single call opens the netCDF dataset, given the dataset name, and returns a netCDF ID that is used to refer to the open netCDF dataset in all subsequent calls. Next, a call to NF90_INQ_DIMID for each dimension of interest gets the dimension ID from the dimension name. Similarly, each required variable ID is determined from its name by a call to NF90_INQ_VARID. Once variable IDs are known, variable attribute values can be retrieved using the netCDF ID, the variable ID, and the desired attribute name as input to NF90_GET_ATT for each desired attribute. Variable data values can be directly accessed from the netCDF dataset with calls to NF90_GET_VAR. Finally, the netCDF dataset is closed with NF90_CLOSE. There is no need to close a dataset open only for reading. 1.3 Reading a netCDF Dataset with Unknown Names =============================================== It is possible to write programs (e.g., generic software) which do such things as processing every variable, without needing to know in advance the names of these variables. Similarly, the names of dimensions and attributes may be unknown. Names and other information about netCDF objects may be obtained from netCDF datasets by calling inquire functions. These return information about a whole netCDF dataset, a dimension, a variable, or an attribute. The following template illustrates how they are used: NF90_OPEN ! open existing netCDF dataset ... NF90_INQUIRE ! find out what is in it ... NF90_INQUIRE_DIMENSION ! get dimension names, lengths ... NF90_INQUIRE_VARIABLE ! get variable names, types, shapes ... NF90_INQ_ATTNAME ! get attribute names ... NF90_INQUIRE_ATTRIBUTE ! get other attribute information ... NF90_GET_ATT ! get attribute values ... NF90_GET_VAR ! get values of variables ... NF90_CLOSE ! close netCDF dataset As in the previous example, a single call opens the existing netCDF dataset, returning a netCDF ID. This netCDF ID is given to the NF90_INQUIRE routine, which returns the number of dimensions, the number of variables, the number of global attributes, and the ID of the unlimited dimension, if there is one. All the inquire functions are inexpensive to use and require no I/O, since the information they provide is stored in memory when a netCDF dataset is first opened. Dimension IDs use consecutive integers, beginning at 1. Also dimensions, once created, cannot be deleted. Therefore, knowing the number of dimension IDs in a netCDF dataset means knowing all the dimension IDs: they are the integers 1, 2, 3, ...up to the number of dimensions. For each dimension ID, a call to the inquire function NF90_INQUIRE_DIMENSION returns the dimension name and length. Variable IDs are also assigned from consecutive integers 1, 2, 3, ... up to the number of variables. These can be used in NF90_INQUIRE_VARIABLE calls to find out the names, types, shapes, and the number of attributes assigned to each variable. Once the number of attributes for a variable is known, successive calls to NF90_INQ_ATTNAME return the name for each attribute given the netCDF ID, variable ID, and attribute number. Armed with the attribute name, a call to NF90_INQUIRE_ATTRIBUTE returns its type and length. Given the type and length, you can allocate enough space to hold the attribute values. Then a call to NF90_GET_ATT returns the attribute values. Once the IDs and shapes of netCDF variables are known, data values can be accessed by calling NF90_GET_VAR. 1.4 Writing Data in an Existing NetCDF Dataset ============================================== With write access to an existing netCDF dataset, you can overwrite data values in existing variables or append more data to record variables along the unlimited (record) dimension. To append more data to non-record variables requires changing the shape of such variables, which means creating a new netCDF dataset, defining new variables with the desired shape, and copying data. The netCDF data model was not designed to make such "schema changes" efficient or easy, so it is best to specify the shapes of variables correctly when you create a netCDF dataset, and to anticipate which variables will later grow by using the unlimited dimension in their definition. The following code template lists a typical sequence of calls to overwrite some existing values and add some new records to record variables in an existing netCDF dataset with known variable names: NF90_OPEN ! open existing netCDF dataset ... NF90_INQ_VARID ! get variable IDs ... NF90_PUT_VAR ! provide new values for variables, if any ... NF90_PUT_ATT ! provide new values for attributes, if any ... NF90_CLOSE ! close netCDF dataset A netCDF dataset is first opened by the NF90_OPEN call. This call puts the open dataset in data mode, which means existing data values can be accessed and changed, existing attributes can be changed, but no new dimensions, variables, or attributes can be added. Next, calls to NF90_INQ_VARID get the variable ID from the name, for each variable you want to write. Then each call to NF90_PUT_VAR writes data into a specified variable, either a single value at a time, or a whole set of values at a time, depending on which variant of the interface is used. The calls used to overwrite values of non-record variables are the same as are used to overwrite values of record variables or append new data to record variables. The difference is that, with record variables, the record dimension is extended by writing values that don't yet exist in the dataset. This extends all record variables at once, writing "fill values" for record variables for which the data has not yet been written (but *note Fill Values:: to specify different behavior). Calls to NF90_PUT_ATT may be used to change the values of existing attributes, although data that changes after a file is created is typically stored in variables rather than attributes. Finally, you should explicitly close any netCDF datasets into which data has been written by calling NF90_CLOSE before program termination. Otherwise, modifications to the dataset may be lost. 1.5 Adding New Dimensions, Variables, Attributes ================================================ An existing netCDF dataset can be extensively altered. New dimensions, variables, and attributes can be added or existing ones renamed, and existing attributes can be deleted. Existing dimensions, variables, and attributes can be renamed. The following code template lists a typical sequence of calls to add new netCDF components to an existing dataset: NF90_OPEN ! open existing netCDF dataset ... NF90_REDEF ! put it into define mode ... NF90_DEF_DIM ! define additional dimensions (if any) ... NF90_DEF_VAR ! define additional variables (if any) ... NF90_PUT_ATT ! define other attributes (if any) ... NF90_ENDDEF ! check definitions, leave define mode ... NF90_PUT_VAR ! provide new variable values ... NF90_CLOSE ! close netCDF dataset A netCDF dataset is first opened by the NF90_OPEN call. This call puts the open dataset in data mode, which means existing data values can be accessed and changed, existing attributes can be changed (so long as they do not grow), but nothing can be added. To add new netCDF dimensions, variables, or attributes you must enter define mode, by calling NF90_REDEF. In define mode, call NF90_DEF_DIM to define new dimensions, NF90_DEF_VAR to define new variables, and NF90_PUT_ATT to assign new attributes to variables or enlarge old attributes. You can leave define mode and reenter data mode, checking all the new definitions for consistency and committing the changes to disk, by calling NF90_ENDDEF. If you do not wish to reenter data mode, just call NF90_CLOSE, which will have the effect of first calling NF90_ENDDEF. Until the NF90_ENDDEF call, you may back out of all the redefinitions made in define mode and restore the previous state of the netCDF dataset by calling NF90_ABORT. You may also use the NF90_ABORT call to restore the netCDF dataset to a consistent state if the call to NF90_ENDDEF fails. If you have called NF90_CLOSE from definition mode and the implied call to NF90_ENDDEF fails, NF90_ABORT will automatically be called to close the netCDF dataset and leave it in its previous consistent state (before you entered define mode). At most one process should have a netCDF dataset open for writing at one time. The library is designed to provide limited support for multiple concurrent readers with one writer, via disciplined use of the NF90_SYNC function and the NF90_SHARE flag. If a writer makes changes in define mode, such as the addition of new variables, dimensions, or attributes, some means external to the library is necessary to prevent readers from making concurrent accesses and to inform readers to call NF90_SYNC before the next access. 1.6 Error Handling ================== The netCDF library provides the facilities needed to handle errors in a flexible way. Each netCDF function returns an integer status value. If the returned status value indicates an error, you may handle it in any way desired, from printing an associated error message and exiting to ignoring the error indication and proceeding (not recommended!). For simplicity, the examples in this guide check the error status and call a separate function to handle any errors. The NF90_STRERROR function is available to convert a returned integer error status into an error message string. Occasionally, low-level I/O errors may occur in a layer below the netCDF library. For example, if a write operation causes you to exceed disk quotas or to attempt to write to a device that is no longer available, you may get an error from a layer below the netCDF library, but the resulting write error will still be reflected in the returned status value. 1.7 Compiling and Linking with the NetCDF Library ================================================= Details of how to compile and link a program that uses the netCDF C or Fortran interfaces differ, depending on the operating system, the available compilers, and where the netCDF library and include files are installed. Every Fortran 90 procedure or module which references netCDF constants or procedures must have access to the module information created when the netCDF module was compiled. The suffix for this file is "MOD" (or sometimes "mod"). Most F90 compilers allow the user to specify the location of .MOD files, usually with the -I flag. (Some compilers, like absoft, use -p instead). f90 -c -I/usr/local/include mymodule.f90 Starting with version 3.6.2, another method of building the netCDF fortran libraries becomes available. With the -enable-separate-fortran option to configure, the user can specify that the C library should not contain the fortran functions. In these cases an additional library, libnetcdff.a (not the extra "f") will be built. This library contains the fortran functions. For more information about configure options, *Note Specifying the Environment for Building: (netcdf-install)Specifying the Environment for Building. Building separate fortran libraries is required for shared library builds, but is not done, by default, for static library builds. When linking fortran programs without a separate fortran library, programs must link to the netCDF library like this: f90 -o myprogram myprogram.o -L/usr/local/netcdf/lib -lnetcdf 2 Datasets ********** 2.1 Datasets Introduction ========================= This chapter presents the interfaces of the netCDF functions that deal with a netCDF dataset or the whole netCDF library. A netCDF dataset that has not yet been opened can only be referred to by its dataset name. Once a netCDF dataset is opened, it is referred to by a netCDF ID, which is a small nonnegative integer returned when you create or open the dataset. A netCDF ID is much like a file descriptor in C or a logical unit number in FORTRAN. In any single program, the netCDF IDs of distinct open netCDF datasets are distinct. A single netCDF dataset may be opened multiple times and will then have multiple distinct netCDF IDs; however at most one of the open instances of a single netCDF dataset should permit writing. When an open netCDF dataset is closed, the ID is no longer associated with a netCDF dataset. Functions that deal with the netCDF library include: * Get version of library. * Get error message corresponding to a returned error code. The operations supported on a netCDF dataset as a single object are: * Create, given dataset name and whether to overwrite or not. * Open for access, given dataset name and read or write intent. * Put into define mode, to add dimensions, variables, or attributes. * Take out of define mode, checking consistency of additions. * Close, writing to disk if required. * Inquire about the number of dimensions, number of variables, number of global attributes, and ID of the unlimited dimension, if any. * Synchronize to disk to make sure it is current. * Set and unset nofill mode for optimized sequential writes. * After a summary of conventions used in describing the netCDF interfaces, the rest of this chapter presents a detailed description of the interfaces for these operations. 2.2 NetCDF Library Interface Descriptions ========================================= Each interface description for a particular netCDF function in this and later chapters contains: * a description of the purpose of the function; * a Fortran 90 interface block that presents the type and order of the formal parameters to the function; * a description of each formal parameter in the C interface; * a list of possible error conditions; and * an example of a Fortran 90 program fragment calling the netCDF function (and perhaps other netCDF functions). The examples follow a simple convention for error handling, always checking the error status returned from each netCDF function call and calling a handle_error function in case an error was detected. For an example of such a function, see Section 5.2 "Get error message corresponding to error status: nf90_strerror". 2.3 NF90_STRERROR ================= The function NF90_STRERROR returns a static reference to an error message string corresponding to an integer netCDF error status or to a system error number, presumably returned by a previous call to some other netCDF function. The list of netCDF error status codes is available in the appropriate include file for each language binding. Usage ===== function nf90_strerror(ncerr) integer, intent( in) :: ncerr character(len = 80) :: nf90_strerror `NCERR' An error status that might have been returned from a previous call to some netCDF function. Errors ====== If you provide an invalid integer error status that does not correspond to any netCDF error message or or to any system error message (as understood by the system strerror function), NF90_STRERROR returns a string indicating that there is no such error status. Example ======= Here is an example of a simple error handling function that uses NF90_STRERROR to print the error message corresponding to the netCDF error status returned from any netCDF function call and then exit: subroutine handle_err(status) integer, intent ( in) :: status if(status /= nf90_noerr) then print *, trim(nf90_strerror(status)) stop "Stopped" end if end subroutine handle_err 2.4 Get netCDF library version: NF90_INQ_LIBVERS ================================================ The function NF90_INQ_LIBVERS returns a string identifying the version of the netCDF library, and when it was built. Usage ===== function nf90_inq_libvers() character(len = 80) :: nf90_inq_libvers Errors ====== This function takes no arguments, and returns no error status. Example ======= Here is an example using nf90_inq_libvers to print the version of the netCDF library with which the program is linked: print *, trim(nf90_inq_libvers()) 2.5 NF90_CREATE =============== This function creates a new netCDF dataset, returning a netCDF ID that can subsequently be used to refer to the netCDF dataset in other netCDF function calls. The new netCDF dataset opened for write access and placed in define mode, ready for you to add dimensions, variables, and attributes. A creation mode flag specifies whether to overwrite any existing dataset with the same name and whether access to the dataset is shared. Usage ===== function nf90_create(path, cmode, ncid) character (len = *), intent(in ) :: path integer, intent(in ) :: cmode integer, optional, intent(in ) :: initialsize integer, optional, intent(inout) :: chunksize integer, intent( out) :: ncid integer :: nf90_create `path' The file name of the new netCDF dataset. `cmode' The creation mode flag. The following flags are available: NF90_NOCLOBBER, NF90_SHARE, NF90_64BIT_OFFSET, NF90_HDF5, and NF90_CLASSIC_MODEL. A zero value (defined for convenience as NF90_CLOBBER) specifies the default behavior: overwrite any existing dataset with the same file name and buffer and cache accesses for efficiency. The dataset will be in netCDF classic format. *Note NetCDF Classic Format Limitations: (netcdf)NetCDF Classic Format Limitations. Setting NF90_NOCLOBBER means you do not want to clobber (overwrite) an existing dataset; an error (NF90_EEXIST) is returned if the specified dataset already exists. The NF90_SHARE flag is appropriate when one process may be writing the dataset and one or more other processes reading the dataset concurrently; it means that dataset accesses are not buffered and caching is limited. Since the buffering scheme is optimized for sequential access, programs that do not access data sequentially may see some performance improvement by setting the NF90_SHARE flag. (This only applies to netCDF-3 classic or 64-bit offset files.) Setting NF90_64BIT_OFFSET causes netCDF to create a 64-bit offset format file, instead of a netCDF classic format file. The 64-bit offset format imposes far fewer restrictions on very large (i.e. over 2 GB) data files. *Note Large File Support: (netcdf)Large File Support. Setting the NF90_HDF5 flag causes netCDF to create a netCDF-4/HDF5 format output file. Oring the NF90_CLASSIC_MODEL flag with the NF90_HDF5 flag causes the resulting netCDF-4/HDF5 file to restrict itself to the classic model - none of the new netCDF-4 data model features, such as groups or user-defined types, are allowed in such a file. `ncid' Returned netCDF ID. The following optional arguments allow additional performance tuning. `initialsize' The initial size of the file (in bytes) at creation time. A value of 0 causes the file size to be computed when nf90_enddef is called. This is ignored for NetCDF-4/HDF5 files. `chunksize' Controls a space versus time trade-off, memory allocated in the netcdf library versus number of system calls. Because of internal requirements, the value may not be set to exactly the value requested. The actual value chosen is returned. The library chooses a system-dependent default value if NF90_SIZEHINT_DEFAULT is supplied as input. If the "preferred I/O block size" is available from the stat() system call as member st_blksize this value is used. Lacking that, twice the system pagesize is used. Lacking a call to discover the system pagesize, the default chunksize is set to 8192 bytes. The chunksize is a property of a given open netcdf descriptor ncid, it is not a persistent property of the netcdf dataset. This is ignored for NetCDF-4/HDF5 files. Errors ====== NF90_CREATE returns the value NF90_NOERR if no errors occurred. Possible causes of errors include: * Passing a dataset name that includes a directory that does not exist. * Specifying a dataset name of a file that exists and also specifying NF90_NOCLOBBER. * Specifying a meaningless value for the creation mode. * Attempting to create a netCDF dataset in a directory where you don't have permission to create files. Example ======= In this example we create a netCDF dataset named foo.nc; we want the dataset to be created in the current directory only if a dataset with that name does not already exist: use netcdf implicit none integer :: ncid, status ... status = nf90_create(path = "foo.nc", cmode = nf90_noclobber, ncid = ncid) if (status /= nf90_noerr) call handle_err(status) 2.6 NF90_OPEN ============= The function NF90_OPEN opens an existing netCDF dataset for access. Usage ===== function nf90_open(path, mode, ncid, chunksize) character (len = *), intent(in ) :: path integer, intent(in ) :: mode integer, intent( out) :: ncid integer, optional, intent(inout) :: chunksize integer :: nf90_open `path' File name for netCDF dataset to be opened. This may be an OPeNDAP URL if DAP support is enabled. `omode' A zero value (or NF90_NOWRITE) specifies the default behavior: open the dataset with read-only access, buffering and caching accesses for efficiency Otherwise, the creation mode is NF90_WRITE, NF90_SHARE, or NF90_WRITE|NF90_SHARE. Setting the NF90_WRITE flag opens the dataset with read-write access. ("Writing" means any kind of change to the dataset, including appending or changing data, adding or renaming dimensions, variables, and attributes, or deleting attributes.) The NF90_SHARE flag is appropriate when one process may be writing the dataset and one or more other processes reading the dataset concurrently; it means that dataset accesses are not buffered and caching is limited. Since the buffering scheme is optimized for sequential access, programs that do not access data sequentially may see some performance improvement by setting the NF90_SHARE flag. `ncid' Returned netCDF ID. The following optional argument allows additional performance tuning. `chunksize' Controls a space versus time trade-off, memory allocated in the netcdf library versus number of system calls. Because of internal requirements, the value may not be set to exactly the value requested. The actual value chosen is returned. The library chooses a system-dependent default value if NF90_SIZEHINT_DEFAULT is supplied as input. If the "preferred I/O block size" is available from the stat() system call as member st_blksize this value is used. Lacking that, twice the system pagesize is used. Lacking a call to discover the system pagesize, the default chunksize is set to 8192 bytes. The chunksize is a property of a given open netcdf descriptor ncid, it is not a persistent property of the netcdf dataset. Errors ====== NF90_OPEN returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible causes of errors include: * The specified netCDF dataset does not exist. * A meaningless mode was specified. Example ======= Here is an example using NF90_OPEN to open an existing netCDF dataset named foo.nc for read-only, non-shared access: use netcdf implicit none integer :: ncid, status ... status = nf90_open(path = "foo.nc", cmode = nf90_nowrite, ncid = ncid) if (status /= nf90_noerr) call handle_err(status) 2.7 NF90_REDEF ============== The function NF90_REDEF puts an open netCDF dataset into define mode, so dimensions, variables, and attributes can be added or renamed and attributes can be deleted. Usage ===== function nf90_redef(ncid) integer, intent( in) :: ncid integer :: nf90_redef `ncid' netCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. Errors ====== NF90_REDEF returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible causes of errors include: * The specified netCDF dataset is already in define mode. * The specified netCDF dataset was opened for read-only. * The specified netCDF ID does not refer to an open netCDF dataset. Example ======= Here is an example using NF90_REDEF to open an existing netCDF dataset named foo.nc and put it into define mode: use netcdf implicit none integer :: ncid, status ... status = nf90_open("foo.nc", nf90_write, ncid) ! Open dataset if (status /= nf90_noerr) call handle_err(status) ... status = nf90_redef(ncid) ! Put the file in define mode if (status /= nf90_noerr) call handle_err(status) 2.8 NF90_ENDDEF =============== The function NF90_ENDDEF takes an open netCDF dataset out of define mode. The changes made to the netCDF dataset while it was in define mode are checked and committed to disk if no problems occurred. Non-record variables may be initialized to a "fill value" as well (*note NF90_SET_FILL::). The netCDF dataset is then placed in data mode, so variable data can be read or written. This call may involve copying data under some circumstances. For a more extensive discussion *Note File Structure and Performance: (netcdf)File Structure and Performance. Usage ===== function nf90_enddef(ncid, h_minfree, v_align, v_minfree, r_align) integer, intent( in) :: ncid integer, optional, intent( in) :: h_minfree, v_align, v_minfree, r_align integer :: nf90_enddef `ncid' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. The following arguments allow additional performance tuning. Note: these arguments expose internals of the netcdf version 1 file format, and may not be available in future netcdf implementations. The current netcdf file format has three sections: the "header" section, the data section for fixed size variables, and the data section for variables which have an unlimited dimension (record variables). The header begins at the beginning of the file. The index (offset) of the beginning of the other two sections is contained in the header. Typically, there is no space between the sections. This causes copying overhead to accrue if one wishes to change the size of the sections, as may happen when changing the names of things, text attribute values, adding attributes or adding variables. Also, for buffered i/o, there may be advantages to aligning sections in certain ways. The minfree parameters allow one to control costs of future calls to nf90_redef or nf90_enddef by requesting that some space be available at the end of the section. The default value for both h_minfree and v_minfree is 0. The align parameters allow one to set the alignment of the beginning of the corresponding sections. The beginning of the section is rounded up to an index which is a multiple of the align parameter. The flag value NF90_ALIGN_CHUNK tells the library to use the chunksize (see above) as the align parameter. The default value for both v_align and r_align is 4 bytes. `h_minfree' Size of the pad (in bytes) at the end of the "header" section. `v_minfree' Size of the pad (in bytes) at the end of the data section for fixed size variables. `v_align' The alignment of the beginning of the data section for fixed size variables. `r_align' The alignment of the beginning of the data section for variables which have an unlimited dimension (record variables). Errors ====== NF90_ENDDEF returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible causes of errors include: * The specified netCDF dataset is not in define mode. * The specified netCDF ID does not refer to an open netCDF dataset. * The size of one or more variables exceed the size constraints for whichever variant of the file format is in use). *Note Large File Support: (netcdf)Large File Support. Example ======= Here is an example using NF90_ENDDEF to finish the definitions of a new netCDF dataset named foo.nc and put it into data mode: use netcdf implicit none integer :: ncid, status ... status = nf90_create("foo.nc", nf90_noclobber, ncid) if (status /= nf90_noerr) call handle_err(status) ... ! create dimensions, variables, attributes status = nf90_enddef(ncid) if (status /= nf90_noerr) call handle_err(status) 2.9 NF90_CLOSE ============== The function NF90_CLOSE closes an open netCDF dataset. If the dataset is in define mode, NF90_ENDDEF will be called before closing. (In this case, if NF90_ENDDEF returns an error, NF90_ABORT will automatically be called to restore the dataset to the consistent state before define mode was last entered.) After an open netCDF dataset is closed, its netCDF ID may be reassigned to the next netCDF dataset that is opened or created. Usage ===== function nf90_close(ncid) integer, intent( in) :: ncid integer :: nf90_close `ncid' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. Errors ====== NF90_CLOSE returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible causes of errors include: * Define mode was entered and the automatic call made to NF90_ENDDEF failed. * The specified netCDF ID does not refer to an open netCDF dataset. Example ======= Here is an example using NF90_CLOSE to finish the definitions of a new netCDF dataset named foo.nc and release its netCDF ID: use netcdf implicit none integer :: ncid, status ... status = nf90_create("foo.nc", nf90_noclobber, ncid) if (status /= nf90_noerr) call handle_err(status) ... ! create dimensions, variables, attributes status = nf90_close(ncid) if (status /= nf90_noerr) call handle_err(status) 2.10 NF90_INQUIRE Family ======================== The NF90_INQUIRE subroutine returns information about an open netCDF dataset, given its netCDF ID. The subroutine can be called from either define mode or data mode, and returns values for any or all of the following: the number of dimensions, the number of variables, the number of global attributes, and the dimension ID of the dimension defined with unlimited length, if any. An additional function, NF90_INQ_FORMAT, returns the (rarely needed) format version. No I/O is performed when NF90_INQUIRE is called, since the required information is available in memory for each open netCDF dataset. Usage ===== function nf90_inquire(ncid, nDimensions, nVariables, nAttributes, & unlimitedDimId, formatNum) integer, intent( in) :: ncid integer, optional, intent(out) :: nDimensions, nVariables, & nAttributes, unlimitedDimId, & formatNum integer :: nf90_inquire `ncid' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `nDimensions' Returned number of dimensions defined for this netCDF dataset. `nVariables' Returned number of variables defined for this netCDF dataset. `nAttributes' Returned number of global attributes defined for this netCDF dataset. `unlimitedDimID' Returned ID of the unlimited dimension, if there is one for this netCDF dataset. If no unlimited length dimension has been defined, -1 is returned. `format' Returned integer indicating format version for this dataset, one of nf90_format_classic, nf90_format_64bit, nf90_format_netcdf4, or nf90_format_netcdf4_classic. These are rarely needed by users or applications, since thhe library recognizes the format of a file it is accessing and handles it accordingly. Errors ====== Function NF90_INQUIRE returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible causes of errors include: * The specified netCDF ID does not refer to an open netCDF dataset. Example ======= Here is an example using NF90_INQUIRE to find out about a netCDF dataset named foo.nc: use netcdf implicit none integer :: ncid, status, nDims, nVars, nGlobalAtts, unlimDimID ... status = nf90_open("foo.nc", nf90_nowrite, ncid) if (status /= nf90_noerr) call handle_err(status) ... status = nf90_inquire(ncid, nDims, nVars, nGlobalAtts, unlimdimid) if (status /= nf90_noerr) call handle_err(status) status = nf90_inquire(ncid, nDimensions = nDims, & unlimitedDimID = unlimdimid) if (status /= nf90_noerr) call handle_err(status) 2.11 NF90_SYNC ============== The function NF90_SYNC offers a way to synchronize the disk copy of a netCDF dataset with in-memory buffers. There are two reasons you might want to synchronize after writes: * To minimize data loss in case of abnormal termination, or * To make data available to other processes for reading immediately after it is written. But note that a process that already had the dataset open for reading would not see the number of records increase when the writing process calls NF90_SYNC; to accomplish this, the reading process must call NF90_SYNC. This function is backward-compatible with previous versions of the netCDF library. The intent was to allow sharing of a netCDF dataset among multiple readers and one writer, by having the writer call NF90_SYNC after writing and the readers call NF90_SYNC before each read. For a writer, this flushes buffers to disk. For a reader, it makes sure that the next read will be from disk rather than from previously cached buffers, so that the reader will see changes made by the writing process (e.g., the number of records written) without having to close and reopen the dataset. If you are only accessing a small amount of data, it can be expensive in computer resources to always synchronize to disk after every write, since you are giving up the benefits of buffering. An easier way to accomplish sharing (and what is now recommended) is to have the writer and readers open the dataset with the NF90_SHARE flag, and then it will not be necessary to call NF90_SYNC at all. However, the NF90_SYNC function still provides finer granularity than the NF90_SHARE flag, if only a few netCDF accesses need to be synchronized among processes. It is important to note that changes to the ancillary data, such as attribute values, are not propagated automatically by use of the NF90_SHARE flag. Use of the NF90_SYNC function is still required for this purpose. Sharing datasets when the writer enters define mode to change the data schema requires extra care. In previous releases, after the writer left define mode, the readers were left looking at an old copy of the dataset, since the changes were made to a new copy. The only way readers could see the changes was by closing and reopening the dataset. Now the changes are made in place, but readers have no knowledge that their internal tables are now inconsistent with the new dataset schema. If netCDF datasets are shared across redefinition, some mechanism external to the netCDF library must be provided that prevents access by readers during redefinition and causes the readers to call NF90_SYNC before any subsequent access. When calling NF90_SYNC, the netCDF dataset must be in data mode. A netCDF dataset in define mode is synchronized to disk only when NF90_ENDDEF is called. A process that is reading a netCDF dataset that another process is writing may call NF90_SYNC to get updated with the changes made to the data by the writing process (e.g., the number of records written), without having to close and reopen the dataset. Data is automatically synchronized to disk when a netCDF dataset is closed, or whenever you leave define mode. Usage ===== function nf90_sync(ncid) integer, intent( in) :: ncid integer :: nf90_sync `ncid' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. Errors ====== NF90_SYNC returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible causes of errors include: * The netCDF dataset is in define mode. * The specified netCDF ID does not refer to an open netCDF dataset. Example ======= Here is an example using NF90_SYNC to synchronize the disk writes of a netCDF dataset named foo.nc: use netcdf implicit none integer :: ncid, status ... status = nf90_open("foo.nc", nf90_write, ncid) if (status /= nf90_noerr) call handle_err(status) ... ! write data or change attributes ... status = NF90_SYNC(ncid) if (status /= nf90_noerr) call handle_err(status) 2.12 NF90_ABORT =============== You no longer need to call this function, since it is called automatically by NF90_CLOSE in case the dataset is in define mode and something goes wrong with committing the changes. The function NF90_ABORT just closes the netCDF dataset, if not in define mode. If the dataset is being created and is still in define mode, the dataset is deleted. If define mode was entered by a call to NF90_REDEF, the netCDF dataset is restored to its state before definition mode was entered and the dataset is closed. Usage ===== function nf90_abort(ncid) integer, intent( in) :: ncid integer :: nf90_abort `ncid' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. Errors ====== NF90_ABORT returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible causes of errors include: * When called from define mode while creating a netCDF dataset, deletion of the dataset failed. * The specified netCDF ID does not refer to an open netCDF dataset. Example ======= Here is an example using NF90_ABORT to back out of redefinitions of a dataset named foo.nc: use netcdf implicit none integer :: ncid, status, LatDimID ... status = nf90_open("foo.nc", nf90_write, ncid) if (status /= nf90_noerr) call handle_err(status) ... status = nf90_redef(ncid) if (status /= nf90_noerr) call handle_err(status) ... status = nf90_def_dim(ncid, "Lat", 18, LatDimID) if (status /= nf90_noerr) then ! Dimension definition failed call handle_err(status) status = nf90_abort(ncid) ! Abort redefinitions if (status /= nf90_noerr) call handle_err(status) end if ... 2.13 NF90_SET_FILL ================== This function is intended for advanced usage, to optimize writes under some circumstances described below. The function NF90_SET_FILL sets the fill mode for a netCDF dataset open for writing and returns the current fill mode in a return parameter. The fill mode can be specified as either NF90_FILL or NF90_NOFILL. The default behavior corresponding to NF90_FILL is that data is pre-filled with fill values, that is fill values are written when you create non-record variables or when you write a value beyond data that has not yet been written. This makes it possible to detect attempts to read data before it was written. *Note Fill Values::, for more information on the use of fill values. *Note Attribute Conventions::, for information about how to define your own fill values. The behavior corresponding to NF90_NOFILL overrides the default behavior of prefilling data with fill values. This can be used to enhance performance, because it avoids the duplicate writes that occur when the netCDF library writes fill values that are later overwritten with data. A value indicating which mode the netCDF dataset was already in is returned. You can use this value to temporarily change the fill mode of an open netCDF dataset and then restore it to the previous mode. After you turn on NF90_NOFILL mode for an open netCDF dataset, you must be certain to write valid data in all the positions that will later be read. Note that nofill mode is only a transient property of a netCDF dataset open for writing: if you close and reopen the dataset, it will revert to the default behavior. You can also revert to the default behavior by calling NF90_SET_FILL again to explicitly set the fill mode to NF90_FILL. There are three situations where it is advantageous to set nofill mode: 1. Creating and initializing a netCDF dataset. In this case, you should set nofill mode before calling NF90_ENDDEF and then write completely all non-record variables and the initial records of all the record variables you want to initialize. 2. Extending an existing record-oriented netCDF dataset. Set nofill mode after opening the dataset for writing, then append the additional records to the dataset completely, leaving no intervening unwritten records. 3. Adding new variables that you are going to initialize to an existing netCDF dataset. Set nofill mode before calling NF90_ENDDEF then write all the new variables completely. If the netCDF dataset has an unlimited dimension and the last record was written while in nofill mode, then the dataset may be shorter than if nofill mode was not set, but this will be completely transparent if you access the data only through the netCDF interfaces. The use of this feature may not be available (or even needed) in future releases. Programmers are cautioned against heavy reliance upon this feature. Usage ===== function nf90_set_fill(ncid, fillmode, old_mode) integer, intent( in) :: ncid, fillmode integer, intent(out) :: old_mode integer :: nf90_set_fill `ncid' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `fillmode' Desired fill mode for the dataset, either NF90_NOFILL or NF90_FILL. `old_mode' Returned current fill mode of the dataset before this call, either NF90_NOFILL or NF90_FILL. Errors ====== NF90_SET_FILL returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible causes of errors include: * The specified netCDF ID does not refer to an open netCDF dataset. * The specified netCDF ID refers to a dataset open for read-only access. * The fill mode argument is neither NF90_NOFILL nor NF90_FILL.. Example ======= Here is an example using NF90_SET_FILL to set nofill mode for subsequent writes of a netCDF dataset named foo.nc: use netcdf implicit none integer :: ncid, status, oldMode ... status = nf90_open("foo.nc", nf90_write, ncid) if (status /= nf90_noerr) call handle_err(status) ... ! Write data with prefilling behavior ... status = nf90_set_fill(ncid, nf90_nofill, oldMode) if (status /= nf90_noerr) call handle_err(status) ... ! Write data with no prefilling ... 3 Groups ******** NetCDF-4 added support for hierarchical groups within netCDF datasets. Groups are identified with a ncid, which identifies both the open file, and the group within that file. When a file is opened with NF90_OPEN or NF90_CREATE, the ncid for the root group of that file is provided. Using that as a starting point, users can add new groups, or list and navigate existing groups. All netCDF calls take a ncid which determines where the call will take its action. For example, the NF90_DEF_VAR function takes a ncid as its first parameter. It will create a variable in whichever group its ncid refers to. Use the root ncid provided by NF90_CREATE or NF90_OPEN to create a variable in the root group. Or use NF90_DEF_GRP to create a group and use its ncid to define a variable in the new group. Variable are only visible in the group in which they are defined. The same applies to attributes. "Global" attributes are defined in whichever group is refered to by the ncid. Dimensions are visible in their groups, and all child groups. Group operations are only permitted on netCDF-4 files - that is, files created with the HDF5 flag in nf90_create. (*note NF90_CREATE::). Groups are not compatible with the netCDF classic data model, so files created with the NF90_CLASSIC_MODEL file cannot contain groups (except the root group). 3.1 Find a Group ID: NF90_INQ_NCID ================================== Given an ncid and group name (NULL or "" gets root group), return ncid of the named group. Usage ===== function nf90_inq_ncid(ncid, name, grp_ncid) integer, intent(in) :: ncid character (len = *), intent(in) :: name integer, intent(out) :: grp_ncid integer :: nf90_inq_ncid `NCID' The group id for this operation. `NAME' A character array that holds the name of the desired group. Must be less then NF90_MAX_NAME. `GRPID' The ID of the group will go here. Errors ====== `NF90_NOERR' No error. `NF90_EBADID' Bad group id. `NF90_ENOTNC4' Attempting a netCDF-4 operation on a netCDF-3 file. NetCDF-4 operations can only be performed on files defined with a create mode which includes flag HDF5. (*note NF90_OPEN::). `NF90_ESTRICTNC3' This file was created with the strict netcdf-3 flag, therefore netcdf-4 operations are not allowed. (*note NF90_OPEN::). `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= This example is from nf90_test/ftst_groups.F. 3.2 Get a List of Groups in a Group: NF90_INQ_GRPS ================================================== Given a location id, return the number of groups it contains, and an array of their ncids. Usage ===== function nf90_inq_grps(ncid, numgrps, ncids) integer, intent(in) :: ncid integer, intent(out) :: numgrps integer, intent(out) :: ncids integer :: nf90_inq_grps `NCID' The group id for this operation. `NUMGRPS' An integer which will get number of groups in this group. `NCIDS' An array of ints which will receive the IDs of all the groups in this group. Errors ====== `NF90_NOERR' No error. `NF90_EBADID' Bad group id. `NF90_ENOTNC4' Attempting a netCDF-4 operation on a netCDF-3 file. NetCDF-4 operations can only be performed on files defined with a create mode which includes flag HDF5. (*note NF90_OPEN::). `NF90_ESTRICTNC3' This file was created with the strict netcdf-3 flag, therefore netcdf-4 operations are not allowed. (*note NF90_OPEN::). `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 3.3 Find all the Variables in a Group: NF90_INQ_VARIDS ====================================================== Find all varids for a location. Usage ===== function nf90_inq_varids(ncid, nvars, varids) integer, intent(in) :: ncid integer, intent(out) :: nvars integer, intent(out) :: varids integer :: nf90_inq_varids `NCID' The group id for this operation. `VARIDS' An already allocated array to store the list of varids. Use nf90_inq_nvars to find out how many variables there are. (*note NF90_INQUIRE_VARIABLE::). Errors ====== `NF90_NOERR' No error. `NF90_EBADID' Bad group id. `NF90_ENOTNC4' Attempting a netCDF-4 operation on a netCDF-3 file. NetCDF-4 operations can only be performed on files defined with a create mode which includes flag HDF5. (*note NF90_OPEN::). `NF90_ESTRICTNC3' This file was created with the strict netcdf-3 flag, therefore netcdf-4 operations are not allowed. (*note NF90_OPEN::). `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 3.4 Find all Dimensions Visible in a Group: NF90_INQ_DIMIDS =========================================================== Find all dimids for a location. This finds all dimensions in a group, or any of its parents. Usage ===== function nf90_inq_dimids(ncid, ndims, dimids, include_parents) integer, intent(in) :: ncid integer, intent(out) :: ndims integer, intent(out) :: dimids integer, intent(out) :: include_parents integer :: nf90_inq_dimids `NCID' The group id for this operation. `DIMIDS' An array of ints when the dimids of the visible dimensions will be stashed. Use nf90_inq_ndims to find out how many dims are visible from this group. (*note NF90_INQUIRE_VARIABLE::). `INCLUDE_PARENTS' If zero, only the group specified by NCID will be searched for dimensions. Otherwise parent groups will be searched too. Errors ====== `NF90_NOERR' No error. `NF90_EBADID' Bad group id. `NF90_ENOTNC4' Attempting a netCDF-4 operation on a netCDF-3 file. NetCDF-4 operations can only be performed on files defined with a create mode which includes flag HDF5. (*note NF90_OPEN::). `NF90_ESTRICTNC3' This file was created with the strict netcdf-3 flag, therefore netcdf-4 operations are not allowed. (*note NF90_OPEN::). `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 3.5 Find the Length of a Group's Full Name: NF90_INQ_GRPNAME_LEN ================================================================ Given ncid, find length of the full name. (Root group is named "/", with length 1.) Usage ===== function nf90_inq_grpname_len(ncid, len) integer, intent(in) :: ncid integer, intent(out) :: len integer :: nf90_inq_grpname_len nf90_inq_grpname_len = nf_inq_grpname_len(ncid, len) end function nf90_inq_grpname_len `NCID' The group id for this operation. `LEN' An integer where the length will be placed. Errors ====== `NF90_NOERR' No error. `NF90_EBADID' Bad group id. `NF90_ENOTNC4' Attempting a netCDF-4 operation on a netCDF-3 file. NetCDF-4 operations can only be performed on files defined with a create mode which includes flag HDF5. (*note NF90_OPEN::). `NF90_ESTRICTNC3' This file was created with the strict netcdf-3 flag, therefore netcdf-4 operations are not allowed. (*note NF90_OPEN::). `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 3.6 Find a Group's Name: NF90_INQ_GRPNAME ========================================= Given ncid, find relative name of group. (Root group is named "/"). The name provided by this function is relative to the parent group. For a full path name for the group is, with all parent groups included, separated with a forward slash (as in Unix directory names) *Note NF90_INQ_GRPNAME_FULL::. Usage ===== function nf90_inq_grpname(ncid, name) integer, intent(in) :: ncid character (len = *), intent(out) :: name integer :: nf90_inq_grpname `NCID' The group id for this operation. `NAME' The name of the group will be copied to this character array. The name will be less than NF90_MAX_NAME in length. `' Errors ====== `NF90_NOERR' No error. `NF90_EBADID' Bad group id. `NF90_ENOTNC4' Attempting a netCDF-4 operation on a netCDF-3 file. NetCDF-4 operations can only be performed on files defined with a create mode which includes flag HDF5. (*note NF90_OPEN::). `NF90_ESTRICTNC3' This file was created with the strict netcdf-3 flag, therefore netcdf-4 operations are not allowed. (*note NF90_OPEN::). `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 3.7 Find a Group's Full Name: NF90_INQ_GRPNAME_FULL =================================================== Given ncid, find complete name of group. (Root group is named "/"). The name provided by this function is a full path name for the group is, with all parent groups included, separated with a forward slash (as in Unix directory names). For a name relative to the parent group *Note NF90_INQ_GRPNAME::. To find the length of the full name *Note NF90_INQ_GRPNAME_LEN::. Usage ===== function nf90_inq_grpname_full(ncid, len, name) integer, intent(in) :: ncid integer, intent(out) :: len character (len = *), intent(out) :: name integer :: nf90_inq_grpname_full `NCID' The group id for this operation. `LEN' The length of the full group name will go here. `NAME' The name of the group will be copied to this character array. Errors ====== `NF90_NOERR' No error. `NF90_EBADID' Bad group id. `NF90_ENOTNC4' Attempting a netCDF-4 operation on a netCDF-3 file. NetCDF-4 operations can only be performed on files defined with a create mode which includes flag HDF5. (*note NF90_OPEN::). `NF90_ESTRICTNC3' This file was created with the strict netcdf-3 flag, therefore netcdf-4 operations are not allowed. (*note NF90_OPEN::). `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 3.8 Find a Group's Parent: NF90_INQ_GRP_PARENT ============================================== Given ncid, find the ncid of the parent group. When used with the root group, this function returns the NF90_ENOGRP error (since the root group has no parent.) Usage ===== function nf90_inq_grp_parent(ncid, parent_ncid) integer, intent(in) :: ncid integer, intent(out) :: parent_ncid integer :: nf90_inq_grp_parent `NCID' The group id. `PARENT_NCID' The ncid of the parent group will be copied here. Errors ====== `NF90_NOERR' No error. `NF90_EBADID' Bad group id. `NF90_ENOGRP' No parent group found (i.e. this is the root group). `NF90_ENOTNC4' Attempting a netCDF-4 operation on a netCDF-3 file. NetCDF-4 operations can only be performed on files defined with a create mode which includes flag HDF5. (*note NF90_OPEN::). `NF90_ESTRICTNC3' This file was created with the strict netcdf-3 flag, therefore netcdf-4 operations are not allowed. (*note NF90_OPEN::). `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 3.9 Create a New Group: NF90_DEF_GRP ==================================== Create a group. Its location id is returned in new_ncid. Usage ===== function nf90_def_grp(parent_ncid, name, new_ncid) integer, intent(in) :: parent_ncid character (len = *), intent(in) :: name integer, intent(out) :: new_ncid integer :: nf90_def_grp `PARENT_NCID' The group id of the parent group. `NAME' The name of the new group. `NEW_NCID' The ncid of the new group will be placed there. Errors ====== `NF90_NOERR' No error. `NF90_EBADID' Bad group id. `NF90_ENAMEINUSE' That name is in use. Group names must be unique within a group. `NF90_EMAXNAME' Name exceed max length NF90_MAX_NAME. `NF90_EBADNAME' Name contains illegal characters. `NF90_ENOTNC4' Attempting a netCDF-4 operation on a netCDF-3 file. NetCDF-4 operations can only be performed on files defined with a create mode which includes flag HDF5. (*note NF90_OPEN::). `NF90_ESTRICTNC3' This file was created with the strict netcdf-3 flag, therefore netcdf-4 operations are not allowed. (*note NF90_OPEN::). `NF90_EHDFERR' An error was reported by the HDF5 layer. `NF90_EPERM' Attempt to write to a read-only file. `NF90_ENOTINDEFINE' Not in define mode. Example ======= C Create the netCDF file. retval = nf90_create(file_name, NF90_NETCDF4, ncid) if (retval .ne. nf90_noerr) call handle_err(retval) C Create a group and a subgroup. retval = nf90_def_grp(ncid, group_name, grpid) if (retval .ne. nf90_noerr) call handle_err(retval) retval = nf90_def_grp(grpid, sub_group_name, sub_grpid) if (retval .ne. nf90_noerr) call handle_err(retval) 4 Dimensions ************ 4.1 Dimensions Introduction =========================== Dimensions for a netCDF dataset are defined when it is created, while the netCDF dataset is in define mode. Additional dimensions may be added later by reentering define mode. A netCDF dimension has a name and a length. At most one dimension in a netCDF dataset can have the unlimited length, which means variables using this dimension can grow along this dimension. There is a suggested limit (512) to the number of dimensions that can be defined in a single netCDF dataset. The limit is the value of the constant NF90_MAX_DIMS. The purpose of the limit is to make writing generic applications simpler. They need only provide an array of NF90_MAX_DIMS dimensions to handle any netCDF dataset. The implementation of the netCDF library does not enforce this advisory maximum, so it is possible to use more dimensions, if necessary, but netCDF utilities that assume the advisory maximums may not be able to handle the resulting netCDF datasets. Ordinarily, the name and length of a dimension are fixed when the dimension is first defined. The name may be changed later, but the length of a dimension (other than the unlimited dimension) cannot be changed without copying all the data to a new netCDF dataset with a redefined dimension length. A netCDF dimension in an open netCDF dataset is referred to by a small integer called a dimension ID. In the Fortran 90 interface, dimension IDs are 1, 2, 3, ..., in the order in which the dimensions were defined. Operations supported on dimensions are: * Create a dimension, given its name and length. * Get a dimension ID from its name. * Get a dimension's name and length from its ID. * Rename a dimension. 4.2 NF90_DEF_DIM ================ The function NF90_DEF_DIM adds a new dimension to an open netCDF dataset in define mode. It returns (as an argument) a dimension ID, given the netCDF ID, the dimension name, and the dimension length. At most one unlimited length dimension, called the record dimension, may be defined for each netCDF dataset. Usage ===== function nf90_def_dim(ncid, name, len, dimid) integer, intent( in) :: ncid character (len = *), intent( in) :: name integer, intent( in) :: len integer, intent(out) :: dimid integer :: nf90_def_dim `ncid' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `name' Dimension name. Must begin with an alphabetic character, followed by zero or more alphanumeric characters including the underscore ('_'). Case is significant. `len' Length of dimension; that is, number of values for this dimension as an index to variables that use it. This should be either a positive integer or the predefined constant NF90_UNLIMITED. `dimid' Returned dimension ID. Errors ====== NF90_DEF_DIM returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible causes of errors include: * The netCDF dataset is not in definition mode. * The specified dimension name is the name of another existing dimension. * The specified length is not greater than zero. * The specified length is unlimited, but there is already an unlimited length dimension defined for this netCDF dataset. * The specified netCDF ID does not refer to an open netCDF dataset. Example ======= Here is an example using NF90_DEF_DIM to create a dimension named lat of length 18 and a unlimited dimension named rec in a new netCDF dataset named foo.nc: use netcdf implicit none integer :: ncid, status, LatDimID, RecordDimID ... status = nf90_create("foo.nc", nf90_noclobber, ncid) if (status /= nf90_noerr) call handle_err(status) ... status = nf90_def_dim(ncid, "Lat", 18, LatDimID) if (status /= nf90_noerr) call handle_err(status) status = nf90_def_dim(ncid, "Record", nf90_unlimited, RecordDimID) if (status /= nf90_noerr) call handle_err(status) 4.3 NF90_INQ_DIMID ================== The function NF90_INQ_DIMID returns (as an argument) the ID of a netCDF dimension, given the name of the dimension. If ndims is the number of dimensions defined for a netCDF dataset, each dimension has an ID between 1 and ndims. Usage ===== function nf90_inq_dimid(ncid, name, dimid) integer, intent( in) :: ncid character (len = *), intent( in) :: name integer, intent(out) :: dimid integer :: nf90_inq_dimid `ncid' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `name' Dimension name, a character string beginning with a letter and followed by any sequence of letters, digits, or underscore ('_') characters. Case is significant in dimension names. `dimid' Returned dimension ID. Errors ====== NF90_INQ_DIMID returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible causes of errors include: * The name that was specified is not the name of a dimension in the netCDF dataset. * The specified netCDF ID does not refer to an open netCDF dataset. Example ======= Here is an example using NF90_INQ_DIMID to determine the dimension ID of a dimension named lat, assumed to have been defined previously in an existing netCDF dataset named foo.nc: use netcdf implicit none integer :: ncid, status, LatDimID ... status = nf90_open("foo.nc", nf90_nowrite, ncid) if (status /= nf90_noerr) call handle_err(status) ... status = nf90_inq_dimid(ncid, "Lat", LatDimID) if (status /= nf90_noerr) call handle_err(status) 4.4 NF90_INQUIRE_DIMENSION ========================== This function information about a netCDF dimension. Information about a dimension includes its name and its length. The length for the unlimited dimension, if any, is the number of records written so far. Usage ===== function nf90_inquire_dimension(ncid, dimid, name, len) integer, intent( in) :: ncid, dimid character (len = *), optional, intent(out) :: name integer, optional, intent(out) :: len integer :: nf90_inquire_dimension `ncid' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `dimid' Dimension ID, from a previous call to NF90_INQ_DIMID or NF90_DEF_DIM. `name' Returned dimension name. The caller must allocate space for the returned name. The maximum possible length, in characters, of a dimension name is given by the predefined constant NF90_MAX_NAME. `len' Returned length of dimension. For the unlimited dimension, this is the current maximum value used for writing any variables with this dimension, that is the maximum record number. Errors ====== These functions return the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible causes of errors include: * The dimension ID is invalid for the specified netCDF dataset. * The specified netCDF ID does not refer to an open netCDF dataset. Example ======= Here is an example using NF90_INQ_DIM to determine the length of a dimension named lat, and the name and current maximum length of the unlimited dimension for an existing netCDF dataset named foo.nc: use netcdf implicit none integer :: ncid, status, LatDimID, RecordDimID integer :: nLats, nRecords character(len = nf90_max_name) :: RecordDimName ... status = nf90_open("foo.nc", nf90_nowrite, ncid) if (status /= nf90_noerr) call handle_err(status) ! Get ID of unlimited dimension status = nf90_inquire(ncid, unlimitedDimId = RecordDimID) if (status /= nf90_noerr) call handle_err(status) ... status = nf90_inq_dimid(ncid, "Lat", LatDimID) if (status /= nf90_noerr) call handle_err(status) ! How many values of "lat" are there? status = nf90_inquire_dimension(ncid, LatDimID, len = nLats) if (status /= nf90_noerr) call handle_err(status) ! What is the name of the unlimited dimension, how many records are there? status = nf90_inquire_dimension(ncid, RecordDimID, & name = RecordDimName, len = Records) if (status /= nf90_noerr) call handle_err(status) 4.5 NF90_RENAME_DIM =================== The function NF90_RENAME_DIM renames an existing dimension in a netCDF dataset open for writing. If the new name is longer than the old name, the netCDF dataset must be in define mode. You cannot rename a dimension to have the same name as another dimension. Usage ===== function nf90_rename_dim(ncid, dimid, name) integer, intent( in) :: ncid character (len = *), intent( in) :: name integer, intent( in) :: dimid integer :: nf90_rename_dim `ncid' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `dimid' Dimension ID, from a previous call to NF90_INQ_DIMID or NF90_DEF_DIM. `name' New dimension name. Errors ====== NF90_RENAME_DIM returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible causes of errors include: * The new name is the name of another dimension. * The dimension ID is invalid for the specified netCDF dataset. * The specified netCDF ID does not refer to an open netCDF dataset. * The new name is longer than the old name and the netCDF dataset is not in define mode. Example ======= Here is an example using NF90_RENAME_DIM to rename the dimension lat to latitude in an existing netCDF dataset named foo.nc: use netcdf implicit none integer :: ncid, status, LatDimID ... status = nf90_open("foo.nc", nf90_write, ncid) if (status /= nf90_noerr) call handle_err(status) ... ! Put in define mode so we can rename the dimension status = nf90_redef(ncid) if (status /= nf90_noerr) call handle_err(status) ! Get the dimension ID for "Lat"... status = nf90_inq_dimid(ncid, "Lat", LatDimID) if (status /= nf90_noerr) call handle_err(status) ! ... and change the name to "Latitude". status = nf90_rename_dim(ncid, LatDimID, "Latitude") if (status /= nf90_noerr) call handle_err(status) ! Leave define mode status = nf90_enddef(ncid) if (status /= nf90_noerr) call handle_err(status) 5 User Defined Data Types ************************* 5.1 User Defined Types Introduction =================================== NetCDF-4 has added support for four different user defined data types. `compound type' Like a C struct, a compound type is a collection of types, including other user defined types, in one package. `variable length array type' The variable length array may be used to store ragged arrays. `opaque type' This type has only a size per element, and no other type information. `enum type' Like an enumeration in C, this type lets you assign text values to integer values, and store the integer values. Users may construct user defined type with the various nf90_def_* functions described in this section. They may learn about user defined types by using the nf90_inq_ functions defined in this section. 5.2 Learn the IDs of All Types in Group: NF90_INQ_TYPEIDS ========================================================= Learn the number of types defined in a group, and their IDs. Usage ===== function nf90_inq_typeids(ncid, ntypes, typeids) integer, intent(in) :: ncid integer, intent(out) :: ntypes integer, intent(out) :: typeids integer :: nf90_inq_typeids `NCID' The group id. `NTYPES' A pointer to int which will get the number of types defined in the group. If NULL, ignored. `TYPEIDS' A pointer to an int array which will get the typeids. If NULL, ignored. Errors ====== `NF90_NOERR' No error. `NF90_BADID' Bad ncid. Example ======= 5.3 Learn About an User Defined Type: NF90_INQ_TYPE =================================================== Given an ncid and a typeid, get the information about a type. This function will work on any type, including atomic and any user defined type, whether compound, opaque, enumeration, or variable length array. For even more information about a user defined type *note NF90_INQ_USER_TYPE::. Usage ===== function nf90_inq_type(ncid, xtype, name, size, nfields) integer, intent(in) :: ncid integer, intent(in) :: xtype character (len = *), intent(out) :: name integer, intent(out) :: size integer, intent(out) :: nfields integer :: nf90_inq_type `NCID' The ncid for the group containing the type (ignored for atomic types). `XTYPE' The typeid for this type, as returned by NF90_DEF_COMPOUND, NF90_DEF_OPAQUE, NF90_DEF_ENUM, NF90_DEF_VLEN, or NF90_INQ_VAR, or as found in netcdf.inc in the list of atomic types (NF90_CHAR, NF90_INT, etc.). `NAME' The name of the user defined type will be copied here. It will be NF90_MAX_NAME bytes or less. For atomic types, the type name from CDL will be given. `SIZEP' The size of the type (in bytes) will be copied here. VLEN type size is the size of one element of the VLEN. String size is returned as zero, since it varies from string to string. Return Codes ============ `NF90_NOERR' No error. `NF90_EBADTYPEID' Bad typeid. `NF90_ENOTNC4' Seeking a user-defined type in a netCDF-3 file. `NF90_ESTRICTNC3' Seeking a user-defined type in a netCDF-4 file for which classic model has been turned on. `NF90_EBADGRPID' Bad group ID in ncid. `NF90_EBADID' Type ID not found. `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 5.4 Learn About an User Defined Type: NF90_INQ_USER_TYPE ======================================================== Given an ncid and a typeid, get the information about a user defined type. This function will work on any user defined type, whether compound, opaque, enumeration, or variable length array. Usage ===== function nf90_inq_user_type(ncid, xtype, name, size, base_typeid, nfields, class) integer, intent(in) :: ncid integer, intent(in) :: xtype character (len = *), intent(out) :: name integer, intent(out) :: size integer, intent(out) :: base_typeid integer, intent(out) :: nfields integer, intent(out) :: class integer :: nf90_inq_user_type `NCID' The ncid for the group containing the user defined type. `XTYPE' The typeid for this type, as returned by NF90_DEF_COMPOUND, NF90_DEF_OPAQUE, NF90_DEF_ENUM, NF90_DEF_VLEN, or NF90_INQ_VAR. `NAME' The name of the user defined type will be copied here. It will be NF90_MAX_NAME bytes or less. `SIZE' The size of the user defined type will be copied here. `BASE_NF90_TYPE' The base typeid will be copied here for vlen and enum types. `NFIELDS' The number of fields will be copied here for enum and compound types. `CLASS' The class of the user defined type, NF90_VLEN, NF90_OPAQUE, NF90_ENUM, or NF90_COMPOUND, will be copied here. Errors ====== `NF90_NOERR' No error. `NF90_EBADTYPEID' Bad typeid. `NF90_EBADFIELDID' Bad fieldid. `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 5.4.1 Set a Variable Length Array with NF90_PUT_VLEN_ELEMENT ------------------------------------------------------------ Use this to set the element of the (potentially) n-dimensional array of VLEN. That is, this sets the data in one variable length array. Usage ===== INTEGER FUNCTION NF90_PUT_VLEN_ELEMENT(INTEGER NCID, INTEGER XTYPE, CHARACTER*(*) VLEN_ELEMENT, INTEGER LEN, DATA) `NCID' The ncid of the file that contains the VLEN type. `XTYPE' The type of the VLEN. `VLEN_ELEMENT' The VLEN element to be set. `LEN' The number of entries in this array. `DATA' The data to be stored. Must match the base type of this VLEN. Errors ====== `NF90_NOERR' No error. `NF90_EBADTYPE' Can't find the typeid. `NF90_EBADID' ncid invalid. `NF90_EBADGRPID' Group ID part of ncid was invalid. Example ======= This example is from nf90_test/ftst_vars4.F. C Set up the vlen with this helper function, since F77 can't deal C with pointers. retval = nf90_put_vlen_element(ncid, vlen_typeid, vlen, & vlen_len, data1) if (retval .ne. nf90_noerr) call handle_err(retval) 5.4.2 Set a Variable Length Array with NF90_GET_VLEN_ELEMENT ------------------------------------------------------------ Use this to set the element of the (potentially) n-dimensional array of VLEN. That is, this sets the data in one variable length array. Usage ===== INTEGER FUNCTION NF90_GET_VLEN_ELEMENT(INTEGER NCID, INTEGER XTYPE, CHARACTER*(*) VLEN_ELEMENT, INTEGER LEN, DATA) `NCID' The ncid of the file that contains the VLEN type. `XTYPE' The type of the VLEN. `VLEN_ELEMENT' The VLEN element to be set. `LEN' This will be set to the number of entries in this array. `DATA' The data will be copied here. Sufficient storage must be available or bad things will happen to you. Errors ====== `NF90_NOERR' No error. `NF90_EBADTYPE' Can't find the typeid. `NF90_EBADID' ncid invalid. `NF90_EBADGRPID' Group ID part of ncid was invalid. Example ======= 5.5 Compound Types Introduction =============================== NetCDF-4 added support for compound types, which allow users to construct a new type - a combination of other types, like a C struct. Compound types are not supported in classic or 64-bit offset format files. To write data in a compound type, first use nf90_def_compound to create the type, multiple calls to nf90_insert_compound to add to the compound type, and then write data with the appropriate nf90_put_var1, nf90_put_vara, nf90_put_vars, or nf90_put_varm call. To read data written in a compound type, you must know its structure. Use the NF90_INQ_COMPOUND functions to learn about the compound type. In Fortran a character buffer must be used for the compound data. The user must read the data from within that buffer in the same way that the C compiler which compiled netCDF would store the structure. The use of compound types introduces challenges and portability issues for Fortran users. 5.5.1 Creating a Compound Type: NF90_DEF_COMPOUND ------------------------------------------------- Create a compound type. Provide an ncid, a name, and a total size (in bytes) of one element of the completed compound type. After calling this function, fill out the type with repeated calls to NF90_INSERT_COMPOUND (*note NF90_INSERT_COMPOUND::). Call NF90_INSERT_COMPOUND once for each field you wish to insert into the compound type. Note that there does not seem to be a fully portable way to read such types into structures in Fortran 90 (and there are no structures in Fortran 77). Dozens of top-notch programmers are swarming over this problem in a sub-basement of Unidata's giant underground bunker in Wyoming. Fortran users may use character buffers to read and write compound types. User are invited to try classic Fortran features such as the equivilence and the common block statment. Usage ===== function nf90_def_compound(ncid, size, name, typeid) integer, intent(in) :: ncid integer, intent(in) :: size character (len = *), intent(in) :: name integer, intent(out) :: typeid integer :: nf90_def_compound `NCID' The groupid where this compound type will be created. `SIZE' The size, in bytes, of the compound type. `NAME' The name of the new compound type. `TYPEIDP' The typeid of the new type will be placed here. Errors ====== `NF90_NOERR' No error. `NF90_EBADID' Bad group id. `NF90_ENAMEINUSE' That name is in use. Compound type names must be unique in the data file. `NF90_EMAXNAME' Name exceeds max length NF90_MAX_NAME. `NF90_EBADNAME' Name contains illegal characters. `NF90_ENOTNC4' Attempting a netCDF-4 operation on a netCDF-3 file. NetCDF-4 operations can only be performed on files defined with a create mode which includes flag NF90_NETCDF4. (*note NF90_OPEN::). `NF90_ESTRICTNC3' This file was created with the strict netcdf-3 flag, therefore netcdf-4 operations are not allowed. (*note NF90_OPEN::). `NF90_EHDFERR' An error was reported by the HDF5 layer. `NF90_EPERM' Attempt to write to a read-only file. `NF90_ENOTINDEFINE' Not in define mode. Example ======= 5.5.2 Inserting a Field into a Compound Type: NF90_INSERT_COMPOUND ------------------------------------------------------------------ Insert a named field into a compound type. Usage ===== function nf90_insert_compound(ncid, xtype, name, offset, field_typeid) integer, intent(in) :: ncid integer, intent(in) :: xtype character (len = *), intent(in) :: name integer, intent(in) :: offset integer, intent(in) :: field_typeid integer :: nf90_insert_compound `TYPEID' The typeid for this compound type, as returned by NF90_DEF_COMPOUND, or NF90_INQ_VAR. `NAME' The name of the new field. `OFFSET' Offset in byte from the beginning of the compound type for this field. `FIELD_TYPEID' The type of the field to be inserted. Errors ====== `NF90_NOERR' No error. `NF90_EBADID' Bad group id. `NF90_ENAMEINUSE' That name is in use. Field names must be unique within a compound type. `NF90_EMAXNAME' Name exceed max length NF90_MAX_NAME. `NF90_EBADNAME' Name contains illegal characters. `NF90_ENOTNC4' Attempting a netCDF-4 operation on a netCDF-3 file. NetCDF-4 operations can only be performed on files defined with a create mode which includes flag NF90_NETCDF4. (*note NF90_OPEN::). `NF90_ESTRICTNC3' This file was created with the strict netcdf-3 flag, therefore netcdf-4 operations are not allowed. (*note NF90_OPEN::). `NF90_EHDFERR' An error was reported by the HDF5 layer. `NF90_ENOTINDEFINE' Not in define mode. Example ======= 5.5.3 Inserting an Array Field into a Compound Type: NF90_INSERT_ARRAY_COMPOUND ------------------------------------------------------------------------------- Insert a named array field into a compound type. Usage ===== function nf90_insert_array_compound(ncid, xtype, name, offset, field_typeid, & ndims, dim_sizes) integer, intent(in) :: ncid integer, intent(in) :: xtype character (len = *), intent(in) :: name integer, intent(in) :: offset integer, intent(in) :: field_typeid integer, intent(in) :: ndims integer, intent(in) :: dim_sizes integer :: nf90_insert_array_compound `NCID' The ID of the file that contains the array type and the compound type. `XTYPE' The typeid for this compound type, as returned by nf90_def_compound, or nf90_inq_var. `NAME' The name of the new field. `OFFSET' Offset in byte from the beginning of the compound type for this field. `FIELD_TYPEID' The base type of the array to be inserted. `NDIMS' The number of dimensions for the array to be inserted. `DIM_SIZES' An array containing the sizes of each dimension. Errors ====== `NF90_NOERR' No error. `NF90_EBADID' Bad group id. `NF90_ENAMEINUSE' That name is in use. Field names must be unique within a compound type. `NF90_EMAXNAME' Name exceed max length NF90_MAX_NAME. `NF90_EBADNAME' Name contains illegal characters. `NF90_ENOTNC4' Attempting a netCDF-4 operation on a netCDF-3 file. NetCDF-4 operations can only be performed on files defined with a create mode which includes flag NF90_NETCDF4. (*note NF90_OPEN::). `NF90_ESTRICTNC3' This file was created with the strict netcdf-3 flag, therefore netcdf-4 operations are not allowed. (*note NF90_OPEN::). `NF90_EHDFERR' An error was reported by the HDF5 layer. `NF90_ENOTINDEFINE' Not in define mode. `NF90_ETYPEDEFINED' Attempt to change type that has already been committed. The first time the file leaves define mode, all defined types are committed, and can't be changed. If you wish to add an array to a compound type, you must do so before the compound type is committed. Example ======= 5.5.4 Learn About a Compound Type: NF90_INQ_COMPOUND ---------------------------------------------------- Get the number of fields, length in bytes, and name of a compound type. In addtion to the NF90_INQ_COMPOUND function, three additional functions are provided which get only the name, size, and number of fields. Usage ===== function nf90_inq_compound(ncid, xtype, name, size, nfields) integer, intent(in) :: ncid integer, intent(in) :: xtype character (len = *), intent(out) :: name integer, intent(out) :: size integer, intent(out) :: nfields integer :: nf90_inq_compound function nf90_inq_compound_name(ncid, xtype, name) integer, intent(in) :: ncid integer, intent(in) :: xtype character (len = *), intent(out) :: name integer :: nf90_inq_compound_name function nf90_inq_compound_size(ncid, xtype, size) integer, intent(in) :: ncid integer, intent(in) :: xtype integer, intent(out) :: size integer :: nf90_inq_compound_size function nf90_inq_compound_nfields(ncid, xtype, nfields) integer, intent(in) :: ncid integer, intent(in) :: xtype integer, intent(out) :: nfields integer :: nf90_inq_compound_nfields `NCID' The ID of any group in the file that contains the compound type. `XTYPE' The typeid for this compound type, as returned by NF90_DEF_COMPOUND, or NF90_INQ_VAR. `NAME' Character array which will get the name of the compound type. It will have a maximum length of NF90_MAX_NAME. `SIZEP' The size of the compound type in bytes will be put here. `NFIELDSP' The number of fields in the compound type will be placed here. Return Codes ============ `NF90_NOERR' No error. `NF90_EBADID' Couldn't find this ncid. `NF90_ENOTNC4' Not a netCDF-4/HDF5 file. `NF90_ESTRICTNC3' A netCDF-4/HDF5 file, but with CLASSIC_MODEL. No user defined types are allowed in the classic model. `NF90_EBADTYPE' This type not a compound type. `NF90_EBADTYPEID' Bad type id. `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 5.5.5 Learn About a Field of a Compound Type: NF90_INQ_COMPOUND_FIELD --------------------------------------------------------------------- Get information about one of the fields of a compound type. Usage ===== function nf90_inq_compound_field(ncid, xtype, fieldid, name, offset, & field_typeid, ndims, dim_sizes) integer, intent(in) :: ncid integer, intent(in) :: xtype integer, intent(in) :: fieldid character (len = *), intent(out) :: name integer, intent(out) :: offset integer, intent(out) :: field_typeid integer, intent(out) :: ndims integer, intent(out) :: dim_sizes integer :: nf90_inq_compound_field function nf90_inq_compound_fieldname(ncid, xtype, fieldid, name) integer, intent(in) :: ncid integer, intent(in) :: xtype integer, intent(in) :: fieldid character (len = *), intent(out) :: name integer :: nf90_inq_compound_fieldname function nf90_inq_compound_fieldindex(ncid, xtype, name, fieldid) integer, intent(in) :: ncid integer, intent(in) :: xtype character (len = *), intent(in) :: name integer, intent(out) :: fieldid integer :: nf90_inq_compound_fieldindex function nf90_inq_compound_fieldoffset(ncid, xtype, fieldid, offset) integer, intent(in) :: ncid integer, intent(in) :: xtype integer, intent(in) :: fieldid integer, intent(out) :: offset integer :: nf90_inq_compound_fieldoffset function nf90_inq_compound_fieldtype(ncid, xtype, fieldid, field_typeid) integer, intent(in) :: ncid integer, intent(in) :: xtype integer, intent(in) :: fieldid integer, intent(out) :: field_typeid integer :: nf90_inq_compound_fieldtype function nf90_inq_compound_fieldndims(ncid, xtype, fieldid, ndims) integer, intent(in) :: ncid integer, intent(in) :: xtype integer, intent(in) :: fieldid integer, intent(out) :: ndims integer :: nf90_inq_compound_fieldndims function nf90_inq_cmp_fielddim_sizes(ncid, xtype, fieldid, dim_sizes) integer, intent(in) :: ncid integer, intent(in) :: xtype integer, intent(in) :: fieldid integer, intent(out) :: dim_sizes integer :: nf90_inq_cmp_fielddim_sizes `NCID' The groupid where this compound type exists. `XTYPE' The typeid for this compound type, as returned by NF90_DEF_COMPOUND, or NF90_INQ_VAR. `FIELDID' A one-based index number specifying a field in the compound type. `NAME' A character array which will get the name of the field. The name will be NF90_MAX_NAME characters, at most. `OFFSETP' An integer which will get the offset of the field. `FIELD_TYPEID' An integer which will get the typeid of the field. `NDIMSP' An integer which will get the number of dimensions of the field. `DIM_SIZESP' An integer array which will get the dimension sizes of the field. Errors ====== `NF90_NOERR' No error. `NF90_EBADTYPEID' Bad type id. `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 5.6 Variable Length Array Introduction ====================================== NetCDF-4 added support for a variable length array type. This is not supported in classic or 64-bit offset files, or in netCDF-4 files which were created with the NF90_CLASSIC_MODEL flag. A variable length array is represented in C as a structure from HDF5, the nf90_vlen_t structure. It contains a len member, which contains the length of that array, and a pointer to the array. So an array of VLEN in C is an array of nc_vlen_t structures. The only way to handle this in Fortran is with a character buffer sized correctly for the platform. VLEN arrays are handled differently with respect to allocation of memory. Generally, when reading data, it is up to the user to malloc (and subsequently free) the memory needed to hold the data. It is up to the user to ensure that enough memory is allocated. With VLENs, this is impossible. The user cannot know the size of an array of VLEN until after reading the array. Therefore when reading VLEN arrays, the netCDF library will allocate the memory for the data within each VLEN. It is up to the user, however, to eventually free this memory. This is not just a matter of one call to free, with the pointer to the array of VLENs; each VLEN contains a pointer which must be freed. 5.6.1 Define a Variable Length Array (VLEN): NF90_DEF_VLEN ---------------------------------------------------------- Use this function to define a variable length array type. Usage ===== function nf90_def_vlen(ncid, name, base_typeid, xtypeid) integer, intent(in) :: ncid character (len = *), intent(in) :: name integer, intent(in) :: base_typeid integer, intent(out) :: xtypeid integer :: nf90_def_vlen `NCID' The ncid of the file to create the VLEN type in. `NAME' A name for the VLEN type. `BASE_TYPEID' The typeid of the base type of the VLEN. For example, for a VLEN of shorts, the base type is NF90_SHORT. This can be a user defined type. `XTYPEP' The typeid of the new VLEN type will be set here. Errors ====== `NF90_NOERR' No error. `NF90_EMAXNAME' NF90_MAX_NAME exceeded. `NF90_ENAMEINUSE' Name is already in use. `NF90_EBADNAME' Attribute or variable name contains illegal characters. `NF90_EBADID' ncid invalid. `NF90_EBADGRPID' Group ID part of ncid was invalid. `NF90_EINVAL' Size is invalid. `NF90_ENOMEM' Out of memory. Example ======= 5.6.2 Learning about a Variable Length Array (VLEN) Type: NF90_INQ_VLEN ----------------------------------------------------------------------- Use this type to learn about a vlen. Usage ===== function nf90_inq_vlen(ncid, xtype, name, datum_size, base_nc_type) integer, intent(in) :: ncid integer, intent(in) :: xtype character (len = *), intent(out) :: name integer, intent(out) :: datum_size integer, intent(out) :: base_nc_type integer :: nf90_inq_vlen `NCID' The ncid of the file that contains the VLEN type. `XTYPE' The type of the VLEN to inquire about. `NAME' The name of the VLEN type. The name will be NF90_MAX_NAME characters or less. `DATUM_SIZEP' A pointer to a size_t, this will get the size of one element of this vlen. `BASE_NF90_TYPEP' An integer that will get the type of the VLEN base type. (In other words, what type is this a VLEN of?) Errors ====== `NF90_NOERR' No error. `NF90_EBADTYPE' Can't find the typeid. `NF90_EBADID' ncid invalid. `NF90_EBADGRPID' Group ID part of ncid was invalid. Example ======= 5.6.3 Releasing Memory for a Variable Length Array (VLEN) Type: NF90_FREE_VLEN ------------------------------------------------------------------------------ When a VLEN is read into user memory from the file, the HDF5 library performs memory allocations for each of the variable length arrays contained within the VLEN structure. This memory must be freed by the user to avoid memory leaks. This violates the normal netCDF expectation that the user is responsible for all memory allocation. But, with VLEN arrays, the underlying HDF5 library allocates the memory for the user, and the user is responsible for deallocating that memory. Usage ===== function nf90_free_vlen(vl) character (len = *), intent(in) :: vlen integer :: nf90_free_vlen nf90_free_vlen = nf_free_vlen(vl) end function nf90_free_vlen `VL' The variable length array structure which is to be freed. Errors ====== `NF90_NOERR' No error. `NF90_EBADTYPE' Can't find the typeid. Example ======= 5.7 Opaque Type Introduction ============================ NetCDF-4 added support for the opaque type. This is not supported in classic or 64-bit offset files. The opaque type is a type which is a collection of objects of a known size. (And each object is the same size). Nothing is known to netCDF about the contents of these blobs of data, except their size in bytes, and the name of the type. To use an opaque type, first define it with *note NF90_DEF_OPAQUE::. If encountering an enum type in a new data file, use *note NF90_INQ_OPAQUE:: to learn it's name and size. 5.7.1 Creating Opaque Types: NF90_DEF_OPAQUE -------------------------------------------- Create an opaque type. Provide a size and a name. Usage ===== function nf90_def_opaque(ncid, size, name, xtype) integer, intent(in) :: ncid integer, intent(in) :: size character (len = *), intent(in) :: name integer, intent(out) :: xtype integer :: nf90_def_opaque `NCID' The groupid where the type will be created. The type may be used anywhere in the file, no matter what group it is in. `NAME' The name for this type. Must be shorter than NF90_MAX_NAME. `SIZE' The size of each opaque object. `TYPEIDP' Pointer where the new typeid for this type is returned. Use this typeid when defining variables of this type with *note NF90_DEF_VAR::. Errors ====== `NF90_NOERR' No error. `NF90_EBADTYPEID' Bad typeid. `NF90_EBADFIELDID' Bad fieldid. `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 5.7.2 Learn About an Opaque Type: NF90_INQ_OPAQUE ------------------------------------------------- Given a typeid, get the information about an opaque type. Usage ===== function nf90_inq_opaque(ncid, xtype, name, size) integer, intent(in) :: ncid integer, intent(in) :: xtype character (len = *), intent(out) :: name integer, intent(out) :: size integer :: nf90_inq_opaque `NCID' The ncid for the group containing the opaque type. `XTYPE' The typeid for this opaque type, as returned by NF90_DEF_COMPOUND, or NF90_INQ_VAR. `NAME' The name of the opaque type will be copied here. It will be NF90_MAX_NAME bytes or less. `SIZEP' The size of the opaque type will be copied here. Errors ====== `NF90_NOERR' No error. `NF90_EBADTYPEID' Bad typeid. `NF90_EBADFIELDID' Bad fieldid. `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 5.8 Enum Type Introduction ========================== NetCDF-4 added support for the enum type. This is not supported in classic or 64-bit offset files. 5.8.1 Creating a Enum Type: NF90_DEF_ENUM ----------------------------------------- Create an enum type. Provide an ncid, a name, and a base integer type. After calling this function, fill out the type with repeated calls to NF90_INSERT_ENUM (*note NF90_INSERT_ENUM::). Call NF90_INSERT_ENUM once for each value you wish to make part of the enumeration. Usage ===== function nf90_def_enum(ncid, base_typeid, name, typeid) integer, intent(in) :: ncid integer, intent(in) :: base_typeid character (len = *), intent(in) :: name integer, intent(out) :: typeid integer :: nf90_def_enum `NCID' The groupid where this compound type will be created. `BASE_TYPEID' The base integer type for this enum. Must be one of: NF90_BYTE, NF90_UBYTE, NF90_SHORT, NF90_USHORT, NF90_INT, NF90_UINT, NF90_INT64, NF90_UINT64. `NAME' The name of the new enum type. `TYPEIDP' The typeid of the new type will be placed here. Errors ====== `NF90_NOERR' No error. `NF90_EBADID' Bad group id. `NF90_ENAMEINUSE' That name is in use. Compound type names must be unique in the data file. `NF90_EMAXNAME' Name exceeds max length NF90_MAX_NAME. `NF90_EBADNAME' Name contains illegal characters. `NF90_ENOTNC4' Attempting a netCDF-4 operation on a netCDF-3 file. NetCDF-4 operations can only be performed on files defined with a create mode which includes flag NF90_NETCDF4. (*note NF90_OPEN::). `NF90_ESTRICTNC3' This file was created with the strict netcdf-3 flag, therefore netcdf-4 operations are not allowed. (*note NF90_OPEN::). `NF90_EHDFERR' An error was reported by the HDF5 layer. `NF90_EPERM' Attempt to write to a read-only file. `NF90_ENOTINDEFINE' Not in define mode. Example ======= 5.8.2 Inserting a Field into a Enum Type: NF90_INSERT_ENUM ---------------------------------------------------------- Insert a named member into a enum type. Usage ===== function nf90_insert_enum(ncid, xtype, name, value) integer, intent(in) :: ncid integer, intent(in) :: xtype character (len = *), intent(in) :: name integer, intent(in) :: value integer :: nf90_insert_enum `NCID' The ncid of the group which contains the type. `TYPEID' The typeid for this enum type, as returned by nf90_def_enum, or nf90_inq_var. `IDENTIFIER' The identifier of the new member. `VALUE' The value that is to be associated with this member. Errors ====== `NF90_NOERR' No error. `NF90_EBADID' Bad group id. `NF90_ENAMEINUSE' That name is in use. Field names must be unique within a enum type. `NF90_EMAXNAME' Name exceed max length NF90_MAX_NAME. `NF90_EBADNAME' Name contains illegal characters. `NF90_ENOTNC4' Attempting a netCDF-4 operation on a netCDF-3 file. NetCDF-4 operations can only be performed on files defined with a create mode which includes flag NF90_NETCDF4. (*note NF90_OPEN::). `NF90_ESTRICTNC3' This file was created with the strict netcdf-3 flag, therefore netcdf-4 operations are not allowed. (*note NF90_OPEN::). `NF90_EHDFERR' An error was reported by the HDF5 layer. `NF90_ENOTINDEFINE' Not in define mode. Example ======= 5.8.3 Learn About a Enum Type: NF90_INQ_ENUM -------------------------------------------- Get information about a user-defined enumeration type. Usage ===== function nf90_inq_enum(ncid, xtype, name, base_nc_type, base_size, num_members) integer, intent(in) :: ncid integer, intent(in) :: xtype character (len = *), intent(out) :: name integer, intent(out) :: base_nc_type integer, intent(out) :: base_size integer, intent(out) :: num_members integer :: nf90_inq_enum `NCID' The group ID of the group which holds the enum type. `XTYPE' The typeid for this enum type, as returned by NF90_DEF_ENUM, or NF90_INQ_VAR. `NAME' Character array which will get the name. It will have a maximum length of NF90_MAX_NAME. `BASE_NF90_TYPE' An integer which will get the base integer type of this enum. `BASE_SIZE' An integer which will get the size (in bytes) of the base integer type of this enum. `NUM_MEMBERS' An integer which will get the number of members defined for this enumeration type. Errors ====== `NF90_NOERR' No error. `NF90_EBADTYPEID' Bad type id. `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 5.8.4 Learn the Name of a Enum Type: nf90_inq_enum_member --------------------------------------------------------- Get information about a member of an enum type. Usage ===== function nf90_inq_enum_member(ncid, xtype, idx, name, value) integer, intent(in) :: ncid integer, intent(in) :: xtype integer, intent(in) :: idx character (len = *), intent(out) :: name integer, intent(in) :: value integer :: nf90_inq_enum_member `NCID' The groupid where this enum type exists. `XTYPE' The typeid for this enum type. `IDX' The one-based index number for the member of interest. `NAME' A character array which will get the name of the member. It will have a maximum length of NF90_MAX_NAME. `VALUE' An integer that will get the value associated with this member. Errors ====== `NF90_NOERR' No error. `NF90_EBADTYPEID' Bad type id. `NF90_EHDFERR' An error was reported by the HDF5 layer. Example ======= 5.8.5 Learn the Name of a Enum Type: NF90_INQ_ENUM_IDENT -------------------------------------------------------- Get the name which is associated with an enum member value. This is similar to NF90_INQ_ENUM_MEMBER, but instead of using the index of the member, you use the value of the member. Usage ===== function nf90_inq_enum_ident(ncid, xtype, value, idx) integer, intent(in) :: ncid integer, intent(in) :: xtype integer, intent(in) :: value integer, intent(out) :: idx integer :: nf90_inq_enum_ident `NCID' The groupid where this enum type exists. `XTYPE' The typeid for this enum type. `VALUE' The value for which an identifier is sought. `IDENTIFIER' A character array that will get the identifier. It will have a maximum length of NF90_MAX_NAME. Return Code =========== `NF90_NOERR' No error. `NF90_EBADTYPEID' Bad type id, or not an enum type. `NF90_EHDFERR' An error was reported by the HDF5 layer. `NF90_EINVAL' The value was not found in the enum. Example ======= 6 Variables *********** 6.1 Variables Introduction ========================== Variables for a netCDF dataset are defined when the dataset is created, while the netCDF dataset is in define mode. Other variables may be added later by reentering define mode. A netCDF variable has a name, a type, and a shape, which are specified when it is defined. A variable may also have values, which are established later in data mode. Ordinarily, the name, type, and shape are fixed when the variable is first defined. The name may be changed, but the type and shape of a variable cannot be changed. However, a variable defined in terms of the unlimited dimension can grow without bound in that dimension. A netCDF variable in an open netCDF dataset is referred to by a small integer called a variable ID. Variable IDs reflect the order in which variables were defined within a netCDF dataset. Variable IDs are 1, 2, 3,..., in the order in which the variables were defined. A function is available for getting the variable ID from the variable name and vice-versa. Attributes (see *note Attributes::) may be associated with a variable to specify such properties as units. Operations supported on variables are: * Create a variable, given its name, data type, and shape. * Get a variable ID from its name. * Get a variable's name, data type, shape, and number of attributes from its ID. * Put a data value into a variable, given variable ID, indices, and value. * Put an array of values into a variable, given variable ID, corner indices, edge lengths, and a block of values. * Put a subsampled or mapped array-section of values into a variable, given variable ID, corner indices, edge lengths, stride vector, index mapping vector, and a block of values. * Get a data value from a variable, given variable ID and indices. * Get an array of values from a variable, given variable ID, corner indices, and edge lengths. * Get a subsampled or mapped array-section of values from a variable, given variable ID, corner indices, edge lengths, stride vector, and index mapping vector. * Rename a variable. 6.2 Language Types Corresponding to netCDF external data types ============================================================== The following table gives the netCDF external data types and the corresponding type constants for defining variables in the FORTRAN interface: Type FORTRAN API Mnemonic Bits byte NF90_BYTE 8 char NF90_CHAR 8 short NF90_SHORT 16 int NF90_INT 32 float NF90_FLOAT 32 double NF90_DOUBLE 64 The first column gives the netCDF external data type, which is the same as the CDL data type. The next column gives the corresponding Fortran 90 parameter for use in netCDF functions (the parameters are defined in the netCDF Fortran 90 module netcdf.f90). The last column gives the number of bits used in the external representation of values of the corresponding type. Note that there are no netCDF types corresponding to 64-bit integers or to characters wider than 8 bits in the current version of the netCDF library. 6.3 Create a Variable: `NF90_DEF_VAR' ===================================== The function NF90_DEF_VAR adds a new variable to an open netCDF dataset in define mode. It returns (as an argument) a variable ID, given the netCDF ID, the variable name, the variable type, the number of dimensions, and a list of the dimension IDs. Usage ===== function nf90_def_var(ncid, name, xtype, dimids, varid) integer, intent( in) :: ncid character (len = *), intent( in) :: name integer, intent( in) :: xtype integer, dimension(:), optional, intent( in) :: dimids integer, intent(out) :: varid integer :: nf90_def_var `ncid' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `name' Variable name. Must begin with an alphabetic character, followed by zero or more alphanumeric characters including the underscore ('_'). Case is significant. `xtype' One of the set of predefined netCDF external data types. The type of this parameter, NF90_TYPE, is defined in the netCDF header file. The valid netCDF external data types are NF90_BYTE, NF90_CHAR, NF90_SHORT, NF90_INT, NF90_FLOAT, and NF90_DOUBLE. If the file is a NetCDF-4/HDF5 file, the additional types NF_UBYTE, NF_USHORT, NF_UINT, NF_INT64, NF_UINT64, and NF_STRING may be used, as well as a user defined type ID. `dimids' Vector of dimension IDs corresponding to the variable dimensions. For example, a vector of 2 dimension IDs specifies a 2-dimensional matrix. If an integer is passed for this parameter, a 1-D variable is created. If this parameter is not passed it means the variable is a scalar with no dimensions. For classic data model files, if the ID of the unlimited dimension is included, it must be first. In expanded model netCDF4/HDF5 files, there may be any number of unlimited dimensions, and they may be used in any element of the dimids array. This argument is optional, and if absent specifies a scalar with no dimensions. `varid' Returned variable ID. Errors ====== NF90_DEF_VAR returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible causes of errors include: * The netCDF dataset is not in define mode. * The specified variable name is the name of another existing variable. * The specified type is not a valid netCDF type. * The specified number of dimensions is negative or more than the constant NF90_MAX_VAR_DIMS, the maximum number of dimensions permitted for a netCDF variable. * One or more of the dimension IDs in the list of dimensions is not a valid dimension ID for the netCDF dataset. * The number of variables would exceed the constant NF90_MAX_VARS, the maximum number of variables permitted in a netCDF dataset. * The specified netCDF ID does not refer to an open netCDF dataset. Example ======= Here is an example using NF90_DEF_VAR to create a variable named rh of type double with three dimensions, time, lat, and lon in a new netCDF dataset named foo.nc: use netcdf implicit none integer :: status, ncid integer :: LonDimId, LatDimId, TimeDimId integer :: RhVarId ... status = nf90_create("foo.nc", nf90_NoClobber, ncid) if(status /= nf90_NoErr) call handle_error(status) ... ! Define the dimensions status = nf90_def_dim(ncid, "lat", 5, LatDimId) if(status /= nf90_NoErr) call handle_error(status) status = nf90_def_dim(ncid, "lon", 10, LonDimId) if(status /= nf90_NoErr) call handle_error(status) status = nf90_def_dim(ncid, "time", nf90_unlimited, TimeDimId) if(status /= nf90_NoErr) call handle_error(status) ... ! Define the variable status = nf90_def_var(ncid, "rh", nf90_double, & (/ LonDimId, LatDimID, TimeDimID /), RhVarId) if(status /= nf90_NoErr) call handle_error(status) 6.4 Define Chunking Parameters for a Variable: `NF90_DEF_VAR_CHUNKING' ====================================================================== The function NF90_DEF_VAR_CHUNKING sets the chunking parameters for a variable in a netCDF-4 file. The total size the chunk must be less than 4 GiB. That is, the product of all chunksizes and the size of the data (or the size of nc_vlen_t for VLEN types) must be less than 4 GiB. This function must be called after the variable is defined, but before nf_enddef is called. Usage ===== function nf90_def_var_chunking(ncid, varid, contiguous, chunksizes) integer, intent(in) :: ncid integer, intent(in) :: varid integer, intent(in) :: contiguous integer, dimension(:), intent(in) :: chunksizes integer :: nf90_def_var_chunking `ncid' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `varid' Variable ID. `contiguous' If non-zero, then contiguous storage is used for this variable. Variables with one or more unlimited dimensions cannot use contiguous storage. If contiguous storage is turned on, the chunksizes parameter is ignored. `chunksizes' An array of chunk sizes. The array must have the one chunksize for each dimension in the variable. If the contiguous parameter is set, then the chunksizes parameter is ignored. Errors ====== NF90_DEF_VAR_CHUNKING returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible return codes include: `NF90_NOERR' No error. `NF90_BADID' Bad ncid. `NF90_ENOTNC4' Not a netCDF-4 file. `NF90_ENOTVAR' Can't find this variable. `NF90_ELATEDEF' This variable has already been the subject of a NF90_ENDDEF call. In netCDF-4 files NF90_ENDDEF will be called automatically for any data read or write. Once enddef has been called, it is impossible to set the chunking for a variable. `NF90_ENOTINDEFINE' Not in define mode. This is returned for netCDF classic or 64-bit offset files, or for netCDF-4 files, when they were been created with NF90_STRICT_NC3 flag. (*note NF90_CREATE::). `NF90_ESTRICTNC3' Trying to create a var some place other than the root group in a netCDF file with NF90_STRICT_NC3 turned on. Example ======= 6.5 Learn About Chunking Parameters for a Variable: `NF90_INQ_VAR_CHUNKING' =========================================================================== The function NF90_INQ_VAR_CHUNKING returns the chunking settings for a variable in a netCDF-4 file. Usage ===== function nf90_inq_var_chunking(ncid, varid, contiguous, chunksizes) integer, intent(in) :: ncid integer, intent(in) :: varid integer, intent(out) :: contiguous integer, dimension(:), intent(out) :: chunksizes integer :: nf90_inq_var_chunking `NCID' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `VARID' Variable ID. `CONTIGUOUS' Set to 1 if this variable uses contiguous storage, 0 if it used chunked storage. `CHUNKSIZES' An array of chunk sizes. The array must have the one element for each dimension in the variable. Errors ====== NF90_INQ_VAR_CHUNKING returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible return codes include: `NF90_NOERR' No error. `NF90_BADID' Bad ncid. `NF90_ENOTNC4' Not a netCDF-4 file. `NF90_ENOTVAR' Can't find this variable. Example ======= 6.6 Define Compression Parameters for a Variable: `NF90_DEF_VAR_DEFLATE' ======================================================================== The function NF90_DEF_VAR_DEFLATE sets the deflate parameters for a variable in a netCDF-4 file. This function must be called after the variable is defined, but before NF90_ENDDEF is called. Usage ===== function nf90_def_var_deflate(ncid, varid, shuffle, deflate, deflate_level) integer, intent(in) :: ncid integer, intent(in) :: varid integer, intent(in) :: shuffle integer, intent(in) :: deflate integer, intent(in) :: deflate_level integer :: nf90_def_var_deflate `NCID' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `VARID' Variable ID. `SHUFFLE' If non-zero, turn on the shuffle filter. `DEFLATE' If non-zero, turn on the deflate filter at the level specified by the deflate_level parameter. `DEFLATE_LEVEL' If the deflate parameter is non-zero, set the deflate level to this value. Must be between 0 and 9. Errors ====== NF90_DEF_VAR_DEFLATE returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible return codes include: `NF90_NOERR' No error. `NF90_BADID' Bad ncid. `NF90_ENOTNC4' Not a netCDF-4 file. `NF90_ENOTVAR' Can't find this variable. `NF90_ELATEDEF' This variable has already been the subject of a NF90_ENDDEF call. In netCDF-4 files NF90_ENDDEF will be called automatically for any data read or write. Once enddef has been called, it is impossible to set the deflate for a variable. `NF90_ENOTINDEFINE' Not in define mode. This is returned for netCDF classic or 64-bit offset files, or for netCDF-4 files, when they were been created with NF90_STRICT_NC3 flag. (*note NF90_CREATE::). `NF90_EPERM' Attempt to create object in read-only file. `NF90_EINVAL' Invalid deflate_level. The deflate level must be between 0 and 9, inclusive. Example ======= 6.7 Learn About Deflate Parameters for a Variable: `NF90_INQ_VAR_DEFLATE' ========================================================================= The function NF90_INQ_VAR_DEFLATE returns the deflate settings for a variable in a netCDF-4 file. Usage ===== function nf90_inq_var_deflate(ncid, varid, shuffle, deflate, deflate_level) integer, intent(in) :: ncid integer, intent(in) :: varid integer, intent(out) :: shuffle integer, intent(out) :: deflate integer, intent(out) :: deflate_level integer :: nf90_inq_var_deflate `NCID' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `VARID' Variable ID. `SHUFFLE' NF90_INQ_VAR_DEFLATE will set this to a 1 if the shuffle filter is turned on for this variable, and a 0 otherwise. `DEFLATE' NF90_INQ_VAR_DEFLATE will set this to a 1 if the deflate filter is turned on for this variable, and a 0 otherwise. `DEFLATE_LEVEL' NF90_INQ_VAR_DEFLATE function will write the deflate_level here, if deflate is in use. Errors ====== NF90_INQ_VAR_DEFLATE returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible return codes include: `NF90_NOERR' No error. `NF90_BADID' Bad ncid. `NF90_ENOTNC4' Not a netCDF-4 file. `NF90_ENOTVAR' Can't find this variable. Example ======= 6.8 Define Fletcher32 Parameters for a Variable: `NF90_DEF_VAR_FLETCHER32' ========================================================================== The function NF90_DEF_VAR_FLETCHER32 sets the fletcher32 parameters for a variable in a netCDF-4 file. This function must be called after the variable is defined, but before NF90_ENDDEF is called. Usage ===== function nf90_def_var_fletcher32(ncid, varid, fletcher32) integer, intent(in) :: ncid integer, intent(in) :: varid integer, intent(in) :: fletcher32 integer :: nf90_def_var_fletcher32 `NCID' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `VARID' Variable ID. `FLETCHER32' If this is non-zero, fletcher32 checksums will be turned on for this variable. Errors ====== NF90_DEF_VAR_FLETCHER32 returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible return codes include: `NF90_NOERR' No error. `NF90_BADID' Bad ncid. `NF90_ENOTNC4' Not a netCDF-4 file. `NF90_ENOTVAR' Can't find this variable. `NF90_ELATEDEF' This variable has already been the subject of a NF90_ENDDEF call. In netCDF-4 files NF90_ENDDEF will be called automatically for any data read or write. Once enddef has been called, it is impossible to set the fletcher32 for a variable. `NF90_ENOTINDEFINE' Not in define mode. This is returned for netCDF classic or 64-bit offset files, or for netCDF-4 files, when they were been created with NF90_STRICT_NC3 flag. (*note NF90_CREATE::). `NF90_EPERM' Attempt to create object in read-only file. Example ======= 6.9 Learn About Fletcher32 Parameters for a Variable: `NF90_INQ_VAR_FLETCHER32' =============================================================================== The function NF90_INQ_VAR_FLETCHER32 returns the fletcher32 settings for a variable in a netCDF-4 file. Usage ===== function nf90_inq_var_fletcher32(ncid, varid, fletcher32) integer, intent(in) :: ncid integer, intent(in) :: varid integer, intent(out) :: fletcher32 integer :: nf90_inq_var_fletcher32 `NCID' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `VARID' Variable ID. `FLETCHER32' NF90_INQ_VAR_FLETCHER32 will set this to 1 if the fletcher32 filter is turned on for this variable, and 0 if it is not. Errors ====== NF90_INQ_VAR_FLETCHER32 returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible return codes include: `NF90_NOERR' No error. `NF90_BADID' Bad ncid. `NF90_ENOTNC4' Not a netCDF-4 file. `NF90_ENOTVAR' Can't find this variable. Example ======= 6.10 Define Endianness of a Variable: `NF90_DEF_VAR_ENDIAN' =========================================================== The function NF90_DEF_VAR_ENDIAN sets the endianness for a variable in a netCDF-4 file. This function must be called after the variable is defined, but before NF90_ENDDEF is called. By default, netCDF-4 variables are in native endianness. That is, they are big-endian on a big-endian machine, and little-endian on a little endian machine. In some cases a user might wish to change from native endianness to either big or little-endianness. This function allows them to do that. Usage ===== function nf90_def_var_endian(ncid, varid, endian) integer, intent(in) :: ncid integer, intent(in) :: varid integer, intent(in) :: endian integer :: nf90_def_var_endian `NCID' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `VARID' Variable ID. `ENDIAN' Set to NF90_ENDIAN_NATIVE for native endianness. (This is the default). Set to NF90_ENDIAN_LITTLE for little endian, or NF90_ENDIAN_BIG for big endian. Errors ====== NF90_DEF_VAR_ENDIAN returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible return codes include: `NF90_NOERR' No error. `NF90_BADID' Bad ncid. `NF90_ENOTNC4' Not a netCDF-4 file. `NF90_ENOTVAR' Can't find this variable. `NF90_ELATEDEF' This variable has already been the subject of a NF90_ENDDEF call. In netCDF-4 files NF90_ENDDEF will be called automatically for any data read or write. Once enddef has been called, it is impossible to set the endianness of a variable. `NF90_ENOTINDEFINE' Not in define mode. This is returned for netCDF classic or 64-bit offset files, or for netCDF-4 files, when they were been created with NF90_STRICT_NC3 flag, and the file is not in define mode. (*note NF90_CREATE::). `NF90_EPERM' Attempt to create object in read-only file. Example ======= 6.11 Learn About Endian Parameters for a Variable: `NF90_INQ_VAR_ENDIAN' ======================================================================== The function NF90_INQ_VAR_ENDIAN returns the endianness settings for a variable in a netCDF-4 file. Usage ===== function nf90_inq_var_endian(ncid, varid, endian) integer, intent(in) :: ncid integer, intent(in) :: varid integer, intent(out) :: endian integer :: nf90_inq_var_endian `NCID' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `VARID' Variable ID. `ENDIAN' NF90_INQ_VAR_ENDIAN will set this to NF90_ENDIAN_LITTLE if this variable is stored in little-endian format, NF90_ENDIAN_BIG if it is stored in big-endian format, and NF90_ENDIAN_NATIVE if the endianness is not set, and the variable is not created yet. Errors ====== NF90_INQ_VAR_ENDIAN returns the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible return codes include: `NF90_NOERR' No error. `NF90_BADID' Bad ncid. `NF90_ENOTNC4' Not a netCDF-4 file. `NF90_ENOTVAR' Can't find this variable. Example ======= 6.12 Get Information about a Variable from Its ID: NF90_INQUIRE_VARIABLE ======================================================================== NF90_INQUIRE_VARIABLE returns information about a netCDF variable given its ID. Information about a variable includes its name, type, number of dimensions, a list of dimension IDs describing the shape of the variable, and the number of variable attributes that have been assigned to the variable. Usage ===== function nf90_inquire_variable(ncid, varid, name, xtype, ndims, dimids, nAtts) integer, intent( in) :: ncid, varid character (len = *), optional, intent(out) :: name integer, optional, intent(out) :: xtype, ndims integer, dimension(*), optional, intent(out) :: dimids integer, optional, intent(out) :: nAtts integer :: nf90_inquire_variable `ncid' NetCDF ID, from a previous call to NF90_OPEN or NF90_CREATE. `varid' Variable ID. `name' Returned variable name. The caller must allocate space for the returned name. The maximum possible length, in characters, of a variable name is given by the predefined constant NF90_MAX_NAME. `xtype' Returned variable type, one of the set of predefined netCDF external data types. The type of this parameter, NF90_TYPE, is defined in the netCDF header file. The valid netCDF external data types are NF90_BYTE, NF90_CHAR, NF90_SHORT, NF90_INT, NF90_FLOAT, AND NF90_DOUBLE. `ndims' Returned number of dimensions the variable was defined as using. For example, 2 indicates a matrix, 1 indicates a vector, and 0 means the variable is a scalar with no dimensions. `dimids' Returned vector of *ndimsp dimension IDs corresponding to the variable dimensions. The caller must allocate enough space for a vector of at least *ndimsp integers to be returned. The maximum possible number of dimensions for a variable is given by the predefined constant NF90_MAX_VAR_DIMS. `natts' Returned number of variable attributes assigned to this variable. These functions return the value NF90_NOERR if no errors occurred. Otherwise, the returned status indicates an error. Possible causes of errors include: * The variable ID is invalid for the specified netCDF dataset. * The specified netCDF ID does not refer to an open netCDF dataset. Example ======= Here is an example using NF90_INQ_VAR to find out about a variable named rh in an existing netCDF dataset named foo.nc: use netcdf implicit none integer :: status, ncid, & RhVarId & numDims, numAtts integer, dimension(nf90_max_var_dims) :: rhDimIds ... status = nf90_open("foo.nc", nf90_NoWrite, ncid) if(status /= nf90_NoErr) call handle_error(status) ... status = nf90_i