성적표 주고 나서 그 자료들을 성적순으로 출력해오라는 과제를 받았는데

어케 하는건가요? 찾아보니까 밑에 내용이 있긴한데 실행하니까 아무것도 안나오네요


      program credit
     
      implicit none
     
      character*15 s_name(15)
      integer num(15),i,j,k
      real math(15),korean(15),english(15)
      real s_avg(15),average(3),variance(3),standard_dev(3)
      real temp

c Read \'data.inp\' file     
      open (10,file="data.inp")
      read(10,*)
      do 20 i=1,15
        read(10,100) num(i),s_name(i),korean(i),english(i),math(i)
 20   continue
      close(10)

c Calculate average,variance,standard devation of each subject 
      call calc_dev(korean,average(1),variance(1),standard_dev(1))
      call calc_dev(english,average(2),variance(2),standard_dev(2))
      call calc_dev(math,average(3),variance(3),standard_dev(3))

c Calculate average of each student
      do 30 i=1,15
       call calc_avg(korean(i),english(i),math(i),s_avg(i))
 30   continue
     
      do 40 j=14,1,-1
      temp=s_avg(j+1)
      do 40 i=1,j
        if (s_avg(i)<=temp) then
        s_avg(j+1)=s_avg(i)
        s_avg(i)=temp
        temp=s_avg(j+1)
        k=num(j+1)
        num(j+1)=num(i)
        num(i)=k
        end if
 40   continue

c Output to file \'cerdit.out\'
      open(11,file="credit.out")
      write(11,*) "-- Korean --"
      write(11,110) average(1),variance(1),standard_dev(1)
      write(11,*) "-- English --"
      write(11,110) average(2),variance(2),standard_dev(2)
      write(11,*) "-- Math --"
      write(11,110) average(3),variance(3),standard_dev(3)
      write(11,*) "------------ Student Data ------------"
      write(11,*) ".Num....Name....Avg..Kor..Eng..Math.."
      do 90 i=1,15
       write(11,120) num(i),s_name(num(i)),s_avg(i),korean(num(i)),
     *  english(num(i)),math(num(i))
 90   end do
      close(11)
     
 100  format(4X,I2,4X,A10,4X,F4.1,7X,F4.1,4X,F4.1)
 110  format("Average : ",F7.2,"  Variance : ",F7.2,
     *"  Standard devation : ",F7.2)
 120  format(1X,I2,2X,A10,F4.1,1X,F4.1,1X,F4.1,1X,F4.1)
      end program credit

      subroutine calc_avg(k,e,m,avg)
      real k,e,m,avg
      avg=(k+e+m)/3.0
      end

      subroutine calc_dev(subj,a,v,sd)
      real subj(15),a,v,sd
      real s_sum,diff
      integer l

      s_sum=0
      do 60 l=1,15
       s_sum=s_sum+subj(l)
 60   continue
      a=s_sum/15.0

      diff=0
      do 70 l=1,15
       diff=diff+(subj(l)-a)**2
 70   continue
      v=diff/15.0

      sd=sqrt(v)   
     
 pause
 end

어디가 잘못된거죠?