C Kazutoshi Ichikawa, University of Cambridge C 26 April 1996 C C Program for getting values of Log likelyfood G(w) for testing data for C best fit and moderated predictions for COMMITTEE MODEL. C C Final version. Thanks Fujii san. C *********************************************************************** C implicit real*8(a-h,k-z), integer(i,j) double precision TARGET(200),OUT(200,200),MO(200,200), &LLFOUT(200),LLFMO(200) C IROWTS is the row number where the test data set begins. C IROWTE is the total number of the row. C MODELMAX is a maxium number of the models for the committee. C C OUT(MODEL No., Row No.) C MO(MODEL No., Row No.) C OPEN(20,FILE='LogLF',STATUS='NEW') IROWTS= 78 IROWTE=154 IMODELMAX=126 CALL GETDATA(IROWTE,TARGET,OUT,MO) IOM=0 DO 10 IMODEL=1,IMODELMAX CALL LOGLF(IMODEL,IROWTS,IROWTE,TARGET,OUT,LLF,I0,I1,IOM) LLFOUT(IMODEL)=LLF IOUT0=I0 IOUT1=I1 10 CONTINUE IOM=1 DO 20 IMODEL=1,IMODELMAX CALL LOGLF(IMODEL,IROWTS,IROWTE,TARGET,MO ,LLF,I0,I1,IOM) LLFMO(IMODEL)=LLF 20 CONTINUE WRITE(20,50) DO 40 IMODEL=1,IMODELMAX WRITE(20,30) IMODEL, LLFOUT(IMODEL), LLFMO(IMODEL),IOUT0,IOUT1 40 CONTINUE 50 FORMAT('No. of Models',' G(w) for _out',' G(w) for mo') 30 FORMAT(I3,' ',F20.6,' ',F20.6,' ',I3,' ',I3) STOP END C *********************************************************************** SUBROUTINE LOGLF(IMODEL,IROWTS,IROWTE,TARGET,Y,LLF,I0,I1,IOM) implicit real*8(a-h,k-z), integer(i,j) double precision TARGET(200),Y(200,200), &YSUM(200,200),YBAR(200),GI(200) C *********************************************************************** C Get the YBAR (Average of the outputs) C DO 30 J=IROWTS,IROWTE DO 10 I=1,IMODEL IF (I .EQ. 1) THEN YSUM(I,J)=Y(I,J) ELSE YSUM(I,J)=YSUM(I-1,J)+Y(I,J) ENDIF 10 CONTINUE DMODEL=DBLE(IMODEL) YBAR(J)=YSUM(IMODEL,J)/DMODEL 30 CONTINUE DO 1 I=IROWTS,IROWTE IF (YBAR(I) .EQ. 0.0) THEN YBAR(I)=0.000001 ENDIF IF (YBAR(I) .EQ. 1.0) THEN YBAR(I)=0.999999 ENDIF IF (TARGET(I) .EQ. 0) THEN GI(I)=DLOG(1.000000-YBAR(I)) ENDIF IF (TARGET(I) .EQ. 1) THEN GI(I)=DLOG(YBAR(I)) ENDIF 1 CONTINUE SUM=0 DO 50 I=IROWTS, IROWTE SUM=SUM+GI(I) 50 CONTINUE LLF=SUM RETURN END C *********************************************************************** SUBROUTINE GETDATA(IROWTE,TARGET,OUT,MO) implicit real*8(a-h,k-z), integer(i,j) double precision TARGET(200),OUT(200,200),MO(200,200) C C *********************************************************************** C Read 'target' file C OPEN(10,FILE='hc5.t',STATUS='OLD') DO 1000 I=1,IROWTE READ(10,*) TARGET(I) 1000 CONTINUE C *********************************************************************** C Read 'outef2/out' files C OPEN(11,FILE='outef2/out1',STATUS='OLD') DO 1 I=1,IROWTE READ(11,*) OUT(1,I) 1 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out2',STATUS='OLD') DO 2 I=1,IROWTE READ(11,*) OUT(2,I) 2 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out3',STATUS='OLD') DO 3 I=1,IROWTE READ(11,*) OUT(3,I) 3 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out4',STATUS='OLD') DO 4 I=1,IROWTE READ(11,*) OUT(4,I) 4 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out5',STATUS='OLD') DO 5 I=1,IROWTE READ(11,*) OUT(5,I) 5 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out6',STATUS='OLD') DO 6 I=1,IROWTE READ(11,*) OUT(6,I) 6 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out7',STATUS='OLD') DO 7 I=1,IROWTE READ(11,*) OUT(7,I) 7 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out8',STATUS='OLD') DO 8 I=1,IROWTE READ(11,*) OUT(8,I) 8 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out9',STATUS='OLD') DO 9 I=1,IROWTE READ(11,*) OUT(9,I) 9 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out10',STATUS='OLD') DO 10 I=1,IROWTE READ(11,*) OUT(10,I) 10 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out11',STATUS='OLD') DO 11 I=1,IROWTE READ(11,*) OUT(11,I) 11 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out12',STATUS='OLD') DO 12 I=1,IROWTE READ(11,*) OUT(12,I) 12 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out13',STATUS='OLD') DO 13 I=1,IROWTE READ(11,*) OUT(13,I) 13 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out14',STATUS='OLD') DO 14 I=1,IROWTE READ(11,*) OUT(14,I) 14 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out15',STATUS='OLD') DO 15 I=1,IROWTE READ(11,*) OUT(15,I) 15 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out16',STATUS='OLD') DO 16 I=1,IROWTE READ(11,*) OUT(16,I) 16 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out17',STATUS='OLD') DO 17 I=1,IROWTE READ(11,*) OUT(17,I) 17 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out18',STATUS='OLD') DO 18 I=1,IROWTE READ(11,*) OUT(18,I) 18 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out19',STATUS='OLD') DO 19 I=1,IROWTE READ(11,*) OUT(19,I) 19 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out20',STATUS='OLD') DO 20 I=1,IROWTE READ(11,*) OUT(20,I) 20 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out21',STATUS='OLD') DO 21 I=1,IROWTE READ(11,*) OUT(21,I) 21 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out22',STATUS='OLD') DO 22 I=1,IROWTE READ(11,*) OUT(22,I) 22 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out23',STATUS='OLD') DO 23 I=1,IROWTE READ(11,*) OUT(23,I) 23 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out24',STATUS='OLD') DO 24 I=1,IROWTE READ(11,*) OUT(24,I) 24 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out25',STATUS='OLD') DO 25 I=1,IROWTE READ(11,*) OUT(25,I) 25 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out26',STATUS='OLD') DO 26 I=1,IROWTE READ(11,*) OUT(26,I) 26 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out27',STATUS='OLD') DO 27 I=1,IROWTE READ(11,*) OUT(27,I) 27 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out28',STATUS='OLD') DO 28 I=1,IROWTE READ(11,*) OUT(28,I) 28 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out29',STATUS='OLD') DO 29 I=1,IROWTE READ(11,*) OUT(29,I) 29 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out30',STATUS='OLD') DO 30 I=1,IROWTE READ(11,*) OUT(30,I) 30 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out31',STATUS='OLD') DO 31 I=1,IROWTE READ(11,*) OUT(31,I) 31 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out32',STATUS='OLD') DO 32 I=1,IROWTE READ(11,*) OUT(32,I) 32 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out33',STATUS='OLD') DO 33 I=1,IROWTE READ(11,*) OUT(33,I) 33 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out34',STATUS='OLD') DO 34 I=1,IROWTE READ(11,*) OUT(34,I) 34 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out35',STATUS='OLD') DO 35 I=1,IROWTE READ(11,*) OUT(35,I) 35 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out36',STATUS='OLD') DO 36 I=1,IROWTE READ(11,*) OUT(36,I) 36 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out37',STATUS='OLD') DO 37 I=1,IROWTE READ(11,*) OUT(37,I) 37 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out38',STATUS='OLD') DO 38 I=1,IROWTE READ(11,*) OUT(38,I) 38 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out39',STATUS='OLD') DO 39 I=1,IROWTE READ(11,*) OUT(39,I) 39 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out40',STATUS='OLD') DO 40 I=1,IROWTE READ(11,*) OUT(40,I) 40 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out41',STATUS='OLD') DO 41 I=1,IROWTE READ(11,*) OUT(41,I) 41 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out42',STATUS='OLD') DO 42 I=1,IROWTE READ(11,*) OUT(42,I) 42 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out43',STATUS='OLD') DO 43 I=1,IROWTE READ(11,*) OUT(43,I) 43 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out44',STATUS='OLD') DO 44 I=1,IROWTE READ(11,*) OUT(44,I) 44 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out45',STATUS='OLD') DO 45 I=1,IROWTE READ(11,*) OUT(45,I) 45 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out46',STATUS='OLD') DO 46 I=1,IROWTE READ(11,*) OUT(46,I) 46 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out47',STATUS='OLD') DO 47 I=1,IROWTE READ(11,*) OUT(47,I) 47 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out48',STATUS='OLD') DO 48 I=1,IROWTE READ(11,*) OUT(48,I) 48 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out49',STATUS='OLD') DO 49 I=1,IROWTE READ(11,*) OUT(49,I) 49 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out50',STATUS='OLD') DO 50 I=1,IROWTE READ(11,*) OUT(50,I) 50 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out51',STATUS='OLD') DO 51 I=1,IROWTE READ(11,*) OUT(51,I) 51 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out52',STATUS='OLD') DO 52 I=1,IROWTE READ(11,*) OUT(52,I) 52 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out53',STATUS='OLD') DO 53 I=1,IROWTE READ(11,*) OUT(53,I) 53 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out54',STATUS='OLD') DO 54 I=1,IROWTE READ(11,*) OUT(54,I) 54 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out55',STATUS='OLD') DO 55 I=1,IROWTE READ(11,*) OUT(55,I) 55 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out56',STATUS='OLD') DO 56 I=1,IROWTE READ(11,*) OUT(56,I) 56 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out57',STATUS='OLD') DO 57 I=1,IROWTE READ(11,*) OUT(57,I) 57 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out58',STATUS='OLD') DO 58 I=1,IROWTE READ(11,*) OUT(58,I) 58 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out59',STATUS='OLD') DO 59 I=1,IROWTE READ(11,*) OUT(59,I) 59 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out60',STATUS='OLD') DO 60 I=1,IROWTE READ(11,*) OUT(60,I) 60 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out71',STATUS='OLD') DO 71 I=1,IROWTE READ(11,*) OUT(71,I) 71 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out72',STATUS='OLD') DO 72 I=1,IROWTE READ(11,*) OUT(72,I) 72 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out73',STATUS='OLD') DO 73 I=1,IROWTE READ(11,*) OUT(73,I) 73 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out74',STATUS='OLD') DO 74 I=1,IROWTE READ(11,*) OUT(74,I) 74 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out75',STATUS='OLD') DO 75 I=1,IROWTE READ(11,*) OUT(75,I) 75 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out76',STATUS='OLD') DO 76 I=1,IROWTE READ(11,*) OUT(76,I) 76 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out77',STATUS='OLD') DO 77 I=1,IROWTE READ(11,*) OUT(77,I) 77 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out78',STATUS='OLD') DO 78 I=1,IROWTE READ(11,*) OUT(78,I) 78 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out79',STATUS='OLD') DO 79 I=1,IROWTE READ(11,*) OUT(79,I) 79 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out80',STATUS='OLD') DO 80 I=1,IROWTE READ(11,*) OUT(80,I) 80 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out81',STATUS='OLD') DO 81 I=1,IROWTE READ(11,*) OUT(81,I) 81 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out82',STATUS='OLD') DO 82 I=1,IROWTE READ(11,*) OUT(82,I) 82 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out83',STATUS='OLD') DO 83 I=1,IROWTE READ(11,*) OUT(83,I) 83 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out84',STATUS='OLD') DO 84 I=1,IROWTE READ(11,*) OUT(84,I) 84 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out85',STATUS='OLD') DO 85 I=1,IROWTE READ(11,*) OUT(85,I) 85 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out86',STATUS='OLD') DO 86 I=1,IROWTE READ(11,*) OUT(86,I) 86 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out87',STATUS='OLD') DO 87 I=1,IROWTE READ(11,*) OUT(87,I) 87 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out88',STATUS='OLD') DO 88 I=1,IROWTE READ(11,*) OUT(88,I) 88 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out89',STATUS='OLD') DO 89 I=1,IROWTE READ(11,*) OUT(89,I) 89 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out90',STATUS='OLD') DO 90 I=1,IROWTE READ(11,*) OUT(90,I) 90 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out91',STATUS='OLD') DO 91 I=1,IROWTE READ(11,*) OUT(91,I) 91 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out92',STATUS='OLD') DO 92 I=1,IROWTE READ(11,*) OUT(92,I) 92 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out93',STATUS='OLD') DO 93 I=1,IROWTE READ(11,*) OUT(93,I) 93 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out94',STATUS='OLD') DO 94 I=1,IROWTE READ(11,*) OUT(94,I) 94 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out95',STATUS='OLD') DO 95 I=1,IROWTE READ(11,*) OUT(95,I) 95 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out96',STATUS='OLD') DO 96 I=1,IROWTE READ(11,*) OUT(96,I) 96 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out97',STATUS='OLD') DO 97 I=1,IROWTE READ(11,*) OUT(97,I) 97 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out98',STATUS='OLD') DO 98 I=1,IROWTE READ(11,*) OUT(98,I) 98 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out99',STATUS='OLD') DO 99 I=1,IROWTE READ(11,*) OUT(99,I) 99 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out100',STATUS='OLD') DO 100 I=1,IROWTE READ(11,*) OUT(100,I) 100 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out101',STATUS='OLD') DO 101 I=1,IROWTE READ(11,*) OUT(101,I) 101 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out102',STATUS='OLD') DO 102 I=1,IROWTE READ(11,*) OUT(102,I) 102 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out103',STATUS='OLD') DO 103 I=1,IROWTE READ(11,*) OUT(103,I) 103 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out104',STATUS='OLD') DO 104 I=1,IROWTE READ(11,*) OUT(104,I) 104 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out105',STATUS='OLD') DO 105 I=1,IROWTE READ(11,*) OUT(105,I) 105 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out106',STATUS='OLD') DO 106 I=1,IROWTE READ(11,*) OUT(106,I) 106 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out107',STATUS='OLD') DO 107 I=1,IROWTE READ(11,*) OUT(107,I) 107 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out108',STATUS='OLD') DO 108 I=1,IROWTE READ(11,*) OUT(108,I) 108 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out109',STATUS='OLD') DO 109 I=1,IROWTE READ(11,*) OUT(109,I) 109 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out110',STATUS='OLD') DO 110 I=1,IROWTE READ(11,*) OUT(110,I) 110 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out111',STATUS='OLD') DO 111 I=1,IROWTE READ(11,*) OUT(111,I) 111 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out112',STATUS='OLD') DO 112 I=1,IROWTE READ(11,*) OUT(112,I) 112 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out113',STATUS='OLD') DO 113 I=1,IROWTE READ(11,*) OUT(113,I) 113 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out114',STATUS='OLD') DO 114 I=1,IROWTE READ(11,*) OUT(114,I) 114 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out115',STATUS='OLD') DO 115 I=1,IROWTE READ(11,*) OUT(115,I) 115 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out116',STATUS='OLD') DO 116 I=1,IROWTE READ(11,*) OUT(116,I) 116 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out117',STATUS='OLD') DO 117 I=1,IROWTE READ(11,*) OUT(117,I) 117 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out118',STATUS='OLD') DO 118 I=1,IROWTE READ(11,*) OUT(118,I) 118 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out119',STATUS='OLD') DO 119 I=1,IROWTE READ(11,*) OUT(119,I) 119 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out120',STATUS='OLD') DO 120 I=1,IROWTE READ(11,*) OUT(120,I) 120 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out121',STATUS='OLD') DO 121 I=1,IROWTE READ(11,*) OUT(121,I) 121 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out122',STATUS='OLD') DO 122 I=1,IROWTE READ(11,*) OUT(122,I) 122 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out123',STATUS='OLD') DO 123 I=1,IROWTE READ(11,*) OUT(123,I) 123 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out124',STATUS='OLD') DO 124 I=1,IROWTE READ(11,*) OUT(124,I) 124 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out125',STATUS='OLD') DO 125 I=1,IROWTE READ(11,*) OUT(125,I) 125 CONTINUE CLOSE(11) OPEN(11,FILE='outef2/out126',STATUS='OLD') DO 126 I=1,IROWTE READ(11,*) OUT(126,I) 126 CONTINUE CLOSE(11) C C *********************************************************************** C Read moef2/mo files C OPEN(12,FILE='moef2/mo1',STATUS='OLD') DO 1001 I=1,IROWTE READ(12,*) MO(1,I) 1001 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo2',STATUS='OLD') DO 1002 I=1,IROWTE READ(12,*) MO(2,I) 1002 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo3',STATUS='OLD') DO 1003 I=1,IROWTE READ(12,*) MO(3,I) 1003 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo4',STATUS='OLD') DO 1004 I=1,IROWTE READ(12,*) MO(4,I) 1004 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo5',STATUS='OLD') DO 1005 I=1,IROWTE READ(12,*) MO(5,I) 1005 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo6',STATUS='OLD') DO 1006 I=1,IROWTE READ(12,*) MO(6,I) 1006 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo7',STATUS='OLD') DO 1007 I=1,IROWTE READ(12,*) MO(7,I) 1007 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo8',STATUS='OLD') DO 1008 I=1,IROWTE READ(12,*) MO(8,I) 1008 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo9',STATUS='OLD') DO 1009 I=1,IROWTE READ(12,*) MO(9,I) 1009 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo10',STATUS='OLD') DO 1010 I=1,IROWTE READ(12,*) MO(10,I) 1010 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo11',STATUS='OLD') DO 1011 I=1,IROWTE READ(12,*) MO(11,I) 1011 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo12',STATUS='OLD') DO 1012 I=1,IROWTE READ(12,*) MO(12,I) 1012 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo13',STATUS='OLD') DO 1013 I=1,IROWTE READ(12,*) MO(13,I) 1013 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo14',STATUS='OLD') DO 1014 I=1,IROWTE READ(12,*) MO(14,I) 1014 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo15',STATUS='OLD') DO 1015 I=1,IROWTE READ(12,*) MO(15,I) 1015 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo16',STATUS='OLD') DO 1016 I=1,IROWTE READ(12,*) MO(16,I) 1016 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo17',STATUS='OLD') DO 1017 I=1,IROWTE READ(12,*) MO(17,I) 1017 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo18',STATUS='OLD') DO 1018 I=1,IROWTE READ(12,*) MO(18,I) 1018 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo19',STATUS='OLD') DO 1019 I=1,IROWTE READ(12,*) MO(19,I) 1019 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo20',STATUS='OLD') DO 1020 I=1,IROWTE READ(12,*) MO(20,I) 1020 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo21',STATUS='OLD') DO 1021 I=1,IROWTE READ(12,*) MO(21,I) 1021 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo22',STATUS='OLD') DO 1022 I=1,IROWTE READ(12,*) MO(22,I) 1022 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo23',STATUS='OLD') DO 1023 I=1,IROWTE READ(12,*) MO(23,I) 1023 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo24',STATUS='OLD') DO 1024 I=1,IROWTE READ(12,*) MO(24,I) 1024 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo25',STATUS='OLD') DO 1025 I=1,IROWTE READ(12,*) MO(25,I) 1025 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo26',STATUS='OLD') DO 1026 I=1,IROWTE READ(12,*) MO(26,I) 1026 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo27',STATUS='OLD') DO 1027 I=1,IROWTE READ(12,*) MO(27,I) 1027 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo28',STATUS='OLD') DO 1028 I=1,IROWTE READ(12,*) MO(28,I) 1028 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo29',STATUS='OLD') DO 1029 I=1,IROWTE READ(12,*) MO(29,I) 1029 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo30',STATUS='OLD') DO 1030 I=1,IROWTE READ(12,*) MO(30,I) 1030 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo31',STATUS='OLD') DO 1031 I=1,IROWTE READ(12,*) MO(31,I) 1031 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo32',STATUS='OLD') DO 1032 I=1,IROWTE READ(12,*) MO(32,I) 1032 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo33',STATUS='OLD') DO 1033 I=1,IROWTE READ(12,*) MO(33,I) 1033 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo34',STATUS='OLD') DO 1034 I=1,IROWTE READ(12,*) MO(34,I) 1034 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo35',STATUS='OLD') DO 1035 I=1,IROWTE READ(12,*) MO(35,I) 1035 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo36',STATUS='OLD') DO 1036 I=1,IROWTE READ(12,*) MO(36,I) 1036 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo37',STATUS='OLD') DO 1037 I=1,IROWTE READ(12,*) MO(37,I) 1037 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo38',STATUS='OLD') DO 1038 I=1,IROWTE READ(12,*) MO(38,I) 1038 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo39',STATUS='OLD') DO 1039 I=1,IROWTE READ(12,*) MO(39,I) 1039 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo40',STATUS='OLD') DO 1040 I=1,IROWTE READ(12,*) MO(40,I) 1040 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo41',STATUS='OLD') DO 1041 I=1,IROWTE READ(12,*) MO(41,I) 1041 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo42',STATUS='OLD') DO 1042 I=1,IROWTE READ(12,*) MO(42,I) 1042 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo43',STATUS='OLD') DO 1043 I=1,IROWTE READ(12,*) MO(43,I) 1043 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo44',STATUS='OLD') DO 1044 I=1,IROWTE READ(12,*) MO(44,I) 1044 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo45',STATUS='OLD') DO 1045 I=1,IROWTE READ(12,*) MO(45,I) 1045 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo46',STATUS='OLD') DO 1046 I=1,IROWTE READ(12,*) MO(46,I) 1046 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo47',STATUS='OLD') DO 1047 I=1,IROWTE READ(12,*) MO(47,I) 1047 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo48',STATUS='OLD') DO 1048 I=1,IROWTE READ(12,*) MO(48,I) 1048 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo49',STATUS='OLD') DO 1049 I=1,IROWTE READ(12,*) MO(49,I) 1049 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo50',STATUS='OLD') DO 1050 I=1,IROWTE READ(12,*) MO(50,I) 1050 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo51',STATUS='OLD') DO 1051 I=1,IROWTE READ(12,*) MO(51,I) 1051 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo52',STATUS='OLD') DO 1052 I=1,IROWTE READ(12,*) MO(52,I) 1052 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo53',STATUS='OLD') DO 1053 I=1,IROWTE READ(12,*) MO(53,I) 1053 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo54',STATUS='OLD') DO 1054 I=1,IROWTE READ(12,*) MO(54,I) 1054 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo55',STATUS='OLD') DO 1055 I=1,IROWTE READ(12,*) MO(55,I) 1055 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo56',STATUS='OLD') DO 1056 I=1,IROWTE READ(12,*) MO(56,I) 1056 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo57',STATUS='OLD') DO 1057 I=1,IROWTE READ(12,*) MO(57,I) 1057 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo58',STATUS='OLD') DO 1058 I=1,IROWTE READ(12,*) MO(58,I) 1058 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo59',STATUS='OLD') DO 1059 I=1,IROWTE READ(12,*) MO(59,I) 1059 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo60',STATUS='OLD') DO 1060 I=1,IROWTE READ(12,*) MO(60,I) 1060 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo71',STATUS='OLD') DO 1071 I=1,IROWTE READ(12,*) MO(71,I) 1071 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo72',STATUS='OLD') DO 1072 I=1,IROWTE READ(12,*) MO(72,I) 1072 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo73',STATUS='OLD') DO 1073 I=1,IROWTE READ(12,*) MO(73,I) 1073 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo74',STATUS='OLD') DO 1074 I=1,IROWTE READ(12,*) MO(74,I) 1074 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo75',STATUS='OLD') DO 1075 I=1,IROWTE READ(12,*) MO(75,I) 1075 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo76',STATUS='OLD') DO 1076 I=1,IROWTE READ(12,*) MO(76,I) 1076 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo77',STATUS='OLD') DO 1077 I=1,IROWTE READ(12,*) MO(77,I) 1077 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo78',STATUS='OLD') DO 1078 I=1,IROWTE READ(12,*) MO(78,I) 1078 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo79',STATUS='OLD') DO 1079 I=1,IROWTE READ(12,*) MO(79,I) 1079 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo80',STATUS='OLD') DO 1080 I=1,IROWTE READ(12,*) MO(80,I) 1080 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo81',STATUS='OLD') DO 1081 I=1,IROWTE READ(12,*) MO(81,I) 1081 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo82',STATUS='OLD') DO 1082 I=1,IROWTE READ(12,*) MO(82,I) 1082 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo83',STATUS='OLD') DO 1083 I=1,IROWTE READ(12,*) MO(83,I) 1083 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo84',STATUS='OLD') DO 1084 I=1,IROWTE READ(12,*) MO(84,I) 1084 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo85',STATUS='OLD') DO 1085 I=1,IROWTE READ(12,*) MO(85,I) 1085 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo86',STATUS='OLD') DO 1086 I=1,IROWTE READ(12,*) MO(86,I) 1086 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo87',STATUS='OLD') DO 1087 I=1,IROWTE READ(12,*) MO(87,I) 1087 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo88',STATUS='OLD') DO 1088 I=1,IROWTE READ(12,*) MO(88,I) 1088 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo89',STATUS='OLD') DO 1089 I=1,IROWTE READ(12,*) MO(89,I) 1089 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo90',STATUS='OLD') DO 1090 I=1,IROWTE READ(12,*) MO(90,I) 1090 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo91',STATUS='OLD') DO 1091 I=1,IROWTE READ(12,*) MO(91,I) 1091 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo92',STATUS='OLD') DO 1092 I=1,IROWTE READ(12,*) MO(92,I) 1092 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo93',STATUS='OLD') DO 1093 I=1,IROWTE READ(12,*) MO(93,I) 1093 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo94',STATUS='OLD') DO 1094 I=1,IROWTE READ(12,*) MO(94,I) 1094 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo95',STATUS='OLD') DO 1095 I=1,IROWTE READ(12,*) MO(95,I) 1095 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo96',STATUS='OLD') DO 1096 I=1,IROWTE READ(12,*) MO(96,I) 1096 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo97',STATUS='OLD') DO 1097 I=1,IROWTE READ(12,*) MO(97,I) 1097 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo98',STATUS='OLD') DO 1098 I=1,IROWTE READ(12,*) MO(98,I) 1098 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo99',STATUS='OLD') DO 1099 I=1,IROWTE READ(12,*) MO(99,I) 1099 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo100',STATUS='OLD') DO 1100 I=1,IROWTE READ(12,*) MO(100,I) 1100 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo101',STATUS='OLD') DO 1101 I=1,IROWTE READ(12,*) MO(101,I) 1101 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo102',STATUS='OLD') DO 1102 I=1,IROWTE READ(12,*) MO(102,I) 1102 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo103',STATUS='OLD') DO 1103 I=1,IROWTE READ(12,*) MO(103,I) 1103 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo104',STATUS='OLD') DO 1104 I=1,IROWTE READ(12,*) MO(104,I) 1104 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo105',STATUS='OLD') DO 1105 I=1,IROWTE READ(12,*) MO(105,I) 1105 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo106',STATUS='OLD') DO 1106 I=1,IROWTE READ(12,*) MO(106,I) 1106 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo107',STATUS='OLD') DO 1107 I=1,IROWTE READ(12,*) MO(107,I) 1107 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo108',STATUS='OLD') DO 1108 I=1,IROWTE READ(12,*) MO(108,I) 1108 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo109',STATUS='OLD') DO 1109 I=1,IROWTE READ(12,*) MO(109,I) 1109 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo110',STATUS='OLD') DO 1110 I=1,IROWTE READ(12,*) MO(110,I) 1110 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo111',STATUS='OLD') DO 1111 I=1,IROWTE READ(12,*) MO(111,I) 1111 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo112',STATUS='OLD') DO 1112 I=1,IROWTE READ(12,*) MO(112,I) 1112 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo113',STATUS='OLD') DO 1113 I=1,IROWTE READ(12,*) MO(113,I) 1113 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo114',STATUS='OLD') DO 1114 I=1,IROWTE READ(12,*) MO(114,I) 1114 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo115',STATUS='OLD') DO 1115 I=1,IROWTE READ(12,*) MO(115,I) 1115 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo116',STATUS='OLD') DO 1116 I=1,IROWTE READ(12,*) MO(116,I) 1116 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo117',STATUS='OLD') DO 1117 I=1,IROWTE READ(12,*) MO(117,I) 1117 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo118',STATUS='OLD') DO 1118 I=1,IROWTE READ(12,*) MO(118,I) 1118 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo119',STATUS='OLD') DO 1119 I=1,IROWTE READ(12,*) MO(119,I) 1119 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo120',STATUS='OLD') DO 1120 I=1,IROWTE READ(12,*) MO(120,I) 1120 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo121',STATUS='OLD') DO 1121 I=1,IROWTE READ(12,*) MO(121,I) 1121 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo122',STATUS='OLD') DO 1122 I=1,IROWTE READ(12,*) MO(122,I) 1122 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo123',STATUS='OLD') DO 1123 I=1,IROWTE READ(12,*) MO(123,I) 1123 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo124',STATUS='OLD') DO 1124 I=1,IROWTE READ(12,*) MO(124,I) 1124 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo125',STATUS='OLD') DO 1125 I=1,IROWTE READ(12,*) MO(125,I) 1125 CONTINUE CLOSE(12) OPEN(12,FILE='moef2/mo126',STATUS='OLD') DO 1126 I=1,IROWTE READ(12,*) MO(126,I) 1126 CONTINUE CLOSE(12) RETURN END C ***********************************************************************