MMCR Data Stream Format Changes

 MMCR Processor Upgrade

mmcrmom netCDF Header


Color Key:  Black = Same ; Blue = Changed ;  Green = Added ;  Purple = Deleted
New MMCR Data Stream

SGP as of 9/9/2003
NSA as of 4/15/2004
 MMCR Data Streams Prior to Upgrade

SGP before 9/9/2003,
NSA before 4/14/2004,
All TWP-C1, TWP-C2

 mmcrmomC1.b1 mmcrcalC1.a1 mmcrmomentsC1.a1
base_time
 base_time
 base_time
time_offset
 time_offset
 time_offset
time

qc_time

ModeDescription (mode)

RxCalTimeStamp (mode)

RxGain (mode)

SkyNoiseLevel (mode)

Rx290KLevel (mode)

CalCheckLevel (mode)

CalCheckTime (mode)

RadarConstant (mode)

MinimumDetectableReflectivity
(hourly, mode, heights)

InterPulsePeriod (mode)
 InterPulsePeriod (time)
 InterPulsePeriod (time)
PulseWidth (mode)
 PulseWidth (time)
 PulseWidth (time)
StartGateDelay (mode)
 StartGateDelay (time)
 StartGateDelay (time)
GateSpacing (mode)
 GateSpacing (time)
 GateSpacing (time)
NumCoherentIntegrations (mode)
 NumCoherentIntegrations (time)
 NumCoherentIntegrations (time)
NumSpectralAverages (mode)
 NumSpectraAveraged (time)
 NumSpectraAveraged (time)
NumFFT (mode)
 NumFft (time)
 NumFft (time)
NumHeights (mode)
 NumHeights (time)
 NumHeights (time)
NumCodeBits (mode)
 NumCodeBits (time)
 NumCodeBits (time)
NumReceivers (mode)


ReceiverMode (mode)

ReceiverNumber (mode)

NyquistVelocity (mode)
NyquistVelocity (time)
DCFilteringONOFF (mode)

WindowingONOFF (mode)

ClutterHeight (mode)

heights (mode, heights)
 Heights (time, height)
 Heights (time, height)
ModeNum (time)

DataQualityStatus (time)

MeanDopplerVelocity (time, heights)
Units: m/s
MeanDopplerValue (time, height)  
Unitless
SignalToNoiseRatio (time, heights)
Units: dB
SignalToNoiseRatio (time, height)
Units: dB * 100
Power (time, heights)

SpectralWidth (time, heights)
Units: m/s
SpectralWidthValue (time, height)
Unitless
NoiseLevel (time, heights)
Units: dB
Noise (time, height)
Units: millibells
Reflectivity (time, heights)
 Reflectivity (time,  height)
RangeCorrectedPower (time, heights)
 RangeCorrectedCalibratedPower
(time, height)
CircularDepolarizationRatio
(time, heights)

AvgNoiseLevel (time)

TimeAvg (hourly)

PeakTxPowerAverage (hourly)
 PeakTransmittedPower (time)

TWTStatus (hourly)

lat
 lat
 lat
lon
 lon
 lon
alt alt
 alt


Delay (time)


Switch (time)

TransmitterPathLoss (time)

ReceiverPathLoss (time)

CirculatorAttenuation1(time)

CodedPulseLoss (time)

PartialDecodingLoss (time, entries)

RxCalibrationNoiseCrv
(time, crventries)

RxCalibrationNoiseCrvSignal
(time, crventries)


Sample sgpmmcrmomC1.b1 netcdf header:

