program rdfneb_csp character*3 mcod(7) data mcod /'m23','m18','m13','m07','m03','p00','p02'/ character*10 gty(13) data gty /'Burst ', 'p1GBurst ', 'p3GBurst ', 'p5GBurst ', : 'p7GBurst ', 'p9GBurst ', 'E ', 'S0 ', 'Sa ', : 'Sb ', 'Sc ', 'Sd ', 'Irr ' / real tb(90),w(1221), fs(1221,90), fneb(1221,90), ftot(1221,90) character cwhat*1,fn21*80 c _________________________________________________________ write(*,*) 'ikind=?, 1: bb; 2:ss' read(*,*) ikind if(ikind.eq.1) cwhat='b' if(ikind.eq.2) cwhat='s' do nz=1,7 do nty=1,13 ii = index(gty(nty), ' ') fn21 = mcod(nz)//cwhat//cwhat//'_csp_'//gty(nty)(1:ii-1)//'. :fneb' open (2,file=fn21, form='formatted') read(2,5,end=26) nn, (tb(i),i=1,nn) read(2,6) inl,(w(kkk),kkk=1,inl) do i=1,nn read(2,7) (fs(kkk,i),kkk=1,inl) read(2,7) (fneb(kkk,i),kkk=1,inl) read(2,7) (ftot(kkk,i),kkk=1,inl) enddo close(2) enddo enddo 26 continue 5 format(i8,90e11.4) 6 format(i8,1221f10.1) 7 format(7e11.4) stop end