제목 : 오염확산론 (Wind Rose analysis)
─────────────────────────────────────────
□ 개요
? 원본 데이터 추출과정
? 풍속, 풍향별 데이터 counting (Fortran90)
? 데이터 수정 및 저장 (Work Sheet)
? 데이터 베이스 입력 후 Wind Rose 출력 (Origin)
? 참고자료
□ 원본 데이터 추출과정
- 기상청 홈페이지(www.kma.go.kr)로 접속 후, 한 달간의 풍향과 풍속분석에 관한
AWS 데이터를 받음. (1)
대상지역 : 전남 순천
기 간 : 2007/08/01 00:00 ~ 08/31 23:00
자료형태 : 1 × 1 형식, 16방위
1분 평균 자료
일일 60분 단위 자료
- Fortran90 프로그래밍을 위해 Work Sheet(MS Exel2007)를 이용하여, 원본 데이터를 ‘공백
으로 분리’한 형태로 저장 후, 텍스트 파일로 전환 후 다시 저장. (2)
□ 풍속, 풍향별 데이터 counting (Fortran90)
- PC로부터 Text 파일을 서버로 다운로드 (1)
- Fortran90을 이용하여 AWS 데이터 중 필요한 풍향과 풍속의 데이터만 추출 (2)
!------------------------------------------------------------------------------------
program wrose_1
!------------------------------------------------------------------------------------
implicit none
…(생략)
|
---------------------------
program wrose_2
!---------------------------------------Declar a Variable-------------------------------------
implicit none
integer :: num
real,dimension(744) :: wind_s
character(3),dimension(744) :: wind_d
!----------------------------------------Executive(1)-----------------------------------------
do num=1,744
open(10,file=`wind_ds.txt`,status=`old`,form=`formatted`)
read(10,`a3,2x,f4.1`)wind_d(num),wind_s(num)
open(20,file=`wind_fin.txt`,status=`unknown`,form=`formatted`)
if(wind_d(num)=="N")then
print`(a,2x,f4.1)`,wind_d(num),wind_s(num)
write(20,`a,2x,f4.1`)wind_d(num),wind_s(num)
else if(wind_d(num)=="NNE")then
print`(a,2x,f4.1)`,wind_d(num),wind_s(num)
write(20,`a,2x,f4.1`)wind_d(num),wind_s(num)
else if(wind_d(num)=="NE")then
print`(a,2x,f4.1)`,wind_d(num),wind_s(num)