|
|
|
|
IAS METHOD FOR 2D DATA ANALYSIS ON PC'sInstitute of Atmospheric Sciences South Dakota School of Mines and Technology 501 E. St. Joseph Street Rapid City, SD 57701-3995
CHAPTER 5: ANALYZE and ANELAPThe rather large archive files created by PRELIM2D or TOUCH2D are not in a format that the typical user would like to play with. Often there is an interest in determining particle concentration statistics. Programs ANALYZE and ANELAP are used to generate much smaller ASCII files which contain particle count information by time, habit, and diameter size. A listing of these statistics often proves satisfactory by itself, or the user may write other programs to further analyze or plot the data. The two programs, ANALYZE and ANELAP, are virtually identical to the user. However, ANALYZE outputs its records whenever a specified number of 2D buffers or clock seconds have elapsed, whereas ANELAP uses as its output criterion whether or not the cumulative active time of the 2D probe has reached a specific value.Running ANALYZE & ANELAP These two programs read an archive file and produce a variable-length ASCII file with information broken down by time, particle type (habit), and particle size. A second file with particle mass information may also be generated. Again, make sure the executable file (.EXE suffix) is in the current directory or path, and enter either ANALYZE or ANELAP The program asks for the 3-digit flight number, flt, from which it creates both input and output filenames. It first looks for an archive file with the name NEWflt.ARC. If not found, it attempts to open PMSflt.ARC. If still unsuccessful, it requests a file pathname from the user. Once a file has been located, ANALYZE prompts for the following information. ANELAP is similar, except where noted below.
CONC_flt.n and MASS_flt.n, whereas those created by ANELAP are nearly the same, namely CONCEflt.n and MASSEflt.n Once again, flt is the 3-digit flight number and n is a sequential number beginning with 1. Note that each time either program is executed, the sequence starts over at 1. Any previously existing file by the same name will be overwritten! The value of n merely increments during the program when more than one time or buffer range (usually corresponding to a particular cloud penetration) is requested. Layout of Mass FileThe output records for the mass files are extremely simple. There is only one record type, which consists of a start time, stop time, and mass value in grams per cubic meter. A Fortran format statement that would match the output record would be something like (217.6,F7.3).Layout of Concentration FileThis ASCII file consists of a number of comment lines which describe the actual data, then an end-of-comments line consisting of asterisks, followed by the actual particle concentration information itself. If the user did not enter "Y" at step 7 of the program instructions, however, the descriptor records (down through the line of asterisks) will be omitted from this file. All of these latter data records are essentially free format; that is, one field is separated from the next by one or more spaces. The sequence of the data records is as follows:Record type 1: Header record consisting of the following fields: flt #, # type 2 records (nsel), # habit records output per time/buffer interval (nout), # particle size categories (bins), bin increment (microns), midpoint value of bin 1 (microns). Only the first three fields are generally of any interest to the user. Record type 2: Names of habits analyzed and/or output during the file's creation. The record consists of an integer habit number followed by the name of the habit. The number of these records in this file equals the value of the 2nd field in record type 1 (nsel). If the habit number is prefixed with a minus sign, it means that there are no records for this individual habit, but the habit was included in the totals line for each output interval. Record types 1 and 2 appear only once in the file. After that comes a sequence of one type 3 record followed by "nout" type 4 records. This sequence repeats, one set for each time/buffer interval that was in the start-to-end range. The description of these records is as follows: Record type 3: Time fields. The interval start and stop times (HHMMSS) are followed by probe active time, pressure, temp, FSSP or JW liquid water, and 2D mass. (The mass values are in grams per cubic meter.) To distinguish the liquid water source and which mass calculation was used (Andy's or spherical), the JW and spherical 2d mass values are output as the negative of their value. The time values are integers; the rest are floating point numbers. Record type 4: Frequency and concentration values. The layout is as follows: Habit #, # diams with non-zero counts (n), {diam 1, freq 1, conc 1}, {diam 2, freq 2, conc 2}, ..., {diam n, freq n, conc n}, {freq for all diams, conc for all diams} (Note: Neither the braces nor the commas appear in the data.) The second field in this record tells how many sets of three values follow it. At the end, there are two totals values representing particle count and concentration values for all diameter sizes. The diameter values are floating point values, whereas the other numbers are integers. Concentration values are in number per cubic meter. The records described above are terminated by a carriage return, line feed combination. Fields are separated with spaces. Sometimes the records are so long (especially All Habits) that when printed or displayed on the screen they wrap around and may appear to be several records. In actuality, they are each just one rather long record.
Chapter 1 - Discussion of Probes Chapter 2 - Particle Shadow Classification |