NCAR-RAF netCDF Conventions

Version 1.3

National Center for Atmospheric Research

Research Aviation Facility

NOAA AOC (GIV)

University of Wyoming (King Air)

Table of Contents

Appendices

Introduction

The goal is to support synchronous time series including the following properties:

RAF's primary datasets consist of synchronous time series data. Data are written to the netCDF file one record per second, even for data at higher rates.

Data acquisition of numerous instruments at high-rate (from 1 to 5K samples per second) for 10 to 11 hours can produce large data files. Typical, raw packed data files from the NCAR C-130 are 600MB, and unpacked with no data reduction would be twice that size. Files reduced to 1 sample per second (sps) are typically 70-100MB.

1.0.1 Vector data

Instruments which measure particle size distributions produce histograms or vector data. These instruments include the suite of PMS-1D probes, the RDMA and quantitative data reduction from the PMS-2D probes.

Note: For all size-distributions (PMS1D, PMS2D, RDMA) the zeroth or first (depending on how you count) bin is undefined, do not use this value. Attributes for the various size distribution variables should include FirstBin and LastBin.

1.0.2 Multi-rate files

RAF's standard data product is a netCDF file where all measurements are output at 1 sps. Certain projects require high-rate data. RAF has standardized at 25 sps for high-rate data. In general we don't want to put out variables sampled at 1 sps at 25 sps, as there is no information in the data, and it would take up excess space.

There is also an occasional need to look at data at the sample rate.

NetCDF File Component Descriptions

1 Filename

RAF uses a naming convention consisting of the 3 digit project number, followed by the flight type, followed by a 2 digit flight number, followed by ".nc".

e.g.

217tf01.nc would be test flight 1 from project 217.
131ff01.nc would be ferry flight 1 from project 131.
803rf17h.nc would be high-rate data from research flight 17 of project 803.

Flight types used by RAF are "ff", "tf", and "rf".

2 Data types

All data are stored using type float, except for base_time which is of type int.

3 Naming conventions

Variable, dimension and attribute names should begin with a letter and be composed of letters, digits, and underscores.

3.1 Variable names

The underscore has special meaning for variable names. Anything following an underscore is considered a suffix for grouping variables from the same instrument or probe. There is no standard for suffix names, but RAF regularly uses about a dozen specific ones for describing an instrument's mounting location on the aircraft.

e.g.,

AFSSP_RPO, CFSSP_RPO, PLWC_RPO, CONCF_RPO all belong to the same intrument.

In our files RPO stands for "Right Pod Outboard".

3.2 Dimension names

Dimension names we use are: although I don't believe it matters what they are called.

4 Dimensions

4.1 First dimension - Time

Time is the unlimited dimension, and every variable (except base_time) has the Time dimension. No more than one record is to be written per second; all high-rate data are to be contained "per record", though data averaged down slower than 1 sps may be written (e.g., when you have 1 sample every 10 seconds).

4.2 Second dimension - data rate

The second dimension, if present, represents the number of samples per second. The size of this dimension will vary according to the output rate of the variable (e.g., RAF uses 1, 5, 10, 25, 50, 250, and 1000). If the dimension does not exist, then there is only one sample for the second.

4.3 Third dimension - vector length

The third dimension, if present, is the vector length of the data. (This will be present only for variables with vector data). Only PMS-1D probes are producing vector data at this time.

A new dimension has been added called Async. Async is a triplet containing the location of the asynchronous data in a separate binary "heap" file. The first value contains the offset in bytes into the file of the first record for that second, the second number tells the record length, and the third is the number of records for that particular second. Async will always be in the 3rd (or Vector) dimension, e.g., float 2D-P_LPI(Time, sps1, Async). Async was added to facilitate future capabilities of being able to provide satellite images, PMS-2D data, and other asynchronous products into one "heap" file. (Note: I have yet to use this feature, so implementation of file name, etc., are incomplete).

5 Variables

http://www.cfconventions.org/.

5.1 Time

