A corrected CIPW program for interactive use

Computers & Geosctences VoI l I No 6 pp 787-797 1985

0098-3004/85 $3 00 + 00 © 1985 Pergamon Press Ltd

Pnnted in the U S A

SHORT NOTE A CORRECTED

CIPW

PROGRAM

FOR

INTERACTIVE

USE

DANIEL FEARS Department of Geology, Eastern Washington Umvermy, Cheney, WA 99004, U S A

(Recetved6 May 1985, accepted26 July 1985) To make an interactive C I P W normaUve calculation program, I selected to modify the C I P W subroutine of the H A R D R O C K package (Till, 1977) To test the program, a variety of analyses, including Tdl's test data, were used Although the program produced the values given by Till for his test data, it proved to be m error for leuclte and kahophdlte normative rocks For example, Table 1 column l shows the analysis o f a leuclte ljollte (Washington, 1917, p 701) C o l u m n 2 shows the C I P W n o r m for the data when computed using the Till program, and column 3 when computed using my corrected versmn of the program My revision is based upon the algorithm of Washlngton (1917), with the clarifications of Kelsey (1965) Changes from Tlll's original program include rewriting the silica undersaturaUon calculations and changes m the way the chnopyroxenes are presented A test condltlOn to prevent dividing by zero when the normative plag~oclase is calculated also has been added (Wheeler, 1978)

As seen in Table 1, the Till program produces erroneous values for leuclte and kahophlhte This is due to the method used in calculating the silica undersaturated minerals My revised program eliminates slhca undersaturatlon In a senes of steps At each step some of the deficiency is removed by converting a speofic silicate mineral The chnopyroxene calculation has been changed from wo-en-fs endmembers to one of dl-hd endmembers Akermamte (2CaO-MgO-2SIO2) similarly has been replaced by calcmm orthoslhcate (2CaO-SIO2) Values in the arrays X F O R M and X N N A M have been changed to accommodate this The molecular weight of zircon [XFORM(35)] also has been changed It is hsted incorrectly m Ttll's program as 123 22, and has been corrected to 183 31 in my revision My revision IS gwen in the Appendix as an rateractive program It can be modified back easdy into a subroutine for insertion into the H A R D R O C K package

Tablel Chemmalanalyslsand CIPWnormsofleuclteOohte COLUMN $I02 TiO2 A1203 Fe203 FeO MgO CaO Na20 K20 P205 CO2 ~O3 H20(+) TOTAL

I

39,59 2,45 17,33 2.53 8,19 3,15 9,79 6.70 7~63 1,10 0,68 0.07 0,08 0,43 99,72

Crbt~MN 2

COLUNN 3

NEPHELINE LEUCITE KALIOPHIL/TE ACNIT? AKERMANITE FnRSTERITF FAYAL/T? ILMENITF NA-MSTLICATE APATITE CALCITE THENARDITF PYRITE TOTAL WATER SUM MTN$+H20

25.29 182.64 ****** 7,32 19.17 5,50 8.36 4,65 0.29 ~.61 1,55 0,12 0,15 0,43 150.34

NEPHELINE LEuCITE KALIuPHILITE ACMITE CA-ORSILCATE FORSTER/TE FAYALITE ILMENITE NA-MSILICATE APATITE CALCITE THENARDITE PYRITE TOTAL WATER SUM NINS+H20

25:,28

DIF?ER,INDEX COLOR INDEX £N % OP OPX F~ % OF OPX FO % OF OLIV FA % OF OLIV Wn % OF DIUP E~ % OF D/UP FS % OF D/OP NORNAT, PLAG

101,t8 25,83 0.00 0,00 39.66 60,34 0.00 0.00 0.00

DIFFER,INDEX COLOR INDEX EN % OF OPX FS % OF OPX FO % OF OLIV FA % OF OLIV DI % OF CPX HD % OF CPX NORNAT. PLAG

~.29 5,83 0,00 0.00 9,~7 0,33

