2014年1月21日火曜日

IDLのプロシージャ(.proってついてるやつ)を追加して使う方法&円形にarrayに数値代入する方法

クソ雑魚でめっちゃ苦労したのでメモしておく。たぶんまた忘れるだろうから。任意の場所にあるプロシージャ(.proとついてるプログラムをこう言うらしい)をIDLで使う方法。

hoge/hogeってところにIDL_codeってディレクトリがあってそこにプロシージャ放り込む。でhomeにある.bashrcに

IDL_PATH=/User/username/hoge/hoge/IDL_code
(export IDL_PATHが必要なら追加)
を追加。export IDL_PATHはそもそもIDL本体のlibへのパス書いてあるのを通すために書いてあるからいらない。
参考:http://shirokumabook.wordpress.com/2011/07/16/idl-ディレクトリへのpathの通し方/
こっちのほうがもっと一般的。

で、円形にarrayに数値代入する方法。というかnan打ち込む方法。
http://astro.berkeley.edu/~johnjohn/idlprocs/make_circ.pro
にあるmake_circ.proを/IDL_codeへ保存。bashやって通しておく。文法は
IDL> result = make_circ(input_name, radius, cent=[X,Y], val=!VALUES.F_NAN)
でnanが円形状に入る。

2013年12月13日金曜日

sshでつないだ時にx11windowが時間が経つとタイムアウトするのを防ぐ方法

つなげる側のPCの/etc/にあるssh_configで
x11Forwarding yes
x11ForwandingTrusted yes
の2行を追加する。
ググって同じようなのあったんだけど勘違いしてて、x11Forwarding Trusted yesとForwardingとTrustedの間に半角スペース入れてて機能してなくて困った困ったしてたけど、半角入れないとできた。

あと公開鍵暗号で設定するとパスワード打つ必要なくなって捗ることを知った。

2013年12月11日水曜日

awkを使ったターミナル上の1行コマンドで.datなどのファイルのデータを計算したりファイル作ったり

.datファイルなどにあるデータの合計や平均値、そのデータを用いた計算をわざわざプログラムファイルを作らずとも、ターミナル上のコマンドでやる方法。個人的学習メモ。awkを使った1行プログラムでdatファイルを処理する。

 1.連続した数字を並べたdatファイルを作る
 awk 'BEGIN{for(i=1; i<=20; i++) print i}' > test.dat
と打つと 1から20までならんだtest.datファイルが作られる
1
2
3
.
.
19
20
これの平均を求める
awk '{sum += $1}END{print sum/NR ;}' test.dat
10.5
とでる。test.datの最後の行にave10.5を追加するには
awk '{sum += $1}END{print "ave""\t" sum/NR ;}' test.dat >> test.dat
と書く。>で上書き保存、>>で後ろに追加。"\t"でタブ。NRはその時点で読み込んだ行数。この行数は#でコメントアウトした行も含まれる。ENDで最後読み終わったあとの処理の意味。BEGINはファイル読み始める前に処理の意味。

2. 1行ごとに処理(計算)して出力
例えば
 awk 'BEGIN{for(i=1; i<=20; i++) print i*3}' > test2.dat
した
3
6
9
.
.
57
60
というtest2.datファイルを作ってtest2.datをtest.datでそれぞれの行を割りたい時。 pasteコマンドを使ってくっつけて|で受け渡してawkで処理して出力。
paste test.dat test2.dat | awk 'BEGIN{for(i=0; i<=20; i++) NR>i}{print $2/$1;}' > test3.dat
とすると
3
3
3
.
.
3
3
というtest3.datができる。ミソはBEGIN{for(i=0; i<=20; i++) NR>i}部分と思われる。for文で1つずつiを増やしていって{}に入らない条件部分でiより大きいNRの行数を見て{print $2/$1;}を処理。forのループでまた1つ行がずれる。(という理解。たぶん)

プログラミング雑魚なので今日やったことをとりあえずメモしておく。printで;つけてなくて出力できなくて30分以上悩んでた雑魚です。これ覚えておけば、解析でのdatファイルの扱いはなんとかなるはず。

2013年11月14日木曜日

研究ノートは研究日記

研究ノート is 何?

誰も教えてくれないというか恥ずかしくて聞けなかったというのが本当のところなんだけど、4月に大学院入学してもう半年以上経ってるのに研究ノートというのがいまいちよくわからなくて困ってた。研究メモだよとか言われてもよくわかんないし、何書いたらいいんだよって思ってて時々メモ帳代わりに使ってただけ。