Time follows the recomendations set forth in the CF conventions section 4.4. 
        int Time(Time) ;
                Time:long_name = "time of measurement" ;
                Time:standard_name = "time" ;
                Time:units = "seconds since 2004-08-02 16:46:11 +0000" ;
                Time:strptime_format = "seconds since %F %T %z" ;

Pre version 1.3 time was represented as follows, it is still included for backwards compatibility but is being deprecated.

Time is stored as a base_time and time_offset pair. base_time is standard Unix time (seconds since Jan 1, 1970), and is typically the time stamp of the first record in the file. base_time has only one value, and does not contain the time dimension. time_offset is added to base_time to give the time of the current record. 

        int base_time ;
                base_time:units = "seconds since 1970-01-01 00:00:00 +0000" ;
                base_time:long_name = "Start time of data recording." ;
        float time_offset(Time) ;
                time_offset:long_name = "Seconds since base_time." ;
                time_offset:units = "seconds since 2004-08-02 16:46:11 +0000" ;

5.2 Missing data

At this time we use the _FillValue attribute.

Prior to 01/27/2003 we used missing_value as the missing value. missing_value (pre 1998 spelled MissingValue) is being deprecated.

6 Attributes

6.1 Global attributes

// global attributes:
                :Source = "NCAR Research Aviation Facility" ;
                :Address = "P.O. Box 3000, Boulder, CO 80307-3000" ;
                :Phone = "(303) 497-1030" ;
                :Conventions = "NCAR-RAF/nimbus" ;
                :ConventionsURL = "http://www.eol.ucar.edu/raf/Software/netCDF.html" ;
                :ConventionsVersion = "1.3" ;
                :ProcessorRevision = "3213" ;
                :ProcessorURL = "http://svn/svn/raf/trunk/nimbus" ;
                :DateProcessed = "2006-02-17 16:57:57 -0700" ;
                :ProjectName = "TREX" ;
                :Aircraft = "N677F" ;
                :ProjectNumber = "503" ;
                :FlightNumber = "tf01" ;
                :FlightDate = "02/17/2006" ;
                :coordinates = "GGLON GGLAT GGALT Time" ;
                :wind_field = "WSC WDC WIC" ;
                :landmarks = "39.9083 -105.116 BJC,40.25 -103.799 FtMrg,43.5838 -96.7346 FSD,46.8392 -92.1897 DLH,32.9048 -97.0459 DFW,36.0847 -115.151 LAS,40.7858 -111.979 SLC,33.4358 -112.013 PHX,45.5894 -122.592 PDX,33.9436 -118.403 LAX" ;
                :TimeInterval = "00:03:16-01:29:43" ;

coordinates is required starting with version 1.3. These are the aircrafts "best" position variables.
TimeInterval is required by ncplot and ncpp, but I hope to get away from that.

6.2 Variable attributes

All variables contain the basic attributes units, long_name, and _FillValue (missing_value is being deprecated). In addition we provide some other attributes; not all are used for every variable:

standard_name we have adopted the CF CF conventions section 3.3 use of standard_name where applicable.

SampledRate is the rate (sps) at which the variable was sampled onboard the aircraft. This does not apply to derived variables.

DataQuality attempts to inform the user of the quality of these data. Some values are Bad, Preliminary and Good

Category Since a file can contain hundreds of variables with not very descriptive names, this is provided to help create sub-lists of variables. Comma separated list.

CalibrationCoefficients are the values used to produce engineering units from a measurement's DC voltage. It is used by the analog/digital group. These values have already been applied! They are present for documentation.

Dependencies are the input variables that were used to produce this derived value. Like CalibrationCoefficients , these are present for documentation.