REFERENCES

Kelsey, C H , 1965, Calculation of the C I P 'd/ norm Mmeralogacal Magazine, v 34, p 276-282 Tdl, R , 1977, The HARDROCK package, a senes of FORTRAN IV computer programs for performing and plotting petrochemical caleulaUons Computers & Geoselences, v 3, no 2, p 185-243

0.00

21:I t1:l|

50 36 65 29 51 ~5

i!

76

~.00 0,00 0.00

Washington, H S, 1917, Chemical analysis of Igneous rocks, pubhshed from 1884 to 1913, reclusive, with a critical discussion of the character and use of analyses U S Geol Survey Prof Paper 99, 1201 p Wheeler, J F 1978, Comments on HARDROCK Computers & Geosc]ences v 4, no 4, p 372

787

788

Sho~ Note APPENDIX C THIS PROGRAM

IS A M O D I F I C A T I O N

OF THE CIPW S U B R O U T I N E FROM

C THE HARDROCK PACKAGE OF ROGER TILL:

TILL,

R,,

1977, THE

C HARDROCK PACKAGE, A SERIES OF FORTRAN I V COMPUTER PROGRAMS C FOR P E R F O R M I N G

AND PLOTTING P E T R O C H E M I C A L

C COMPUTERS AND GEOSCIENCES, V.

3,

P.

CALCULATIONSl

185-243,

C C MODIFIED BY DANIEL FFARS C C ORDER OF STORAGE FOR THE ARRAYS OXNAM AND XFORM C

CELL

VALUE

CELL

VALUE

C

1

$102

I

QUARTZ

C

2

TIO2

2

ORTHOCLASE

C

3

Zr02

3

ALBITE

C

4

A1203

4

ANORTHITE

C

5

Fe203

5

NFPNELINE

C

6

Cr203

6

LEUCITE

C

7

V203

7

KALIOPHIL~TE

C

8

FeD

R

CORUNDUM

C

9

MnO

9

ACMITE

C

10

MID

10

WOLLASTONITE

C

11

CoO

11

DIOPSZDE

C

12

MgO

12

HEDENBERGITE

C

13

CaO

13

.....

C

14

SrO

14

ENSTATITE

C

15

BaD

15

FERROSILITE

C

16

LI20

16

Ca O R T H O S I L I C A T E

C

17

Na20

17

FORSTERITE

C

18

K20

18

FAYALITE

C

19

Rb20

19

MAGNETITE

C

20

Cs20

20

HEMATITE

C

21

P205

21

ILMENITE

C

22

C02

22

RUTILE

C

23

SO3

23

TITANITE

C

24

S

24

PEROVSKITE

C

25

F

25

Na M E T A S I L I C A T E

C

26

CI

26

K METASlLICATE

C

27

H20(÷)

27

APATITE

C

28

H20(-)

28

CALCITE

C

29

TOTAL

29

CHROMITE

C

30

FLUORITE

C

31

HALITE

C

32

THENARDITF

C

33

Na CARBONATE

C

34

PYRITE

C

35

ZIRCON

C

36

Ha FLUORIDE

C

37

TOTAL WATER

C

38

SUM MINERALS + WATER

C

39

DIFFERENTIATION

C

40

COLOR INDEX

C

41

E N S ? A T I T E % OF OPX

C

42

FERRO$ILITE

INDFX

% OF OPX

Short Note C

789

43

FORSTERITE

% OF OLIVINE

C

44

FAYALITE

% OF OLIVINE

C C

45

DIOPSIDE

% OF CPX

46

HEDENSERGITE

C

47

N O R M A T I V EPLAGIOCLA~E

C

48

.....

C

49

.....

C

50

.....

% OF CPX

C REAL wTMOL(28),PNOL(28) REAL OXIDE(2g),xFORM(37) REAL XNORM(50)

N,IO/2/

INTEGER

CHARACTER*6

XNNAM(2,50),OXNAM(29)

