% This is an example metadata run script for a Nortek Aquadopp. % Reqired NetCDF attributes are noted and must not be removed. % Users may remove any non-required attributes or add any % additional attributes they choose % % K.Rosenberger, 1/4/07, version 1.0, modified from Charlene Sullivan's metaRDI %get the global attributes - need the metadata=read_globalatts('../gatt_OBX10S21.txt'); baseFile = '063401'; instnum = '01'; %the instrument number on the tripod metadata.LatLonDatum = 'NAD83'; %if there was a significant clock error enter it here metadata.ClockError = 0; metadata.transducer_offset_from_bottom = 0.14; % height of transducer heads off the bottom metadata.orientation = 'UP'; % use this to identify orientation of profiler metadata.DATA_CMNT = 'Upward looking Aquadopp Profiler'; metadata.instrument_type = 'Nortek Aquadopp Profiler'; % type of instrument and instrument manufacturer metadata.initial_instrument_height = metadata.transducer_offset_from_bottom; metadata.instrument_depth = metadata.WATER_DEPTH - metadata.initial_instrument_height; metadata.head_rotation = 'horizontal'; % metadata.pred_accuracy = 2.0; % horizontal predicted accuracy metadata.magnetic_variation_applied = metadata.magnetic_variation; % degrees between magnetic and true north at data location metadata.inst_type = metadata.instrument_type; % type of instrument and instrument manufacturer %if you dont have a global attributes text file to read them in, fill in %these below... % metadata.MOORING = gatt.MOORING; % mooring identifier number % metadata.longitude = gatt.longitude; % (required) instrument longitude, negative is West % metadata.latitude = gatt.latitude; % (required) instrument latitude, positive is North % metadata.WATER_DEPTH = gatt.WATER_DEPTH; % water depth in meters % metadata.WATER_DEPTH_NOTE = gatt.WATER_DEPTH_NOTE; % source of water depth information % metadata.PROGRAM = gatt.PROGRAM; % program name % metadata.PROJECT = gatt.PROJECT; % project name % metadata.EXPERIMENT = gatt.EXPERIMENT; % experiment name % metadata.Region = gatt.Region; % project name % metadata.DESCRIPTION = gatt.DESCRIPTION; % location and/or site of data collection % metadata.magnetic_variation = gatt.MagneticVariation; % degrees between magnetic and true north at data location % metadata.SciPi = gatt.SciPi; % % N.B. - Negative is westerly variation % metadata.Site = gatt.Site; % metadata.MooringType = gatt.MooringType; % metadata.CruiseID = gatt.CruiseID; % metadata.Deployment_date = '26-Aug-2008 19:50:00'; % This needs to be actual deployment date and time on the Bottom % metadata.Recovery_date = '27-Aug-2008 17:48:00'; % This needs to be actual recovery date and time on the Bottom %these attributes are part of the EPIC standard and don't change fromnc %deployment to deployment, so keep them for the next deployment metadata.DATA_SUBTYPE = 'MOORED'; metadata.DATA_ORIGIN = 'USGS WRCMG Sed Trans Group'; metadata.COORD_SYSTEM = 'GEOGRAPHIC'; metadata.Conventions = 'PMEL/EPIC'; metadata.WATER_MASS = '?'; % water mass flag used for EPIC contouring programs metadata.POS_CONST = 0; % consistent position flag (=1 not consistent) metadata.DEPTH_CONST = 0; % consistent depth flag (=1 not consistent) metadata.DRIFTER = 0; % drifter flag (=1 if drifter) metadata.COMPOSITE = 0; % number of pieces in composite series %use these attributes if you have external sensors, otherwise comment them %out. For the OBS3+, Campbell lists the coefficients in the wrong order - %the 'Low' range coefficients are to be applied to the second channel, and %the 'High' range are to be applied to teh first channel. % metadata.AnalogInput1.sensor_type = 'OBS3+'; % metadata.AnalogInput1.manufacturer = 'D+A Instruments'; % metadata.AnalogInput1.serial_number = 'T8318'; % metadata.AnalogInput1.sensor_height = 0.34; %in meters % metadata.AnalogInput1.range = 250; %in NTU's - get from calibration certificate % metadata.AnalogInput1.cals.NTUcoef = [-0.15873 177.24707 19.0]; %[x^0 x^1 x^2] - get from calibration certificate % metadata.AnalogInput2.sensor_type = 'OBS3+'; % metadata.AnalogInput2.manufacturer = 'D+A Instruments'; % metadata.AnalogInput2.serial_number = 'T8318'; % metadata.AnalogInput2.sensor_height = 0.34; %in meters % metadata.AnalogInput2.range = 1000; %in NTU's - get from calibration certificate % metadata.AnalogInput2.cals.NTUcoef = [0.03757 47.30507 0.48189]; %[x^0 x^1 x^2] - get from calibration certificate % if you have done a tank calibration - leave empty if not %from calibration certificate % metadata.AnalogInput1.cals.SEDcoef = [0.2 30 3.5]; %[x^0 x^1 x^2] % metadata.AnalogInput1.cals.SEDcoefUnits = 'g/L'; % or whatever your units are % metadata.AnalogInput2.cals.SEDcoef = [0.2 30 3.5]; %[x^0 x^1 x^2] % metadata.AnalogInput2.cals.SEDcoefUnits = 'g/L'; % or whatever your units are %------------------Some Options for Correcting Data------------------------ %use this if you need to correct the time, for instance if the instrument %was not set in GMT - a positive value adds time, thus moving it forward %(e.g., for Hawaii, use 10 hours to correct to GMT). %comment out if not necessary - use hours % metadata.TimeCorrection = []; %if you want to crop bad data near the surface or bottom; % TODO want to work in backscatter at some point. % 'Water Level' if the instrument is pointing up and you want to use the water surface % 'Bin Number' where Number is the actual bin number, e.g. 'Bin 7', if instrument is pointing down... metadata.trim_method = 'Water Level'; %if you want to set your own amplitude cutoff for cutting out bad data %comment out if unnecessary % metadata.cutoff_ampl = 25; %use this to fill in any potential gaps in the record, e.g. if the wave %bursts were longer than the current interval. Comment out if not. metadata.fixgaps = 1; %-------------------------------------------------------------------------- % version of the metafile metadata.metadatafile_name = mfilename; % name of metadatafile metadata.metadatafile_author = 'Kurt Rosenberger, USGS MP'; % author of the metadatafile metadata.metadatafile_version = '2.0'; % version of the metadatafile %good ensembles for this deployment are 218-35684 % %now run the scripts diary(sprintf('run%s',datestr(now, 30))) % cdfFile = [metadata.MOORING instnum 'aq']; cdfFile = [metadata.MOORING instnum 'aq']; metadata.filename = cdfFile; if 1 aqdprf2cdf(baseFile,cdfFile,metadata); end %convert raw data to netcdf CDFFile = [cdfFile '-raw'];outFile = cdfFile; if 1 VEL = aqdprf2nc(CDFFile,outFile,metadata); end diary off