cf.read

cf.read(files, verbose=False, ignore_read_error=False, aggregate=True, umversion=4.5, squeeze=False, unsqueeze=False, fmt=None, select=None, select_options={})[source]

Read fields from files into cf.Field objects.

A file may be on disk or on a OPeNDAP server.

Any amount of any combination of CF-netCDF, CFA-netCDF and Met Office (UK) PP format files may be read.

netCDF files
  • The netCDF variable names of the field and its components are stored in their ncvar attributes.
  • Fields referenced within coordinate reference or ancillary variables objects are not included in the returned list of fields.
PP files
  • If any PP files are read then the aggregate option 'relaxed_units' is set to True for all input files.
  • STASH code to standard conversion uses the table in cf/etc/STASH_to_CF.txt.
Files on OPeNDAP servers
  • All files on OPeNDAP servers are assumed to be netCDF files.

See also

cf.write

Parameters:

files : (arbitrarily nested sequence of) str

A string or arbitrarily nested sequence of strings giving the file names or OPenDAP URLs from which to read fields. Various type of expansion are applied to the file names:

Expansion Description
Tilde An initial component of ~ or ~user is replaced by that user‘s home directory.
Environment variable Substrings of the form $name or ${name} are replaced by the value of environment variable name.
Pathname A string containing UNIX file name metacharacters as understood by the glob module is replaced by the list of matching file names. This type of expansion is ignored for OPenDAP URLs.

Where more than one type of expansion is used in the same string, they are applied in the order given in the above table.

Example: If the environment variable MYSELF has been set to the “david”, then '~$MYSELF/*.nc' is equivalent to '~david/*.nc', which will read all netCDF files in the user david’s home directory.

verbose : bool, optional

If True then print information to stdout.

umversion : int or float, optional

For PP format files only, the Unified Model (UM) version to be used when decoding the PP header. Valid versions are, for example, 4.2, '6.6.3' and '8.2'. The default version is 4.5. The version is ignored if it can be inferred from the PP headers, which will generally be the case for files created at versions 5.3 and later. Note that the PP header can not encode tertiary version elements (such as the 3 in '6.6.3'), so it may be necessary to provide a UM version in such cases.

Ignored for any non-PP input files.

ignore_read_error : bool, optional

If True then ignore any file which raises an IOError whilst being read, as would be the case for an empty file, unknown file format, etc. By default the IOError is raised.

fmt : str, optional

Only read files of the given format, ignoring all other files. Valid formats are 'NETCDF' for CF-netCDF files, 'CFA' for CFA-netCDF files and 'PP' for PP files and ‘FF’ for UM fields files. By default files of any of these formats are read. default files of any of these formats are read.

aggregate : bool or dict, optional

If True (the default) or a (possibly empty) dictionary then aggregate the fields read in from all input files into as few fields as possible using the CF aggregation rules. If aggregate is a dictionary then it is passed as keyword arguments to the cf.aggregate function. If False then the fields are not aggregated.

squeeze : bool, optional

If True then remove size 1 axes from each field’s data array.

unsqueeze : bool, optional

If True then insert size 1 axes from each field’s domain into its data array.

select, select_options : optional

Only return fields which satisfy the given conditions on their property values. Only fields which, prior to any aggregation, satisfy f.match(select, **select_options) == True are returned. See cf.Field.match for details.

Returns:
out : cf.FieldList

A list of fields.

Examples:
>>> f = cf.read('file*.nc')
>>> f
[<CF Field: pmsl(30, 24)>,
 <CF Field: z-squared(17, 30, 24)>,
 <CF Field: temperature(17, 30, 24)>,
 <CF Field: temperature_wind(17, 29, 24)>]
>>> cf.read('file*.nc')[0:2]
[<CF Field: pmsl(30, 24)>,
 <CF Field: z-squared(17, 30, 24)>]
>>> cf.read('file*.nc')[-1]
<CF Field: temperature_wind(17, 29, 24)>
>>> cf.read('file*.nc', prop={'units': 'K'})
[<CF Field: temperature(17, 30, 24)>,
 <CF Field: temperature_wind(17, 29, 24)>]
>>> cf.read('file*.nc', attr={'ncvar': 'ta'})
[<CF Field: temperature(17, 30, 24)>]
>>> cf.read('file*.nc', prop={'standard_name': '.*pmsl*', 'units':'K|Pa'})[0]
<CF Field: pmsl(30, 24)>
>>> cf.read('file*.nc', prop={'units':['K', 'Pa']})
[<CF Field: pmsl(30, 24)>,
 <CF Field: temperature(17, 30, 24)>,
 <CF Field: temperature_wind(17, 29, 24)>]

Previous topic

cf.pickle

Next topic

cf.write

This Page