CHARACTER~40 SAMPL MOLECULAR

WEIGHTS

OF OXIDES

DATA W T N O L / 6 0 , O 9 , 7 9 , 9 0 , 1 2 3 , 2 2 , 1 0 1 , 9 6 , 1 5 9 • 6 9 , I 5 1 , 9 9 e 1 4 9 • 8 5 , 7 1 , 8 5 , I?oeg4,74~70,74,92,40,30,56,08,103,62,153•34,29,$SeSL,90,94,20, 2186,93,281,80,141,g5,44,01eSO,06,32,06,19•O0,35,45,18.01,18•O$/ MOLECULAR

WEIGHTS

OF ~INEPAL

SPECIES

DATA XFORM/60,og,556,70,524,48,278,22,284,12,43b.52,316,34, 1101,96,462•03,116,]7,216,56,248,11,0,0,I00,39,131,94,172,25, 2140,69,203,79,231,54,159,6g,151e75,79,go,Igs,0T,135,98,122,07, 3154,29,336,22,100,09,223•84,78,08,58,44,142,04,105•99,I19e97,

4183,31,41,99,18.01/ NAMES OF MINERAL SPECIES

DATA XNNAM/•OUARTZ't "

I'ANORTH*e•IT[

2*HILITE','CORUND*,*UM 3"DIOPSI*,WDE 4•lIE

5"FAYALI*,•TE

,'ACMITE

•,'CHROMIe,eTE

1"

*,*HEMATIe~eTE *TE

* , ,?LUORI•

•TE

• INDEX,

•COLOR

• ,•INDEX

",•HALITE',•

• •e

",•ZIRCON',

",'SUM M I ' e e N S ÷ H 2 0 * ,

" , • E N % OeeeF OPX "perS % 0",

3"F OPX •,'?0 % O*,•F OLIV e ,*FA % O*,•F OLIVee'DI

PLAG*, •

% Oe,'F CPX •

","

'e"

",

OF OXIDE SPECIES

DATA OXNAM/'SI02

1"V203

*KITE

",'CALCIT*,

* , * N A - C A R e ,•BONATE ° *PYRITE',"

4•HD % O*eeF CPX *,*NOR~AT" , * , 5e •, e e e/ NAMES

*•'~LMENI',

*e•PEROVS•

* *APATIT*,•E

• , ' N A FLU'e*ORIDE " ,'TOTAL *,•WATER

2,DIFFER,

" , ' K A L I O : e,

",*WOLLASe,'IONITE',

"," *,'ENSTAT'. elLCATE*e'FORSTE*,*RITE ,

• ,*TITANI*

7*NAwNSIe,eLICATE•e*K-MSIL • ,*ICATE 9"THENAR',eDITE

"

• ,*CA-ORS*

",°MAGNET*w•ITE

",'RUTILE*, •

8"E

",•ALBIT[',"

°,'LEUCZT•,*E

*,•HEDFNB" ,*ERGITE • *

e,*FERROS°,'ILITE

6"TE

*,'ORTHOC*,°LASE

°weNFPHEL* ,*INE

e,eFeO

"ceil02 e,eMBO

e,eZr02

e,•A1203

•,eNIO

",'COO

*,'re20]

",'Cr20$

*,*MgO

*eeCaO

* •,

2"SFO

*,*~IO

•,*LI20

•,*Na20

• eeK20

','Rb20

"ewC$20

•,

3"P205

",*CO2

",'$O3

",*S

• ,*F

",'C1

"e'H20(÷)

• ,

4"H20(-~e,'TOTAL PRINT*,"

•/

•

P R I N T ~ , * T H I 6 PROGRAM COMPUTES A CIPW NORMATIVE CALCULATION ' PRINT*,'WHEN GIVEN A CHEMICAL ANALYSIS, DATA MUST BE INPUT e PRINTee'AS REAL NUMBERS ( £ e e , THEY MUST CONTAIN A DECIMAL • PRINT*e'POINT) IN TH~ FORM OF WEIGHT PERCENT OXIDES, e PRINTVw'ACCURACY SHOULD BE TO TWO DECIMAL PLACES,"

