포트란 : 파일로 부터 읽어들인 값들을 벡터로 저장하는 방법 질문요!

포트란으로 데이터 파일을 읽어들여서 계산을 하고는 다시 저장 시키는데요.

'파일 - 읽어들이기 및 계산하기 - 새로운 파일 생성/저장' 은 할 수 있는데

계산된 새로운 데이터를 행렬로 저장했으면 좋겠어서요.

현재 do loop를 돌려서 이렇게 하고 있습니다.

  1 program HW
  2   implicit none
  3
  4   real :: wavleng
  5   real :: anfq
  6   real, parameter :: c = 2.998E+8
  7   integer :: s
  8
  9   open(unit=1, file='HWdata.dat')
10
11   do
12    read(1,*,IOSTAT=s) wavleng
13    if(s /= 0)  exit
14
15    anfq = 2*4*atan(1.0)*(c/wavleng)
16
17
18
19    open(unit=2, file='angular_fq')
20    write(2,*) anfq
21   end do
22
23
24 end program HW
25

위에서 보시는 봐와 같이 wavleng 이라는게 do loop 한번 돌때마다 바뀌는데요.

이걸 행렬로 저장해서 계산을 loop 없이 한꺼번에 할 방법 없을까요?

감사합니다.

잡담 | 1501명이 읽었어요. 52.14.84.187

비공개 손님 … 2016-01-20 14:28:44

integer length
real, allocatable :: wavleng(:)
open(unit=1, file='HWdata.dat')
read(1,*) length
allocate(wavleng(length))
do i=1, length
read(1,*) wavleng(i)
end do
deallocate(wavleng)

간단하게 이런식으로 하면 wavleng을 행렬로 읽을수 있는데 질문을 정확하게 이해한지 모르겠네요.
이건 편의상 인풋 파일 맨앞에 wavleng행렬의 사이즈를 넣어줄때구요.

비공개 손님 … 2016-01-20 15:30:49

답변 감사합니다.
3번째 줄에서 HWdata.dat을 열고 4번째 줄에서 length 라고 지정된 변수를 읽어들이는데
이렇게 되면 한줄만 읽는거 아닌가요? do loop 이용해서 1열에 있는 데이터를 행렬로 읽어들이려 하거든요..

비공개 손님 … 2016-01-20 15:46:56

예를 들어 1x1000 인 실수의 데이터를 행렬로서 읽어드리려는 건데 do loop을 사용해서 데이터가 1000 번째 행까지 읽게 해야 하지 않나요?

비공개 손님 … 2016-01-20 15:50:20

1000x1 이네요...

비공개 손님 … 2016-01-20 16:03:44

음... 지금 인풋파일이 3x1행렬형태로 되어 있고 예를들어
100
200
300
이렇게 저장되어 있는거면 wavleng(1)=100, wavleng(2)=200, wavleng(3)=300 이렇게 읽으면서 저장하면 되는게 맞나요? 그럼 답변단게 얼추 통하긴 할텐데요

비공개 손님 … 2016-01-20 16:20:24

네. 맞아요. 근데 또 궁금한게 length를 integer로 두셨는데 만약 읽어 들이는 데이터가 real 이면 allocate 할때는 어떻게 해야하나요?

비공개 손님 … 2016-01-20 16:29:43

read(1,*) length 하면
100 -> 이거 한줄만 읽고 땡이어서 length 값은 100이 되고...

wavleng(length) 의 크기는 100x1 이 되는거 아닌가요?

그럼 행렬 데이터의 사이즈에 따라 wavleng의 사이즈가 결정되는게 아니라 데이터 값에 따라 wavleng의 사이즈가 결정되는거 같은데요...

질문만 너무해서 죄송합니다..

비공개 손님 … 2016-01-20 16:30:32

1 답변에서도 사실 length만 integer이고 데이터들은 real로 할당되어 있어요. 1답변 보면 첫번째 줄은 integer값을 읽고 그 다음줄부터는 do loop로 real값을 읽고 있어요. 그래서 인풋파일 맨 윗줄에 데이터 벡터 사이즈만 추가로 넣어주면 되요 예를들어
4
5.14112
6.1512
7.121231
8.1192
이런식으로요. 인풋파일 수정하는게 싫으면 그냥 read(1,*) length 를 빼고 코드에서 행렬 사이즈를 정해주어도 괜찮구요

비공개 손님 … 2016-01-20 16:32:05

8// 아 그럼 인풋 데이터가 얼마나 사이즈가 되는지는 안정하고 하고싶다는 말이죠?

비공개 손님 … 2016-01-20 16:42:06

10// 네 맞습니다. 데이터가 예를들어 셀수도 없이 많거나 세기 귀찮거나 하면 그 갯수를 do loop 으로 돌려서 파악하면서 데이터를 행렬로 읽어들이고 싶어서요.

비공개 손님 … 2016-01-20 16:45:51

이걸 하려고 read 의 IOSTAT 이라는 옵션을 사용하고 변수 i 를 선언해서 카운터를 만들어서 파악한 i 값으로 행렬 사이즈를 allocate 했었거든요.

근데 그 후에 어떻게 해야할지를 모르겠네요...

비공개 손님 … 2016-01-20 16:47:43

데이터의 행렬 사이즈는 파악이 되는데 행렬로서 데이터를 읽어들이지를 못해서 데이터를 가지고 계산을 못하고 있답니다;;;;;

비공개 손님 … 2016-01-20 16:48:09

integer length , s, i
real temp
real, allocatable :: wavleng(:)

open(unit=1, file='input.txt')

length = 0
do
read(1,*,IOSTAT=s) temp
if(s /= 0) exit
length = length + 1
end do

close(unit=1,status='keep')
open(unit=1, file='input.txt')

allocate(wavleng(length))

do i=1, length
read(1,*) wavleng(i)
end do

deallocate(wavleng)
close(unit=1,status='keep')
end

input파일을 두번 읽기는 하지만 이게 젤 나을거 같네요

비공개 손님 … 2016-01-20 16:50:29

처음 input을 읽을때는 전체 데이터의 사이즈만 파악하고
파악된 데이터 사이즈로 allocate하고 나서
데이터를 행렬에 읽으면서 저장하는 방식이에요

비공개 손님 … 2016-01-20 16:51:59

오. 답변 감사합니다. 근데.. deallocate를 하는 이유는 왜 그런가요?
wavleng을 써먹으려면 가만 둬야하는거 아닌가요?

비공개 손님 … 2016-01-20 16:55:28

아 deallocate는 그냥 프로그램 끝날때 해주시면 되고 그 전에 맘대로 써먹으시면 됩니다.ㅋ

비공개 손님 … 2016-01-20 16:57:07

아 ㅎㅎㅎ. 메모리 덤프 시키는 역할인가요?
여튼 너무나도 감사합니다 ㅎㅎㅎ.

댓글을 작성하실 수 없습니다.
(권한이 없는 회원레벨)

목록으로

△ 이전글: 안드로이드 도와주세요 ㅠㅠ [3]

▽ 다음글: pyinstaller 아시는분... [2]

TOP arrow_upward