ABSTRACT SIMULTANEOUS ANALYSIS OF THE FIRST POSITIVE BAND SYSTEM OF N2 IN THE NEAR INFRARED EXCITED BY A HOLLOW CATHODE DISCHARGE BY Lamar E. Bullock A flowing gas, double-ended, water-cooled hollow cathode discharge tube has been designed and constructed to operate in a longitudinal magnetic field. This source was used to excite the first positive system of N2, and high resolution spectra of the bands between 1 and 1.5 microns were recorded. A non-linear least squares analysis using the existing theory of the excited electronic A32: and Bsflg states was made to determine a set of simultaneous constants for the first positive system of N2. A systematic difference between the theory and experimental data for the spin-orbit inter- action energy was found. The effect of this discrepancy on this analysis and its relationships to previous analyses are discussed. Two sets of simultaneous constants and their 95% simultaneous confidence intervals are given and discussed. One set came from the analysis of the near infrared data presented here and the other from the analysis of energy Lamar E. Bullock level data found in the literature. The vibrational dependence of the spin-orbit coupling constant A and the spin-splitting constant A are calculated for the first time. Computer programs were written to calculate and identify the transition frequencies and to perform a least squares fit of the theory to the experimental data. A complete description of their function and use is included. The appendices contain listings of the computer programs and a complete list of the identified frequencies with their weights and deviations from the calculated values. SIMULTANEOUS ANALYSIS OF THE FIRST POSITIVE BAND SYSTEM OF N2 IN THE NEAR INFRARED EXCITED BY A HOLLOW CATHODE DISCHARGE BY Lamar E} Bullock A THESIS Submitted to Michigan State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Department of Physics 1969 To Patricia ii ACKNOWLEDGMENTS I am thankful to Professor C. D. Hause for his guidance and encouragement throughout the course of this research projecc. The freedom he has given me to make decisions on mv own, coupled with his patience and understanding, has provided a very educational research experience. Professor J. A. Cowen has been most helpful in supplying equipment and advice at various times. Discussion with Professors T. H. Edwards and P. M. Parker concerning various aspects of this work have been very useful. The cooperation and guidance of former graduate student and close friend Dr. Donald Keck, on projects of mutual interest and responsibility, have been helpful and is greatly appreciated. I also wish to express my gratitude to fellow graduate student Richard Peterson for his as— sistance and advice on many projects. The cordial re- lationships of all members of the Molecular Spectrosc0py group have been a great asset. The staff of the Michigan State Computer Center have been very c00perative in providing assistance and facilities necessary for the completion of this work. Members of the staff of the Glass Shop and Physics machine Shop, as well as Mr. Mercer of the Cyclotron Shop,have iii provided the skill and patience necessary to construct the equipment described herein and their willing cooperation has been helpful. Several undergraduate students, paid on an hourly basis, have contributed greatly to this project by assisting in a wide range of tasks from measuring, tabulating and checking data; to designing, constructing and operating equipment; also in helping to write, check out and use computer programs. I am indebted to these students for their willing and skillful assistance and to the MSU Breakthrough Research Committee as well as the Work Study program for providing funds for their support. The Breakthrough grant was also helpful in providing funds for the purchase of the necessary supplies and equipment. The support of my graduate studies through a National Aeronautics and Space Administration Traineeship is sincerely appreciated. The dedication of this thesis to my wife is appropriate because without her understanding, cooperation and assistance throughout my graduate education, this thesis would not exist. Also her patience, care and diligence in the preparation of the many drafts of this work was invaluable. iv TABLE OF CONTENTS LIST OF TABLES LIST OF FIGURES LIST OF APPENDICES INTRODUCTION Chapter I. II. III. DOUBLE ENDED HOLLOW CATHODE DISCHARGE TUBE Gas Handling Electrical Power Magnetic Field Initial Discharge Tube Redesigned Discharge Tube THEORY FOR THE FIRST POSITIVE SYSTEM OF N2 Electronic Term Vibrational Term Rotational Term A32: Energy Levels B3Hg State of N2 Selection Rules Frequency Determination NEAR INFRARED SPECTRUM OF N2 SHAFT Page vii viii 10 14 19 20 20 20 22 23 25 27 29 33 Chapter Page III. Cont. Structure of Data Deck 42 IV. ANALYSIS OF N2 SPECTRUM 54 Linearized Least Squares 55 Analysis Programs 56 Structure of Data Deck 57 Observed Data Cards 62 D3SIG 63 Simultaneous Constants from the Near Infrared 64 Simultaneous Constants from the Data of Dieke and Heath 65 SUMMARY AND CONCLUSIONS 77 LIST OF REFERENCES 79 APPENDICES 81 vi Table LIST OF TABLES Experimental Characteristics of the Near Infrared Spectra Sample Tabulated Output from SHAFT 1 Simultaneous Constants in cm- for the A32 State of N2 _1 Simultaneous Constants in cm for the B3n State of N2 _1 Molecular Constants in cm for the A52 State of N2 from the Literature _1 Molecular Constants in cm for the BJH State of N2 from the Literature P Branches of the (0,0) Band Q R "U ”0 'U 2310 W O Branches Branches Branches Branches Branches Branches Branches Branches Branches Branches Branches of of of of of of of of of of of the the the the the the the the the the the (0,0) (0,0) (0,1) (0,1) (0,1) (1,2) (1,2) (1,2) (0,2) (0,2) (0,2) vii Band Band Band Band Band Band Band Band Band Band Band Page 32 4O 66 67 69 70 140 142 144 146 148 150 152 154 156 158 160 162 Table 19. 20. 21. 22. 23. 24. ’U FUD IO Branches Branches Branches Branches Branches Branches of of of of of of the the the the the (1,3) (1,3) (1,3) (2.4) (2,4) (2,4) viii Band Band Band Band Band Band Page 164 166 168 170 172 174 LIST OF FIGURES Figure 1. Gas Distribution System 2. Electrical System 3. Magnet Characteristics 4. Water Cooled Stainless Steel Hollow Cathode 5. End Window and Anode Assemblies 6. End Window-Anode and Gas-Coolent Feedthrough Assemblies 7. Water Cooled Hollow Cathode 8. Allowed Rotational Transitions from a 3Hg State to a Eu State 9. (0,0) Band of First Positive System of N2 10. Block Diagram of SHAFT 11. SHAFT Input for Averaging Two Groups 12. SHAFT Input for Calculating and Using Calibration Constants 13. SHAFT Input for Calculating and Using Calibration Constants and for Identifying the Resulting Frequencies 14. Ban Energy Levels for v=4 Lower A-Doublet 15. Ban Energy Deviations for v=4 Lower A-Doublet 16. Deviations from the 13400 Band of PH Reported by Gilbert ix Page 11 12 13 16 17 21 30 35 51 52 53 72 73 75 INTRODUCTION Hollow cathode discharge tubes have been used as spectral sources for many years and recently as lasers (1,2). This thesis describes the design, construction, and testing of a double ended hollow cathode discharge tube and the analysis of the molecular spectra excited using this dis- charge tube. The discharge tube was designed such that it could be used as an excitation source for gaseous laser studies. The hollow cathode source consisted of a cylindrical stainless steel cathode with stainless steel anodes at each end. The gas handling system was designed to allow the use of any gas in either a flowing or non-flowing configuration. The electrical system contained separate water cooled ballast resistors for each end and a meter system for moni- toring the electrical characteristics of the discharge. End window assemblies were mounted coxial with the cathode to allow the use of Brewster angle windows for laser appli- cation. The first positive system of N2 was excited by the hollow cathode source and high resolution spectra in the near infrared were recorded. An analysis of the bands between 1. and 1.5 microns was performed. The first positive 2 3 system of N2 consists of transitions from the B Hg electronic state to the A32: electronic state. Extensive analysis of bands at wave lengths below 1 micron appear in the literature (3, 4, 5, 6, 7). However, all previous analyses have been done using combination difference methods. This study determines the simultaneous values of the molecular constants and their 95% simultaneous confidence intervals (SCI). Simultaneous values of the molecular constants means that set of estimators of the molecular constants for which the weighted sum of the squares of the deviations is a minimum. All the vibration bands were fit simultaneously so that the least squares techniques used gives the best set of constants for all the bands studied. The meaning of the 95% SCI is: the probability that all the true constants are within a range of i SCI of the estimated values is .95, assuming they were determined from one sampling of a normally distributed infinite population of similar sets of data. The molecular model of the first positive system of N2 is described in Chapter II. This model is non-linear in the molecular constants; thus, a linearization technique was employed so that linear regression analysis combined with an iteration procedure could be used to determine the molecular constants. Determination of simultaneous molecular constants in a model with 26 variables using in excess of 2000 transitions requires a tremendous number of calculations. One of the reasons this was not done in the earlier works was that 3 high speed digital computers were not available at that time. This study required extensive development of computer pro- grams both for determining and identifying the transition frequencies and for determining the molecular constants. Descriptions of these programs and instructions for their use are included. CHAPTER I DOUBLE ENDED HOLLOW CATHODE DISCHARGE TUBE A double ended, hollow cathode discharge-tube was designed and constructed so that it could be used both as a spectral source and as a laser. The laser application placed several design restrictions on the discharge tube which might not be included in the design of a hollow cathode for use as a spectral source; however, the two applications were not found to be incompatible. The hollow cathode was designed to have a continuously flowing gas system and to operate in a longitudinal magnetic field. Two different designs were tried. The first was a very compact, stainless steel, water cooled device with rapidly interchangeable discharge tubes and side mounted anodes. Heating problems restricted this design to rela- tively low current levels. The second design used a cy- lindrical pyrex vacuum chamber with an axially mounted, water cooled, stainless steel cathode and end mounted anode assemblies. Again the system was designed to allow for rapid interchangeability of the discharge tube. Gas Handling The gas was distributed through a manifold system as 5 illustrated in Figure 1. When mixtures were used, the gas was mixed and stored in a #1A nitrogen cylinder, which was cleaned by pumping to 10-“ torr for several days. Also during the process of testing the discharge tube, the cylinder was flushed several times with ultra high purity helium at a pressure of approximately 500 torr. During operation of the discharge, the gas flowed back into the manifold from the cylinder through a Gilmont 3230-10 flowmeter. The gas flow from the manifold was controlled by a Gilmont 3235 micrometer capillary valve. The reservoir gas pressure was measured by means of a mercury manometer attached to the manifold. The pressure on the input and output side of the discharge tube was measured using an oil manometer. From the capillary valve, the gas flowed through the anode, into the discharge area, and out through the cathode into a Welch 140 liter/minute Duo-Seal mechanical pump. The establishment of a stable discharge proved to be quite difficult. In general, a new discharge tube was very unstable until it had operated for several hours. It was also unstable for some time after being exposed to air at atmospheric pressure. The stability problem was attributed, at least in part, to the release of trapped gas by sputtering action in the cathode. This would also explain the decrease in stability when the discharge characteristics were changed, as the discharge would then be acting on a different pro- portion of the surface. Examination of the discharge tube after several hours of operation revealed a considerable Emumhm coaunnwuumwo mom .a mun—mew 1 _ 1 1 _ 1 “ acme _ 1 09.332 1. a 1 3 J. .. .1 mm mz " . Md moofimo L.“ 1 . m M 1 1 BOHHOM 1. 11.. . _ . 1.1. How 080:0: 1 1 1 . V Shauna: . M , r . _ r v . .11. .1 fix “mum: . a . 1. Sam . Levee 1mmmw ,o a 1. 17. 1 Ar (\ i r X /1\ ax odomeemz mmo mmmao I .II'IIIIO; . 7 redistribution of metal, indicating that a great deal of sputtering occurred. This problem with stability made a study of discharge characteristics, with respect to gas parameters, difficult. However, some general trends were observed. There was a very broad range from about .5 to 7 torr, over which gas pressure had very little effect on the output energy. At lower pressures the energy dropped off considerably; and at higher pressures, the discharge became unstable. It was also found that the percent concentration of N2 had very little effect on the energy output between about 5 percent and 100 percent; however, the discharge current was slightly higher at large concentrations, and the stability was not as good. Argon was also used as a carrier gas. It seemed to provide a stable operation with a somewhat lower current and voltage than with helium. However, the excitation of the N2 first positive system was not as pronounced with argon. Mixing a small percentage of argon with the helium had no apparent effect on the discharge characteristics or the nitrogen excitation; thus it was useful as an internal calibration standard. An increase in gas flow caused a momentary increase in output; but this returned to near the previous value when equilibrium was reached. Electrical Power Electrical power-for the discharge was obtained from five 0-150 volt D. C. generators connected in series. The output voltage of each generator was controlled by a 8 rheostat in series with the field windings of the generator. The series output from these generators provided the input to the electrical system shown in Figure 2. With this system it was possible to control the discharge current, both by variations in the input voltage, and by adjustment of the ballast resistors. The meter system monitored the total current, total voltage, and the voltage drop across either end of the discharge tube. Again the instability of the discharge made it impos- sible to get reliable plots of output energy, as a function of discharge conditions. However, it was necessary to have a total voltage of about 300 volts in order to maintain a continuous discharge.9 This corresponded to a current of about .6 amperes. With this and most voltages used, it was necessary to start the discharge using a hand held Tesla coil. By increasing the voltage to 650 volts, the current increased to about 2 amperes. Above this current, excessive heating of the end window assemblies occurred, and the O- ring seals failed. There was a slight increase of out- put with current. The energy output and current also in- creased when the ballast resistance was decreased from 148 to 111 ohms. Decreasing the ballast also had a detrimental effect on stability. Magnetic Field The magnetic field was provided by two air core solenoids. They were made by winding six layers of Anaconda fiber? glass coated rectangular copper tubing on an aluminum mandrel “mum uao> Ewumhm HMOHHuomHm .N mudmwm wponumo 1 v Amuououocmuv Magnum Hm3om >omn+.o 1 1 muoumammm ago em use mMU 10 and then potting them in thermo-setting epoxy. The six windings were connected electrically in series, mainly because of power supply characteristics. However, the cooling water was fed into each winding in parallel for maximum cooling efficiency. The size of the tubing was dictated by what was available, and the length and diameter of the coil was determined by the size of the discharge tube with which it was to be used. The additional re- quirements established in designing the magnet were: a field of at least 2000 gauss with a maximum current of 600 amps, maximum voltage of 160 volts, and maximum temperature rise of the cooling water of 60 degrees centigrade. The requirements indicated we would need the six layers, and could obtain a 2000 gauss field, with 200 amps and 73 volts, with the rise in water temperature being about 50 degrees centigrade. As can be seen from the data in Figure 3, the actual characteristics of the magnet, especially the electrical resistance, varied somewhat from the calculated values. However, the magnets met the nominal requirements indicated above. Initial Discharge Tube Components of the hollow cathode discharge tube are shown in Figures 4 and 5. The stainless steel discharge tube was sealed into the stainless steel water jacket with O-ring seals. The four discharge tubes were tested under various discharge conditions, to find a combination which ‘would produce a stable high intensity hollow cathode Field (Kilo Gauss) 11 _ — —Temperature / _ A . / 6) 3.0 - Field _ 60 \, m _ /' _ g / 15’ 2.0 _ /’// L *0 El / 5 — / _ E-I / x’ u 1'0 " / / h 20 J3 __.—~ " 3 ' ' u '2 0 0 a: l I T I 0 100 200 Current (Amperes) Size Tubing 3/15" sq. l/u" dia. hole Coil Length 12" Coil Diameter 81/» " O.D. by 43/. " I.D. Test results at 200 amps with coils in series Field 2600 gauss Voltage Dr0p 96 volts Inner Wall temperature 45 oC Outer Wall temperature 36 oC Exit water temperature 44 oC Water flow rate 2 gal./min. Figure 3. Magnet Characteristics 12 oposuoo onHom Hmwum mmmacwmum pwaooo mono: .v musmam maze mmumaomwo . { sew _ 411.14“: \ umxomn Hmums J L l3 mmwanammmfl moosd cam sopcfla cam .m mudmflm zma> meem moose _ pwaooo “mums v ammum mmmacflmum MM .1 1.11. uwacw Hmumz F.2b! umauso nouns x- manfimmmm scone: cam Nuumso 3mw> mos / \I/ ®\ m. k\ (I 7.\ 1 _. 1 :NWIV 14 discharge. With this setup, the 10 mm diameter tube was found to be the most useful. The 3 mm and 5 mm tubes became unstable when the discharge was increased above about 1 amp. In the 3 mm tube, the silver solder joint, which bonded the end flange to the discharge tube, failed with a current of about 1 amp. The instability of these small tubes was at- tributed, in part, to the large stainless steel surface area at the end of the discharge tube. The 15 mm tube provided a very stable discharge; but when the current was increased to obtain a higher current density in the plasma, the end window assemblies overheated and the seals failed. The end window assemblies were made from 19 mm 0. D. quartz tubing, with a 15 mm Salv-Seal quartz joint sealed into the tubing as a side port. A window was attached to one end of this assembly using sealing wax. Silicone 0- rings were used in the 3/4 inch Veeco quick couplings, which sealed the end window assemblies to the discharge tube. Silicone was used because the discharge heated the quartz tubing above the melting point of butyl rubber. The anodes were water cooled stainless steel cylinders machined to seal into the 15 mm Salv-Seal joints. The anodes were symmetrical so that the other end of the cylinder would fit into another Salv-Seal joint for connection to the gas handling system. Redesigned Discharge Tube In order to eliminate some of the problems with stability and heating, which were encountered with the original 15 design, a new discharge tube was designed and constructed. The vacuum chamber was a 4% in. O. D. by 4 in. I. D. by 28 in. long pyrex cylinder. Two holes were drilled in the walls of this cylinder for entrance and exit of the gas. A 18/7 stainless ball joint soldered into a brass flange was epoxied to the wall around one of these holes for connection to the vacuum line. A more complicated assembly for handling cooling water and gas inlet was epoxied to the wall around the other hole. This assembly, shown in Figure 6, consisted of a double walled cylinder. The outer section was Split axially into two sections, to handle both water inlet and outlet. The outer section was again terminated by a 18/7 stainless steel ball joint for connection to the gas handling system. The water cooled discharge tube, shown in Figure 7, was attached to this assembly by means of a double O-ring seal. The fine ground ends of the pyrex cylinder were sealed with water cooled stainless steel flanges, which also acted as anodes. These anode assemblies had an axially mounted stainless cylinder, terminated at Brewster's angle, for mounting an exit window by means of an O—ring seal. The magnets were designed to fit over this assembly, pro- viding an axial magnetic field. Two water cooled discharge tubes were available for this apparatus. They consisted of a 25 mm 0. D. outer cylinder, and either a 8 mm or 15 mm I. D. inner cylinder. The end flanges connecting these two stainless steel tubes 16 mwflansmmmm nmsounuwmmm ucwHooolmmo can wvocdlzopcwz cam manEmmm< nmsousupmmm ucmHooUImmw 1.1.1 limT hanEmmmfi donousupmmm mo wmcmam Eouuom mm>oouw maomno .w ousmwm OUO§I30©CH3 Ufim l7 moozumo 3oHHom pmfiooo “mums .b mnsmflm uuoe amusemo mo 3ww> mpflm “Mom Hmuucmu mo 30w> mos hanfimmmd owonumo mo soapomm mmouo .7? //////////////// //,\\ Egg/Z1. \\ as 11.1.1 1 18 were sealed with high temperature silver solder. However, the intense discharge vaporized the solder and produced a metal vapor arc when the joint was left unprotected. Fuzed quartz shields were designed to slip over each end of the discharge tube. These quartz pieces were about 18 inches long and protected the stainless steel tube inside and out, for this distance, on each end. With this design, it was found that a magnetic field of about 2000 gauss was required to maintain a stable discharge. Stable discharges in argon, helium, N2, and mixtures were obtained with currents up to 10 amps. Instability was a problem at higher currents. Other characteristics of the discharge were much the same with this setup as with the former one, except no heating problems were encountered with the present operating currents. CHAPTER I I THEORY FOR THE FIRST POSITIVE SYSTEM OF N2 Born and Oppenheimer (8) have shown that, to a good approximation, the energy of a diatomic molecule can be expressed as the sum of electronic, vibrational, and roe tational energies, viz., E = E + E + E . e v r Even when this is not sufficiently accurate, additional terms can be added to account for small interactions between the electronic, vibrational, and rotational motions. The first positive system of N2 consists of transitions between two excited electronic states, the upper state being the Bang state and the lower the A32: state. A theory for these states was first presented by Hebb (6) in 1936. Several other authors (9,10,11,12) have reported theoretical models of these states. The model used in this study, to determine numerical values and 95% confidence intervals of all the molecular constants, was based on the work of Schlapp (12) for the 32 state and that of Hebb and Gilbert (10) for 3 the H state. 19 \’ ‘ 20 Electronic Term The electronic energy of an excited state is given by Te1 = TO + AAZ where To is the electronic energy measured from the ground state of the molecule and AAZ represents the spin-orbit interaction energy for multiplet states. A is the component of the electronic orbital angular momentum along the inter- nuclear axis, and Z is the component of the spin angular momentum along the internuclear axis. Since the spin- orbit interaction energy is not completely independent of the rotational state, it will be included in the rotational energy expression. Thus, the electronic energy To will be considered simply as an additive constant for each state. Vibrational Term The vibrational energy is given by (13): 2 3 u G(v) - we(v+%) - wexe(v+k) + weYe(v+k) - weze(v+k) where wexe’ weYe’ and weze are small constants to allow for the anharmonic character of the molecular vibrations. Rotational Term 3 3 The A 2: and B Hg states have different rotational 3 energy expressions. Each rotational state of the 2 state is split by the spin-rotation interaction into three closely spaced but not equally spaced levels as shown in Figure 8. 21 .muMUm + 5 w m m on mumum m = m m Scum mcowuwmcmua Hocownmuom pascaafl .m «Human .e «m Fm a NI:M .7...“ 3M ~+:M ~+..M .H 22 3 In the H state the spin-orbit interaction splits the state into three states designated by 0 = 0, l, 2 separated by _1 approximately 40 cm . In addition, each of these states is split into two levels by A type doubling, also shown in . 3 Figure 8. Thus, for each value of J in the H state there are six energy levels. 3 A 2: Energy Levels The formula for determining the spin-split energy levels of a 32 state as developed by Kramers (ll),and later by Hebb (5) appear to fit observed data quite well for large K values. Schlapp (12) reported a more accurate formula, which takes into account the fact that the rotation does not completely decouple the spin from the axis of the molecule. Schlapp showed that this gave excellent agreement with experimental data. The energy levels are given by: I >1 I F1(K) = BVK(K+1) + (2K+3)Bv 2 2 2 g _ + +1 [(2K+3) BV + AV ZAVBV] Y(K ) F2(K) = BVK(K+1) xv + 2 2 2 1% K [(ZK-l) BV + AV - ZAVBV Y F3(K) = BVK(K+1) - (2K-1)BV where F1, F2, and F3 refer to the levels with J = K+l, K * and K—l respectively, and where A and y are constants. Assuming that the effect of centrifugal distortion on the * For K = 0 the sign in front of the square root of F3 must be reversed. 23 spin splitting will be negligible, we modified the above equations by adding a centrifugal distortion term DvK (K+1)2. A rotation-vibration interaction is implied in the above equations by the subscript v. The explicit vibrational dependence of the rotational and spin splitting constants is given by: 2 EV = Be - aB(v+%) + BB(V+%) Av = Xe - ax(v+k). Combining the above equations into a single equation (where the possible subscripts of F are indicated by the letter N) we obtain: 2 2 } FN(K) — BVK(K+1) - DVK (K+l) - (N-2){(2K+(s-2N) BV + , 2 2 2 1% + [(N-2)lv- (2K+(5-2N) BV + 1V - 213v (K+l.5 - .5)yV as the complete description of the rotational energy levels 3 of the A 2 state of N2. 3 B H State of N2 ___2 The B3H state of N2 is a state with coupling intermediate (15) between Hunds cases (a) and (b). Hill and Van Vleck (16) have considered the general equations for-intermediate coupling of multiplet states. Hebb (6) used their work to study 3H - 32 transitions and developed a formula for the A- type doubling in 3H states. This work was extended by Gilbert (12) to include centrifugal distortion terms. He specifies 24 the apprOpriate Hamiltonian as the real symmetric matrix * with elements H(H0;H0) = BV[J(J+1)+1] - AV - D[J2(J+1)2 + 6J(J+1) + 5/3] H(H1;H1) = BV[J(J+1)+1] - D[J2(J+1)2 + 8J(J+1) + 2/3] H(H2;H2) = BV[J(J+1)-3] + AV - D[J2(J+l)2 - 2J(J+l) + 17/3] H(H1;H0) = H(H0;H1) = [av - 2D{J(J+1) + 1}]{2J(J+1)}is H(H2;H1) = H(H1;H2) = [BV - 2D{J(J+1)—l}]{2J(J+1) - 4}’5 H(H2;H0) = H(H0;H2) = -ZD{J(J+1)-2}l’{J(J+1)}l5 A similarity transformation on H is then performed to obtain the diagonalized Hamiltonian HD' Here the diagonal elements of H represent the energy levels D 3 0 of the H state not including the A-type doubling. Let us again designate them by F1, F2, F3 corresponding to the Ho, H1,H2 states. Hebb shows that the A-type doubling shift of energy level F can be expressed in terms of the matrix I elements of the transformation matrix as: 2 % AVI = C0 S(HO;HI) + C1[2J(J+1)] S(H1;HI)S(H0;HI) + 2 C2{J(J+1)S(H1;HI) + 2[(Jsl)J(J+l)(J+2)]%S(H2;HI)S(H0;Hi)} ‘where C0, C1, C2 may be treated as empirical constants. * The sign of the centrifugal distortion terms has been changed to conform to modern convention. 25 3 I I I 0 Again as in the E state the v subscript implies a rotation- vibration interaction which is explicitly given by: B 2 v Be - aB(V+%) + BB(V+%) A 2 v Ae - aA(V+%) + BA(V+%) . 3 The complete rotational energy of the B H state is then given in terms of the eigenvalues of the energy matrix and the A-type doubling as: FI(J,V) i %AVI(J,V). Selection Rules In order to determine which transitions from the 3H state to the 32 state will be observed, we must consider the applicable selection rules (17). In the absence of A- type doubling, transitions from any of the 3H states to any of the 32 states are allowed as long as the selection rule on total angular momentum AJ = 0, i1 is not violated. This results in a total of twenty-seven allowed branches. These will be designated as R and PI corresponding IN' QIN’ N th 3 to AJ = +1, 0, -1 transitions from the I H level to the 3 th 2 level. The question still remains, which of the A- N 3 doublet levels of the H state is involved in a given trans- ition. The two components of the A-type doublet have different symmetry with respect to reflection at the origin. The lower energy component of the A—type doublet is +(-) when 26 J' is even (odd). For the A32: lower state, the even numbered (K) rotational levels are positive and the odd are negative. Since transitions between states of the same symmetry are not allowed, it is possible to determine which component of the A-type doublet is to be used in each tran- sition. If, for example, K is even, then we have a positive 32 state and need a negative 3H state; thus, we use the upper (lower) component of the A-type doublet, if J' is even (odd). The reverse is true for K odd. In addition, since N2 is a homonuclear molecule, a rotational state is either symmetric or antisymmetric with respect to interchange of the identical nuclei. For a given electronic state, either the positive rotational levels are symmetric and the negative are antisymmetric throughout, as is the case for the BSHg state of N2; or the positive are antisymmetric, and the negative are symmetric, as is the case for A32: state of N2. Since the nitrogen atom has a nuclear spin I = 1, the total nuclear spin quantum number of the molecule (T) can be 2, l or 0, i.e., the corresponding statistical weights are 5, 3 and 1. The even T values correspond to symmetrical states and the odd T values to anti- symmetrical states, thus, the probability of the occurrence of the two states is in the ratio 6:3. Therefore, alternate rotational lines appear with an intensity ratio of 2:1. (Herzberg (18) gives a more complete description of intensity alternation. 27 Frequency Determination Using the energy expressions and selection rules established above, it is possible to obtain a single expression for the observed transitions in the first positive system of N2. Standard spectrographic notation will be used to identify the quantum dependence of each transition. Let X denote either R, Q or P for AJ = +1, 0, -1 respectively. Let I = 1, 2, 3 represent the three triplet states for the upper (3H) level as discussed above, and N = l, 2, 3 represent the triplet states of the lower (32) level. K =-K" is the rotational angular momentum quantum number of the molecule in the final state. Thus, using single prime to represent the initial state, and double prime to represent the final state, a complete expression for the transition frequency in terms of the energy relationships develOped above is given by: I f) = FI'flufi mam NZ «0 ET» 53 a comm _ a o. co m 00 ed .4 .4 a co m 3.01 o1 31 (0,0) band with the grating turning at one degree/minute. At this speed it is possible to see the overall structure of the band and the positioning of some of the calibration lines. The high resolution infrared spectra were recorded in the first order with the grating turning at .01 degrees/ minute and the paper moving at 4 inches/minute. Thus, some of the charts were as long as 40 yards. These spectra were photographed in 15 inch sections using a 35 mm Nikkon F camera with a special minimum distortion copying lens. The photography was done on a specially constructed copying stand where the camera was rigidly mounted in such a way that the plane of the chart paper and that of the film were parallel. Bulk loaded high contrast copy film was used. In order to accurately determine the position of the infrared lines with respect to that of the fringes, the position of the fringes and infrared lines were measured on a Hydel semi-automatic digitized film reader, connected to an IBM 526 printing summary card punch. This instrument punches x—y coordinates of any desired point on the film, to an accuracy of one micron. The line positions measured in this way were then converted to frequencies by the computer program SHAFT, described below. The region from 1 to 1.5 microns contains the (0,0), (0,1), (0,2), (1,2), (1,3), (2,4) bands. This was broken up into three subregions as shown in Table 1. Three calibrated runs of each subregion were made so that the 32 phone mo nausea Home evenness xuflmcoucfi mHo: as :upa3 Mamme ATEOV eee. see. see. see. eee. nee. wee. eee. «He. A) sense .eem Neemmm. Neeemm. eHeemm. eaemee. ememee. eeemee. eeemee. meemee. meemee. 11-aoe e 1~.eeee m~.eeee mm.eeee em.e1e1 me.eeea ~e.e1ee e1.11ee m1.11ee No.11ee . 11-5u1 < em em eN 11 11 11 e N1 N1 moses e zeHe1 eee eee ewe emo eee eom eHe eeo eeo omeeeo> eo\m\~1 ee\em\11 ee\e1\11 eo\em\11 ee\em\11 eo\e1\11 ee\~\me ee\1m\11 ee\11\11 me ‘ zee H-50emm1.ueeee H-eueeee.ueeee 1-5ueeeeuneeee zeHeme same $.11 3.51 3.3 2.8 8.8 82% . .mnuowmm UTHMHMGH Hmmz Tau mo moepmflnmuomumsu HmucmEHHmmxm .H wanes 33 observed frequencies could be averaged to reduce experimental errors. Because of difficulities with the Hydel system, and to guard against large random errors by the Hydel Oper— ators, each record was measured twice, usually by different operators. Thus, there were from two to six measurements of each frequency depending on whether or not it appeared on all three charts. The tabulation and averaging of these multiple measurements was also performed by SHAFT as described below. TEES. This program was specifically designed to be applicable to the work performed in this laboratory. Its purpose is to perform, as completely as possible, the task of reducing raw Hydel data to a tabulated output of the identification average fringe number, weighted average frequency, weighted standard deviation, individual measurements and coded weights, for up to eight measurements of each line in a band. In addition, the program generates the frequency and quantum dependence inputs for any of the many analysis programs used in this laboratory. The general concept for forming this program was to minimize Operator inputs, but at the same time, to cover all possible Options the operator might desire with the maximum data capability and the minimum time requirement. Quite obviously, the final outcome was a compromise among these goals. The versatility of tasks was accomplished by means of a list of option para- meters. Overlays and large multi-purpose variable dimensioned 34 data blocks increased the data capabilities. The run time was minimized by ordering of the data, so that repeat searching was kept to a minimum, and also by careful planning of the sequence of operations. Figure 10 shows the overall structure of the program. The main program SHAFT is a driver program which contains the required common storage blocks, and controls the overall operation of the program. Program IDENTIFY reads in and stores information for identifying the measured frequencies. -ThlS information can be used for calibration purposes, as well as for identification. SCAN is a program (written by M. D. Olman (23)) that uses line positions measured on the Hydel, to determine the positions of the spectral lines relative to the fringes. AVFRQ is a tabulating and averaging program. It takes lists of information provided by the other programs or read in directly, and performs the necessary weighted averaging. Subroutine IIDENT compares the average frequencies to the identification list. Function CWT converts the code weights to actual weights. Program SEARCH reads in two or more initially identified calibration lines, and using fringe numbers generated in SCAN along with calibration data read in by IDENTIFY, calculates initial values of the calibration constants. These constants are then used to convert the remaining fringe numbers to frequencies. As each frequency is formed, it is compared to the calibration data to locate additional calibration lines, which are then used to update the calibration constants. When all the Rheum mo smumswo xooam .oa musmflm 1» 112.25 _ 111.201 - _ m 18:26 N 18:25 _ 112.26 .. 110111131 om..._>< 2481 11.1.1219 3 1 1 1 1 5960.5 522 .1....1mo cvq'a ooo enee.o cmeo.o enee.e 0.0... 0.0... ocea.o COO... 0000.0 0.0... 0.000. 000... 000.00 000... .0000. 00.... mueo.e 0...... pcmo.o 00.000 00.... .00... 0.0... 000... c.0006 race.e o~m~.o wcoe.o so 223. >um t‘ «Hat. DO Oinqfl ifll‘ChflC) (\C’ficsvcd sen page.e «nee.eu nuac.e anon.e coa>.o cacao. oece.e cocoa. naro.e «ano.o ace-co unc\.e oonn.o anvv.o coco.: coda.» coo... cave.e seso.e oooo.a coo... COO..- se.e.« sauv.e «cso.e canv.o oms~.o conv.o he (a: e—wax v‘fll' dandn ' vi... IUCMHIDvOd amaze scum useuso ememaenms madamm nee oeae.e eeee.e coco.o earn.o eeeo.o 0.000. «116.: .00... eeeo.o no.9.o .00... ~ecs.o nxvo.o ca...e «neo.o noea.e COO... nona.o ammo.e anec.e 00000. 0C0... eoao.e enav.o oaco.o eonv.o «uc~.o eecv.o so 132 )wa itarOO 2 : 'QO o1; . ~.1e.e P033910 us~¢.e Saar.e eeoe.e COO..- coo... oovu.e neno.e nxomas ruwone cc...- c.0000 0CD... '00... CO...- aaaaau COO... mound: eye... sans.) ween.e 30.9.9 .00... [OI-l. smut"; aan~.e Oates. xc 2:: >134 3 q 2 c '0 a a A "£02.31...sz me, 1oae.e snea.ns myac.e reen.o o~¢t.o .0000. one... sewn.e mono.c cauw.n weav.o 000... any... .09... IOU... COO... roeu.o FCC... swan.o Clo... (3v=.o watn.n tur>.u scene. 90.... mnfiv.o anwxuu one... No :5: sin .zu::. Tnam 11 year 1:» mm a» rwxamma mac .m means .7 .mcz macs ree.e 1:12_ urn" u: 1:: ach fleas. omee.e ”wafloi fuse-z .vsae.: 0591.: 0H=:.1 amae.. vans. 003:. are... m\:;? naae.. oveu. K$ULI ease. n:a1.. ”Laue :oade newe.1 mea;.1 menu.1 arose. «Hae.e mefl..g «sage; «swoon -- - - S. b- - ‘v . --b-u§- .sm: .apm ”(or 320: Stat 91.: fic.c 92.: ac.; one; ec.e ne.1 aw.n Wee; oe. ec. 0e. 0e. ac.e or.» oe.c an.e .143 nenv.snpm nucu.on.a vrmo.mfler swouofinrm mfintevfiem tcmv.nnbm «encownrm \emneflnem ecmn.o~rm ommhohwfm vacovarn cenv.m~2e acttoflflfim OCVOoOHrm xuea.cacm eeoo.edne eaea.efl:e meec.canm aece.ma.o o:mw.v«ze nmoo.na:n memv.~«hn ownc.rorc meec.veum oaon.mono ~wxv.mo~o >uzmacwum woa v— nm~u_FZHe. nwxrmcmr :a1ea1swr u~ao2 Hum wmm mmmmw mozmm ommmmmzH mamz «Baa msmHo Nz mo mnmum m < on» MOM :80 CH muamumcoo mSOOGOUHOEwm .m OHQDB m H 67 o NH.on mso.on «H.on «H.on maooo.on mmoo.on moo.on -osxmn.on .uoaxv.mn «mooo.on 5000.0” smo.on ma.on mmo.oH Hum wmm mDde e m was new n w h «woo.o- mo~.s oe~.ma hm.smha maooo.o mmoo.ou mm.~ nonxom.m- noaxmm.flu mamao.o ommm.n moo.o ssoo.o Hmm.mv me¢> QmmH mooo.o mmao.o HH.N -oaxmm.m 0 nl oaXmm.ml namao.o ommm.a mmoo.o mmma.o mmv.m¢ mmmmw «Eda mmmHQ .v magma 68 concerning the relative accuracy of these energy levels was not given, and therefore they were all weighted equally. Using the 12 molecular constants for the model Of the A32 state, presented in Chapter II, 436 energy levels were fit to a standard deviation of .0212 cm-1, and all Observed values were within .06 cm.1 of the calculated values. Dieke and Heath indicate their frequency determinations are accurate to approximately .01 cm"1 thus the fit is good to the same order of magnitude as the experimental error of the data. In the final iteration for the constants the predicted corrections were all at least an order of magnitude smaller than their respective standard errors. The resultant constants and corresponding 95% SCI's are given in Table 3. These constants include a the vibrational dependence 1’ of the spin-splitting constant A. Since this vibrational dependency is only about three times its 95% SCI, it would be difficult to Observe in a study using a limited number Of vibrational levels and/or a less accurate molecular model. Dieke and Heath determined the constants given in Table 5 using combination difference methods on one band at a time. This simultaneous analysis includes D and 7, which were not determined by Dieke and Heath. After correcting the few obvious typographical errors in the data of Dieke and Heath and eliminating the F1 (17) level for which the Observed minus calculated values was _1 -8.7 cm , there remained 1138 energy levels to fit to the 15 3 molecular constants of the B H state. As was the case with 69 Nbaoo.o mm.ma om.om¢a em.al muoaxve.w- mtho.o mvmv.a flammav mMmHo omoo.o- moo.ou mm.H- Hm.H- noaxm.m -oaxm.m m m mmao.o mao.o omme.a osv.a Ammaac Aammav qqommmo mosmz .Oudumumuflq map Eoum NZ «0 mumum w « may now ISO ca mucmumcou HMHDOOHOE .m OHQMB m H 70 ml mmao.o th.o mmH.mH Nv.mmwa Ioaxm.m o oaxmm.h vmhao.o mmao.o wvnmm.a wmmw.a Ammmav Ammmav mmmH Q Aqommdu mumnm e m may now m H mmooo.o mmoo.o n~.m suoaxm.m suoaxm.m mwao.o mmo.o memm.a om.H N.Ne m.~e remade immmav flammai mmmm comm mosmz mudumnwqu map Eoum Nz mo :80 ca mucmumcoo HMHDOOHOS .m OHQOB 71 the infrared data, a systematic difference between the Observed and calculated results was found. The average deviation of the Q = 1 level is +.155 cm-1 while that of the Q = O and Q = 2 levels is -.157 cm“1 and -.001 cm"1 respectively. The overall standard deviation Of the fit is .198 cm"1 and the molecular constants are shown in Table 4. To a first approximation the spin-orbit interaction constant A gives the separation Of the F1, F2 and F3 levels. In Figure 14 the solid line represents the theoretical energy for the F1, F2 and F3 levels. The broken line shows the deviation of the observed data from this theoretical line with the scale increased by a factor of 100. This shows a systematic deviation in the separation of the levels which decreases with increasing rotational energy. It thus appears that the asymmetry of the F1 to F2 and F2 to F3 separation is not adequately described by the theory. Figure 15 shows the Observed minus calculated deviation for the three levels with V = 4. All other vibration levels showed the same general trend. The larger number Of degrees Of freedom in the simul- taneous fits resulted in a different distribution Of the deviations in each vibrational level. However, for each vibrational level at least one of the following two state- ments was found to be true: 1) The observed energy sepa- ration of the F1 and F2 levels was greater than predicted. 2) The observed energy separation of the F2 and F3 levels was less than expected. It is not possible to improve the fit by varying A. Energy Above F2(0) of A 2 (cm 72 16800 ‘ Theoretical Energy --— Observed Energy Deviation (Scale increased 100 times) 16600 ' 16400 ‘ 16200 ' 2 4 6 8 10 12 14 16 18 Rotational Quantum Number (J) 3 Figure 14. B H Energy Levels for v=4 lower A-Doublet 73 umandonl< HUBOH ¢n> How wGOflumH>ma mmumcm m m .ma m any Honesz Esucmso accoflumuom musmflm «a ma NH HA OH m m b m m v N F _ _ _ L _ t _ u _ _ l + + + + + L... + + + + + 4 i + O O 0 1H. AU 0 Hm>ma mm + o 0 V IN. Hw>0H .m o O O Hm>ma 1m < o O mo) Abieug penetnotea-peAxesqo (T- 74 3 Gilbert (10)has calculated the energies of the H state of PH and compared these to the values determined from the A 3400 band (3H + 32 transition) data reported by Pearse (26). He used a combination difference method to determine these energies. The deviation of the experi- mental and theoretical results are shown in Figure 16. As with the N2 data, the experimental data indicate a larger asymmetry in the spin-orbit interaction than the theory predicted. Professor I. Kovacs (27) has indicated that these deviations might be due to either a spin-spin interaction or a spin-rotation interaction. The above discussion would seem to indicate that it would be difficult to obtain a reliable value for Av and therefore to calculate the vibrational dependence of the constant Av‘ However, the change in Av as v increases is quite large, and although GA and BA are probably not as accurate as their 95% SCI's indicate, a decrease in A as v increases is certainly apparent. BudO (9) also determined AV for several vibrational levels between v = 0 and v = 12 and although his values are not monotonically decreasing, his result A0 = 42.3 and A12 = 41.3 also suggests a vibra- tional dependence of A. When theory is incapable of producing agreement to within the experimental error Of observed data, the results of a least square fit based on this theory must be interpreted ‘with caution. As stated earlier the 95% SCI's are calcu- lated assuming one sampling of a normally distributed 75 .uumnawo an pmuuommu mm no team oovm< Dav Eoum mcowumw>mo .OH Ousmwm any umnasz Enhance Hmcofiumuom NH Ha ca m m b m m e m N H ¢. *1 _ P t m A) _ _ _ L a a. Q. Q Q G. 0 no 0 O o Cos no no q + .+ + + q. .+ + + Hm>mH mm + Hm>wa mm o +. Hm>ma is < +. m.l v.1Mw S e . M m I e m. D N I w o n . a H I 1 e p . 3 o u e 1 5 H. .A m, m N. . m. v. 76 infinite pOpulation of similar sets of data. When this is not the case the 95% SCI's give sgmg_indication of the match between the theory and experiment but it does not mean that there is a probability of .95 that the estimators are within a range of i SCI of the true values of the constants. Tables 5 and 6 Show the values of the molecular constants of the A3Z and B3H states of N2 obtained by previous workers. These constants were found using combi- nation difference techniques On one band at a time. For a simultaneous fit different results might be expected even- when the theory is adequate. When, as in this case, the theory has shown to be inadequate, one might expect even larger variations in the results of a single band and simultaneous fits. Tables 3 and 4 show for the first time in a single work all the constants describing the first positive system Of N2. S UFDLARY AND CONCLUS 1 ONS A newly designed double-ended hollow cathode has been develOped and stabilized for use as a spectral source. It has been used to excite the Np first positive system in the near infrared. Computer programs for efficiently handling large quantities of spectrographic data have been develOped, tested and described in detail. These include a program for calculating and identifying transition frequencies and programs for non—linear least squares fits of the experimental data to the existing theory Of the A32 and B3H excited electronic states comprising the first positive system of N2. A method of using punched cards to communi- cate between the fitting and identifying programs has been described to make the analysis semi-automatic. This process could be improved and further automated by elimination of the punched card interfacing. The reported analysis of the first positive system of N2 represent the first case in which all the molecular constants have been determined from a single set of data. This work also constitutes the first reported simultaneous analysis of this system. The analysis has indicated an inadequacy in the theory of 3H states. The 77 78 major portion of the discrepancy appears to be in the description of the asymmetry of the energy separation of the F1, F2 and F3 levels. LI ST OF REFERENCES 14. 15. 16. 17. 18. 19. 20. LIST OF REFERENCES I. Agarbiceanu, A. Agafitei, A. Preda and V. Vasilio, Rev. Roum. Phys. 11, 649 (1966). D. A. Huchital, and J. Dane Rigden, IEEE J. Q.E. QE-3, 378 (1967). P. K. Carroll, Proc. Phys. Soc. 54, 369 (1962). G. H. Dieke and D. F. Heath, John Hopkins Spec. Report #17, (1959). M. W. Feast, Proc. Phys. Soc. 63, 568 (1951). M. H. Hebb, Phys. Rev. 42, 610 (1936). S. M. Naude, Proc. Roy. Soc. 136, 114 (1932). M. Born and R. Oppenheimer, Ann. Physik. 84, 457 (1927). A. Budo, Zeits. f. Physik. 26, 219 (1935). Cecil Gilbert, Phys. Rev. 42, 619 (1936). H. A. Kramers, Zeits. f. Physik. 53, 422 (1929). Robert Schlapp, Phys. Rev. 51, 342 (1937). G. Herzberg, S ectra of Diatomic Molecules, 2nd Ed, D. Van Nostrand, (I950). Ibid., p. 106. Ibid., p. 231. E. Hill and J. H. VanVleck, Phys. Rev. 32, 250 (1928). G. Herzberg, Op. cit., pp. 240-245. G. Herzberg, Op. cit., pp. 133-135. J. L. Aubel, Thesis. M. S. U. (1964). D. B. Keck, J. L. Aubel, T. H. Edwards, and C. D. Hause, "Symposium on Molecular Structure and SpectroscOpy", Paper H-2, Columbus, Ohio, (1966). 79 21. 22. 23. 24. 25. 26. 27. 80 D. B. Keck, Thesis. M. S. U. (1967). K. N. Rao, C. J. Humphreys, and D. H. Rank, Wavelength Standards in the Infrared. Academic Press, New York. pp. 160-1627 W TI'_)"966 .— M. D. Olman, Thesis. M. S. U. (1967). C. D. Coleman, R. Boxman, W. F. Meggers, Table gf Wavenumbers, National Bureau of Standards Monograph 3, p. III (1960). J. Mandel, The Statistical Analysis of Ex erimental Data, Interscience Publishers, New York, pp. 147-150 (1964). R. W. B. Pearse, Proc. Roy. Soc. A129, 328 (1930). I. Kovacs, Dept. of Atomic Physics, Polytechnical Univ. Budapest, Hungary. Private Communication. APPENDICES MAIN C S 10 C C C 15 20 C C C 25 APPENDIX A FOPTPAM LISTING OF PROGRAM SHAFT ppOGQAM SHAFT SCAN HYDFL AVFWAGF FPEQUFNCIES IARULAIE COWMON /l/ IOP(3)9IL§(R)QTOLIQTOLgoTOL-‘IGAORQLDOIDNTQICAL. IIIONJQJLNWD COMMON /2/ NoVISAQLF COMMON /3/ NQMFI(2000)CNAMEZ(POOO)9NU(?000)9NMAX9TOL4 COMMON /S/ MUM COMMON /R/ ICHEADIQ)QIHFAD(992)CIDN(4) COMMON /Q/ CCWI(RQQ) DIMENSION FPQ(1910000)0 WTS(1010000) INTEGFQ WTI9WT9WTS¢VISAQLE QEAD c509 NGPSOIOPOIDNTQICALQ IDONQVISARLFQIOLI cTnL70T0L39 l TOL4-IFDONE IF (IFDONF.FO.8H STOP) STOP IDN(1):RH IDN(?)=RH IDN(3)=RH IDN(A)=8H JUMP:0 NGPSZ:HGPS IF (ICAL.FQ.O.OQ.ICAL.FQ.3) GO TO 10 NGPSP=NGDS+2 IDEM:10000./NGPSP READS IN CALIPQATION OP IDFNTIFICATION DECKS IF (IDNI.FO.2.OP.ICAL.FQ.2) CALL OVFPLAY (1.9109) IF (IDNT) 90915920 Lp=l NAMF?(1)=RH NAME1(])=NAMEZ(I) CONTINUE N21 D0 25 I=I9NGPS GIVES TOLFQANCFS PPFSET VALUES IF THEY ARE NOT PFAD IN IL5(I)=1 IF (TOL4.FQ.O.) TOL4=1. IF (TOL1.F0.0.) TOLI=.050 RI 000 000 8? IF (TOL?.FU.0.) TOLP=TOL1 IF (TOL3.NF.O.) GO TO 30 IF (ICAL.NF.0) TOL3=.03 TOL3=.04-TOL3 CALLS SCAN 30 IF (IOP(I).NF.0) CALL OVERLAY (29 9P9J9FRQ9WIS9NGDSo 1 NGDSP-IDFM) CALLS AVFRQ 35 IF (IOP(2).NF.0) CALL OVERLAY (39 939J9FPQ9WT89NGPS9 1 NGPSP9IDEM) IF (ICAL) 40945940 CALLS SEARCH AND CALFII 40 CALL OVFDLAY (49949J9FQ09NGPS?9IDEM) JUMpzl ICAL=0 IONI=1 GO I0 35 45 RFAD 559 IFDONE SKIPS NEXT CASF IF SFI UP WPONG IF (IFDONF.FQ.RH ) GO TO 45 IF (IFDONF.FQ.RHEND HEAD) GO TO 45 IF (IFDONE.FQ.RHNFW CASF) GO TO 5 IF (IFDONF.FQ.8H STOP) STOP 50 FORMAT (9II9?X94F10.R9??X9AR) 55 FORMAT (72X9AR) END OVERLAY 1 000 PPOGQAM IDENTIFY PFADS IN CALIRHATIOM OP IDFNIIFICAIION DFCKS COMMON /I/ 1013(3) 9ILS(8)9T0L19IOL29IOL39A9ROLP9IIINI9ICAL0 1 IDON9JUMP COMMON /3/ NAMEI(2000)9NAMF2(?000)9NU(?000)9NMAX9IOL4 DEAL NU INTEGFQ OQDFQ N=1 NTzl ASSUMFS OPDFPS APF I UNLFSS GIVFN 83 IICAL=1 IP=IICAL S DEAD 3§9 INAMFIOINAMCPQXNUQURDFD IF (OROFQ.FO.RH OQOEP) GO TO 10 GO TO 15 10 IR=IMAMFl-8H 0 IICAL=INAMF?-RH 0 GO T0 5 15 CONTINUE C C CONVFQTS FQFQUENCIFS TO PPOPFR VALUF DEPENDING 0N OPDFPS USED C XNU=XMU*IQ/IICAL IF (N.EO.1) GO TO 25 C C OPUFPS FREQUENCIES LOW TO HIGH C 20 IF (NU(NT-l).LF.XNU) GO TO 25 NU) FORMAT ( 379 AVFPAGF DFVIATION OF LINES FROM F1=9FR.49 1 4H F2:9FR.49 4H F3=9FR.4/9 21H NUMBER OF LINES FROM9 2 13X9 3HF1=9IQ9 4H F2=9I99 4H F3=9I81 FORMAT ( 9H Y.7Il7./) FORMAT (X97El7.69F16.6.//) 103 FORMAT (11097Il79/1 104 FORMAT(*OSTAT FROM LINFS WHTO LT 10.0“/*O*2SX*STANOARO* 123 1 *OFVIATION* 2F11.4.19XF11.4) 10$ FORMAT ( 16H pAw SSCP MATRIX9/5126) 106 FORMAT (X.1?.SE?6.101 107 FORMAT ( 3H Y.SE?6.101 10R FORMAT ( 8H Y VS Y.E?6.10) 109 FORMAT ( 21H-PARTIAL CORR. COFFF.9//1X.ISIR) 110 FORMAT (I3.16FR.2) 111 FORMAT ( 3H Y.16F9.2) 11? FORMAT (9 F LEVEL OF x-*I?.F10.2.QX*STO OEV OF (O-P)* 1 F12.4//10X *VARIARLF COFFFICIFNT STO FRROR*/) 113 FORMAT (Ill.A0.?Fl7.101 11467 FORMAT(*1LFAST SOUARFS FIT OF ”EAR /*FIT “14* OATIHt 1 ”FOINTS TO * I?” VARIARLFS*/* DELETES UP TO *1?“ POINTS* 2 *IF (O-P) IS GT “F6.3/* WHT NORM * F6.P/§ STD DEV OF“ 3 *(O-C)*FR.4/* COMPLFTFO *12* STEPS*) 11S FORMAT ( RHOSTFP N0913/9 ISH VAR. RFMOVEO.I3) 116 FORMAT RHOSTFP NO9I3/9 lSH VAR. ENTEREO9I3) 117 FORMAT 19H FRROR QESIO SO VAQ.IB9 7H IS NFG) 119 FORMAT 4POVAR9IS9 9H 18 CONST) 119 FORMAT IOHOFRPOR. VMIM POS) 120 FORMAT ISHOFRROR NOIN NEG) 121 FORMAT 13HOY QOUAQF NFG STEP.IS) 12? FORMAT P2H07FRO OFG FRFEOOM STFP.13) 123 FORMAT IOH SQUARF X-.T2. 14H NFGATIVF STFP9I3) 124 FORMAT QH Kzfl STFP9I3) 12S FORMAT *0 OIAG FLEMENT§*/* VAR NO VALUF*//(14.Dl4.4)) 126 FORMAT RH LINF NO.I4. 23H BFING OELFTED FROM FIT) ENO SURROUTINF SOP C THIS ROUTINF NFFOEO TO KILL OVER OPTIMZATION FNO AAAAAAAAAAA SURROUTINE RRNTOUT SURROUTINE PHNTOUT (DATA9FOMC9FCALC9WT9NUSE9NCONP19INOATA) COMMON /17/ IRPS(SOO)9IS(ROO)9IN(SOOIQIK(SOOIQIHEAOIQIO VU(15)9VL11819NVSV115)9FORSISOOO)OITQYS9ITQCN19 IQAQPAYI1PO)9ISAQRAY(120191NARPAYI1ZO)9 IKARRAY(1?O) COMMON /1?3/ NK12610NDELMAX9XOEVMAX COMMON /23/ NOOATA9NOVMI9AVEWHT9STDY9NOSTEDOIWOCNT’ IWTN9FQFQDRO9YDPEO9DEV9N9MS DIMENSION DATA(NCONDI)9 WT(INDATA)9 NUSEIINOATA) ENTRY DRINT? TFOONanH RRINT 30 IHEAD PQINT 3 RETURN ENTRY PRINT3 FDRPO:FORS(N)-FQFQPDO FRPMczFomc—FRFODRO FRCALC=EOQS(N)’FOMC .JNF‘ 124 IF (ITRCNT.LT.ITRYS) GO TO 1 IF (N.FO.IHOATA) IFOONF=SHFNO OATA IF (IKAQRAY(IP).FO.4S) TKARHAY(IP)=O WRITE (9294) VUIMS).VLIMS).IHARRAY(IP)9ISARRAYIIP). IINARRA((IP191KARRAYIIP).FORS(N)9WT(N).FOMC9IFOONE 1 CONTINUE PRINT S9 N9VU(MS)9VL(MS)9IRARRAY(IR)9ISARRAY(IP). lINARRAY(Ip).IKARRAY(IP).WT(N).FORS(N).FRCALC9FOMC. P FRFOPROoFRPQO9FPPMC RFT-JRN 2 FORMAT ( 1H199AS/1 3 FORMAT ( 65” N RANO TRANS. WT ORS FQFO CALC FREQ“ 1* ORS‘CALC OPS-RQF0916X9 ZOHPRED FREQ PRED-CALC9/1 4 FORMAT (2X91191X91192X991971191X91296X9E13o494x9FR939 14X9F12.49I?X9A8) S FORMAT (1X9149 ?H (9119 1H99119 2H) 99192119 1H(9 1129 1H)9F6.392(FQ.39?X)92(Fgo492X)914X92(FQ.492X)1 EMT) SURROUTINF E I GEN SURROUTINE F IOEN (A9VALU9N9N9IA9‘3) FIGFNVALUES ANO EIGFNVFCTORS OF A RFAL SYMMETRIC MATRIX 001’.) OIMENSION A(IA91)9 R(IA.1)9 VALU(3)9 DIAG(3)9 SUPFRO(2)9 1 O(2)9 VALL(3)9 S(2)9 C(219 O(3)9 IND(3)9 U(3) CALCULATE NORM OF MATRIX F300 1 ANORM2=O.O OO 2 I=I9N OO 2 J=19N P ANORM22ANOWMP+A(I9J)**? ANORMzSQRTIANORM?) C GFNERATE IOFNTITY MATRIX IE (“I 39793 3 OO 6 1:19N OO 6 J=19N IF (I‘J) 59493 4 RII9J1=190 GO TO 6 S 9(I9J)=O.O 6 CONTINUE C ppppnpm QOTATTONS TO RFOUCE MATRIX TO JACOHI FOPM 7 IEXIT=1 Mszxl-P IF (NN) 6S91598 (‘30 10 11 12 IS 16 17 19 19 20 21 23 24 PS 26 12S OO 14 I=I9NN 1131+? OO 14 J=II9N T1=A(I9I*1) T2=A(10J) GO TO 66 O0 IO K=I9N T?=COS*A(K91+11+SIN*A(K9J) A(K9J)=COS*A(K9J)-SIN*A(K9I+1) A(K9I+1)=T2 OO 11 K=I9N T2=COS*A(I+19K)+SIN*A(J9K) A(J9K)=COS*A(J9K)’SIN*A(I*19K) A(I+19K)=T2 ‘ IF (”1 12914912 OO 13 (=19N T2=COS*Q(K9I*1)*SIN*S(K9J) Q(K9J)=COS*S(K9J1‘SIN*R(K91*1) RIK9I+1)=T2 CONTINUE MOVE JACORI FORM FLFMFNTS ANO INITIALIZF EIGENVALUF ROUNOS OO 16 I=I9N OIAG(T)=A(I9I) VALUIII=ANORM VALL(I)=-ANO9M OO 17 I=P9N SHDEQOII-1)=A(I’19I) O(1’1)=(9UPFQO(I-I))**? TAU=0.0 I=1 MATCH=O T2=0.0 T1=I.O OO 33 J=19N D=OIAGIJ1‘TAU IE (T2) 19921919 IF (T1) 2O924920 T=P*T1-O(J-1)*T2 GO TO 2% IF (T1) 22923923 71=-1.o Tz-p GO TO 23 T1=1.0 sz GO TO 23 IF (OIJ’111 2§97397S 1F (T?) 27926926 T=-1.0 (705') 29 P9 30 31 32 33 34 3S 36 37 38 39 40 41 42 43 44 45 46 47 49 49 60 SI S2 S3 126 GO TO 29 T=1.0 COUNT AGREEMENTS 1N SIGN IF (T1) 90.29.29 IF IT) 32.31.31 IF (T) 31.32.32 MATCH=AATCH+1 T2=TI T1=T FSTARLISH TIGHTFR gnumns ON EIGFNVALUES OO 33 K=I9N IF (K—MATCH) 34.34.36 IF (TAU-VALL(K)) 3R.3R.3S VALL(K)=TAU GO TO 3Q IF (TAU-VALU(K)) 37.38.38 VALUIK1=TAU CONTINUE IF (VALUII)-VALL(I)-S.OF—R) 42942.40 IF (VALU(I)) 41.43941 1F (ARS(VALL(I)/VALU(I)-l.0)-S.0E-R) 42.42.43 I=I+1 IF (I-N) 39939.44 TAU=IVALL(I)+VALU(I))/P.0 GO TO IR JACORI EIGENVECTORS RY ROTATIONAL TRIANGULARIZATION IF (M) 4§9GC§94S IEXIT=3 OO 46 1:19N OO 46 J=19N A(I9J13090 OO 6? 1:19N IF (1'1) 50990947 IF (VALUII’I)-VALU(I)-q.0E-7) 56956949 IE (VALUII-111 499S0949 IE (ARSIVALUII)/VALU(1-1)-190)-§.0E*71 56996950 COS=1.0 RETA=0 SIN=0.0 OO 54 J=19N IE (J-I) S39S39S1 GO TO GO S(J-1)=SIN C(J“1)=COS O(J*1)=T1*COS+T2*SIN TI=(OIAG(J)-VALU(I))“COS‘RETA*SIN T2=SUPEQOIJ1 £39 94 57 SR S9 60 61 62 C C C 63 64 6S C C C 66 67 69 99059 127 DETA=OHPFRQ(J)*EOS O(N)=T1 OO SS J=19N 10U3(\I)=(1 SMALLO=ANORM OO GO J=19H 1E (II‘JO(J)-1) 579599S9 IE (ARS(SMALLO)’AOS(O(J))) 59959958 SMALLO=O(J) MM3J CONTINUE INOU‘JN) =1 DROOS=1.O IE (NN-I) h2962960 OO 61 K2794N II=NN+1-K A(11+19I)=C(TI)*PROOS PROOSz-PQOOS*S(II) A(19I)=DQOOS FORM MATRIX L’ROOUCT OF ROTATION MATRIX WITH JAFORI MATRIX Or) 64 11:1951 OO 63 K=I9N U(K)=A(K9J) OO 64 I=l9N A(IQ.J)=OOO OO G4 K=I9N A(I9J):R(I9K)*U(K)+A(I9J) QFTURM CALCULATE SINF ANO COSINE OF ANGLE OF ROTATION IF (T2) 67969967 T=SORT(T1**2+T2**2) COSle/T SIN=T2/T GO TO (99SP)9 IFXIT GO TO (1496719 TFXIT ENO AM 0391 DROGRAM 0391 (INPUT9OUTRUT9TAPE1) COMMON 41(3.3)9H2I3o3)oFPII3).NUSFI130019CONIIF).FOMC. 1WT(SIO)9CONSTSTNT(1S)9FCONS(15)9SCON(15) COMMON /A/ NAMFIIS) COMMON /TFMP/ ONU(130019FONU(1300) COMMON /1P/ INFOI9INF029INFO39INFO4.NUMCON9JUMP9CONST(1S)9 IEEIN9EEOUT9TOL9IOENT1219TOLF COMMON /17/ 18(13001.INI1300).IJ(1300)o1HEAOI9).VU(IS). INVSV(16)9RS(13OO)9ITRYS9ITRCNT 0(3F) IPR COVMQN /?l/ NT(16)9VFCTOP(17017)9INDEX(16)oIDFX(15)- TSIGMA(16)~CUFN(16)9910MC0(16)9NOTIN(16)9XCONST(16)oNOIN COMMON /123/ NK(lq)oNDFLMAXoXDFVMAX DIMENSION FUT(R)- INFIT(2)9 COD(3) DIMFNSION AFHN(3)~ DFUN(3) VQUIVALENCE (COHQTNT( 1)9 API)9 l (COM§TmT( 3). UPI). 2 (CKNTSTNT( ‘3). C1)~ 3 (CONgTNTT 7)! 2590?). 4 (CONSTNT( 9). XU)’ G (COM§INT(ll)o 2U). 7 (com§TNT(lS)o RAHPI) TNTEGFQ VU INITALIZES NAMFQ AND VAWIARLFS TNFTT(1)=3H M0 IMFIT(9)=3HYFS NAMF(1):QHFA91 NAMF(2)=4HFRDT NA”F(3)=4HFUDT NAMF(4)=3HFCO NAMF(§)=3HFC1 NAMF(6)=3HECB NAMF(7)=6HF2FPOF NAME(R)=3HFWU MAMF(Q):}HEXH MAV€(10)=3HFYU NAVF(11)=3HF2U NA”E(12)=6HFALADT MAMF(13)=6HFHAAPT NAMF(14)=6HEALRDT MAWF(15)=6HFRARDI MUM-CON: 1‘5 NC0N01=JUMCOH+1 MVSV(1)=0 JK=1 KI=O JHMD=0 TRAnzo FFRTzo. CONTINUF ITPCNTz0 DFAD§ TV AND DRINTS OUT HFADING DFAD 26o IHFADqIFDOMF IF (IFDONE.FQ.HH STOP) QTOP DPINT 269 IHFAH IF (IFDUNF,ME.RHFMH HEAD) GO TO 1 QFAD§ ODTION CAQD (CONQTNT( 2). RpI)9 (CONqTNT( 4). C0)9 (CONgTNT( 6). C2)9 (CONQTNT( R). WU). (CONQTNT(10)¢ YU). (CON§TNT(I?)o ALAPI)9 000 F)O:O 129 DFAO 279 IOFMT9ITWY§9INFO19INFO29TNFO39INFOQ9(MK(T)9 1T=191C5) 9"“)FL‘J‘AX0XflFVMAXOTUL'TOLFOEFINQFFOUT’ICASF ASSIGN A TO MGO IF (ICASF.F9.RHNFW CASE) GO TO 2 IF (ICASF.F0.QHHFW CONS) ASQIGN 13 TO “GO IF (ICASF.FO.RHADU CONS) GO TO 13 CONTINUE QFAOQ In CUNSTANTS AND FRACTIOMAL CHANGFS QFAD 299 COM§THT9OFLAP19DELRDI DFAU 789 FCOWS DQINT 999 CONgTNT TF (FCOVSTI).FQ.0.0) FCON5(1)=1.0 O0 3 TI=791§ §FTS FRACTIOMAL CHAMGF IF ANY OR ALL NOT DEFTNFO TF (FCON§(TT).F0.0.0) FCONSTII)=FCONS(1) C(VJTIAMJF an TO “809 (1394) CONTINUE TJCT=0 IT=IJCT J=IT COMTTNUF DFAOS IN EMFRGTFS QFAU 309 FUT CFVJTTKNHT TF (FUT(R).NF.RRRR.QRRH) GO TO 7 VU(JK)=FUT(1) NVSV(JK)=KI DPTNT 319 JK9VU(JK)9NVSV(JK) JK=JK+1 K120 J=0 TJCT=O QFAL) 309 FWJT TF (FUT(9).FO.QQQQ.QQQQ) GO TO 11 GO TO 6 CONTTNUH O0 10 KRzloR CALCULATES QUANTUM OFPFNOENCIES KT=KI+1 IT=IT+1 IJCT=IJCT+I 1F (IJCT.L£.8) GO T0 8 IJCTzl J=J+1 10 11 19 13 14 19 130 TF (J.NE.0) GO TO Q TT=IT-1 KI=O GO TO 10 COMTTNHF FO*%S(IT)=FUT(KR) T§(TT)=(TJCT+1)/? TJ(TT)=J L‘IT(IT):1. IF (FUT(KR).F0.0.) WT(IT)=O. C(WJTT’WJE CO TO 5 CONTTNWF ITM=TT OO 12 M31015 §C(‘)T\T(M):Oo CONTINUF PQIVTS INITIAL CONSTANTQ ANn VAPIARLES DPIMT 329 ICA§F9ITQYS9FFOUT9FFTN9TOLF9(NAME(M)9CONSTNT(M)9 lFCON§(W)9TNFTT(NK(M)+1)9M=l915) CONTINUE JK=0 ”FWIND l DDITP=DPI*2. MVSVTzO no 24 JD:191TM ASSIGN IR TO NGO T=IS(JP) J=IJ(JD) X=J*(J+l) IMP=I-? TF (NVSVT.E0.0) GO TO 19 CALCULATES VIRPATTONAL OEPFNOFNT VAPIARLES NVQVT=NV§VT-l GO TO 16 JKzJK+l MVSVT=NVSV(J<)-1 IViN3=VLI(JK) VIQU=IVUD+.S VIEU?=VI%H%VIQ1 VIHU3=VIQU*VIRU? VTRUA=VIRU9VIRU3 APIV:API-ALADI%VIRU+RAAPT*VIRU? QDIVzRQI-ALHQT*VI%U+RARPI*VIHU? VFUZWU*V[RU-XH*VIRU7+YU*VIQU3-2U*VIPU4 DRINT 339 VIKH-APIV9RPIV9VFU CON(Q)=VIPJ FON(Q)=-VI%U2 cnu(1n)=v{qua 7705') 16 17 19 COW(1])=-VIFUA COMTINUF POOPI=(J+I-1.R)**4 §OQ2X=SQPT(2.*X) COZP=-2. 1F (I.FQ.?) COZ?=4. 1KK=(-l)**JP DARO?:.5*IKK XM32X-3 QORX12=SUQT(X*(J-1)*(J+P)) X91=X91 SOHX?M4=SORT(?*x-4)+1.5-RO XQO=X*X COOOFJ=(XSO+A.*X+1.666686666666) COO(1)=CODOF)+2.#X COD(?)=COOHFJ+A.*X-4. COO(3)=COOUFJ-6.*X+4. CONTIMuE Hl(lol)=QPIV*XPI-APTV H1(?91):4PIV*SOQZX H1(391)=O. H1(29?)=RDIV*X91 ”1(392)=RPIV*SODX2MA H1(393)=QPTV*XM7+APTV H1(19?)=H1(291) H1(193)=Oo H1(293)=H1(39?) H1(lo])=%l(191)-OPI*COO(1) H1(?.2)=Hl(29?)-OPI*COO(?) H1(393)2H1(J93)-ODI*COO(3) H1(291)=H1(?9])-OPIT2*SOQEX*XP1 H1(19?.)‘-'Hl(291) H1(1.2)=H1(392)-ODIT2*SQRK2M4*(X-l.) H1(293)=H1(39?) H1(3.1)=H1(3.1)-ODTT2*SORT(X*(X-2.)) H1(193)=H1(391) 1NV=1 IF (T.FO.?) INV22 IF (T.EQ.1) INV=3 CALL FTGQN (H19F919391939H2) F1=FDI(IMV)*VFU COCfiszQflP*H1(IMV9l)**? COC1:DA&O?*9OQ2X*H](INV9?)*H1(INV91) COC?:PAHO?*(X*H](INV92)**2+2*SOPX12*H1(INV93)*H1(INV91)) DFLMH2COC0*C0+COC1*C1+COC2*C? FHC2FI+DELMu+ZEQOF GO TO N909 (199209219??o?3) CONTINUE FOMC=FO?§(Jp)-FHC TF LTNF IS OUT OF TOLFQANCF THFN WFIGHT IS QET TO ZFRO TF (AQ§(FOMC).LF.TOLF) CO TO 19 013') 10 20 .72 23 24 132 IF (1-IT(JD).F_O.O) GO TO IQ WT(JD)=O. DPTNT 149 JP9FOR§(JD)9FOMC FF9T=FE9T+FOWC IQAO=IRAO+I CONTTNOF THIS SFCTION CALCULATES PARTIAL DFQIVATIVES ONU(1P)=AG€(2*OFLNU) FOOU(JD):AOS(FOR§(JD)-FORS(JP-TKK))-DNU(JP) CON(3)=-COO(I) COM<41=COCO CON(§)=COC1 COM(A)=COC? COH(7)=1. CON(1?)=—VIPO*CON(1) CON(13)=VIOUP*CON(1) CON(14)=-VIWU*COM(2) CON(15)=VTOU2*CON(2) COM(16)=FOMC FOCL=FUC RFHM(1)=FUC AFUM(1)=FUC RDIVzRPIV+OFLRPI ASSIGN 20 TO NGO GO TO 17 CONTINUE RFUN(P)=FUC qDIV=RDIV+OFLBPT ASSTGN 21 TO MGO GO TO 17 CONTINUE RFUN(3)=FUC P°1V=RDIV-8*OFLPDI APIV=APIV+OFLAPI ASSIGN 2? TO NGO GO TO 17 CONTINUF APIV=APTV+OFLAPI AFUM(2)=FHC ASSIGN 23 TO NGO GO TO 17 PONTIAMHZ AFUN(3)=FHC ADIV=APIV-2*OELAPI CON(11=DFRIV(AFUN9OFLAPI) CON(?)=OFQIV(RFUN9OFLRPI) WRITF (1) COM CONTINUE 1F (TRAO.NF.O) pRINT 3%. IRAncFERT CALL STFPFIT (CON9FOMC9FCALC9NUSE9WT9NCONPI9ITM) C C P5 C C C C 26 27 ER 29 3O 31 3? 33 3A 35 36 DPOGQ 133 CALCULATES NFW CONSTANTG F904 STEPFIT CORRECTIONS OO 29 JD=I9NOTN ICONG=INOFX1JP) ICONGO=IOFX(ICONG) CONSTNT(ICONGO)=CON§TNT(ICONGO)+COFNIJD)*FCONS(ICONGO) GCON(ICONGO)=SIGMCO(JP) CONTINUE TTQCNT=ITPCNT+1 POINTS CONSTANTS AFTFQ FACH FIT PQINT 3G9 IHEAD9ITRCNT9ITRYS9(NAMEIM)9CONSTNT(M)9 1 SCON(M)9M=1915) IF (ITQCNToLEolTRYg) GO TO 14 GO TO 1 FORMAT (1011“) FORMAT (2AR920I199X9IP99F5.09AR) VOQMAT (OF10.2) FORMAT (RF15.7) FORMAT (QFIO.4) FORMAT ( 7H INDEX=9I39§X9 3HVU=9I39SX9 QHNV§V=9I41 FORMAT ( 7HOTHI§ I9A89 16H] CASE WILL MAKE9I29 3RH TTERATIONS FOQ MOLFCULAQ CONSTANT§.9/9 41H CONSTANTS MILL NOT RE USEO IN THE FIT IF ?QH THFIR F-LFVFL IS GQEATFR9/9 5H THAN9ER.29 43H. OATA POINTS FOP WOICH THF ORSFRFVED MINU9 11H CALCULATFD9/9 19H VALUE IS LFSS THAN9F6.29 14H WILL HF USFO.9/9 QHOCONSTANT 9RX9SHVALUE9GX9 17HFPACTIONAL CHANGF93X9 4HUSED9/2X99993X9E12.6 9HX9F7.497X9A3)) FOPMAT ( 6H VIHU=9F5.I9‘3X9 SHAPIV=9E1‘5.79‘3X9 IQHQDIV=9F15.79§X9 4HVEU=9Elso71 FOQMAT ( QHOLIN? NO.9I49 77H IS OUT OF TOLFRANCF. FORS=9 IF10.49 7H FOMC=9FIO.41 FORMAT ( 17H TOTAL FDPOQ FOR 9149 19H DEJECTED LINES IS I F11.4) FODMAT (QAR//9 3HFTT9119 3H OF9I39 51H ITFRATIONS I 4?H IQ COWULFTFD WITH THE FOLLOWING CONSTANTS9//9 2 9H CONQTANT9RX9 RHVALHE912X9 QHSTn £99099I/(IX999 3 2F?O.l?)) FNO DNJ‘JTT-‘JVH AM O3§IG PPOGQAM 03316 (INPUT9OUTPUT9TAPE11 COMMON HI(393)9H21393)9FpII319NUSEISIO)9C0NIIO)9FOMC9 IWTISIO)9COM§TNT(12)9FC0MS(12)99CON(12) COMMON /A/ NAMF(12) COMMON /1?/ INFUI9INFOP9INFO39INF049NUMCON9JUM99CONST(I?)9 1FFIN9FFOUT9TOL9IOFNT(2)9TOLF 000 134 IFOH§(GIU)9IT4YS9ITPCNT COMMON /17/ IN(G10)9IK(510)9IHFAO(919VL(1519NV§V(1519 COMMON /?1/ NI(13)9VFCTOR(1491419INOEX(I319TOFX(1319 ISIGMA(13)9COFN(13)9§IGMCO(1319NOTIN(13)9XCONST(1319NOIN JTL‘QJVH COMMON /1?T/ NK(12)9NDFLMAX9XDFVMAX OIMFWWSIOM FW1T(819 IAN7IT(81 FOUIVALFNCF (COMSTNTT 119 GSIGI9 (CONGTNT( (COMSTNT( 319 LAMOA19 (CONQTNT( (CONSTMT( 919 ZEPOF19 (CONGTNTI (CONSTNT( 719 XL19 (CONOTNT( (CONSTNT( Q19 7L19 (CONSTNTI1119HAHSIG19 PFAL LAWOA9LAMDAV IMTFGFQ VL INITALIZFS NAMFG ANO VAQIARLFS INFITII1=3ONO INFIT(?1=3HYES NAMF11125HFGSIG MAMF121=SHFO§IG NAMF(3)=6HFLAMOA NAMF(A):GHFGAMMA NA%E(G)=6HF7FQOF NAMF(A):1HCWL NAME(7)=3HEXL NAMF191=3PFYL NAMF(Q)=3HF7L MAMF1101=7HFALHSIG MAWF1111=7HFRAHSIG NA‘FII71=6HFALLAM 1J'WCON-Tl? NCOlezmqucOw+1 JK=1 KI=0 IQAO=0 FFPT:0.0 JUWD=O CONTINUE TTWCMTIO UFAO 209 IHEAO9TFOONF pranq 1H A40 DHTNTS OUT HEAOING 1F (1FOOMF.FO.HH STOD) STOP POINT P09 IHFAO 1F (TFOONF.NF.RHFMO HFAO) GO TO QFAOO ODTION CAQO 1 ?19 419 619 819 DSIGI9 GAMMAI9 WL19 YL19 (CONSTNT(1019ALRSIG19 (CONGTNT(1?19 ALLAM) QFAO 219 IOFMT9ITQY99INFOI9INFO29IHFO39INFO49(NK(I19 11:]91219MDFLMAX9XOFVMAX9TOL9TOLF9EFIN9FFOUT9ICASF A§§IGN 19 TO WIGO 07)?) 7100 135 IF (1CA$F.FO.RHMFw CASF) GO TO 2 IF (ICASF.EC.RHOFN CONQ) ASSTGM 13 TO MGO TF (ICASF.LO.RHAOO CONS) GO TO 13 C(HJTTPMIE ”FAOG IN CONGTANTS AND FRACTIONAL CHANGFS QFAO 229 COHSTWT DEAD 7P9 FCCWg qug FQACTIONAL CHANGE IF ANY OQ ALL NOT DEFINFO IF (FCOIS(11.FO.O.O1 FCONSII1=19O OO 3 11:2916 IF (FCONSIII19FO90901 FCON§(II1=FCONS(11 GO TO ”40(19 (119-591 CONTIKHFT IJCT=0 IT=IJCT K=IT CONTINUE PFAO 239 FUT PFAOQ 1w ENEQGIFQ CONTIANNE PPINT 949 FJT IF (FUTIR1.VF.ORQQ.QRRO1 GO TO 7 VL(JK1=FOT(I1 NVSV(JK1=KI DQINT 239 J‘9VL(JK19NVSV(JK1 JK:JK91 KIZO QFAO P39 FUT IF (FUT(R1.FO.QQQQ.QQQQ1 GO TO 10 K=O TJCT=O GO TO 6 CONTINUF 00 Q KH:19R KI=KT+1 CALCULATFG OOANTHM DFDFHDFNCIES 1T=IT+1 IJCT=IJCT+1 IF (TJCT.LF.T) GO TO 8 IJCT=1 K:K+] FOOS(IT)=FUT(KP) TN(IT)=IJCT 1K(IT)=K WT(IT1=1. 0531') 000 OOUO no 0 10 11 1? 13 1A 1% 136 IF (FUT1KQ),FO.O.) wT(IT)=0. CONTINUE GO TO 9 CONTTKOHT ITO=1TI3+1 quTCHFs F1 nun F2 OO 11 I=19TTO TQJR1=3*1-2 TOUR9=3*I-1 TFM92F04§(ISO91) FORS(TQJQI)=FOR§(ISHRP) FOHS(ISUR9)=TFMD CONTINUE ITM:TT OO 1? M=191? CCO\J(M)=O.O DOINTS INITIAL CONQTANTQ ANO VARIARLES CFHHTTFHJE PRINT P69 ICA§F9ITPYS9FFOUT9FFIN9TOLF9(NAMF(M19CONSTNT(M19 1FCONS(M19INFIT(NK(H1+1)9M=19121 CONTINHF JKIO QFWINO 1 NVSVT=O CALCULATFS ENFRGY VALUFQ AND PARTIAL PFPIVATIVFS WITH DFGDFCT TO THF CONSTANTO OO 19 JD=19ITM M=IN C CO . '9- I ,_4 ,-. S'\ r\ [x If I." r\ L .J 4\ 9.x I (.2 C r. C _ ft 4L (CC Q C r. : C .1. rCC ... r. f. ccc 141? OK If u? (L P. (C ft VJ fugu J: 1\ 9x (X C: Ct :94 C C MULL7 . 9 .\ r \ . «I . . iii a / x , ,:Z r t ' I ,..\ :1: Kl mu nep. mfl.m 0co.9¢r: gcc.c odoc CCCoCI mph. NH.» cnc.c¢us ((5.: (H. CFC. ,rrxmuu QflCol WK. (UF.OEJJ C.r,.‘.o (Co 'JOoxrmurJ "Fro Fwow NJJoruUfO ,rco fin. mc<.mwwo $5.». (No (V0.3UUO arc.“ CO.C Ccc.cl ARC UJP ‘,.0 CBC r: d¢:.~r¢0 FQH m~.c manuofl.flr&3 film (.9. , mam rcn CFF INC {NV CZC ccc. Nuxm. bnm. Ucfl CFC ¢:m. rwq. r:c. ccc. kc meUZ0<.JOWO AoFJCUV an? .x CCC Cru LCL (No bra if C C. C C .C F. C C(C (CC Ccc Unlk (GK 7. .2 C .LLC (CC CCC c;c ccc 6;” br\ UJTCP C QC C- 0 cc C- Q: VUNU ; cc cl : cc cl 3 CC CI c cc Cl C' rufl‘ Jpn: m¢ hero a 1: cl 1 A? U340 C CC C- C .C: Cl C Of c: {did HNO 1...; ON» ; 1., a- ' NI_. ‘ JVLU L CC CI I k? {CNS r (1 LI .\ »\.» LI r. 3.: (I .... LC CI ..L c: at .l r; A ... AVF— I .h r N .IA+I\.\ C. C Z I f..'. 3.3 CU g Cr .~II «I AI. .P .7 fl 0 I C: C¢r.~1qo FPC 4‘ 4.x r CI Fl~ C (CC r. K. I. ,. (K I: V a . . ‘ ‘1 3 f. If CCC CCC C C C CL.C (CC ccc CLL ICC Ccc C(C CCC {UM MUM v. .... ..u A. v. _u r. Cl CC C «U. 0 a». 0 U5. C CC.C C(oc cc.c cc.c n».~ O cc.c cc.c (Cot cc.c cc.c (co: cnw. uwc. hue. (pa. :Ln c::.II nu;.;I :cg.: c;.c rcc.; cc.c czc.c c¢.c gcc. cm.~ rmr. rc.r (xa.I cm. 2:2.I Ir. nwr.s :c.; at.. e.. CCCocI CCCoCI cccocl OFQ..~.:W_O .¢mma mcu.VCrd ccc.CI nUc.:.m~¢o 52“.. Km. Mcc.wwea mxpo cu.F dmfoGL¢o C:;.: CCoC Cccocl . o 0 OJ, (U » ?wC.bc¢a C;C.t .CCoc CCCoCI d¢w. cu.r Cr3.thJ c;c.c cc.c C:CoCl fi~f\*\0 CrCcC CCoC ccc.:l CC”. FCCC 11¢.nficn Occ.l ow. bxo..\,cc CCC..: CCoC :::.cl C‘Coc CCoC CCCoCI WAC CZCQ Amory r. g r, r,\ ,r\ mac.m¢¢1 4» F C. I O a\ (m p (6.» (n. r p. rC U' 1%, floc ..LCC CC.( CC (CC C .uI» Caxfiu rx \‘<\.01 N, rd. 3 C C C C Q C “42.. ff #0 003...] C- A. R C, C .r. c. C c C bf mutcrl .tc __(. C; r. C 7.: .... C C C CCU- ' r'ld' 3Cd1 . CC. _.,- I “A...\ mzcw V- \\ ....K 9.1m C .1: \| ' rf‘ Awr\4\.1. CC .C' r. 0px C- 5 CC C. (C «.....n .Qn... .» .xr. Cl. C .2: :Il ~C Nae} . Cr C- PR? 06' FIN . . Cr. /C_C.C .1. I :;C _(..\._ . . .. 1 C:: 0 L \(C RCA .c bxr {.10 ac (my: rm w! 0 .. _, Ix..~(~ or r h . “x. C o; . 3 ..2.. h. r. C . .2... C I . bk 7...._ (Ctr u .1}... Cr: (1‘ ,V ‘ C ~ .N...,n\~| p» »C.C n... ,1 v .CII a. _\ ll . . II C ‘L ,n y I f C . C‘. Al x . a <. y . I. \ \ 1.. ... |l \r r! .r | I» p w ,v 0| r \ r all n. r K . F C C: Cl FPC z, \/\ C C .C C. p. C c'\ V L: Lf C C CC c- C C: C- C CC CZ.- g‘ CC. CI c cc C' I CC Ll C CC LII! C. .CC C‘ I CC C- C CC \' r; C (c L- C CC C- I hr {Cauq .l .CU .KGGZ. a CC C. I (C GNCG RFC C C G C .\ ... (Lt CC .L C(C cc; ccc fr.... (:0 r .CC. CC: C C .x‘. Ct; CC; (CC JJJC erC Ccc CCC C (1.. ccc C CC c- C LC C- c cc Cl a CC C- L CC C- .k. CC (I ; cc C- C CC CI- I CC dc¢d 6 CC C- C (C CI g C: Cl I O{ CJU: I {w OdUu mc JFUQ I CC coma hr? 2 C C ccc (CC CCC C .L C ccc .C r. C C L C C C C C C c .Cfc (CC C;C CC< :rc Ccc Ccc CCC CT: CC. C cw: Gab (Lg CCC C (C (_| c c: C' C CG Cl C CC C- 2 CC I CC Cl. C c: Cl C LC Cl C C( C- C CC C' 0 cc C- .. ,1 Jn¢q C CC cl L C: Cl - chla- ”H.350 1 CC Cl. wnc h\~ dQU K\.(._r.¥ CCC JIM $00 3,: ca.» C C C no: C1: 3“: non C C ,1 CLC R7. r fico ,N. . C UR K LF C IV r r60 QJC owh Kl PAM flu NRC C r\.Cv rxmt uQC.Tan o C ~.O.I.O‘ IF“ mfid.¢0¢1 {C(0: .Jo. ccc.:I :r.:0: A‘:O Etc-(l »\..\’\._§ C7,. :CCOC' In, ... cc. cflc.~JUg ....CoQ 2C. cra.cl LCCQC :Co CCCOCI 1(c. «c. 3dhgmnvc (cc.c c:. cc:.cI ccc.c cc. cccocl CCLOC CC. (cc-cl -.Ccoc CC. ccc.cI 3(1ol Tho cm».rq.TI rJC.I WC. n\\«.n\0¢le4~ .. (.HC. .l.C. FCC.h(c0 an: ukkcoy .dw .C CC: C.~C ('\. r C? a. ../\ f‘ .(_ .a r‘ CT f. I.-. A n .. c c m cro CCC (C; P}. 4% .fl v FTF uqsqr R C CU P rx .3 f\ C (\ C C ... C C C _II _ ..k' P»m.~l ’I dua... .w Lv,‘ inl- r mun ~Ur~fll ,._' m.' , , IL. L I CC I a .l ,. I I _ r \ I ,r , Is an: . \ . C I. r» 0.3 r _. , C. I V C .4\ <\ A .F ”13 . t. .l .., P :17 x;» . C ..I rC ccr ‘I ~- K cc: ccc WIVLC m ‘\ Fc CCC CCC . , . LC C amcd ' r\c. (ace I A(cfl URJQ II '\ P U . .Q .C (r UPCC I mc.n dmdc I (C (ace C (Co: C- O I C» flwdc I cc.c hficd c cc.c C' C CCoC Cl C cc.c C‘ c cc.c C- mro (LC ccc UOC Own ccc CCC qfic a f I ... r r I x x | I ‘I r I ma. C A. II. . — c. c ,x I I ,_ \ 35rd .. L. l .I ‘ \r c; r .... /. (_. pm- r: l[[lllllllllal I ‘ [III-ll I ll llil II l!| I'lll a S l r V f F , k P WP; gcn CJC QCN ~. 1‘ a\ nu FJJO I Cw 03¢: I (w Crc~ Orb rqw Old Cd: OUC l m: (WC: UCPU PCFU {GNU fijva I L 6: ~ [\ (v G F C P h F '( Ck n: nunc has C CfC JFC CFC CrC CCC (¢< :3» (hr ff: CC. CC» Mac 6:: FRF C(r CPO KP: r~|r MO: hco I LC n. J. on .v, r. 3 “d C CT UFdJ {PCT .U r 4 f.».r,C CTC CC: CCC CCC cw: hnb quck C C :{J \1, I . JWUT Onwu C R U I \CUi 33d? ..C ,\ C K r. L C C lllllllllllill‘l F m C «.m‘rx. q.: 0 C "v .C C (CC (Cc OraJ 4 ‘- 4 f‘ \{ cc ruflsfirfl .... C C Cl Ffia ”I .c CC: an. CIA cu: cc: .5 P. C <¢n C x.’ x! C. O (\ cbc C..Uo\ mrw #0.? (CC (CC F3». QJN .u rrC hem hi: 1. 4\ n d C ..7 ,U3 FQC rcm .\‘ ...» p»¥ (cc CCC CCC ‘ C CC C .2 I .JR- (hh( I C(oc Gihf | JG. drxuurx L ccoc c:c cc: Kc” Nu. .p C C cc. C C C (CC ficw fife v .x r.“ rx CCN 6UQ< ~( (q.. n«. c~3.cmn< flam- (00C c:o.«.~¢ CCCoC ..kCoC cccocl CC... 0 ... C... C ccc.CI Czcoc :c.c ccc.ci c.g.c (:.c ccc.c| FRO (....dn. Amoc. (MP. p FJd .. :C r A "Q Ru... «nib. (I . hmn _\¢fi <( 1c¢.~cxr we..nvr¢ an.;0:< ctr. qu doc. UF. n< cg: no»; ch.~ch< rnr.~n~; n m.nu\¢ : {.rr;u :\_.a:;; .a;.i (» V/;\ l .4 . r o e ( L.\.'ol h'o Cfirui “Co CCCoL C(oC (pic: TCoC Ad.o PC. PX. a.. \X. .I \.o VV.» ... ,,c r, I QwJ‘oFQbru Cflfiorxfhrx CCCoCI .CCCoCI 11:;aner firmwcrkchm. -(.....{V. r.,..V.V_..;r./; ..V.,.x.......,,. .1 l O «G , O I O 0 , O O c I _ 0 V hcc.l :J , am; I << an, Uh 0J0 mp psr. rm , ::,.; c: c 5,. I V r V V .. rV: a . O V. o Vac.wonx v,m.m¢s¢ ~¢a.ac~¢ u., ¢cr¢ rcu.oakr gzc.rn ,;~.,;;a VV. pg; nV;.V :c V . O O . 1 3.”.I :Ccc Jmc.l ?r. (F: (m an; I am “2‘.I ~C. rw\.I ::.r m ..l ... )P..! VL. a, .I V ((rom‘nfirx CkuéoCUbl (UFoUFbVJ. Oblortcfiru (arm-thrx CFC-UFFJ (11.5.11: r. .1 .V.rf. ,on .r . ffpol PC. CC'ol (CoC Cero “New (.meol {Co EUR. CCoc. _,.:V_a«. CCoC ....fnl »; to if ..I ..(c r.,:o.l r (aroflcxyrx (RhonuUfiJ pCcujrxhrx Chronckflu {Jfioqrfif CCC..CI “71.04;_\. F: .-\r;.f ,;.-.4 JV Vic. Ac. ccc.c cc.c (or... 0‘... +\KV..I cc.c 43c. 3...; 5:. 1:. V1,. n:,. V: .I r... .rVi ? ;(.::r( thccl h: » .PJNIW .JVCC.3:\.< (.\.Uo~. fifii .mfi M .LfihJ rV..h;<.fJ.J¢ fir... av.“ «.1 761;. . .. 1:th CCIVL \lkaO Ml“. ”ORCO' CCOC Cn_COC CCOC nmker NJN. (Jr:I “M. Mr. 1L. . 0.. (“\Or ,VVTVVVOC (CO: r V. r £ ”1 Lgcocl Ran.nUb( nooofiUhq ccc.cl Mboorrhw »!~.nn>{ Vig.,l :::.cl 1H,..V; ¥:C.c Cc.c ccc.c cc.c ccgoc cc.c c:c.c (Coc :cc.c reoc g:c.c cc.o ,1V.: r:.V 2c;.V 6L.: :c<.p r ccc.:0 Cccocl ccc.c| ccc.:l ccccct ccc.CI ctc.\l ccrccl ccccct FPC amt WmC FRO RAG ”AC PFC x»? r,5 CIcU Am.c. bl» hC Umlu}cfi; 0 .h. ugscr IC! *\ rx. 0 r“ ’\ J I r CC I Obi ¢~ CCC. CC.; C C C (CC CC: CC( ...r. L C C C C.a‘_.\h ccc C3: \.. nxrp I CC flnOé : cc. C' C CC .Ll C CC CI C CC C' C CC C' L cc Cl C; C€1( I CC LFC¢ Arc CCC CCC. (cc CCC CCC Jaw CNVP C C 1 (CC CCC CCC CCC 7;, (V C ccc C- CC CC C- C m n \ h c I in dmhc NI. cnn< (C 0, $30 (N: or» «w; (DU IR ”A. r . '3? C C P ”Uu CL: (cc (RC 6. C L ccc an arc Ka< I} 4N fl \_ ~¢u CCC crc (CC CCC ' II. C T J K- J s cc I CC rqu (r «Ifhrx I n_ €UF< I NC RUN¢ C C: L. C cc C- oC P-I {two an. ¢1qowki¢ (my. cm. bnw.m>¢c C~Ir0 r\CO Uc~.-(¢ an». an. an.<~¢q c;;.c cc. CCCocI bx.. Uh. snx.1~¢¢ C(;o AC. hcmopdqq Vac. hr. brr.cox< puflCO Cr. n¢a.num< Che. (Co gcy.w3¢¢ CCVa CC. ¢u0.(3¢¢ ccc.c cc. C1CoCI cg:.c cr. ccc.cl CCfoC C .t (crocl baker cc. ccc.cl rm: mtzfis “mot. \ 0ft Cr< JRC NRC CFC W's‘U t CK» AFC OCC ulk .JMN .flwhw . wc onch{ . wc ohrmw . «c .CCuI o ...‘C oTLFJ . Tc .Jchq . 7c . rm. .1 k. f. o G: .mc~< .' m(a .ICN¢ .I :0 .chc 7‘ ADD hf UbleGQ 4\~.I_fl r—h F»; U.€ < ~I.NI . U7; 3 4\44I.. wi\'§k w A” : Ric Jhc T1 I r (‘I .\ VJ» fa» wwi .1: rFV ( Chc ..L 4Wf.pw 1)_.f.... (4m\ \.I..I\ fly f béV I. L .V 3: C ml. C p. '— ,_ 3;“ .w q p «NSC hwn a: ncn r,‘ C C C(‘C smc Fro Jncc daflx! (JO NP; fiVfi C» C (x CLC :ic bruc hco Q .U‘P KRJJ To Hnrw rah c 1: UV .6 r r.‘ ‘.~_. Wu: ray .; 33"“ C I: (if .(.f rs C L c. 3.... f p 3. A ‘ x -c.’ ncflxw. Tww PPN fa C, 5 HA..\ hIi 3.9; '.(...\u .r ' y F f\’ ’7“ .w N. 4.x. r‘ r\.\f i. .. n rnl JIC flkm 4... . c p f 1‘ If 57.5. CJC u_Gc» f‘ .I‘ C. C (”I r. firm gua.| on. ¢w,.| cr. gcc.c cc.c rpm. cc.c :cr. ;.. pcc. :c.¢ x1:.| r< » U . \. _,. . . . pr1~.xnxr ;;_.c cc.: c0:.| ac.c ;c;.c c:.c rg:.c cc.c ug:.v c:.c ptw. "c. c _.: .V :;r.ct undovrwh cc;.c| ccc.cl :cc.ct {zxomew c;c.cl dcg. PW. has. rc- fizcot NF. dCro ;C.C UM,. om. dlm. :C.n :r:.\ .: .7r.maFn ban.¢<,h (ch.hcnn <7,.c~c~ m ,em.cccn nq,.rccp NJCoI (Lo FLO-l V0. JTCoI (mu-w RIG... 0.». (QFo {Co hair 0 PM. __ V... TC V V V V .4 V . .V V , . V o o (Lronvck. créocrrh ..Jofifirk. FPhochh Oh'oPrMCP .ru. (Lab (T: :l V16... rvom «....vmol Um. .rcgfo croc and”. CH. chol Uh. {Scot 1:3“. 3:7, . u~:.curb ctJ.Ncrh ccm.¢l qpuoxccn m:£.w ”1&4er. C7. rmno6.~ :(,b CYCo: CCoC (Cicc CCcC :cgaC (Co: on... JR. (in. c“. u» ..I. hf. Tr ,.L.C.CI cccocl ccc.cl ...:C-Ccm; .r.,...?:c\r_ .,\.,.r..r...:k 1,: .,.,..:. cc. 5:1,. ”dc. ..onl h”. 9:. Cr. 14?. 1.x... 73;.I «(om : ._, O . I C I . . SCI \IJih. RFJ. (Crk. .Unur. (\wa 4\\I;\0CZCN Rruror.o.\r.~. UUaxoPJQL. C. C Fro «#0 pro «an was 'mfi , »:CCC .Hm u.Cgr f 1Y0 Q C C (ac (\ ... U r- r» r C x- k. Vb: 11¢ ox ’V C Q CC h K. (C J. C . QCJ hhh 4 . Vb (‘L_\. (a: r‘ L' k d 4 r (ON. C r\ .L. (V C (’&.. C Q C cc: CC; .I CC Ctr. :cc c:c ccc (CC Ccc C r. L c:c ”(H Fwd _f. P. C ccc CW.» 6N0 (r (cc :CC CC; Cc; (CL 0 C CC C- c (C C- C CC C- c cc C- 9 cc CI I J: ruth l WC Foch I. CC fiaon I «H c~¢h c: deb c C? C- I CC 02qu 2 CC C- c cc Cl : c: Q- wra lac. cc.c Can.r~qn KMC. Wk.-. 30(0NC4VF ”FR. FC. .axrflofi‘cdh C..,Co ...Co OfinonTLIF md;.I cccr n1”.momn O 4. rxc (r d¢c.>:mn Criom cco: cccocl nwfio 6:. QJrocxrr THC. CCOC [caobhfih Cac. cw. m.(.».\orI.\IF.\. $.po (:- UJ.P.~FNI..CI hxw. rr. UJvowwmh :c;.c cc.c crc.:l C, (0: (CoC Ccc.cl c. Co: CC...<1 c<:.cI was Clflf ’30Rv IJr. Pdc C. CR Uer (C r CC... C CC C _..._. C .r. C r. Llr .A It I CH nqqk C .r._ ~wd> on ¢¢c~ C: mdqm anN C? ”raw EC fnm~ I mp (OWN l :C ERG f. qr uJ 0 5.. a: rc 3..” r. CC. C. UL1C3<1C C 3 J C 0 ' ... ,C C Awfh . 0C C L s. C Cl C II .... L Nada. . | VIfi. ruh4uh V C. (”fix . . ~ __ ’ I.fl\‘\ . ....I‘ .(« ;.l I. KC. swap Ceca r \ch. ' r‘. p O (pcx C;- L {CJN . C | C r. C Ccowh O J Cr. A. r»\r\. . . ... A» (xrh . ..C C: C cl ,A0 .A\ uwdu (. \ r\ (. .,.. <\ *. I?! .... O I. CC 0 6 Ci. .0 ..l. CC 0 : ...... ._ O .z\. C,\ O C L i. I O .... (:0 Fr . I z... 0 7.. :C.‘ C ..v...‘ 0 C. _ o, o, . .0 I 0 1C: F ' C C C C II C C ..V C - C C C C | : C C I r... f C r. ' (ccuo CLoC ((Co: CCo: CLCoC CCoC (......C (cor. C;Cor Cro; (1...; CCCoCl LCCcCI CCCoCI CCCoCI CCCuCI (CC-Cl :c_,._.c Cc..c (Cgor CCoC ccc.c (Co: (Ceca Cc.c :CC.: c;.c :.;.-..; (crocl cccocl ccc.cI ccc.cl cccocl canal o 0 .0 . 0.. 0 . . o .... (I... c Cc ; cc: r. (Coc r...c.c c:.c .....c . cc : "7,... cc. c (r... : _ 0.. 0 O . O o 1 .cc c (L: (I cc: :I (Cr. (I (Cc (I cc: cl ....u (Lac {CL}; ccoc cacol. ccoc (an. (C... fir... :Coc :...:.: (......Cl .LCLoCI CCQofl um‘CoflRUF CCCoCl C...C.Cl CC... CQoC (fig-C CC.C Cicoc LCoC (.5. RC. ;.Cog 19.6 (that choCI CCC..C| C(CcCl [pans—pmvb C;C.Cl C..._o;I cgr.c (Coc :;c.c cc.c c:c.c cc.c Pq;.I :M.q gcc.c :c.c c::.c _._:,:.cI :_::.CI gczocI hchobn. cgc.;n (Cpol pk:- CLC-C (Cor... rgCo: C100 rukfol FL. H1.11.3 4.:- r, .r..,.. www.mCUh ccc.cI c I CC 1h<~ CUQN C (C C- C- c C: Cl I AH flack I PC nu<~ C CC C- ; or C Cr 1‘. ”re Iv:.I cg.o UJUon¢fi> ,UQPO CUOIP 0¢O.m¢m~ OJ“. CCoC £O£orwrr (aw. FL. OthoflJnr fl¢Co WC. aprocqrn (rwr. ....C0rs ,r,.er~ {tiol RH. ”Lu‘yofififirn flu». n». rq¢.cur~ (Izoc cc.c TCCOr' v.s.UrO NICO c¢u.¢urn CCCoC CCoC ccc.cn 0mg... (CCC. hmr.n(rb Y}. O . C C~IC PC ¢;:oh(rh CrVCoC :(o.( Crc.rl Crroc ccoc CCCoCI NRC FIX/NS Adam.» u i r. 3» JR. C Trfih b 1n. I CCJN .C . r CC C J CI . sh. r QUFN : r‘ ... a L NI ..U NI \ 5 AI; r O .1» thr I x o r. Torr h a O 9 [I n ‘ ndflh u . 0 I CC C ( n. \r P n F 1. O f. r r {hrfi N (:0 p Phflh C -o \ RF n Abra r I \flc C rKI 9..» C C .I .r r. A (I. C \.\.. (r, ~ N R J _ . LI‘ L ’7 R e C: I. #w. . ~ .~ C C ,C o. .I C nIrWC. CAF . (... Cm. C o C L r gr: 0 .. f ., (. r\{, c r -. r.rfl ax . t r. a . .rxm o : J , 4\. r; H: V. Rah . \ C Thf Cwfi O C a r ..\ .\...,I, o T I - ill: 0 x P «...:~ 0 fl .. .I Cr _, 0 ‘ .r . r,r,(. b.~rr F rwy (Y , \ I'll-[I‘llIIIIIIIlII‘ ,. .I., III CC C ZCC .2 CC : (CC h. (C r\ II\» (.C _.( .7: .7 .7; : f, .. .. :1 Q .x, . II LI . I , I (”I C r ( ,.. .I, . I. O .0 - . I. O ' OI _ O I .I I Ii 0 17:. CI QCL L (.C( C (n..- “Ibex: ((6. (I Cfc a. .. T-I . :, ,I..I (fI. ..I LQCIIKI. (CDC ,..L:OCI. CCOC CCCOC CCOC JNIAKO CLON f..(I( CCOC ....u....l.,. CCIC ... .I: .50.. I., .011 ...: I_.A. .. a. l. I. O . I O I. I I (gr CI (CC (I (CG CI CJC 5!.1. Crc.CI (Cc :I (cc :I .;.. (<4? c,..:..CI I . ..o o .I o. I o .. o . . I _ , cacac. cc c (5: c cc c c,Ic.c cc c 4:: I cJ.d :_c<.c stor (kc.: cc.: c: I... ;:.c \II .I :C._: :3. :I...,, O .0 I O. O . I . I t (CC CI (CC CI. CLC C- Phat fkm‘dkn (CCOC' (1: Cal C. N. 0:. rvfo. r \\I g \( II a O o 0 CI 0 _ I O . o O _ I0, I o .0 . o .c c o I ...,I L cc c c: . :c L (cc c c: c (cc r. cc C ..I, c c. t: c .-I. _: CC C (III I : . .I I , , .I I g ,I . .. I, 0, II _ . I III I. 0. II .. O. .I. .01, f. O. C .I : _ I. It: 2 (CC .2 {CC 1 (CC C C(L cl (fr. Cl Id.» .I._n\ H.Ir. g I . I c O O ‘ . .. I . .. I I I . . l Tn: I Kr L;:.c CCoC :ccac (Soc Jr». CCoC (.croC (Coc (1.....C Cf... rtkcl c...,.; .7. ... ..:.; I, I .. I o _ o I o. .o . I I0. I I . . o I . . ....rIp nIAJb (CC CI CCC CI FIG :PQN. (CC (I C_I:o:I .er. .Isfivf 7 Soil. ... r K .0 o In. 0 o ,,O O o O O o o . o . o o a .. (“Jr I CC C (C C C: C (r; L c: c hf I WC Tc; T: c : C ; ... _ T 1., I .I _ . h I ( : x, 1 .I O . \ O O C ... I , . cnr 4,ar ccc cl (cc cI q§r cn¢n :wv.m.rk cc:.c| wcr.ntnn :... I ,1I.:;:: NE... v.0. (22.; cc... 9:37: cc.c ....JH. Irc. ,I,.__.I.< (Co: :.._-_._ _;Ic.: : ...I .(I .I p II. J.,.. , . ..nx...I¢., ,anxI CCCICI CCCOCI («MICKOQVMI CCICICI CCILOCI TFIZILI I.I\. 0.....a., .LCZ.‘ ,7.I..\ . . .I acoc rc:., c:.: cccor cc.; :....c “in: JITI.I ...Iroc ,I:_c.: :i.: ....I.; _;I,.,. IIItI. “2’2. .1 I. Ii}, 0. . O . O I. I .. . ( ”PfioFCCF (Cr. . I CCC CI (CC CI IdioCC£P CQC CI fr .1. LI r TI c.u..._: ,. ,1. . a . .. O O I I O O U 0 O O O l. I O I O 1 nix. I (C (CC ; CC C c.C C Cr c Frxc FL .(I fir I: C .,( \ 2.. a (, . a . I c : . I, I I .. o o a I o ... I. o v . . , .. I. (In I:«I.I (CC LI (CL CI PJC Nfdb :Cwortdh C: ..I .II I h . o‘,...,. I \ 4 ,.\ Y I. J,....,I,.I cc..c. :c...c Lao: 6.1.6.: (6.6 nit. hI:.C .....(o Ir. c. csn.np.a J «.quu pr.I (XI. III muio .\...I\..OII .LCOC QCW‘I‘C (CDC ..‘Ccr CCOC dur.\' .1an Avrflo :70: «I. . QII 750$. . . v; r..( z o .a.<,I.1.rc\ .....33713. cc:.cI ccg.cI 0{<.« a?» 231. :4. III. a}... u ,:.. L.:I . ,.II , .I \ o ... .I .0 . O I. .\0 0 .0 O 0 .10 o .0 II. I o o a 9 (It? hr hr. I {c cc: c c: c ”.1; «r rcc N" 3n: r: w: _,: I :A :. .. .I _ , o . . 0 I o o I. .0. I a o o I I. kc: QP¢N ucf fia¢h (cc CI Ucc nrrh O§C .ndr bfu rst J1; c,In :JI nv.\ rfg ::Ib . Ix, ¢ - . i n\ I ~- . I \ PFC RFC rPC PAC On? TYC Pfir r: ”7‘ (V .\ ICQ C. C be C C (C 4 an“?! -. III : pr wxfi.t cc.c Chg. cc.c :;:.( cc.c arc. mw. a” nafir amm.~md> mm<.qun ccc.cI mac.rrmw v. (a mnfl.t «a. Uxm.l «w. Jcc. LC. qua. mc. ;m garb roo.~aqn cm5.¢¢c> «(n.xurn pw¢.~¢nn or I {c Cgc.r cc.c (29.: rg.c CU”. Cflo (:c.I :c.c 14 affix: CCLoCI (ICC-CI (GnoCRfiF (”Co‘VJFN. Qr I ¢q on..l C, ,cm.n ow. wyc. mc. Crr. wc. fix bhrr 3¢1.nw m3>orr