790

Sho~ N~e PRINT*,eOUTPUT IS IN THE FORM OF A DATA r i L E | PRINT*,' ' PRINTSpeENTER NUMBER Or SAMPLES" READ*eN DO I J#I,N

PRINTSe~ENTER SAMPLE NAME (40 SPACES MAXIMUM) w ~EAD ( * , 5 0 0 ) 500

SAMPL

FORMAT (A40) PRINTS,'ENTER $102" READ*,OXIDE(1) PRINTS,'ENTER T102' READS,OXID~(~) PRINT*,'ENTER ZrO2" READS,OXIDE(3) PRINT*,eEhTER AI20]" READSeOXIDE(4) PRINTV,'ENTER Fe203* READSeOXIDE(5) PRINTS,'ENTER Cr203" READS,OXIDE(6) PRINTSe'ENTER V203" READS,OXIDE(7) PRINT*w'ENTER FeO" READ*,OXIDE(8)

PRINT*e'ENTER Mno" READSeOXIDE(9)

PRINTSweENTER NIO w READS,OXIDE(10) PRINTSe'ENTER COO" READ*eOXIDE(11) PRINT*e'ENTER MgO* R£AD*eOXIDE(12) PRINT*,'ENTER CaO" READ**OXIDE(13) PRINT*,eENTER SrO" READ*,0XIDE(14) PRINTS,'ENT~R BaO ° R/AD*,OXIDE(|5) PRINT*,WENT£R L120" READ* OXIDE(16) PRINTS,'ENTER Na20 ° READ* OXIDE(|7)

PRINT*e'ENTER K20' READ*,OXIDE(18) PRINT*weENTER Rb20* READ*,OXIDE(19) PRINT*w'ENTER CS20* READ* OXIDE(20) PRINTS,*ENTER P205" READS OXIDE(21) PRINT**'ENTER CO?" READ* OXIDE(22) PRINTS,'ENTER S03"

FOROO2,DA?/|"

Short Note

READS,OXIDE(23) PRINTIeeENTER S' READS,OXIDE(24) PRINTS,CENTER F" READ*,OXIDE(25) PRINTS,CENTER CI" READ*,OXIDE(26) PRINT~,'ENTER H20(+)" READ~,OXIDE(27) RRINT~,'ENTER H20(-)" READ~,OXIDE(28) DO 501 I=1,28

PMOL(I)=OXIDE(1)/WTMOL{Z) 501

CONTINUE IF(PMOL(S].EO.0.0) GOTO 504 IF(PMOL(B).EO.0.0) GOTO 504 SI-PNOL(1) TI-PMOL(2) ZR=pMOL(3) AL=PMOL(4) F3.PMOL(5)+PMOL(7) CRwPNOL(6) F2aPNOL(8)+PMOL(9)÷PMOL(10)+PMOL(II) FMGfPNOL(I~)÷PMOL(16)

CAmPMOLCI3)+PMOL(14)÷PMOL(15) ZNAmPMOL(17) yK.PNOL(IB)ePMOL(19)+PMOL(20) P=PMOL(21) C-PNOL(22) 83=PMOL(23) 3aPMOL(24) FsPNOL(25)

CLsPMOL(26) DO 502 I = 1 , 5 0 502

50]

