- This topic has 4개 답변, 2명 참여, and was last updated 9 years, 3 months 전에 by 익명.
-
글쓴이글
-
-
익명비활성
아래와 같이 짜려고 생각해보았지만 잘 되질 않습니다.. 8개의 지역은 각각 txt파일로 월별로 쪼개져 있습니다.
7개지역 12개월 해서 84개이고요 각 지역의 데이터를 data1~data7으로 지정하여 일평균,월평균 1년 전체 변화도 등을 그려보고 싶습니다. data1~data7 로 지정하는 과정에서 문제가 생기는데 도움 부탁드립니다.
매번 도움주시는 모든분들 감사드립니다Site_A=[‘BG’,’GC’,’GH’,’GW’,’KI’,’SM’,’TB’]
;—————————————————————————————
data=[[data1],[data2],[data3],[data4],[data5],[data6],[data7]]for k = 1, 1 do begin
spc=s[k-1]
for month = 1, 1 do begin
for j = 1, 7 do begin
site = Site_A[j-1]
;—————————————————————————————
if (month le 9) then month_s=’0’+strtrim(month,2) else month_s=strtrim(month,2)
input_file=dir_in+month_s+’_’+site+’.txt’
FMT = ‘L,I’readcol, input_file, F=FMT, Date, Time, v1,v2,v3,v4,v5,v6,v7,v8,v9,v10,$
v11,v12,v13,v14,v15,V16,v17,v18,v19,v20,v21,v22,v23,v24,v25,v26,v27,$
v28,v29,v30,v31,v32,v33,v34,v35,v36,v37,v38,v39,v40,v41,v42,v43,v44,$
v45,v46,v47,v48,v49,v50,v51,v52,v53,v54,v55,/silentV=[[v1],[v2],[v3],[v4],[v5],[v6],[v7],[v8],[v9],[v10],[v11],[v12], $
[v13],[v14],[v15],[v16],[v17],[v18],[v19],[v20],[v21],[v22],[v23], $
[v24],[v25],[v26],[v27],[v28],[v29],[v30],[v31],[v32],[v33],[v34], $
[v35],[v36],[v37],[v38],[v39],[v40],[v41],[v42],[v43],[v44],[v45], $
[v46],[v47],[v48],[v49],[v50],[v51],[v52],[v53],[v54],[v55]]
;——————————————————————————————
whe=where(s eq spc)data[*,j]=v[*,whe]
print, data[*,j-1]
;——————————————————————————————
endfor
endfor
endfor
end -
Jonghyuk회원
꽤 복잡한 코드입니다.
일단, 아래 부분을 IDL이 인식하는 상황인지요?data=[[data1],[data2],[data3],[data4],[data5],[data6],[data7]]
아래 부분은 IDL이 잘 처리할 것 같긴 한데,
if (month le 9) then month_s=’0’+strtrim(month,2) else month_s=strtrim(month,2)
다음과 같이 처리하는 것이 더 일반적입니다.
month_s=string(month, format='(I02)’) -
익명비활성
data=[[data1],[data2],[data3],[data4],[data5],[data6],[data7]]
이부분의 경우 data1~7이 정의되지 않았다면서 인식을 못하더군요…
그문제만 해결하면 진행이 될것 같은데 해결이 안되서 막혔습니다..아래 지적해주신 부분은 더 간단한 방법이 있었는데 몰랐네요; 수정했습니다
-
익명비활성
나중에도 편의를 위해
data1~7 까지 각각 지정해주고 싶은데 지정하는게 어렵네요… -
익명비활성
입력 데이터의 포맷을 바꾸어서 해결하였습니다
감사합니다
-
-
글쓴이글