で、ここ最近というか今日なんだけど、観測実習の準備と来日してるCrossfieldからのアドバイスをメモったりそれ参考にしてやったりしててで気づいたんだけど、研究ノートというのは研究「日記」だということ。今日自分が何をしたのか、何のコマンドを使ったのか、明日のタスクは?というのを書くものだということだと理解しました。やっとかもしれんが。

いやまあそれを気づくに至った理由なんですが、至ってシンプルであれこれ前やった時どの数値、コマンド使ったっけというのが全くを持って忘れててまたグーグル先生に聞き直したり計算しなおしたりしてだるかったからです。

じゃあこれまでどうしてたのかというとgeneral educationが重要ということで普通に勉強なわけでメモする重要性を感じてなかったんだよね。ゼミの準備とかも普通にゼミ資料とかプレゼン資料で残ってるわけだし、わざわざノート書く理由がなかった。

研究ノートって実験ノートと同じじゃんって思ったけど、実験ノートと違うのは論文の存在。読んだ論文とか参照する論文とか感想とかって実験ノートじゃ書かないじゃん。でも、日記なら書くでしょ。ログなわけだ。忘れやすいから。

だから、研究ノートは研究日記。些細なことでもあとで読み返してわかるように日記のように書こう。
というので自分では腑に落ちました。

というのを忘れないうちに書いておきます。

好きなものリスト

研究で疲れて帰ってきてブログ書く気力がなかった上に艦これやりまくってたからブログ書いてなかった。

なんていうかこの間の京都伊勢高山弾丸旅行の時にも話したんだけど、自分が好きだったものが好きじゃなくなるってすげー悲しい。なんていうかな…あれほど情熱的に熱狂的に好きだったものに対する魔法が解けたように自分が変わってしまったということがわかってしまう。いやまあしょうがないというか変化は仕方のないことでそうなんだけど、あの時の気持ちも本物で今の気持ちも本物だとするといつかは全てのものに飽きが来ちゃってたぶん自殺すると思うんだよね。怖いわけじゃないけど、あの時の気持ちを思い出せなくなっちゃうのは悲しいから今好きなものをちょっとリスト化しようかなと思いました。2013年に買ったものとかかな。twitterでいろいろとぼやいてるからtwilog漁るのでもいいんだけど、それだと探すのがだるい。

・輪るピングドラム
BDついに買いました。全話見直すのは時間なくて無理だけど1ヶ月ぐらい前に1話と10話見て、18話見たらやっぱり泣きました。愛の話なんだよなあ。

・やはり俺の青春ラブコメはまちがっている。
アニメの3話あたりでお前が俺を交通事故に遭わせたことに負い目があって優しくしているならやめろというヒッキーのセリフで一気に好きになりました。BDは買ってないけど原作は全巻買った。

・GJ部
フェチアニメ。かわいい。BD買った。2軍のファンタジスタTシャツ欲しかった…。2期楽しみ。チカラメシとのコラボにも行ったけど楽しい。ダメダカラネー

・ガールズアンドパンツァー
12年の11月と13年の5月に大洗行った。結構劇中の場所回る聖地巡礼楽しい。あと5月行った時はガルパンかくれんぼやってて等身大パネルが街中にあって捗った。また折を見て行きたい。あ、BDも買いました。OVAは買ってない。それ以外のグッズは買ってないな―なんかアニメで完結してる。

・艦これ
7月に舞鶴に着任して今、提督Lv.99。先輩も某スタッフもやってる。イベントで武蔵取ったけど、最近ゲームバランス崩れすぎてゴミ化しそう。5-3くそだる。自分の気に入った艦娘は集中的に育てて装備とかも揃えるけど、嫌いな朝潮型とかどうでもいいのは弾除けとかストレス溜まった時に轟沈させたりしてるある意味ブラック鎮守府です☆ 五十鈴と那珂ちゃんは割りと理不尽にいじめたい。ストレスたまってるんですかね。運ゲーにシフトしてゴミゲーになりつつあるので、どうにかしてください。とりあえずまだやってるし、気に入った子がいるからやめないけど、アニメ化ってどうなるんでしょ。

こんなもんですかね。アイマスは赤羽のシークレットイベントの発表部分だけ中継みたけど、割りと新作発表でもときめかないで今頃かよって感じになってた。たぶん完全に情熱が失せてる。春香さんアイコンは替えないけど。春香さんに対してはやっぱり好きとか愛してるとかじゃなくて完全に刺さってるんだと思います。