XNORM(I)=0.0 XNORM(37)mPMOL(27)+PMOL(28) WRITE (I0,503} J,$AMPL FORNAT(IHO,SX,'CIPW NORM FOR SAMPLE NO. ° f I 2 e e l e e 2 X e A 4 0 / I H ) GOTO 508

504

DO 505 1:I,50

XNORN(I)=O,O

506

CONTINUE WRITE(IO,506) J,SAMPL FORMAT(IHOe'NO CIPW NORM FOR SAMPLE NO, ' , I 2 , ' t e , 2 X e A 4 0 )

507

WRITE(IO,507) FORMAT(IHO,'BECAUSE Fe 2+ AND Fe 3+ WERE NOT OUOTED SEPARATELY')

505

508 C

GOTO 1 CONTINUE APATITE, FLUORITE, NA-FLUORIDE CALCULATION XNORMC27)=P

2

CAmCA-3.33*P IF(F-0.66#P) 5,2,2 F=F-O,66~F

791

792

Short Note

3

XFCF-2.08CA) 3 , 3 , 4 XNORH(30)=F*0.5 CAuCA-F*0.5 GOTO 5

4

XNORM(30)aCA F:F-2,0~CA CA:0.0 XNORM(36)=F ZNA:ZNA-F*0.5

5

C 6

TF(CL,EQ,O,O) GOTO 8 HALITE CALCULATION IF(CL~0.5*ZNA) 6,6,7 XNORM(31):CL ZNA=ZNA-CL*0.5 CL:O.0 GOTO 8

7

XNORM(31)=ZNA CL:CL-ZNA

ZNA:O,O 8 C 9

IFCS3.EG.0e0.OR.ZNA.EG.0.0) THENARDITE CALCULATION

GOTO 11

IF(S3-ZNA) 9,9,10 XNORM(32)mS3 ZNA=ZNA-S3 $3:0,0 GOTO 11

10

XNORM(32)=ZNA S3:S3-ZNA ZNA=0.O

11 C

I F [ S , E O , 0 . 0 ) GOTO 14 PYRITE CALCULATION

12

IF(~-2.0~r2) 12,12,13 XNOR~(34):O,5~S r2:r2-o.5~s $:0,0

13

GOTO 14 XNORM(34)uF2 Sm~q2,0~F2 F2:0.0

14

C 15

IF(C.EO.0.0) GOTO 19 CALCITE, CANCRINITE CALCULATION IFCC-CA) 15,15,16 XNORM(28)=C CA:CA-C C:O,O

16

GOTO 19 XNORMC28)=CA C:C'CA

CA:O.0

17

IF(ZNA.EO.0.0) GOTO 19 IF(C-ZNA) 17,17,18 XNORM(33)-C ZNA=ZNA-C

Short Note Cg0.0 GOTO 19 18

XNORM(33~sZNA CsC-ZNA

ZNA=0,0 C

CHROHITECALCULATION 19

Zr(cR-r2)

20,20,21

20

F2~F2"CR XNORH(29)zCR GOTO 22

21

XNORM(2g)=F2 F2=0,0 ILMENITE CALCULATION

C 22

IF(F2-TI) 23,23,24

23

XNORM(21)=F2

TI=TI-?2 P2=O,O GOTO 25 24

XNORH(21)=TI

F2=F2-TI 25 C

TIm010 IF(ZR.EO.0.0] GOYO 26 Z I R C O NCALCULATION XNORN(35)aZR

SIISI-ZR ZR=0.0 PROVISIONAL ORTHOCLASE CALCULATION

C 26

IF(ALoYK) 28,28,27

27

PORmYK

AL:AL-YK YK=0.0 COTO 29 28

POR=AL

YK=¥K-AL AL=0.0 XNORH(26)=YK ¥K=0.0 PROVISIONALALSIYE CALCULATION

C 2g

IF(AL-ZNA)

30

PAB=ZNA

31,31,3n

AL=AL-ZNA ZNA=O.O

GOTO 32 31

PAB=AL ZNA=ZNA-AL

C 32 33

AL=0.0 ACMITE CALCULATInN IF(ZNA-F3) 34,34,33 XNOR~(g)=F3 ZNAzZNA-F3 F3=O.0 XNORM(25)uZNA ZNAm0.0

793

Short Note

794

GOTO 35 34

XNORM(9)sZNA F3mF3-ZNA ZNA=0.0

C

ANORTHITErCORUNDUM CALCULATION

35

IrCAL-CA)

