How to convert the original Lana DMS climatology (in the CSV format) into the netcdf format

I used this script (http://stackoverflow.com/q/22933855) as a reference to create my Python script to convert the CSV formatted Lana DMS climatology (http://www.bodc.ac.uk/solas_integration/implementation_products/group1/dms/) into the netCDF format. Here is my script:

import netCDF4
from pylab import *
from pandas import read_csv

ncout = netCDF4.Dataset('test.nc','w')

data=array(read_csv('../dmsclimatology/DMSclim_APR.csv',header=None))
dim0 = size(data,0)
dim1 = size(data,1)
lats_out = 89.5 - 1.0*arange(dim0,dtype='float32')
lons_out = -179.5 + 1.0*arange(dim1,dtype='float32')

ncout.createDimension('latitude',dim0)
ncout.createDimension('longitude',dim1)

lats = ncout.createVariable('latitude',dtype('float32').char,('latitude',))
lons = ncout.createVariable('longitude',dtype('float32').char,('longitude',))

lats.units = 'degrees_north'
lons.units = 'degrees_east'

lats[:] = lats_out
lons[:] = lons_out

dms = ncout.createVariable('DMS',dtype('float32').char,('latitude','longitude'))

dms.units =  'nmol L-1'

dms[:] = data
ncout.close()

If you want to use this script, you would need to make sure you have those three modules and change the file name associated with the read_csv function to your own.

The output looks like this (on the right) and you can confirm that this is consistent with the published figure (on the left). NOTE: the color bar scales are different, so they don't look exactly the same.

k_dmspd

今の所k_dmspdはバクテリアが行うので光が存在する時はほぼ0になるようにしているがこの場合、k_dmspの数値を上げても氷が解けた後の北極の場合光が24時間存在するためdmspdの濃度が思ったよりも下がらない。

dmspd濃度.k_dmsp=3（左）とk_dmspd=13（右）の場合

光によるk_dmspdの影響をもう一度考え直す必要がありそう。細菌は光に対しそこまで柔じゃないのだろう。

冬のバンクーバー

学会で月曜から金曜までバンクーバーへ行ってきた。今回は車で行ったのでリッチモンドへ寄れた。まさにアジアの町という感じで、ダイソーがあるAberdeen Centreでは自分がいた時間帯は客のアジア（東洋）人率100パーセントだった。セキュリティガードのおじさんたちが唯一南または中東系だった（⇦それでもアジアw）。まぁ店はダイソー以外は特にいい店はなかった。Saba Streetにあるハラール中華料理店Silkwayはあいにく閉まっていた。

宿泊はダウンタウンのウエストエンドでロブソン通りというかなりいい立地のBarclay hotelに泊まった。エレベーターがないのでベビーカーや車椅子の人にはお勧めできない。ダウンタウンで安く泊まりたいけどホステルはちょっとって人にお勧め。冬のオフシーズンなので一泊$50だった。

以下ダウンタウンで行ったレストラン＆感想：

Banana leaf：マレーシア料理。お店の人が子供の対応に親切で好印象。料理はAuthenticというよりはカナダ人向けに味付けられた感じ。食べたのはRoti canaiとSeafoodカレー＆炒め物。

Pappa Roti:これまたマレーシア発のカフェ。Roti自体はメキシコ発祥らしい。Rotiひとつ$3.5と高い。店の雰囲気はいいけど。

Gyoza bar：トマトスープベースの海鮮ラーメンがあるということで期待してランチで行ってみたが期待はずれだった。スープは美味しいけど拉麺とは合わない。餃子もエビとベジタリアン両方食べてみたけどどっちもフュージョン。Authenticだと思っていたので余計残念。

Guu izakaya Garlic：Authentic！普通に日本の居酒屋という感じで賑やかな店内でお店の人の叫び声が飛び交う。子供も連れて行ったが丁寧に対応してくれた。座敷に座れた。おでんが美味しい。値段も高くないのでお勧め。

Shuraku restaurant：これまたAuthenticな日本食レストラン。値段もそんなに高くない。 鰻のひつまぶしをいただいた。

Konbiniya：おにぎりが２個で$3なのが嬉しい。賞味期限切れ間近のカルピスが$1で売られていたのも良かった。

excelやspreadsheetで各列（または行）で同じ定数で割り算する

雲のデータの値を半減したかったのでLibreofficeのspreadsheetでデータを開いて各列２で割ろうとすると一列目を

=A1/2

とした後に計算されたセルを選択して右下の四角を以下の列にドラッグすれば

=A2/2
=A3/2
・・・

と計算してくれるんですがデータが多いとこれは面倒くさい。

しかし全ての列で割り算をしてくれる方法がありました。

それは

「ドラッグの代わりにダブルクリックする」

という簡単な方法でした。

以下参照：
http://superuser.com/a/180357

余談：計算が終わったら元のデータはいらないけど消すと半減されたデータも無くってしまう（式を使っているので元のデータに依存している、例えば=A1/2）のでこれを防ぐためには半減されたデータをコピーしてペーストするときにPaste Specialを選択してOperationsのDivideのところをNoneに変えれば消えずにすむ。

インターネットの使用容量

これは自分の過去１年間のインターネットの使用容量を表したグラフ。

カナダのインターネットプロバイダー各社が今年に入ってから容量の制限を始めたので使用容量を気にし始めなければいけないなった。

２月を除けば月平均は50GBといったところか。これならShaw、Telus各社の一番低いプランでも制限内なので特に気にすることなさそう。