testMadrigalHdf5File.m
81 lines
| 2.9 KiB
| text/octave
|
MatlabLexer
| r0 | % test/example script to exercise MadrigalHdf5File | |||
| % | ||||
| % $Id: testMadrigalHdf5File.m 6539 2018-07-05 19:27:38Z brideout $ | ||||
| filename = '/Users/brideout/Documents/workspace/mad3_trunk/source/madmatlab/example.h5'; | ||||
| % Uncomment the statement below to test just creating a *.mat file | ||||
| % filename = '/Users/brideout/Documents/workspace/mad3_trunk/source/madmatlab/example.mat'; | ||||
| oneDParms = cellstr(char('azm', 'elm', 'beamid')); | ||||
| independent2DParms = cellstr(char('range')); | ||||
| twoDParms = cellstr(char('ti', 'dti', 'po+')); | ||||
| % Use {} for independent2DParms and twoDParms if all scalar parameters | ||||
| arraySplittingParms = cellstr(char('beamid')); | ||||
| % Use arraySplittingParms = {}; for no splitting | ||||
| % some hard-coded fake data | ||||
| ut1_unix = 1.0E9; | ||||
| ut2_unix = 1.0E9 + 100; | ||||
| kindat = 3410; | ||||
| kinst= 30; | ||||
| numRows = 5; | ||||
| azm = 100.0; | ||||
| elm = 45.0; | ||||
| poplus = [0.5, 0.5, 0.6, 0.6, 0.7]; | ||||
| range = [100.0, 150.0, 200.0, 250.0, 300.0]; | ||||
| ti = [1000.0, 1100.0, 1200.0, 1300.0, 1400.0]; | ||||
| dti = [100.0, 150.0, 200.0, 250.0, 300.0]; | ||||
| beamids = [1,2,1,2,1,2,1,2,1,2]; | ||||
| % new file must not already exist | ||||
| delete(filename); | ||||
| madFile = MadrigalHdf5File(filename, oneDParms, ... | ||||
| independent2DParms, twoDParms, arraySplittingParms); | ||||
| for rec = 1:10 | ||||
| madFile = madFile.appendRecord(ut1_unix, ut2_unix, kindat, ... | ||||
| kinst, numRows); | ||||
| % set 1D values | ||||
| madFile = madFile.set1DParm('azm', azm, madFile.lastRecord); | ||||
| madFile = madFile.set1DParm('elm', elm, madFile.lastRecord); | ||||
| madFile = madFile.set1DParm('beamid', beamids(rec), madFile.lastRecord); | ||||
| % set 2D and independent variables | ||||
| madFile = madFile.set2DParm('range', range, madFile.lastRecord); | ||||
| madFile = madFile.set2DParm('ti', ti, madFile.lastRecord); | ||||
| madFile = madFile.set2DParm('dti', dti, madFile.lastRecord); | ||||
| madFile = madFile.set2DParm('po+', poplus, madFile.lastRecord); | ||||
| % time must always increase in Madrigal 3 | ||||
| ut1_unix = ut1_unix + 100; | ||||
| ut2_unix = ut1_unix + 200; | ||||
| end | ||||
| % add catalog and header info | ||||
| % Note: this method is optional. Even if this method is not | ||||
| % called, the catalog record will contain a description of the | ||||
| % instrument (kinst code and name) and kind of data brief | ||||
| % description, along with a list and descriptions of the | ||||
| % parameters in the file, and the first and last times of the | ||||
| % measurements. | ||||
| principleInvestigator = 'Bill Rideout'; | ||||
| expPurpose = 'Measure the ionosphere'; | ||||
| expMode = 'Vector measurements'; | ||||
| cycleTime = 20.0; | ||||
| correlativeExp = ''; | ||||
| sciRemarks = 'Big solar storm'; | ||||
| instRemarks = 'Working well!!!'; | ||||
| madFile = madFile.setCatalog(principleInvestigator, expPurpose, expMode, ... | ||||
| cycleTime, correlativeExp, sciRemarks, instRemarks); | ||||
| kindatDesc = 'Regular processing'; | ||||
| analyst = 'Phil Erickson'; | ||||
| comments = 'Include unusual parameter description here'; | ||||
| history = 'Reprocessed four times'; | ||||
| madFile = madFile.setHeader(kindatDesc, analyst, comments, history); | ||||
| write(madFile); | ||||