36

XNORM(4)uCA

37,37e36

ALuAL-CA CA:O.0 XNORN(B):AL Ab:0t0 GOTO 38 37

XNORMC4)uAL CAmCA-AL ALs0.O

C

PROVISIONALTITANITEe RUTILE CALCULATION

38

IFCCAoTI) 40,40p39

39

XTN=TI CA:CA-TI TZmO.O GOTO 41

40

XTNzCA TI=TI-CA

CA:0.O XNORN(22)-TI TI:0,0 C

MAGNETITE, H[MATYTE CALCULATION

41

IF(F3-F2)

42

XNORMC19)mF2

43,43e42

F3:F3-F2 F2:0o0 XNDRM(20)=F3 ~3:0.0 GOTO 44 43

XNDRM(19)=g3 r2:r2-F3 V3:o.o

C 44

PROVISIONAL PYROXENES CALCULATION XPYRO:FMG+?2

IF(XPYRO.EG.0.0) XRAM:0.0 IF(XPYRO.EOeO.O) GOTO 45 XRAM:FMGI(FNG+V2) 45

XRAI:lt0"XRAM

46

XDIsXPYRO~XRAM

IF(CA-XPYRO) 47e47,46 XHDsXPYRO~XRA1 CA=CA-CXDI+XHD) XNO=CA CA:0,0 GOTO 48 47

XPYRO=CA XDIuXPYRO~XRAM XHD:XPYRO~XRAI

795

Short Note

CAIO,O FHG=FHG-XDI P2=F2-XHD PEN=FNG C 48

PFS=F2 THIS IS THE HOMENT OF TRUTH SA=b,O$(PAB÷POR)÷2,0$XNORH(4)+XTN÷XWO+2,0$(XDI~XHD) l~PEN+PFS+4,0*XNORH(g)+XNORM(25)+XNORM(26) IF

C

( S I ' S A ) 50~49,49

SILICA DISTRIBUTION

49

OVERSATURATEO

XNORM(I):SI-SA XNORH(2):POR XNORM(3)~PAB XNORM(IO)=XWO XNORH(ll)=XDI XNORM(12)=XHD XNORM(14~=PEN XNORH(Ib)=P?~ XNORM(23)=XTH GOTO 600

50

SD:SA-SI SB:O.5*(PEN+PFS) SILICA DISTRIBUTION OLIVINE

C

IF(SD-SB)

51

51,52,52

XNORH(2)sPOR XNORM(3)=PAB XNORM(IO)=XWO XNORH(ll)-XDI XNORH(12)=XHD OL:SD XNORH(17):OL*XRAM XNORH(18)=OLSXRA1 PH¥=P£N+PFS-2,0*OL XNORN(14)=PHYSXRAH XNORM(15)=PHY~XRA1 GOTO 600

52

OL:SB XNORM(2)mPOR XNORM(3)wPAB XNORM(10)sXWO XNORM(II)mXDI XNORM(12)=XHD XNORM(17)mOL*XRAM XNORH(1B)=OL*XRA1 SD=SDoOL IF(XTN,£O,O,O)

GOTO 55

SB=XTN C

SILICA DISTRIBUTION IF(SD-SB)

53

54

53,54,54

XNORM(23)=XTN-SD XNORM(24):SD COTO 600 XNORY(23)=O.O

PEROV~KITE

ROCKS

796

Short Note XNORM(24)=XTN SD=SD-SB 55 C

SBm4.0*PAB SILICA DISTRIBUTION IF(SD-SB)

56

NEPHEbINE

56,57,~7

XNORM(5)=O.25SSD XNORM(3)=PAB-XNORM(5) COTO 600

57

XNORM(S)=PAB XNORM(3)=O.O SDmSD-SB

SB=2.0VPOR C

SILICA DISTRIBUTION I?(SP-SB) 58

LEUCITE

58,59,59

XNORM(6)=O.5S$~

XNORM(2)=POR'XNORM(6) GOTO 600

59

XLC=POR XNORM(2)=O.O SD=SD-SB IFCXWO.EO.O.O)

GOTO 62

SBwO.5*XWO

C

CA ORTHOSILICATE

SILICA DISTRIBUTION IF (SD-SB) 60,61,61

60

XNORM(16)=SD XNORM(IO)wXWO-2,0~SD GOTO 600

61

62 C

XNORMCJ6)~SB XNORM(IO)mO.O SD-$D-SS SBuXDI÷XHD SILICA DISTRIBUTION IF($DmSS)

