廖董不懂的部落格

分享一下去年底寫的一隻程式，雖然是去年寫好的，但是因為工作所需今年才上線，就在這邊分享一下了
源起：會做這工作主要是因為工作所需，這邊的博士後把去年的TEDS分析好了，接下來就是要更新我們那用了十年的TEDS排放資料了，但是基本上程式資料都是人家給的，根本沒有自己做的，而且裡面的物種也是不完整，這位博士後很認真的幫我們分析資料，並加入了亞洲的排放，點線面都完成了，把汙染物從26種提升到50種，所以排放資料有了，接下來就是預報了。所以我的工作就是把他處理完的資料(TEDS)跟自然排放源(MEGAEN)加在一起，寫成一組新的netcdf檔案，提供給模式(CMAQ、WRFCHEM)使用，廢話不多說先來開始說程式吧。
其實一開始我並不想用NCL來寫的，我一開始用fortran去寫nc檔案，但是寫了一兩天就發生問題了，寫入nc矩陣很簡單，也一下就寫好了，但是寫矩陣說明跟global_attribute的時候格式都不能被模式接受，有很多小地方的錯誤，後來我就放棄了，使用NCL是因為我之前做氣象模擬的時候，有利用NCL去改變初始資料，測試溫度的敏感度，所以想說有了基礎，那就用這個好了，所以就去查了一下官網的說明https://www.ncl.ucar.edu/Applications/netcdf4.shtml，就從基礎開始了。

這一年實驗室變化很大(因為之前很久沒有變化了)，老師的第一位直系博士後高升到中研院，又畢業了兩位博士生，所以很多事情都不一樣了，最近來的博士後一位是大陸人一位是韓國人，基本上差一位日本人就可以完成日本帝國的大東亞共榮圈了(哈哈哈，這是開玩笑的題外話啦)，最近韓國的金小姐跟我說他需要使用這套軟體GMT，但是我對這套軟體還真的是一問三不知，所以我大概找了一個下午查詢如何安裝。
先來介紹何謂GMT(Generic mapping tools)，它是一個開放性軟體，基本為卡氏座標，可以提供地型特徵，可以畫x-y 圖形到3D圖形以PostScript圖型輸出，共有81個指令可以提供使用。要如何使用我就不清楚了，在這邊做的簡單的安裝記錄，因為本實驗室的cluster(叢集電腦)為suse server 10.1版本核心版本為4.1.2算是6~7年前的產品了，所以對照OS查詢的結果應該是沒有辦法安裝這套軟體，所以我後來用去年購買的伺服器(OS: centos 6)來安裝這套軟體。
首先要先介紹一下安裝時需要的套件，一般我們安裝好linux的時候，基本上它會給你他官方的yum server，但是還有另外一個EPEL(Extra Package for Enterprise Linux)維護的套件軟體庫，要在自己的/etc/yum.repos.d下面把這資料庫加進去，或者是直接用rpm到網路安裝EPEL套件
#>rpm -Uvh http://mirror01.idc.hinet.net/EPEL/6/x86_64/epel-release-6-8.noarch.rpm

這次事件是linux使用者成長的路上一定會失誤或者是不小心犯的錯誤也就是"rm -rf *"，話說我們這邊有一位年輕小姐因為放假前太開心不小心下了這條指令，把自己辛辛苦苦寫了兩三個禮拜的程式刪掉了，於是急急忙忙的求救，我在幾年前好像也有遇過一樣的問題，不過那時候是用dd這個指令去把block的空間寫到另一磁區內，不過年代久遠也忘記該怎樣操作，所以就google了一下救回硬碟資料的方法，有很多方法，不過主要的方式就是先找到被刪除檔案的block並利用block的位置來把資料寫回來，我大概試了三個方法提供大家做參考，在執行這些方法的時候有幾個重要的要提醒一下1. 不能重開機 2. 把掛載的硬碟先掛載成只可以讀取(mount -o remount, ro $mountdirectory)。
1. extundelete 這個程式使用方法很簡單(不過要先安裝其他lib，如果有缺啥lib用yum裝起來吧),
#>./configure --prefix=$path

要如何更改linux時區呢??在這邊寫出來提供大家做一個參考。
為什麼我需要更改linux時區呢，這是因為我們家系統去年做了更新，但是系統的時間一直都很奇怪，我直接用date去更改時間，每次改完都還是會跳回來，最近因為要用時間去判斷迴圈，所以這個問題一定要解決，一開始我並不知道為什麼設定好了(local time)時間會跳回來，後來又再去查了UTC時間才發現問題在這裡，所以卡了一下子。
現在就跟大家介紹一下!!
我遇到的問題其實就是這樣

 在大部分的linux主機裡面對於顯示卡的支援好像都不是太完善(不是說電動方面的)，我的伺服器一開始本來可以使用，但是後來不知道更新什麼東西導致google earth、vapor這些需要吃顯示卡的軟體都不能使用，但是內建的gun圖片處理器還有koliurPaint等繪圖工具反而是正常，於是乎這個軟體就休息了好一陣子，最近看到了一篇文章是解決
