Rayinvr installation on Linux

Rayinvr是一隻能夠計算二維走時以及振幅modeling的小程式，相信許多玩RayGUI的朋友都知道他。只是許多人只是想要用這隻程式，並不想要真正用到這麼大型的程式。你可以下載Rayinvr原始碼，或是直接連到Rayinvr的網頁。

不過因為大部分的使用者並不能成功編譯這隻程式。主要原因當然是萬惡的Fortran77，而作者也只靠fortran77來撰寫程式，這使得許多人沒有如此老舊的編譯器及古董級函式庫，根本無法成功編譯或執行。有些人退而求其次試用intel的fortran編譯器，可是這會變得讓程式只能在自己的系統上跑，我非常不贊成用這類型的編譯器。幸好Scott Pearse救了作者 Colin Zelt一命 -- Scott寫了一個patch讓新一代的地球物理學家也能使用這隻程式，本文的目的就是紀錄這隻程式如何使用，並幫助各位成功安裝rayinvr。

1. 準備好編譯器以及所需的libraries

• For Debian/Ubuntu users:

aptitude install fort77 libf2c -y
ln -s /usr/bin/f77 /usr/bin/g77

• For CentOS/Fedora users:
以我自己的Fedora14為例，套件名稱會像下面這樣：

yum install compat-libf2c-34 compat-gcc-34-g77-3.4.6-20.fc14.i686 -y
ln -s /usr/bin/g77 /usr/bin/f77
ln -s /usr/lib/libg2c.so.0 /usr/lib/libf2c.so.2

如果你的版本不同，你需要用yum search找確切的套件名稱。你或許會發現兩套linux對於fortran77的稱呼剛好相反。

2. 下載rayinvr以及patch

cd ; wget -c http://terra.rice.edu/department/faculty/zelt/rayinvr.tar.gz
wget -c http://terra.rice.edu/department/faculty/zelt/rayinvr_g77_patch.tgz
tar zxvf rayinvr.tar.gz && tar zxvf rayinvr_g77_patch.tgz
cp rayinvr_g77_patch/rayinvr_g77.diff rayinvr_g77_patch/build_all.bash rayinvr/

3. 開始編譯:

cd rayinvr/
sed -i 's/usr\/X11R6\/lib/usr\/lib\/X11/g' rayinvr_g77.diff
sed -i 's/usr\/X11R6\/include\/X11/usr\/include\/X11/g' rayinvr_g77.diff
patch -p0 < rayinvr_g77.diff

打完patch應該顯示下列結果：

patching file misc/iseis.f
patching file misc/order_picks.f
patching file misc/phercount.f
patching file misc/reciprocity2.f
patching file misc/reciprocity.f
patching file misc/resample.f
patching file misc/smooth_picks.f
patching file misc/tx2pois.f
patching file misc/tx2rec.f
patching file misc/txmatch.f
patching file misc/txoffset2.f
patching file misc/txoffset.f
patching file misc/txphase.f
patching file misc/txshift.f
patching file misc/vex.f
patching file misc/xshot.f
patching file pltlib/pltlib.f
patching file pltsyn/Makefile
patching file pltsyn/pltsec.f
patching file rayinvr/main.f
patching file rayinvr/Makefile
patching file rayinvr/plt.f
patching file rayinvr/rngkta.f
patching file rayinvr/trc.f
patching file rayplot/Makefile
patching file rayplot/rayplot.f
patching file tramp/amp.f
patching file tramp/Makefile
patching file tramp/plt.f
patching file tramp/rngkta.f
patching file tramp/trc.f
patching file vmodel/Makefile
patching file vmodel/plt.f

另外建議大家額外在補個增加解析度的patch上來:

wget -c http://seismic.ocean.dal.ca/pwp_wiki/static/upload/rayinvr_hiprecision.diff
patch -p0 < rayinvr_hiprecision.diff

打完patch正常結果應該如下：

patching file read.me.precision
patching file vmodel/vmodel.par
patching file vmodel/main.f
patching file rayinvr/main.f
patching file rayinvr/rayinvr.par
patching file misc/v2c.com
patching file misc/v2c.par
patching file misc/rayinvr.par
patching file misc/v2c.f
patching file misc/c2v.f
patching file misc/README
patching file misc/dmplstsqr.f
patching file examples/setup9
patching file examples/d9.in
patching file examples/tx9.in
patching file examples/r9.in
patching file examples/v9.in
patching file examples/vm9.in

開始進行自動編譯：

mkdir bin
./build_all.bash | tee rayinvr.log 2>&1

你可以參考rayinvr.log來檢視是否有任何錯誤發生。如果有，請把這個log寄給Colin Zelt，然後他會告訴你：『I do not help people install the code on their computers, since I don't know about such things.』。也千萬不要寄給我，我很討厭fortran。

4. 設定環境變數：
至此你已經完成編譯。現在你只需要把這東西放到/opt裡，再設定環境變數即可：

cd ..; rsync -av rayinvr /opt/

新增(append)幾行到/etc/profile裡：

export Rayinvr="/opt/rayinvr"
export PATH="$PATH:$Rayinvr/bin"
			

然後重新登入(re-login)，就可以快樂的使用這些指令了。關於rayinvr的操作，請參考射線追縱法與時間反演法的使用

後記：
1. I HATE fortran...想退休了，不想在寫文章描述如何debug該死的fortran。拜託新一代的地球物理學家，麻煩好好學fortran以外的語言...別葬送你們的青春在fortran上！
2. Fedora的使用者跑到patch -p0 < rayinvr_g77.diff可能會出現下列錯誤：

patch: **** rejecting target file name with ".." component: ../rayinvr_g77/misc/iseis.f

原因或許出在先前有cracker利用patch漏洞可任意使用../..來當作檔案路徑，理論上這應該是所有linux的問題，但只有Fedora出現，恐怕真的是Fedora比較安全＠＠...
解決的方法很簡單，從非Fedora的linux上複製/usr/bin/patch給Fedora服用即可。
3. 至於另一個程式zplot，請參考Wanna use zplot? Give up!!一文。

Share