For processing purposes, variables in our files are organized into 3 sorted lists. The first sorted list consists of Analog or Digitally-sampled variables. These only have calibration coefficients applied to them. The second group is raw data from "block" probes (e.g., PMS1D, IRS, GPS, and a few from group 1 that need special processing). The third group consists of purely derived variables. I have included one example from each of the 3 groups (which also show different dimension schemes): 

        float ADIFR(Time) ;
                ADIFR:units = "mbar" ;
                ADIFR:long_name = "Vertical Differential Pressure, Radome" ;
                ADIFR:Category = "Uncorr\'d Raw" ;
                ADIFR:_FillValue = -32767.f ;
                ADIFR:SampledRate = 25 ;
                ADIFR:DataQuality = "Preliminary" ;
                ADIFR:CalibrationCoefficients = -0.0115f, 5.1637f, 0.001f ;

        float AFSSP_RPC(Time, sps1, Vector16) ;
                AFSSP_RPC:units = "count" ;
                AFSSP_RPC:long_name = "FSSP-100 Raw Accumulation (per cell)" ;
                AFSSP_RPC:Category = "PMS Probe" ;
                AFSSP_RPC:_FillValue = -32767.f ;
                AFSSP_RPC:SampledRate = 10 ;
                AFSSP_RPC:DataQuality = "Preliminary" ;
                AFSSP_RPC:SerialNumber = "FSSP076" ;
                AFSSP_RPC:FirstBin = 1 ;
                AFSSP_RPC:LastBin = 15 ;
                AFSSP_RPC:CellSizes = 2.f, 5.f, 8.f, 11.f, 14.f, 17.f, 20.f, 23.f, 26.f, 29.f, 32.f, 35.f, 38.f, 41.f, 44.f, 47.f, 2.f, 4.f, 6.f, 8.f, 10.f, 12.f, 14.f, 16.f, 18.f, 20.f, 22.f, 24.f, 26.f, 28.f, 30.f, 32.f, 1.f, 2.f, 3.f, 4.f, 5.f, 6.f, 7.f, 8.f, 9.f, 10.f, 11.f, 12.f, 13.f, 14.f, 15.f, 16.f, 0.5f, 1.f, 1.5f, 2.f, 2.5f, 3.f, 3.5f, 4.f, 4.5f, 5.f, 5.5f, 6.f, 6.5f, 7.f, 7.5f, 8.f ;
                AFSSP_RPC:CellSizeUnits = "micrometers" ;
                AFSSP_RPC:DepthOfField = 2.7f ;
                AFSSP_RPC:BeamDiameter = 0.18f ;
                AFSSP_RPC:Density = 1.f ;
                AFSSP_RPC:PLWfactor = 1.e-06f ;
                AFSSP_RPC:DBZfactor = 1000000.f ;

        float WIC(Time, sps25) ;
                WIC:units = "m/s" ;
                WIC:long_name = "GPS-Corrected Wind Vector, Vertical Gust Component" ;
                WIC:standard_name = "wind_from_direction" ;
                WIC:Category = "Wind" ;
                WIC:_FillValue = -32767.f ;
                WIC:DataQuality = "Good" ;
                WIC:Dependencies = "9 TASX VEWC VNSC PITCH ROLL THDG ATTACK SSLIP VSPD" ;

Appendices

A  References

[NetCDF]
NetCDF Software Package, from the UNIDATA Program Center of the University Corporation for Atmospheric Research.
[UDUNITS]
UDUNITS Software Package, from the UNIDATA Program Center of the University Corporation for Atmospheric Research.

B  PMS1D Probe Variables Description

B.1 Variable names

Each PMS probe has two vectors (histograms) stored in the netCDF file. First is the raw counts and this variable always starts with the letter 'A' (for "Actual"). Second is the concentration array which always starts with the letter 'C'. Below you will see the variable names for the orginal PMS probe and for the DMT converted electronics.

PMS1D# binsCounts nameConcentration nameDMT Converted# binsCounts nameConcentration name
FSSP16AFSSP_xxxCFSSP_xxxS10031AS100_xxxCS100_xxx
PCASP16APCAS_xxxCPCAS_xxxS20031AS200_xxxCS200_xxx
F30032AF300_xxxCF300_xxxS30031AS300_xxxCS300_xxx
260X64A260X_xxxC260X_xxx------
2D-C32A1DC_xxx1C1DC_xxx1------
2D-C64A2DC_xxx2C2DC_xxx2------

