http://w.atwiki.jp/nsota/ nsota @ ウィキ ja 2011-04-27T16:27:21+09:00 1303889241 https://w.atwiki.jp/nsota/pages/25.html パラメトリック手法 母集団の特性を規定する母数についてある仮説を設ける 平均値の差の検定（t 検定）や分散分析（F 検定）・・・母集団の正規性や等分散性が仮定される。 ノンパラメトリック手法（分布によらない手法） 母集団の分布型（母数）について一切の仮定を設けない 標本サイズが小さい場合には，それから求められた統計量の分布型は不正確なことが多く，パラメトリックな手法を適用することは不適切になりやすい パラメトリックな手法 vs ノンパラメトリックな手法 平均値 vs 代表値 分散 vs 散布度 積率相関係数 vs 関連性係数，順位相関係数 ? vs 度数 間隔尺度，比例尺度 vs 名義尺度, 順序尺度, 間隔尺度, 比例尺度 正規分布, 等分散性を仮定 vs 不問 標本サイズ小さすぎてはいけない vs 不問 2011-04-27T16:27:21+09:00 1303889241 https://w.atwiki.jp/nsota/pages/24.html VISTA memo 2010/12/15 「送る」に追加 windowsボタン+R 「shell:sendto」と記入 追加したいプログラムのショートカットを入れる 2010-12-15T16:44:05+09:00 1292399045 https://w.atwiki.jp/nsota/pages/23.html 大阪湾リアルタイムデータ http://kouwan.pa.kkr.mlit.go.jp/kankyo-db/ の大阪湾水質定点自動観測データ配信システムバナー 2010-05-11T15:41:11+09:00 1273560071 https://w.atwiki.jp/nsota/pages/22.html 20km格子全球大気モデルでは強い熱帯低気圧の出現頻度は現実よりも少ない傾向 弱めの低気圧もカウントすることにより、現実の熱帯低気圧（最大風速17.2 m/s 以上）の全球発生数と同程度の数になるように条件調整 下層（850hPa面）の最大風速基準を上下させると、カウントされる低気圧の数は大きく変わる tcstep2.f90を書き換え 発生時の条件(iflag)と、継続時の条件(iflag2)にわけて、継続時の条件に関しては条件を緩和させる 日本列島付近は対流圏ジェットが強いため、850hPaの風速は300hPaの風速に比べて小さくなりやすいことから、w8mw3crit2を大きめに取る これまでは発生数でチューニングしていたが，現在はトラックに合わせたチューニングを試験中 このまま使うとtcstep3で36時間以上継続するものも増えてしまい、台風の個数も増えるので，発生数の過小評価も合わせて調整 GPIは主に気候モデルの台風発生数過小評価の原因解明のため 気候モデルは北西太平洋で台風発生数が少ないのは確かだが，GPIがそこそこ大きいのに発生数が少ないとなるとトラック抽出法にも欠陥があるのではないか 2010-04-26T17:27:19+09:00 1272270439 https://w.atwiki.jp/nsota/pages/21.html Improved subgrid-scale models for large-eddy simulation Bardina, J.; Ferziger, J. H.; Reynolds, W. C. Taylor, Karl E. Summarizing multiple aspects of model performance in a single diagram J. Geophys. Res., Vol. 106, No. D7, pp. 7183–7192 2010-04-16T19:08:23+09:00 1271412503 https://w.atwiki.jp/nsota/pages/20.html CMIP3 >Yokoi(2009) -3rd phase of Coupled Model Intercomparison Project designed for IPCC AR4 24 CGCMs（24モデル） 12 experiments（12実験ケース） 比較的長期（各実験ケースで数十年） 水平解像度は粗い（最大でT106（スペクトル法の解像度）150km程度） 台風強度の再現性は疑問？頻度はそれなりに合う -北西太平洋を対象に（Yokoi） Model List (12) 1. BCCR-BCM2.0, Norway, T63(T42) 2. CGCM3.1(T63), Canada, T63(T42) 3. CNRM-CM3.0, France, T63(T42) 4. CSIRO-Mk3.0, Australia, T63 5. CSIRO-Mk3.5, Australia, T63 6. ECHAM5/MPI-OM, Germany, T63 7. FGOALS-g1.0, China, 2.8deg 8. GFDL-CM2.0, USA, 2.0(lat)x2.5(lon)deg 9. GFDL-CM2.1, USA, 2.0(lat)x2.5(lon)deg 10. INGV-SXG, Italy, T106 11. MIROC3.2(hires), Japan, T106 12. MIR-CGCM2.3.2, Japan, T42 T42~300km アウトプットは日平均値 熱帯低気圧の定義 (1)850hPa相対渦度の局所最大値（中心）が閾値より大きい = 回転渦 (2)中心における300hPa温度と環境の300hPa温度の差が閾値よりも大きい = 温かい核 (3)条件1と2は少なくとも２時刻は維持される (4)発生地点は海である (5)発生時刻に850hPaの最大風速は300hPaのそれより大きい = 追加の台風を除くため -20世紀実験からみた精度について 年の発生個数（平均値）...CGCM3.1(T63)，CSIRO-Mk3.0，CSIRO-Mk3.5，ECHAM5/MPI-OM，INGV-SXGがよさげ 年変動 2010-04-16T16:45:28+09:00 1271403928 https://w.atwiki.jp/nsota/pages/19.html >台風の運動を記述するモデル（２次元台風モデル） Myersモデル 藤田モデル -陸地形状は考慮されていないので必ずしも風速場の再現性はよろしくない >気象モデル 陸上地形の影響が導入される 気象研究所による非静力学モデル 平尾（1999） 　ARPS（メソスケールモデル） 大澤（2001） 　MM5（第５世代メソスケールモデル） -初期気象場の情報不足 台風ボーガス（擬似渦度）の導入によって初期条件をつくる >台風の経路や中心気圧などを確率的に算定するモデル（確率台風モデル） 高橋（1964） 長友（1980） 端野・桑田（1987） 山口 加藤 橋本 -台風経路や台風諸元を確率的に求めた後，風速や気圧分布を求める際には２次元台風モデル（や気象モデル）を使用 -台風のモデル化「台風経路」「中心気圧」「進行速度」 -台風経路「位置」「進行方向」「進行速度」 -中心気圧（時系列の連続性を考慮すべき） 国富（2005）では1950-1995気象庁資料（1228個）に基づいて構築（緯度経度は１度刻み，取得個数が少ない場合は周囲1セルに取得範囲を拡大） -取得データは「中心気圧」「進行速度」「進行方位」「諸元変動率（当該セルを通過するまでに変化した諸元の単位時間変動量）」 2010-04-12T18:42:09+09:00 1271065329 https://w.atwiki.jp/nsota/pages/18.html ENSO　：　エルニーニョ・南方振動 ロスビー波って何？ http://kobam.hp.infoseek.co.jp/meteor/rossby.html まだよくわからん． Broccoli, A. J., and S. Manabe (1990), Can existing climate models be used to study anthropogenic changes in tropical cyclone climate?, Geophys. Res. Lett., 17(11), 1917-1920. 2010-04-21T18:30:53+09:00 1271842253 https://w.atwiki.jp/nsota/pages/17.html track type データがすべて-127 異常値が入っている コメントではタイプintのはずがdoubleで入っているもの lat lon wind pres pos_stdev wind_stdev wind_min wind_med wind_max pres_stdev orig_wind orig_pres これらはdoubleのほうがいいと思うが．．． pres_maxとpres_minはintでいいのか？ distInMeters = distance(latL, lonL,... latK, lonK, earthRadiusInMeters) 現在の格子データへの落とし方は国富（p.45）とやや異なる 2010-04-12T18:29:14+09:00 1271064554 https://w.atwiki.jp/nsota/pages/16.html -[[戻る>http://www37.atwiki.jp/nsota/pages/1.html]] -ファイル操作系 cd c:\matlab\ ディレクトリ移動 cd .. 上部ディレクトリへ移動 dir ディレクトリ内表示 ls ディレクトリ内表示 whos 現在読み込んでいる変数の詳細表示 path パスの表示 memory 使用可能なメモリーについて -便利ツール h1 = waitbar(0,'Please wait ...'); for i = ... waitbar( i / imax ) end 経過時間表示 -NetCDF系 test_snctools NetCDFの読み込みtoolboxであるsnctoolsの動作確認コマンド netcdf NetCDFデータへの低レベル読み込み（古い） nc_dump('C:\Users\Nakajo\Desktop\Allstorms.ibtracs.v02r01.nc') 変数の一覧を表示（書き方は色々ある） genesis_basin=nc_varget('C:\netcdf_test\Allstorms\Allstorms.ibtracs.v02r01.nc','genesis_basin'); 変数の取得（genesis_basinという名前で変数genesis_basinを読み込む） -Mapping toolbox -距離算定 earthRadiusInMeters = unitsratio('earthRadius','meter') （登録）単位換算 ellipsoid = [1.56,0.3]; 回転楕円体の特性パラメター（長軸半径，離心率） distance([0,0],[180,0], ellipsoid) 球面上の二点間距離（大円最短距離？）を求める．[緯度1，経度1]，[緯度2，経度2]（回転楕円体も指定可能） 以下のように単位換算指定もできるが，回転楕円体パラメターと同時指定はできない？ので長軸半径に実際の単位付き距離を入れるのがよい distance([0,0],[180,0] 2010-04-12T17:38:27+09:00 1271061507