netcdf sgpmmcrmomC1.b1.20031230.000016 {
dimensions:
        time = UNLIMITED ; // (28863 currently)
        mode = 10 ;
        namelength = 32 ;
        heights = 167 ;
        hourly = 24 ;
variables:
        int base_time ;
                base_time:string = "30-Dec-2003,0:00:16 GMT" ;
                base_time:long_name = "Base time in Epoch" ;
                base_time:units = "seconds since 1970-1-1 0:00:00 0:00" ;
        double time_offset(time) ;
                time_offset:long_name = "Time offset from base_time" ;
                time_offset:units = "seconds since 2003-12-30 00:00:16 0:00" ;
        double time(time) ;
                time:long_name = "Time offset from midnight" ;
                time:units = "seconds since 2003-12-30 00:00:00 0:00" ;
        int qc_time(time) ;
                qc_time:long_name = "Results of quality checks on sample time" ;
                qc_time:units = "unitless" ;
                qc_time:description = "The qc_time values are calculated by comparing each sample\n",
    "time with the previous time (i.e. delta_t = t[n] - t[n-1]).\n",
    "If the \'qc_check_prior\' flag is set the first sample time\n",
    "from a new raw file will be compared against the time just\n",
    "previous to it in the stored data. If the \'qc_check_prior\'\n",
    "flag is not set the qc_time value for the first sample time\n",
    "will be set to 0\n",
    "\n",
    "The qc_time bit values are as follows:\n",
    "=========================================================\n",
    "0x0 = delta time is within the specified range\n",
    "0x1 = delta time is equal to 0, duplicate sample times\n",
    "0x2 = delta time is less than the \'delta_t_lower_limit\'\n",
    "0x4 = delta time is greater than the \'delta_t_upper_limit\'\n",
    "\n",
    "" ;
                qc_time:delta_t_lower_limit = 0 ;
                qc_time:delta_t_upper_limit = 60 ;
                qc_time:prior_sample_flag = 1 ;
        char ModeDescription(mode, namelength) ;
                ModeDescription:long_name = "radar mode char identifier" ;
                ModeDescription:units = "unitless" ;
                ModeDescription:missing_value = "0" ;
        int RxCalTimeStamp(mode) ;
                RxCalTimeStamp:long_name = "Receiver Cal Time Stamp" ;
                RxCalTimeStamp:units = "s" ;
                RxCalTimeStamp:missing_value = -9999 ;
        float RxGain(mode) ;
                RxGain:long_name = "Receiver Gain" ;
                RxGain:units = "dB" ;
                RxGain:missing_value = -9999.f ;
        float SkyNoiseLevel(mode) ;
                SkyNoiseLevel:long_name = "Receiver Sky Noise" ;
                SkyNoiseLevel:units = "dB" ;
                SkyNoiseLevel:missing_value = -99999.f ;
        float Rx290KLevel(mode) ;
                Rx290KLevel:long_name = "Receiver 290K Level" ;
                Rx290KLevel:units = "dB" ;
                Rx290KLevel:missing_value = -9999.f ;
        float CalCheckLevel(mode) ;
                CalCheckLevel:long_name = "Receiver Cal Check Level" ;
                CalCheckLevel:units = "dB" ;
                CalCheckLevel:missing_value = -9999.f ;
        int CalCheckTime(mode) ;
                CalCheckTime:long_name = "Receiver Cal Check Time Stamp" ;
                CalCheckTime:units = "s" ;
                CalCheckTime:missing_value = -9999 ;
        float RadarConstant(hourly, mode) ;
                RadarConstant:long_name = "Radar Constant" ;
                RadarConstant:units = "dB" ;
                RadarConstant:missing_value = -9999.f ;
        float MinimumDetectableReflectivity(hourly, mode, heights) ;
                MinimumDetectableReflectivity:long_name = "Minimum detectable reflectivity" ;
                MinimumDetectableReflectivity:units = "dBZ" ;
                MinimumDetectableReflectivity:missing_value = -9999.f ;
        int InterPulsePeriod(mode) ;
                InterPulsePeriod:long_name = "Inter-Pulse Period" ;
                InterPulsePeriod:units = "ns" ;
                InterPulsePeriod:missing_value = -9999 ;
        int PulseWidth(mode) ;
                PulseWidth:long_name = "Pulse Width" ;
                PulseWidth:units = "ns" ;
                PulseWidth:missing_value = -9999 ;
        int StartGateDelay(mode) ;
                StartGateDelay:long_name = "Delay to First Range Gate" ;
                StartGateDelay:units = "ns" ;
                StartGateDelay:missing_value = -9999 ;
        int GateSpacing(mode) ;
                GateSpacing:long_name = "Gate Spacing" ;
                GateSpacing:units = "ns" ;
                GateSpacing:missing_value = -9999 ;
        int NumCoherentIntegrations(mode) ;
                NumCoherentIntegrations:long_name = "Number of Coherent Integrations" ;
                NumCoherentIntegrations:units = "count" ;
                NumCoherentIntegrations:missing_value = -9999 ;
        short NumSpectralAverages(mode) ;
                NumSpectralAverages:long_name = "Number of Spectral Averages" ;
                NumSpectralAverages:units = "count" ;
                NumSpectralAverages:missing_value = -9999s ;
        short NumFFT(mode) ;
                NumFFT:long_name = "Number of Points in FFT" ;
                NumFFT:units = "count" ;
                NumFFT:missing_value = -9999s ;
        short NumHeights(mode) ;
                NumHeights:long_name = "Number of Range Gates" ;
                NumHeights:units = "count" ;
                NumHeights:missing_value = -9999s ;
        short NumCodeBits(mode) ;
                NumCodeBits:long_name = "Number of Code Bits" ;
                NumCodeBits:units = "count" ;
                NumCodeBits:missing_value = -9999s ;
        short NumReceivers(mode) ;
                NumReceivers:long_name = "Number of Receiver" ;
                NumReceivers:units = "count" ;
                NumReceivers:missing_value = -9999s ;
        short ReceiverMode(mode) ;
                ReceiverMode:long_name = "Receiver Mode" ;
                ReceiverMode:units = "count" ;
                ReceiverMode:missing_value = -9999s ;
        short ReceiverNumber(mode) ;
                ReceiverNumber:long_name = "Current Receiver Number" ;
                ReceiverNumber:units = "count" ;
                ReceiverNumber:missing_value = -9999s ;
        float NyquistVelocity(mode) ;
                NyquistVelocity:long_name = "Nyquist Velocity" ;
                NyquistVelocity:units = "m/s" ;
                NyquistVelocity:missing_value = -9999.f ;
        short DCFilterONOFF(mode) ;
                DCFilterONOFF:long_name = "DC Filtering ON-OFF Status" ;
                DCFilterONOFF:units = "count" ;
                DCFilterONOFF:missing_value = -9999s ;
        short WindowingONOFF(mode) ;
                WindowingONOFF:long_name = "Windowing ON-OFF Status" ;
                WindowingONOFF:units = "count" ;
                WindowingONOFF:missing_value = -9999s ;
        float ClutterHeight(mode) ;
                ClutterHeight:long_name = "Max. height of clutter removal" ;
                ClutterHeight:units = "m" ;
                ClutterHeight:missing_value = -9999.f ;
        float heights(mode, heights) ;
                heights:long_name = "Range Heights (center of radar sample volume)" ;
                heights:units = "m MSL" ;
                heights:missing_value = -9999.f ;
        short ModeNum(time) ;
                ModeNum:long_name = "Operating Set for this Record" ;
                ModeNum:units = "count" ;
                ModeNum:missing_value = -9999s ;
        int DataQualityStatus(time) ;
                DataQualityStatus:long_name = "Data Quality Status" ;
                DataQualityStatus:units = "code" ;
                DataQualityStatus:missing_value = -9999 ;
        float MeanDopplerVelocity(time, heights) ;
                MeanDopplerVelocity:long_name = "Mean Doppler Velocity" ;
                MeanDopplerVelocity:units = "m/s" ;
                MeanDopplerVelocity:resolution = 0.001f ;
                MeanDopplerVelocity:missing_value = -9999.f ;
        float SignalToNoiseRatio(time, heights) ;
                SignalToNoiseRatio:long_name = "Signal to Noise Ratio" ;
                SignalToNoiseRatio:units = "dB" ;
                SignalToNoiseRatio:resolution = 0.001f ;
                SignalToNoiseRatio:missing_value = -9999.f ;
        float Power(time, heights) ;
                Power:long_name = "Power (uncalibrated)" ;
                Power:units = "dB" ;
                Power:resolution = 0.001f ;
                Power:missing_value = -9999.f ;
        float SpectralWidth(time, heights) ;
                SpectralWidth:long_name = "Spectral Width" ;
                SpectralWidth:units = "m/s" ;
                SpectralWidth:resolution = 0.001f ;
                SpectralWidth:missing_value = -9999.f ;
        float NoiseLevel(time, heights) ;
                NoiseLevel:long_name = "Mean Noise Level" ;
                NoiseLevel:units = "dB" ;
                NoiseLevel:resolution = 0.001f ;
                NoiseLevel:missing_value = -99999.f ;
        float Reflectivity(time, heights) ;
                Reflectivity:long_name = "Reflectivity" ;
                Reflectivity:units = "dBZ" ;
                Reflectivity:resolution = 0.001f ;
                Reflectivity:missing_value = -9999.f ;
                Reflectivity:comment1 = "Nearfield correction has been applied" ;
                Reflectivity:comment2 = "Reference: Sekelsky, S. (2002), JAOT, v19, no. 4" ;
        float RangeCorrectedPower(time, heights) ;
                RangeCorrectedPower:long_name = "Range Corrected Calibrated Power" ;
                RangeCorrectedPower:units = "dBm" ;
                RangeCorrectedPower:resolution = 0.001f ;
                RangeCorrectedPower:missing_value = -9999.f ;
        float CircularDepolarizationRatio(time, heights) ;
                CircularDepolarizationRatio:long_name = "Circular Depolarization Ratio" ;
                CircularDepolarizationRatio:units = "dB" ;
                CircularDepolarizationRatio:resolution = 0.001f ;
                CircularDepolarizationRatio:missing_value = -9999.f ;
        float AvgNoiseLevel(time) ;
                AvgNoiseLevel:long_name = "Average Noise Level (S/N<0)" ;
                AvgNoiseLevel:units = "dB" ;
                AvgNoiseLevel:resolution = 0.001f ;
                AvgNoiseLevel:missing_value = -9999.f ;
        int TimeAvg(hourly) ;
                TimeAvg:long_name = "Time associated with hourly averaged values" ;
                TimeAvg:units = "seconds since January 1, 1970" ;
                TimeAvg:missing_value = -9999 ;
        float PeakTransmittedPowerAvg(hourly) ;
                PeakTransmittedPowerAvg:long_name = "Peak tranmitted power, averaged over the course of the hour" ;
                PeakTransmittedPowerAvg:units = "dBm" ;
                PeakTransmittedPowerAvg:missing_value = -9999.f ;
        int TWTStatusCode(hourly) ;
                TWTStatusCode:long_name = "TWT Status Code, details forthcoming from MMCR vendor" ;
                TWTStatusCode:units = "unknown" ;
                TWTStatusCode:missing_value = -9999 ;
        float lat ;
                lat:long_name = "north latitude" ;
                lat:units = "degrees" ;
                lat:valid_min = -90.f ;
                lat:valid_max = 90.f ;
        float lon ;
                lon:long_name = "east longitude" ;
                lon:units = "degrees" ;
                lon:valid_min = -180.f ;
                lon:valid_max = 180.f ;
        float alt ;
                alt:long_name = "altitude" ;
                alt:units = "meters above Mean Sea Level" ;

// global attributes:
                :ingest_software = " mmcrmom_ingest.c,v 1.1 2003/11/07 14:35:14 koontz process-ingest-mmcrmom_ingest-7.8-0 $" ;
                :input_source = "mmcrmom1:/data/collection/sgp/sgpmmcrmomC1.00/20033640000MMCRMom.nc" ;
                :site_id = "sgp" ;
                :facility_id = "C1 : Central_Facility" ;
                :sample_int = "variable" ;
                :averaging_int = "None." ;
                :serial_number = " " ;
                :comment = " " ;
                :resolution_description = "The resolution field attributes refer to the \n",
    "number of significant digits relative to the decimal point that should\n",
    "be used in calculations.  Using fewer didgits might result in greater uncertainty,\n",
    "using a larger number of digits should have no effect and thus is unnecessary.\n",
    "However, analyses based on differences in values with a larger number of scientific \n",
    "digits than indicated could lead to erroneous results or misleading scientific\n",
    "conclusions, resolution for lat = 0.001, resolution for lon = 0.001,\n",
    "resultion for alt = 1" ;
                :comment_on_time = "The time stamp comes at the beginning and end of the sample period.  Sample period is made up of data collection and processing time." ;
                :radar_wavelength = "8.600115e-003 m" ;
                :antenna_diameter = "3.05 m" ;
                :radar_operating_frequency = "34.86 GHz" ;
                :nearfield_correction_description = "The nearfield correction to Reflectivity has been applied. Reference: Sekelsky, S. (2002), JAOT, v19, no. 4" ;
                :proc_level = "b1" ;
                :configuration_file = "/data/conf/sgp/sgpmmcrmomC1/RCS/sgpmmcrmomC1.conf,v" ;
                :zeb_platform = "sgpmmcrmomC1.b1" ;
                :history = "created by user dsmgr on machine left at 30-Dec-2003,2:19:31, using $State: ds-zebra-zeblib-4.12-0 $" ;




Last Updated: 5/5/2004
kjohnson@bnl.gov