2010年10月29日 星期五

t_tide(Matlab)

以往都將站名與年份先打在commond window裡,
經過昨天與今天的小小努力,
已經改成可以將站名、年份與緯度分開在兩個檔案裡直接用load就可以了~
方便多了~


-------------------(code)--------------
lat=load('./data/nlsc2010/lat.txt'); %相對應測站所在的緯度
fsite=fopen('./data/nlsc2010/fileyr2.dat','r'); %打開放站名和年份的檔案,並命名為fsite。r是指開啟檔案供讀取,此檔案須先存在磁碟機內,若不存在則開檔失數。開檔成功,傳回大於0的整數,開檔失數,傳回-1。

for nof=1:10 %10個站,或是1個站10年(一年跑一次),or others
    fn=fgetl(fsite); %fgetl:讀取一行字串,但不讀換行字元,若改成可讀換行字元(指令:fgets),就會一次讀兩行的樣子…

    da1=load(strcat('./data/nlsc2010/',fn,'.txt'));  %load data...
    .
    .
    .
end

fclose(fsite); %打開檔案,就要關掉檔案,才是好孩子哦! 不過我不是孩子,所以… 嘿嘿嘿~

matlab logo

沒有留言:

張貼留言