program rdLha_fneb_sdb 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), nphot(90), lha(90), lhb(90), lbol(90) character fn21*80 c _________________________________________________________ do nz=1,7 do nty=1,13 ii = index(gty(nty), ' ') fn21 = mcod(nz)//'bb_csp_'//gty(nty)(1:ii-1)//'-sdb.quan' open (2,file=fn21, form='formatted') read(2,*) do i=1,90 read(2,59) tb(i), nphot(i), lha(i), lhb(i), lbol(i) enddo close(2) enddo enddo 59 format(e12.5, 4f10.5) stop end