63

KALIOPHILITE

63,64,64

XNORMCI6}=XNORM(16)+OeS*SD XOL=O,5~SD XNDRM(17}=XNORM(t7)+XOLSXRAM XNORM(18)=XNORM(18)+XOL~XRA1 XNORM(II)sXDZ-SD~XRAM XNORM(12)=XHD-SDSXRA1

64

GOTO 600 XNDRM(16)=XNORM(16)+O.SWSB XOLsO,5*SB XNORM(JT)aXNORM(17)÷XOLSXRAM XNORM(IB)=XNORM(18)+XOL*XRAI XNORM(11)=O.O XNORM(12)=O.O SDsSD-SB XNORM(7)mO.5*6D XNORM(6)=XLC'XNORH(7~

C

CALCULATE NORMATIV? MINERALS

600

DO 601 I=I,37

XNORMCI)mXNORM(I)*XFnRM(I) XNnRM(38)zXNDRM(]8)+XNORM(I)

Short Note

601

CONTINUE OXIDE(29)=O,O DO 602 I=1w28 OXIDE(2g)=OXIDE(29)÷OXTDECI)

602

CONTINUE DO 604 I = 1 , 2 9

Ir(OXIDEtI)eEOeOoO) GOTO 604 WRITE(IO,603)

(OXNAMCI),OXIDE(I))

603

FORNAT(IH elSXoA6pgXwF~,2)

604

CONTINUE WRITE(IO,700) DO 606 Imi,38 IF(XNORMCI)oEQ.OtO) GOTO 606 WRITE (10,605) (XNNAH(L,I),L=I,2),XNORM(1)

605

FORMAT(IH ,15X,2A6,3X,?6o2)

606

CONTINUE

XNORN(3g)=XNORM(1)+XNORM(2)÷XNORM(3)+XNORH(5)÷XNORN(6)+XNORN(7~ HOLmXNORM(14)+XNORM(15) HOLImXNORM(I?)+XNORM(18) HOL2~XNORM(II)+XNORM(12)*XNORM(13) HOL3mXNORNC4)+XNORM(3)+O,541985SXNORM(5) XNORMC40)=HOL+HOLI÷HOL2÷XNORM(9)+XHORR(19)+XNORM(20)÷XNORN(21) IF(HOL.EO.O.O)

GOTO 607

HOL=IOO.O/HOL XNORMC41):XNOR~C14)mHOL XNORM(42):XNORM(15)~HOL 607

IFCHOLI.EO.O.O) GOTO 608 HOLIIlOO.O/HOL1 XNORM(43)aXNORMCIT)~HOLI XNORM(44):XNORM(18)*HOL1

60B

IF(HOL2eEO.O.O) GOTO 609 HOL2slOO.O/HOL2 XNORM(45)=XNORM(ll)*HOL2 XNORM(46)=XNORM(12)~HOL2

609

IFCHOL3,EO.O,O) GOTO 610 XNORM(47)=IOO.O~XNORMC4)/HOL3

610 700

WRITE(IO,700) FORNAT(1HO) DO 701 1=39,47

701

WRITE(IO,605) WRITE(IO,700)

1

CONTINUE STOP END

(XNNAHCL,I)eLmle~)

,XNORM(I)

797