Lawson 1995について

海洋生態系モデルでデータ同化を試みた論文。

自分もやってみようと思いpythonでコードを書いてやってみた。

いくつか注意点が ↓

1. シミュレーションは１５日間。タイムステップは0.001d（毎日1000ステップ）。
2. 同化するデータ（あるパラメーターで走らせたモデルの解）は 最初の５００ステップ（半日）の一つ飛ばし→２５０データポイント。

 これをpythonのoptimize.leastsqを使ってやってみようとするんだが初めに予測するパラメーターが解と程遠く離れている場合以下のようなエラーが発生する：

  lsoda--  warning..internal t (=r1) and h (=r2) are
       such that in the machine, t + h = t on the next step
       (h = step size). solver will continue anyway
      In above,  R1 =   .1482457577093E+06   R2 =   .1149053549906E-10
 lsoda--  at current t (=r1), mxstep (=i1) steps 
       taken on this call before reaching tout   
      In above message,  I1 =       500
      In above message,  R1 =   .1482457577093E+06
Excess work done on this call (perhaps wrong Dfun type).
Run with full_output = 1 to get quantitative information.

 leastsqの限界なのか 。。。

↑例えば以下のパラメーターの場合：
p_data = [1.,1.,4./spd,-2./spd,-4./spd,-6./spd,2./spd,4./spd]
p_model = [1.,1.,4/spd,-2./spd,-4./spd,-6./spd,2./spd,4./spd]


 三番目の4/spdというパラメーターの初期のguessを0.1までは下げることはできたが0またはそれ以下にするとさっきのエラー文が発生する。同様に、６番目の-6/spdは-5は大丈夫だが-4にするとダメ。

よくわからん。。。

子供が生まれてからすることリスト

ソース→http://www.islamcan.com/raising-children-in-islam/general-instructions-concerning-a-newborn-for-muslim-parents.shtml#.VIUOJKbCc7w

• 初めて入浴させた後すぐに、アザーン（右耳）とイカーマ（左耳）を赤ちゃんの耳元で唱える。
• Tahnik　デイツを食べさせる（スコーラーに頼んだ方がベター）
• ７日目にすると良いこと
• 命名。
• アキーカ。
• 頭剃。髪の毛の重さの分のシルバー（銀）を貧しい人に与える。
•  
• アキーカ：７日目が一番よい（肉）。遅くても思春期を迎える頃までにすること。坊主にすること。goat is best, if dont have then sheep, then cow. 2 for boy, 1 for girl. 1 is ok for boy, but 2 is recommended. slaughter with intention for the baby.
• BC州に出生届(https://www.vs.gov.bc.ca/births/）して同時に以下の書類・証明書登録・発行:
• 出生証明書（$27)
• 社会保障カード
• 健康保険証
• 子供手当
• 日本に出生届（バンクーバー領事館に出向き３ヶ月以内にすること→http://www.vancouver.ca.emb-japan.go.jp/jp/consular_j/consular_services/koseki_shomei_j/shussei_todoke_j.htm)

Telusに電話した件

１１＆１２月分のインターネット料金が$3上がっていた（$92→$95）のでTelusに問い合わせてみるとどうやら１１月からネット料金がそれだけ上がったよう。９月から６ヶ月間、インターネット１５と同じ料金でインターネット２５が使えるプロモーションをもらったんだけど。。。ってことを話してみたら、そしたら１１＆１２月分は$3ずつクレジットを入れてくれて、さらにこれからまた６ヶ月は$5のディスカウントを継続してくれるとのこと。

結果、来月（１月１１日請求分）は$88.69で、それ以降は６ヶ月間、$91.84になる。オファーが消える７月２７日に再度連絡する必要あり。