#>glxinfo   (利用這個指令來查看顯示卡是否在系統上有作用) 
xlib : extension "GLX" missing on display ":0.0"
的解決方法，所以照本宣科的做了一遍，目前倒是很順利的可以執行(https://duxyng.wordpress.com/2012/01/26/finally-working-nvidia-optimus-on-fedora/)。不過我有做了一些解釋跟修改應該不算盜用吧！

這篇文章主要是給自己做一個紀錄，如果有遇到相同問題的人也可以參考一下，之前九月份有跟廠商購買infiniband軟體跟一台光纖diskarray，因為我們一般使用者不可能會用得這麼高級的機器，所以我有很多東西不敢碰，想說教她廠商就好了 ~~ 廠商也是很爽快跟我說:"我們經驗很多，沒有問題"，於是硬體規劃跟配置我就交給他們了，但是這才是我這幾個月噩夢的開始，這家廠商我從在中央就跟他們合作了，我跟他們第一任的關係算是滿愉快的，我也從中學到了不少，但是後來的工程師素質是每況愈下，最近更是差到極點。
先來說說他們的失誤好了，首先是硬體方面，當初跟我說安裝infinband原有的兩個diskarray可以無縫接軌，結果後來計算插槽錯誤，每一台pcie只有兩個!! 所以插光纖+infiniband就花了兩個，所以兩台login server有四片插槽，總共有兩片disk array轉接卡，一片光纖網卡，兩張infiniband網卡，五張卡要插四個槽。(因為我一開始以為是用infiniband來做io，所以就相信他了，這是我的錯誤，不過我當初有問過兩三次這些問題，業務都跟我再三保證，他絕對不會損失我們的硬體空間。)後來我想說我們自己把空間放到別邊使用好。
後來就是軟體安裝了，叫他們裝一下infinfband的驅動跟library這件事就可以花我兩個月，這也是我最受不了的地方，期間還給我刪掉我的/opt檔，更讓我痛不欲生。

這篇文章主要是給自己做個紀錄，如果有問題的朋友可以來信，我會仔細地回答我可以回答的問題。
因為最近有幾位朋友跟我說他們windws密碼被改了，所以沒有辦法登入，不過沒有人去動你的密碼，密碼怎麼可能變更呢？這是我心中的疑問，但是事實就是發生了，所以我就上網找了一些方法，主要的方法有三種。
第一種：
    這是最簡單的方法，首先在開機的時候"F8"-->" 選擇安全模式"-->利用administrator登入-->控制台-->使用者帳戶，去更改密碼就可以了，但是如果你不知道administrator密碼，這就不行了。

要讓PUTTY顯示中文，除了自己要安裝中文的編譯器之外，就是要有良好設定了 ~~~
參考網路文章
http://blog.yam.com/cola1028/article/23040078
Server 端設定
#vi /etc/profile
LC_CTYPE=en_US.ISO8859-1; export LC_CTYPE
LC_CTYPE=zh_TW.UTF-8; export LC_ALL
LANG=zh_TW.UTF-8; export LANG
SSH 連線軟體 Putty 設定：
Change Settings ->> Translation ->> Character set translation on received data ->> 選UTF-8

最近因為我大姊要報稅，所以下載了報稅軟體 !!安裝好之後跑出來看到的是一大堆亂碼，於是做為弟弟的我只好幫他解決 ~~根據我對window的認識我一開始是認為因為是字型不支援，所以本來以為可以去安裝補丁就好，但是好像不是如此，這時候啥都不知道的我只能求助GOOGLE大神，Google大神說：http://tax.nat.gov.tw/a_info_ans_7_isTop.html於是我就照者裡面的步驟去找，國稅局網站說：[控制台]-->[時鐘、語言和區域]-->[地區及語言]可是我找了好久~~都找不到 [時鐘、語言和區域] 這個選項，不過有找到[地區及語言]這個選項，設定好台灣的地區就可以了，記得要重開機。後來大概了解，因為地區的不同會導致unicode轉換(UTF)出現問題，所以才會有亂碼出現，所以只要設定好就可以了，但是這跟regedit應該也可以改設定 ~後來查到可以利用在cmd下面regedit-->HKEY_CURRENT_USER-->Contrl Panel-->International有一個Locale的選項 ~ 按右鍵-->修改成00000404就可以了。

說到Netcdf格式，做大氣跟海洋的人應該都不陌生吧~~應該是跟我一樣有種又愛又恨的Feel，最近本來想說來把的語法好好學一學，因為未來可能要更改某些default設定，所以想說先來學一學netcdf的讀取等等，一開始直覺就是用fortran直接include netcdf的library(http://www.unidata.ucar.edu/software/netcdf/examples/programs/index.html，可以參考這個網站)，基本上是這樣寫啦：
讀檔
include "netcdf.inc"