c Program to explore secondary structures of protein sequences c LSTMOD 17-3-2008 c scope: ? helices,Sheets & Turns 72% character*1 AACID(1000),RA(21),secs(1000,5) character*1 cnatv(1000),cpsip(1000),csabl(1000),csecsc(1000) character*1000 SEQ,seq0,seq1,seq2,seq3 character*200 flname1,flname2,label real amh,ams,amt,amh1,ams1,amt1 real sumh(1000),sums(1000),sumt(1000),coeff(3) real ha(280),sa(280),ta(280),DA(84) real ha1(160),sa1(160),ta1(160) integer mht,mst,mtt,IR(1000) integer mhtrue(1000),mstrue(1000),mttrue(1000) real sd1(1000),sd2(1000),sd3(1000),sd4(1000),avg(4) real sumhe,sumsh,sumtu,nh(1000),ns(1000),nc(1000) c------------------------------------------------------------------------- data RA/'A','C','D','E','F','G','H','I','K','L','M', 1'N','P','Q','R','S','T','V','W','Y','Z'/ data coeff/1.00,1.00,1.00/ ! 65.2% c data coeff/1.50,1.45,1.536/ ! 65.2% c data coeff/1.000,1.,0./ ! 80.6% c data coeff/0.00,1.015,1./ ! 80.6% data avg/3.562,2.068,3.460,2.910/ c data DA/0,0,2,1,0,1,2,0,0,0,1,2,1,0,0,2,0,0,0,3, c 1 0,0,3,0,0,2,0,1,2,0,1,0,1,1,0,1,3,0,0,0, c 2 1,0,1,1,0,2,1,0,2,0,0,1,1,2,0,0,2,1,0,0, c 3 0,1,1,1,1,1,1,0,1,0,2,0,0,3,0,0,0,1,0,2, c 4 0,0,0,0/ data DA/0,0,1,1,0,1,1,1,0,0,1,1,1,0,0,1,0,0,0,1, 1 0,0,1,0,0,1,0,1,1,0,1,0,1,1,0,1,1,0,0,0, 2 1,0,0,1,0,1,1,0,1,0,1,1,1,1,0,0,1,1,0,0, 3 0,1,1,0,1,1,1,0,1,0,1,0,0,1,0,0,0,1,0,1, 4 0,0,0,0/ c--------------------------------------------------------------------------------- print *,' ' write(6,2007) print *,'WELCOME TO PROTEIN SECONDARY STRUCTURE CHARACTERIZATIO 1N (PROSECSC)' 50 format(a1000) 60 format(a200) 1000 format(A80) 2000 format(80A1) write(6,2007) c write(6,2003) c write(6,2007) c read(5,60)flname1 c flname1='test.dat' flname2='prosecsc.out' open(unit=10,file=''//flname1//'',form='formatted',status= 1 'old') open(unit=11,file=''//flname2//'',form='formatted',status= 1 'unknown') c c nfile=9999 nfile=1 pr=0 do 4000 ifile=1,nfile mhelix=0 msheet=0 mturn=0 lfile=0 mshort=0 tp=0 c do k=1,1000 nh(k)=0 ns(k)=0 cnatv(k)='-' cpsip(k)='-' csabl(k)='-' csecsc(k)='-' end do c do 4000 ifile=1,nfile c do 4000 ifile=1,5 read(10,1000)label if (label(1:1) .eq. '>')goto 19 if (label(1:1) .eq. '$')goto 6000 if (label(1:1) .eq. ' ')goto 6000 19 read(10,50,end=6000)SEQ c if (SEQ(1:1) .eq. ' ')goto 17 c read(10,1000)SEQ0 17 read(10,1000)label if (label(1:1) .eq. '>')goto 21 if (label(1:1) .ne. '>')goto 6000 21 read(10,50)SEQ1 read(10,1000)label if (label(1:1) .eq. '>')goto 23 if (label(1:1) .ne. '>')goto 6000 23 read(10,50)SEQ2 c print *,'type sequence' c read(5,1000,end=6000)SEQ NAA=0 do 10 i=1,1000 if (seq(i:i) .eq. ' ')goto 10 NAA=NAA+1 AACID(NAA)=SEQ(i:i) c cnatv(NAA)=seq0(i:i) cpsip(NAA)=seq1(i:i) csabl(NAA)=seq2(i:i) do 20 l=1,20 if (AACID(NAA) .eq. RA(l))IR(NAA)=l 20 continue 10 continue c if (NAA .lt. 13)then c mshort=mshort+1 c goto 4000 c end if c lfile=lfile+1 if (lfile .gt. 5000) then lfile=5000 goto 4001 end if c 4004 continue c write(11,2005)(AACID(k),k=1,NAA) do 4003 ifl=1,NAA c do 4010 ij=ifl,ifl+9 do k=1,100 secs(k,1)='-' secs(k,2)='-' secs(k,3)='-' secs(k,4)='-' mhtrue(k)=0 mstrue(k)=0 mttrue(k)=0 sumh(k)=0 sums(k)=0 sumt(k)=0 enddo c do 4010 ij=ifl,ifl+10 sumhe=0 sumsh=0 sumtu=0 c do 4010 ij=ifl,ifl+10 mht=0 mst=0 mtt=0 c if (NAA .gt. 30)NAA=30 do 4010 ij=1,NAA-10 c do 4010 ij=ifl,ifl+11 ib=ij-3 do 310 m=ij+1,ij+7 AMH=0 AMS=0 AMT=0 AMH1=0 AMS1=0 AMT1=0 ib=ib+1 if (ib .ge. 11)ib=6 l=IR(m) j=(l-1)*14+ib AMH=HA(j) AMS=SA(j) AMT=TA(j) do 311 iu=1,8 l=IR(m) j=(l-1)*8+iu AMH1=AMH1+HA1(j) AMS1=AMS1+SA1(j) AMT1=AMT1+TA1(j) 311 continue c sumh(ij)=sumh(ij)+amh*coeff(1)+amh1*12 sums(ij)=sums(ij)+ams*coeff(2)+ams1*12 sumt(ij)=sumt(ij)+amt*coeff(3)+amt1*12 310 continue sumhe=sumhe+sumh(ij) sumsh=sumsh+sums(ij) sumtu=sumtu+sumt(ij) if (sumhe .gt. sumsh )mht=mht+1 if (sumsh .gt. sumhe )mst=mst+1 if (sumh(ij) .gt. sumh(ij-5))mht=mht+1 if (sums(ij) .gt. sums(ij-5))mst=mst+1 if (sumtu .gt. sumhe)mtt=mtt+1 if (sumtu .gt. sumsh)mtt=mtt+1 if (sumt(ij) .gt. sumt(ij-1))mtt=mtt+1 ijp=ij+5 if (ijp .gt. NAA-4)goto 4010 if (mht .gt. mst )secs(ijp,4)='H' if (mst .gt. mht )secs(ijp,4)='E' c if (mht .gt. mst )secs(ij,4)='H' c if (mst .ge. mht )secs(ij,4)='E' if (mtt .gt. mht .and. mtt .gt. mst)secs(ijp,4)='C' c if (mtt .gt. mht .and. mtt .gt. mst)secs(ij,4)='C' c goto 4010 sumcnt=sumcnt+1 c if (ij .lt. 4 )goto 4010 d1b=0 d2b=0 d3b=0 d4b=0 d1f=0 d2f=0 d3f=0 d4f=0 delg=0 deld=0 dels=0 dell=0 if (ij .lt. 4)goto 7300 do mi=1,NAA sd1(mi)=0 sd2(mi)=0 sd3(mi)=0 sd4(mi)=0 end do do 320 mi=ij-3,ij l=IR(mi) j=(l-1)*4+1 d1b=d1b+DA(j) d2b=d2b+DA(j+1) d3b=d3b+DA(j+2) d4b=d4b+DA(j+3) sd1(mi)=DA(j) sd2(mi)=DA(j+1) sd3(mi)=DA(j+2) sd4(mi)=DA(j+3) 320 continue do 321 mi=ij+1,ij+4 c do 321 mi=ij-1,ij+6 l=IR(mi) j=(l-1)*4+1 d1f=d1f+DA(j) d2f=d2f+DA(j+1) d3f=d3f+DA(j+2) d4f=d4f+DA(j+3) sd1(mi)=DA(j) sd2(mi)=DA(j+1) sd3(mi)=DA(j+2) sd4(mi)=DA(j+3) 321 continue c if ((d2f+d4f)*3.0 .gt. d1f+d3f)then c mht=mht+1 c else c mst=mst+1 c end if c goto 4010 sd12=0 sd22=0 sd32=0 sd42=0 do mi=ij+1,ij+4 sd12=sd12+(sd1(mi)-sd1(mi-1))**2 sd22=sd22+(sd2(mi)-sd2(mi-1))**2 sd32=sd32+(sd3(mi)-sd3(mi-1))**2 sd42=sd42+(sd4(mi)-sd4(mi-1))**2 end do sd12=sqrt(sd12) sd22=sqrt(sd22) sd32=sqrt(sd32) sd42=sqrt(sd42) c delg=d1f-d1b deld=d2f-d2b dels=d3f-d3b dell=d4f-d4b sumd1b=sumd1b+d1b sumd2b=sumd2b+d2b sumd3b=sumd3b+d3b sumd4b=sumd4b+d4b sumd1f=sumd1f+d1f sumd2f=sumd2f+d2f sumd3f=sumd3f+d3f sumd4f=sumd4f+d4f sumdelg=sumdelg+delg sumdeld=sumdeld+deld sumdels=sumdels+dels sumdell=sumdell+dell sumd12=sumd12+sd12 sumd22=sumd22+sd22 sumd32=sumd32+sd32 sumd42=sumd42+sd42 d1423=(d1f*d4f-d2f*d3f) d1324=(d1f*d3f-d2f*d4f) d1234=(d1f*d2f-d3f*d4f) sd1234=sd1234+d1234 sd1324=sd1324+d1324 sd1423=sd1423+d1423 xs=d1423+d1234 sumxs=sumxs+xs c 7200 continue c------------------------------------------------------------ c if (d1b .le. 3)mhtrue(ij)=mhtrue(ij)-1 c if (d2b .lt. 3)mhtrue(ij)=mhtrue(ij)-2 c if (d3b .ge. 3)mhtrue(ij)=mhtrue(ij)-2 c if (d4b .gt. 3)mhtrue(ij)=mhtrue(ij)-1 c if (d1f .lt. 4. .and. d3f .lt. 4)mhtrue(ij)=mhtrue(ij)+3 c if (d1f .gt. 4. .and. d3f .gt. 4)mstrue(ij)=mstrue(ij)+3 c if (d2f .gt. 2.)mhtrue(ij)=mhtrue(ij)+2 c if (d3f .lt. 4.)mhtrue(ij)=mhtrue(ij)+3 c if (d4f .gt. 2.)mhtrue(ij)=mhtrue(ij)+3 if (delg .le. 0 .and. deld .ge. 0)mhtrue(ij)=mhtrue(ij)+1 if (delg .le. 0 .and. dell .ge. 0)mhtrue(ij)=mhtrue(ij)+1 if (delg .gt. 0 .and. dell .le. 0)mstrue(ij)=mstrue(ij)+1 if (delg .le. 0 .and. dels .le. 0)mhtrue(ij)=mhtrue(ij)+1 if (delg .gt. 0 .and. dels .gt. 0)mstrue(ij)=mstrue(ij)+1 c if (sd12 .gt. 3.5)mhtrue(ij)=mhtrue(ij)-1 c if (sd12 .lt. 3.0)mstrue(ij)=mstrue(ij)-1 if (delg .gt. 0 .and. dels .ge. 0)mstrue(ij)=mstrue(ij)+1 c if (delg .gt. 0)mhtrue(ij)=mhtrue(ij)-1 c if (deld .lt. 0)mhtrue(ij)=mhtrue(ij)-1 c if (dels .gt. 0)mhtrue(ij)=mhtrue(ij)-2 c if (dell .lt. 0)mhtrue(ij)=mhtrue(ij)-2 c if (d1324 .lt. 3.0)mhtrue(ij)=mhtrue(ij)-1 c if (d1324 .gt. 5.0)mstrue(ij)=mstrue(ij)-1 c if (d1423 .lt. 2.0)mhtrue(ij)=mhtrue(ij)-1 c if (d1234 .gt. -1.)mhtrue(ij)=mhtrue(ij)-2 c if (sd12 .gt. 3.)mhtrue(ij)=mhtrue(ij)-1 c if (xs .le. 0)mhtrue(ij)=mhtrue(ij)+1 c if (xs .ge. 0)mstrue(ij)=mstrue(ij)+1 c------------------------------------------------------------ c if (d1b .ge. 3)mstrue(ij)=mstrue(ij)-1 c if (d2b .ge. 4)mstrue(ij)=mstrue(ij)-1 c if (d3b .lt. 2)mstrue(ij)=mstrue(ij)-2 c if (d4b .le. 2)mstrue(ij)=mstrue(ij)-1 c if (d1f .gt. 2.)mstrue(ij)=mstrue(ij)+2 c if (d2f .lt. 3.)mstrue(ij)=mstrue(ij)+2 c if (d3f .gt. 3.)mstrue(ij)=mstrue(ij)+2 c if (d4f .lt. 4.)mstrue(ij)=mstrue(ij)+3 c if (delg .lt. 0)mstrue(ij)=mstrue(ij)-2 c if (deld .gt. 0)mstrue(ij)=mstrue(ij)-1 c if (dels .lt. 0)mstrue(ij)=mstrue(ij)-1 c if (dell .gt. 0)mstrue(ij)=mstrue(ij)-2 c if (d1324 .lt. 3)mstrue(ij)=mstrue(ij)-2 c if (d1423 .gt. 4)mstrue(ij)=mstrue(ij)-1 c if (d1234 .lt.-3)mstrue(ij)=mstrue(ij)-2 c if (sd12 .lt. 2.)mstrue(ij)=mstrue(ij)-1 c if (xs .gt. 0.)mstrue(ij)=mstrue(ij)+1 c------------------------------------------------------------ c if (d1b .le. 3)mttrue(ij)=mttrue(ij)-1 c if (d2b .le. 1)mttrue(ij)=mttrue(ij)+1 c if (d3b .le. 2)mttrue(ij)=mttrue(ij)+3 c if (d4b .le. 2)mttrue(ij)=mttrue(ij)+1 if (d1f .lt. 1 .and. d3f .gt. 4)mttrue(ij)=mttrue(ij)+1 if (d2f .gt. 4)mttrue(ij)=mttrue(ij)+1 c if (d3f .gt. 4)mttrue(ij)=mttrue(ij)+1 c if (d4f .lt. 2)mttrue(ij)=mttrue(ij)+1 if (delg .lt. 0)mttrue(ij)=mttrue(ij)+1 if (deld .lt. 0)mttrue(ij)=mttrue(ij)-1 c if (dels .lt. 0)mttrue(ij)=mttrue(ij)-1 c if (dell .lt. 0)mttrue(ij)=mttrue(ij)-1 c if (sd12 .le. 2.)mttrue(ij)=mttrue(ij)-2 c if (sd22 .ge. 1.5)mttrue(ij)=mttrue(ij)-2 c if (sd32 .le. 2.)mttrue(ij)=mttrue(ij)-2 c if (sd42 .ge. 2.)mttrue(ij)=mttrue(ij)-2 c if (xs .lt. -1.)mttrue(ij)=mttrue(ij)+1 c------------------------------------------------------------ c7300 continue c if (mhtrue(ij-1) .le. mhtrue(ij)) mht=mht+1 if (mstrue(ij-1) .lt. mstrue(ij)) mst=mst+1 if (mhtrue(ij) .ge. mstrue(ij-4)) mht=mht+1 if (mstrue(ij) .gt. mhtrue(ij-4)) mst=mst+1 if (mttrue(ij) .gt. mhtrue(ij).and. mttrue(ij).gt. mstrue(ij)) 1 mtt=mtt+1 if (mttrue(ij) .gt. mttrue(ij-4)) mtt=mtt+1 7300 continue c---------------------------------------------------- 4010 continue c---------------------------------------------------- c mh=0 c ms=0 c mt=0 c do 1005 k=5,NAA-5 c do 1005 k=1,NAA c if (secs(k,4) .eq. 'H')mh=mh+1 c if (secs(k,4) .eq. 'E')ms=ms+1 c if (secs(k,4) .eq. 'C')mt=mt+1 c if (secs(k,4) .eq. '-')mh=mh+1 c1005 continue c do 1006 k=6,NAA-5 c if (mh.gt.ms)secs(k,4)='H' c if (ms.ge.mh)secs(k,4)='E' c1006 continue c do 1007 k=1,5 c secs(k,4)='-' c secs(NAA-k+1,4)='-' c1007 continue c c c if (mh.gt.ms .and. mh.gt.mt)then c mhelix=mhelix+1 c end if c if (ms.gt.mh .and. ms.gt.mt)then c msheet=msheet+1 c end if c if (mt.gt.mh .and. mt.gt.ms)then c mturn=mturn+1 c end if c3000 format(i10) c c write(11,2005)(AACID(k),k=1,NAA) c write(11,2005)(secs(k,4),k=1,NAA) c write(6,2005)(AACID(k),k=1,NAA) c write(6,2005)(secs(k,4),k=1,NAA) c sumcnt=sumcnt+1 do k=1,NAA nh(k)=0 ns(k)=0 if (secs(k,4) .eq. 'H')nh(k)=nh(k)+1 if (secs(k,4) .eq. 'E')ns(k)=ns(k)+1 end do 4003 continue do k=1,NAA secs(k,5)='C' if (nh(k) .gt. ns(k))secs(k,5)='H' if (ns(k) .gt. nh(k))secs(k,5)='E' end do do j=1,NAA-2 if (secs(j-1,5) .eq. 'C' .and. secs(j,5) .eq. 'H' .and. 1 secs(j+1,5) .eq. 'E' .and. secs(j+2,5) .eq. 'E')secs(j,5)='E' if (secs(j-1,5) .eq. 'C' .and. secs(j,5) .eq. 'E' .and. 1 secs(j+1,5) .eq. 'H' .and. secs(j+2,5) .eq. 'H')secs(j,5)='H' if (secs(j-2,5) .eq. 'C' .and. secs(j-1,5) .eq. 'H' .and. 1 secs(j,5) .eq. 'C' .and. secs(j+1,5) .eq. 'E' .and. 1 secs(j+2,5) .eq. 'E')secs(j,5)='E' if (secs(j-2,5) .eq. 'C' .and. secs(j-1,5) .eq. 'E' .and. 1 secs(j,5) .eq. 'C' .and. secs(j+1,5) .eq. 'H' .and. 1 secs(j+2,5) .eq. 'H')secs(j,5)='H' end do ne=0 do 6001 k=1,NAA-2 if (cpsip(k) .eq. 'E' .and. cpsip(k+1) .eq. 'E')then do 6002 j=k+3,NAA-2 if (cpsip(j) .eq. 'E' .and. cpsip(j+1) .eq. 'E')ne=ne+1 6002 continue if (ne .eq. 0)then cpsip(k)='C' cpsip(k+1)='C' end if end if 6001 continue do j=2,NAA-3 if (secs(j-1,5) .eq. 'C' .and. secs(j,5) .eq. 'E' .and. 1 secs(j+1,5) .eq. 'C' .and. secs(j+2,5) .eq. 'C')secs(j-1,5)='E' end do do j=2,NAA-2 if (secs(j-1,5) .eq. 'H' .and. secs(j,5) .eq. 'H' .and. 1 secs(j+1,5) .eq. 'H' .and. cpsip(j-1) .eq. 'H' .and. 1 cpsip(j) .eq. 'H' .and. cpsip(j+1) .eq. 'E')cpsip(j+1)='H' end do do k=1,NAA nh(k)=0 ns(k)=0 nc(k)=0 csecsc(k)='C' if (cpsip(k) .eq. 'H')nh(k)=nh(k)+1 if (csabl(k) .eq. 'H')nh(k)=nh(k)+1 if (secs(k,5) .eq. 'H')nh(k)=nh(k)+1 if (cpsip(k) .eq. 'E')ns(k)=ns(k)+1 if (csabl(k) .eq. 'E')ns(k)=ns(k)+1 if (secs(k,5) .eq. 'E')ns(k)=ns(k)+1 if (cpsip(k) .eq. 'C')nc(k)=nc(k)+1 if (csabl(k) .eq. 'C')nc(k)=nc(k)+1 if (nh(k) .gt. ns(k) .and. nh(k) .gt. nc(k))csecsc(k)='H' if (ns(k) .gt. nh(k) .and. ns(k) .gt. nc(k))csecsc(k)='E' end do do j=1,NAA-2 if (csecsc(j-1) .eq. 'C' .and. csecsc(j) .eq. 'H' .and. 1 csecsc(j+1) .eq. 'E' .and. csecsc(j+2) .eq. 'E')csecsc(j)='E' end do do j=2,NAA-3 if (csecsc(j-1) .eq. 'C' .and. csecsc(j) .eq. 'E' .and. 1 csecsc(j+1) .eq. 'C' .and. csecsc(j+2) .eq. 'C')csecsc(j+1)='E' end do do j=2,NAA-3 if (csecsc(j-1) .eq. 'C' .and. csecsc(j) .eq. 'E' .and. 1 csecsc(j+1) .eq. 'H' .and. csecsc(j+2) .eq. 'C')csecsc(j+1)='E' end do do j=2,NAA-3 if (csecsc(j-2) .eq. 'C' .and. csecsc(j-1) .eq. 'C' .and. 1 csecsc(j) .eq. 'H' .and. csecsc(j+1) .eq. 'H' .and. 1 csecsc(j+2) .eq. 'C' .and. csecsc(j+3) .eq. 'C')then csecsc(j)='C' csecsc(j+1)='C' end if end do write(11,1000)label write(6,1000)label write(6,2005)(AACID(k),k=1,NAA) write(11,2005)(AACID(k),k=1,NAA) write(11,1000)label c write(6,2005)(cnatv(k),k=1,NAA) c write(6,2005)(cpsip(k),k=1,NAA) c write(6,2005)(csabl(k),k=1,NAA) c write(6,2005)(secs(k,5),k=1,NAA) write(6,1000)label write(6,2005)(csecsc(k),k=1,NAA) c write(11,2005)(cnatv(k),k=1,NAA) c write(11,2005)(cpsip(k),k=1,NAA) c write(11,2005)(csabl(k),k=1,NAA) c write(11,2005)(secs(k,5),k=1,NAA) write(11,2005)(csecsc(k),k=1,NAA) c mtrue=0 c do i=1,NAA c if (cnatv(i) .eq. secs(i,5))mtrue=mtrue+1 c if (cnatv(i) .eq. cpsip(i))mtrue=mtrue+1 c if (cnatv(i) .eq. csabl(i))mtrue=mtrue+1 c if (cnatv(i) .eq. csecsc(i))mtrue=mtrue+1 c end do c tp= float(mtrue)*100/float(NAA) 4001 continue 4002 continue c print *,' Number of helices=',mhelix c print *,' Number of Sheets=',msheet c print *,' Number of Turns=',mturn c print *,' Total Number of Sequences=',lfile c print *,' Number of Sequences Skipped=',mshort c print *,' Prediction success rate =',tp c write(11,*)' Prediction success rate =',tp c pr=pr+tp 4000 continue 6000 continue print *,' Total Number of Sequences=',nfile c pr=pr/float(nfile) print *,'*********' c print *,' overall Prediction success rate =',pr c print *,'*********' c write(11,*)'*********' write(11,*)' Total Number of Sequences=',nfile c write(11,*)' Overall Prediction success rate =',pr write(11,*)'*********' c 2005 format (1000(a1)) 2003 format(' Protein Secondary Structure Characterization: ',/, 1 ' This Program identifies a query sequence', 2 ' as a Helix or a Sheet or a Turn. ',/, 3 ' The program requires 5 flanking residues on either side', 4 ' of the query sequence',/, 4 ' in the input for better results. ',/, 5 ' The accuracy is 72% (q3) are higher if the query sequence is', 6 ' 4 residues or longer.',/) 2007 format('********************************************************* 1**********************') write(6,2007) close(10) close(11) c stop#1 c--------------------------------------------------------------------------- data HA/144.225, 144.650, 126.442, 109.657, 190.575, 189.759, 1 188.944, 243.093, 207.139, 210.949, 177.151, 111.613, 1 130.653, 133.601, 1 28.240, 29.131, 25.087, 23.138, 17.801, 20.080, 1 24.120, 29.131, 24.132, 25.113, 32.209, 23.127, 1 22.110, 21.094, 1 115.985, 116.524, 156.547, 238.430, 129.842, 176.706, 1 149.748, 98.442, 100.553, 102.461, 99.647, 127.702, 1 135.678, 142.641, 1 126.071, 123.556, 105.368, 106.639, 124.607, 303.212, 1 203.015, 159.718, 196.078, 169.763, 121.791, 100.553, 1 125.628, 138.623, 1 84.720, 82.370, 77.270, 51.307, 73.298, 65.261, 1 91.457, 84.379, 70.387, 89.402, 86.562, 63.348, 1 77.386, 76.343, 1 178.517, 207.935, 180.632, 172.032, 114.136, 107.429, 1 89.447, 75.339, 60.331, 64.289, 241.570, 326.797, 1 183.919, 147.664, 1 51.437, 51.230, 46.161, 52.313, 41.884, 48.192, 1 46.231, 40.180, 49.270, 54.244, 59.386, 53.293, 1 49.246, 46.207, 1 108.925, 104.470, 99.347, 49.295, 89.005, 66.265, 1 102.512, 114.515, 112.619, 80.361, 76.497, 71.392, 1 104.522, 104.470, 1 108.925, 101.456, 93.326, 86.519, 124.607, 116.465, 1 107.537, 128.578, 165.912, 159.718, 132.863, 130.718, 1 133.668, 127.574, 1 171.457, 167.754, 181.635, 107.645, 176.963, 129.518, 1 207.035, 265.193, 238.310, 255.148, 206.341, 139.768, 1 165.829, 154.696, 1 38.325, 36.162, 42.147, 22.132, 33.507, 30.120, 1 43.216, 52.235, 45.248, 47.212, 45.294, 29.160, 1 31.155, 29.131, 1 86.737, 92.415, 112.393, 149.899, 54.450, 78.313, 1 72.361, 61.275, 66.365, 92.415, 123.804, 119.658, 1 95.477, 91.411, 1 108.925, 122.551, 124.435, 179.074, 230.366, 81.325, 1 42.211, 74.334, 50.276, 20.090, 3.019, 157.868, 1 160.804, 170.768, 1 69.591, 65.293, 58.203, 54.325, 76.439, 91.365, 1 107.537, 75.339, 91.503, 92.415, 88.575, 74.409, 1 67.336, 63.284, 1 101.865, 87.393, 82.288, 74.446, 98.429, 92.369, 1 93.467, 114.515, 131.724, 123.556, 111.726, 100.553, 1 95.477, 91.411, 1 125.063, 126.569, 152.533, 221.327, 123.560, 136.546, 1 96.482, 99.447, 111.613, 126.569, 132.863, 120.663, 1 113.567, 119.537, 1 113.968, 114.515, 131.460, 169.014, 85.863, 92.369, 1 99.497, 69.311, 77.425, 94.424, 82.536, 92.508, 1 94.472, 109.492, 1 139.183, 130.587, 119.417, 63.380, 108.900, 89.357, 1 132.663, 110.497, 111.613, 86.388, 82.536, 83.459, 1 120.603, 137.619, 1 28.240, 24.108, 25.087, 20.120, 38.743, 29.116, 1 30.150, 37.167, 28.154, 26.117, 20.130, 17.094, 1 25.125, 25.113, 1 69.591, 71.320, 60.210, 49.295, 67.015, 56.224, 1 72.361, 67.302, 61.337, 79.357, 75.490, 56.309, 1 67.336, 69.311/ data SA/124.623, 130.784, 121.730, 116.582, 118.712, 127.510, 1 138.763, 131.592, 130.718, 120.603, 112.506, 142.713, 1 134.605, 130.653, 1 23.115, 22.132, 24.144, 32.160, 39.235, 39.156, 1 40.221, 36.162, 41.226, 41.206, 40.180, 36.180, 1 26.117, 22.110, 1 159.798, 163.983, 182.092, 46.231, 46.277, 60.240, 1 82.453, 68.307, 46.254, 40.201, 139.628, 175.879, 1 194.876, 189.949, 1 144.723, 128.772, 108.651, 106.532, 84.507, 85.341, 1 87.481, 92.415, 81.447, 100.502, 97.438, 98.492, 1 135.610, 154.773, 1 55.276, 49.295, 57.344, 107.537, 116.700, 118.473, 1 108.597, 102.461, 125.691, 127.638, 103.465, 62.311, 1 43.194, 55.276, 1 233.165, 384.305, 274.647, 128.643, 62.374, 95.381, 1 125.691, 149.673, 111.613, 84.422, 93.420, 236.180, 1 276.243, 240.201, 1 49.246, 44.265, 47.283, 48.241, 38.229, 44.176, 1 41.226, 38.171, 39.215, 43.216, 54.244, 51.256, 1 44.198, 54.271, 1 63.316, 60.362, 66.398, 176.884, 229.376, 212.851, 1 180.995, 184.831, 215.183, 219.095, 156.705, 77.386, 1 46.207, 59.296, 1 145.728, 124.748, 134.808, 132.663, 83.500, 74.297, 1 75.414, 94.424, 74.409, 91.457, 86.388, 95.477, 1 104.470, 121.608, 1 111.557, 85.513, 94.567, 173.869, 234.406, 234.939, 1 213.172, 201.908, 223.227, 217.085, 173.782, 149.748, 1 79.357, 97.487, 1 24.120, 20.120, 25.150, 40.201, 40.241, 42.168, 1 42.232, 38.171, 42.232, 45.226, 34.153, 31.155, 1 19.085, 24.120, 1 122.613, 141.851, 128.772, 50.251, 40.241, 50.200, 1 65.359, 65.293, 41.226, 42.211, 85.384, 137.688, 1 129.583, 120.603, 1 134.673, 136.820, 161.971, 27.135, 45.271, 39.156, 1 67.370, 54.244, 40.221, 22.110, 91.411, 96.482, 1 206.931, 125.628, 1 75.376, 67.404, 69.416, 72.361, 50.301, 51.204, 1 49.270, 60.271, 50.276, 59.296, 57.257, 50.251, 1 58.262, 72.361, 1 96.482, 83.500, 106.639, 109.547, 85.513, 81.325, 1 76.420, 83.375, 78.431, 94.472, 83.375, 81.407, 1 84.379, 87.437, 1 150.753, 120.724, 114.688, 101.507, 85.513, 96.385, 1 111.613, 96.433, 88.486, 88.442, 134.605, 147.738, 1 169.763, 151.758, 1 118.592, 88.531, 107.645, 146.733, 126.760, 118.473, 1 127.702, 122.551, 117.647, 125.628, 134.605, 137.688, 1 116.524, 124.623, 1 93.467, 86.519, 95.573, 251.256, 331.991, 281.124, 1 248.366, 259.166, 305.681, 290.452, 201.908, 116.582, 1 69.311, 90.452, 1 16.080, 13.078, 24.144, 35.175, 38.229, 35.140, 1 35.193, 33.149, 42.232, 39.195, 32.144, 21.105, 1 17.076, 23.115, 1 57.286, 47.283, 54.325, 96.482, 102.615, 112.449, 1 82.453, 87.393, 104.575, 107.537, 87.393, 54.271, 1 44.198, 54.271/ data TA/188.944, 186.934, 156.547, 120.663, 122.551, 124.560, 1 129.778, 148.669, 134.673, 132.927, 117.647, 106.639, 1 162.244, 169.678, 1 30.150, 30.150, 32.112, 24.132, 21.094, 19.085, 1 21.126, 27.122, 23.115, 24.353, 21.116, 24.144, 1 21.428, 24.096, 1 67.336, 70.351, 120.421, 148.818, 176.795, 167.754, 1 135.814, 121.546, 145.728, 148.148, 173.956, 171.026, 1 85.714, 113.453, 1 146.733, 135.678, 116.407, 105.580, 126.569, 134.605, 1 145.875, 138.623, 134.673, 129.883, 110.608, 110.663, 1 152.040, 216.867, 1 97.487, 105.527, 89.312, 57.315, 58.262, 62.280, 1 69.416, 77.348, 65.326, 61.897, 62.342, 68.410, 1 90.816, 85.341, 1 68.341, 69.346, 193.677, 327.802, 251.130, 201.908, 1 195.171, 169.763, 228.140, 342.973, 230.266, 151.911, 1 86.734, 93.373, 1 41.206, 48.241, 59.207, 51.282, 48.216, 48.216, 1 45.271, 50.226, 52.261, 46.676, 45.248, 50.301, 1 35.714, 45.180, 1 170.854, 135.678, 99.347, 66.365, 72.325, 81.366, 1 94.567, 94.424, 72.361, 74.074, 83.459, 96.579, 1 162.244, 141.566, 1 124.623, 133.668, 120.421, 126.696, 122.551, 128.578, 1 112.676, 135.610, 131.658, 117.706, 119.658, 109.657, 1 103.061, 97.389, 1 256.281, 250.251, 182.639, 131.724, 122.551, 124.560, 1 136.820, 169.763, 135.678, 108.574, 133.735, 123.742, 1 209.183, 179.718, 1 46.231, 48.241, 39.136, 30.165, 26.117, 26.117, 1 29.175, 35.158, 29.145, 24.353, 35.193, 29.175, 1 7.142, 37.148, 1 51.256, 70.351, 109.382, 129.713, 120.542, 112.506, 1 99.597, 96.433, 110.552, 120.750, 119.658, 115.694, 1 41.836, 57.228, 1 23.115, 13.065, 33.115, 130.718, 170.768, 158.714, 1 145.875, 109.492, 124.623, 142.059, 161.890, 107.645, 1 147.959, 52.208, 1 74.371, 78.391, 80.280, 65.359, 62.280, 67.302, 1 65.392, 75.339, 73.366, 64.941, 62.342, 61.368, 1 61.224, 71.285, 1 109.547, 115.577, 105.368, 95.525, 92.415, 91.411, 1 89.537, 106.479, 101.507, 88.280, 93.514, 89.537, 1 89.795, 87.349, 1 97.487, 111.557, 135.474, 126.696, 139.628, 138.623, 1 131.790, 133.601, 136.683, 124.809, 129.713, 178.068, 1 95.918, 115.461, 1 88.442, 107.537, 107.375, 105.580, 107.483, 123.556, 1 120.724, 110.497, 117.587, 99.441, 116.641, 172.032, 1 107.142, 106.425, 1 204.020, 171.859, 122.428, 89.492, 87.393, 109.492, 1 144.869, 116.524, 100.502, 100.456, 107.591, 141.851, 1 225.510, 193.775, 1 32.160, 29.145, 22.077, 16.088, 17.076, 22.099, 1 25.150, 21.094, 19.095, 15.220, 22.121, 25.150, 1 36.734, 33.132, 1 81.407,88.442, 75.263, 50.276, 54.244, 57.257, 1 61.368, 62.280, 63.316, 32.470, 53.293, 66.398, 1 77.551, 79.317/ data HA1/1.04,1.28,1.28,1.28,1.39,1.32,1.32,1.28, 2 0.36,0.26,0.30,0.38,0.46,0.38,0.40,0.51, 3 1.38,1.11,1.25,1.18,0.99,1.00,1.01,1.00, 4 1.03,1.10,1.48,1.31,1.20,1.29,1.23,1.09, 5 0.71,0.86,0.81,0.96,0.92,0.85,0.95,0.94, 6 1.24,1.06,1.03,0.95,0.88,0.78,0.81,1.38, 7 0.72,0.62,0.68,0.66,0.60,0.69,0.73,0.77, 8 0.69,0.95,0.82,1.01,1.06,1.05,0.90,0.88, 9 0.93,1.10,1.07,1.03,1.11,1.22,1.20,1.12, 1 1.03,1.25,1.11,1.32,1.42,1.38,1.41,1.31, 1 0.34,0.53,0.48,0.63,0.72,0.65,0.67,0.65, 2 1.18,0.73,0.89,0.86,0.79,0.82,0.96,1.09, 3 1.26,1.36,0.91,0.62,0.87,0.70,0.30,0.00, 4 0.73,0.88,0.96,1.03,0.88,0.96,0.96,0.95, 5 0.87,0.99,0.96,0.97,1.06,1.12,1.09,1.05, 6 1.35,1.09,1.14,0.98,1.00,1.05,1.10,1.12, 7 1.23,0.93,0.96,1.00,0.84,0.89,0.98,0.91, 8 0.80,1.04,0.95,1.12,1.05,1.05,0.93,0.92, 9 0.30,0.59,0.46,0.48,0.57,0.45,0.41,0.30, 1 0.69,0.83,0.75,0.86,0.83,0.79,0.90,0.88/ data SA1/1.07,1.08,1.11,1.14,1.12,1.12,1.08,1.05, 2 0.51,0.59,0.59,0.60,0.56,0.61,0.61,0.60, 3 0.66,0.66,0.78,0.91,0.83,0.66,0.60,1.15, 4 1.03,0.92,0.93,0.94,0.96,0.91,1.00,0.99, 5 1.03,1.07,1.08,1.04,1.01,1.10,1.11,1.02, 6 1.11,0.79,0.98,1.10,1.18,1.05,0.92,0.97, 7 0.68,0.58,0.64,0.61,0.58,0.59,0.63,0.73, 8 1.25,1.36,1.33,1.26,1.27,1.33,1.34,1.20, 9 1.12,0.92,0.87,0.88,0.97,0.87,0.96,0.93, 1 1.24,1.37,1.37,1.33,1.31,1.35,1.34,1.24, 1 0.60,0.60,0.62,0.62,0.58,0.62,0.65,0.53, 2 0.70,0.60,0.70,0.81,0.81,0.61,0.62,0.93, 3 0.43,0.65,0.59,0.83,0.73,0.60,0.34,0.96, 4 0.86,0.70,0.71,0.70,0.78,0.70,0.77,0.76, 5 1.04,0.93,0.91,0.88,0.92,0.89,0.97,0.92, 6 1.00,0.93,0.98,1.05,0.98,0.94,0.94,1.13, 7 1.17,1.10,1.08,1.11,1.09,1.07,1.10,1.13, 8 1.40,1.52,1.45,1.40,1.41,1.49,1.46,1.31, 9 0.54,0.58,0.54,0.54,0.52,0.62,0.59,0.51, 1 0.98,1.01,1.05,0.91,0.94,1.02,1.03,0.94/ data TA1/1.08,1.09,1.10,1.11,1.17,1.13,1.12,1.07, 2 0.38,0.32,0.28,0.32,0.43,0.36,0.38,0.32, 3 1.17,1.25,1.23,1.13,1.09,1.16,1.17,1.24, 4 1.03,1.10,1.13,1.16,1.14,1.13,1.11,1.05, 5 0.76,0.76,0.79,0.84,0.89,0.81,0.79,0.79, 6 1.51,1.40,1.31,1.29,1.23,1.36,1.54,1.36, 7 0.71,0.68,0.68,0.65,0.70,0.72,0.67,0.65, 8 0.82,0.86,0.91,0.98,0.97,0.86,0.87,0.92, 9 1.10,1.09,1.11,1.05,1.13,1.12,1.07,1.08, 1 1.12,1.09,1.10,1.14,1.23,1.13,1.04,1.13, 1 0.48,0.41,0.41,0.46,0.54,0.46,0.38,0.54, 2 1.11,1.08,1.05,1.00,0.98,1.05,1.08,1.08, 3 1.12,1.23,1.20,1.16,1.04,1.10,1.15,1.21, 4 0.81,0.79,0.83,0.81,0.88,0.86,0.81,0.79, 5 0.98,0.96,0.96,0.95,1.03,1.00,0.94,0.97, 6 1.10,1.15,1.14,1.12,1.13,1.14,1.10,1.11, 7 1.03,1.03,1.09,1.08,1.05,1.07,1.00,1.07, 8 0.95,0.94,1.04,1.16,1.07,1.00,1.00,1.03, 9 0.20,0.23,0.34,0.40,0.32,0.28,0.18,0.34, 1 0.70,0.73,0.76,0.79,0.79,0.80,0.52,0.72/ c--------------------------------------------------------------------------- stop end c---------------------------------------------------------------------------