Two sets of arrays are produced for each 2D probe. The 1DC pair is the produced from the 2D probe data using an "entire-in" algorithm. The 2DC pair is produced using a "center-in" algorithm.

B.2 Number of bins

The original PMS probes always had the first (zeroth) bin set to zero or a value which was not part of the size-distribution histogram (e.g. the FSSP-100 recorded 16 bins of which the last 15 contained data). However we wanted all raw data available to an investigator. For consistency this required the concentration array to be the same size, so it's first (zeroth) bin is also set to zero.

DMT converted probes can record either 10, 20, 30 or 40 bins of data. To maintain backwards compatability, we pad the front of the array with a zero and increase the number of bins by one (11, 21, 31 or 41).

NOTE: Always ignore the very first bin of any size-distribution histogram in our netCDF files.

B.3 Cell sizes

Cell sizes are stored as an attribute of the counts array. They are end-points. In the example below, there are 31 cell sizes to go with the 30 used values of the size- distribution. Also note the attributes FirstBin and LastBin, these should be used to determine which bins are considered valid, both are inclusive. 
        float AS100_LPC(Time, sps1, Vector31) ;
                AS100_LPC:_FillValue = -32767.f ;
                AS100_LPC:units = "count" ;
                AS100_LPC:long_name = "FSSP-100 Raw Accumulation (per cell) - DMT" ;
                AS100_LPC:Category = "PMS Probe" ;
                AS100_LPC:missing_value = -32767.f ;
                AS100_LPC:SampledRate = 10 ;
                AS100_LPC:DataQuality = "Preliminary" ;
                AS100_LPC:SerialNumber = "FSSP122" ;
                AS100_LPC:FirstBin = 3 ;
                AS100_LPC:LastBin = 30 ;
                AS100_LPC:CellSizes = 0.7f, 1.4f, 2.35f, 3.9f, 5.45f, 7.f, 8.55f, 10.1f, 11.65f, 13.2f, 14.75f, 16.3f, 17.85f, 19.4f, 20.95f, 22.5f, 24.05f, 25.6f, 27.15f, 28.7f, 30.25f, 31.8f, 33.35f, 34.9f, 36.45f, 38.f, 39.55f, 41.1f, 42.65f, 44.2f, 45.75f ;
                AS100_LPC:CellSizeUnits = "micrometers" ;
                AS100_LPC:DepthOfField = 2.3f ;
                AS100_LPC:BeamDiameter = 0.2f ;
                AS100_LPC:Density = 1.f ;
                AS100_LPC:PLWfactor = 1.e-06f ;
                AS100_LPC:DBZfactor = 1000000.f ;

        float CS100_LPC(Time, sps1, Vector31) ;
                CS100_LPC:_FillValue = -32767.f ;
                CS100_LPC:units = "#/cm3" ;
                CS100_LPC:long_name = "FSSP-100 Concentration (per cell) - DMT" ;
                CS100_LPC:Category = "PMS Probe" ;
                CS100_LPC:missing_value = -32767.f ;
                CS100_LPC:DataQuality = "Preliminary" ;
                CS100_LPC:Dependencies = "5 AS100_LPC TASX REJAT_LPC OVFLW_LPC FRNG_LPC" ;

        float CONCF_LPC(Time) ;
                CONCF_LPC:_FillValue = -32767.f ;
                CONCF_LPC:units = "#/cm3" ;
                CONCF_LPC:long_name = "FSSP-100 Concentration (all cells)" ;
                CONCF_LPC:Category = "PMS Probe" ;
                CONCF_LPC:missing_value = -32767.f ;
                CONCF_LPC:DataQuality = "Preliminary" ;
                CONCF_LPC:Dependencies = "1 CS100_LPC" ;

bin#0123....2930
n/aFirstBinLastBin

For further information contact:

Chris Webster
National Center for Atmospheric Research
Research Aviation Facility
(303) 497-1044

email
Last update: Sat Feb 18 15:06:44 MST 2006