여러 .dat 파일에 대해 동일한 계산을 수행하는 방법

Question

여러 .dat 파일에 대해 고정 된 계산 단계를 수행하려고합니다.

dimension t(128716),x(128716) 

    open (unit=88,file='ALFA-gua-100m-2.dat',status='unknown') 

    do i=1,128716 
     read(88,*)t(i),x(i) 

    enddo 
    sum=0 
    do j=1,128716 
     sum=sum+x(j) 
    enddo 



    write(*,*)sum/128716 
    close(88) 
    stop 
    end 

가 어떻게이 가야합니까 : 여기

은 내가 한 .dat 파일 즉 계산과 함께하고 싶은 내 코드? 제발 제안!

dimension t(128716),x(128716) 

    open (unit=11,file='ALFA-gua-100m-2.dat',status='unknown') 
    open (unit=12,file='ALFA-gua-100m-5.dat',status='unknown') 

    do i=1,2 
     ii = i + 10 

     do j=1,128716 
     read(ii,*)t(j),x(j) 
     enddo 

     sum=0 
     do k=1,128716 
     sum=sum+x(k) 
     enddo 

    enddo 

    do l=1,2 
     ll = l + 10 
     write(ll,*)sum/128716.0 

     close(ll) 
    enddo 

    stop 
    end 

하지만 그 작동하지 :

여기에 여러 파일에 대한 내 코드입니다.

Answer 1

루프가 두 개인 경우 sum 루프를 i 루프의 두 반복에 사용할 수 없습니다. 두 번째 파일을 처리 할 때 덮어 씁니다.

1. 루프를 하나로 결합 할 수 있습니다.

   do i=1,2 
       ii = i + 10 
   
       do j=1,128716 
       read(ii,*)t(j),x(j) 
       enddo 
   
       sum=0 
       do k=1,128716 
       sum=sum+x(k) 
       enddo 
   
       write(ii,*)sum/128716.0 
   
       close(ii) 
   enddo 
   

2. sums(i)의 배열을 사용할 수 있습니다.

   do i=1,2 
       ii = i + 10 
   
       do j=1,128716 
       read(ii,*)t(j),x(j) 
       enddo 
   
       sum=0 
       do k=1,128716 
       sums(i)=sums(i)+x(k) 
       enddo 
   
   enddo 
   
   do l=1,2 
       ll = l + 10 
       write(ll,*)sums(l)/128716.0 
   
       close(ll) 
   enddo 
   

Answer 2

@ VladimirF의 대답에 대한 추가.

우리가 단순히 모든 루프를 작성하는 프로그래머에 대한 필요가 없습니다

sumx = sum(x) 

쓸 수 x라는 배열의 모든 요소를 ​​요약합니다. 합계 배열을 사용하는 경우

sums(1) = sum(x) 

과 같은 것이 적합합니다.

그럼 내가 쓰고있는 동안 나는

meanx = sum(x)/size(x) 

를 써서 배열의 평균을 계산합니다 :이 변수 sum를 호출하는 것은 좋은 생각이 아니다. 그 이름의 기존 내장 함수가 있는데 독자들 (컴파일러가 아니더라도)은 그 이름의 변수도 혼동 할뿐입니다.