--- ungrib 失敗 log ---
失敗 1：
Vtable.GFSに比湿の項を加えてVtable.JRAにする。
先にNCEPをVtable.GFSでungribして、そのあとJRAをungrib。
その後、
NCEP, JRA
としてmetgridすると、
NCEPの鉛直格子数が27層、JRAの鉛直格子数が100hPaより上層に6層あって、
その関係で、全体のRHが33層になって上層部分が足りないとおこられる。
おそらく、100hPaより上層のSpecfic fumidityが足りないと考えられる。
失敗 2：
さっき作ったVtable.JRAからRH、SPECHUMDを消し、
Vtable.GFSからはSPECHUMDを消し、
再度NCEP, JRA
としてmetgridすると、
WRF_DEBUG: Warning DIM 4, NAME num_metgrid_levels REDIFINED by var
RH 24 23 in wrf_io.F90 line 2424
ERROR: Error in ext_pkg_write_field
とErrorを吐いて、置き換わるはずのRHデータ数が合わなくて怒られる。
失敗 3：
ungrib_jra.shを作成して、MSM方式でやってみる。
この処理中に、サブルーチンとして、比湿→相対湿度の変換プログラムを入れて、
ungribと同様の処理を行う。
結果、ungrib
metgridともに成功。
(i) NCEP, JRAで結合
→ 何故かwrf.exeでエラーがでる。
(ii) NCEP, JRA, MSMで結合
→ きちんとwrf.exeまで成功！！
ただ・・・これならNCEP, MSMってやったほうがいいような・・・
しかし、10hPaまで作成した分のRHが置き換わってて、走らせる分には問題ないかと。
構想 1：
GMS使って、
NCEP、GSM、MSM
でmetgridする。
NCEPはVtable.GFSでungrib
GSMは、*Rgl*のデータ, Vtable.JMAGSMを使用してungrib
MSMは、ungrib_msm.shを使用。
--- ungrib 試行錯誤 log ---
NCEP, JRAの組み合わせ
で初期値を作ろうとすると、
Domain 1: Current date being processed: 2008-05-13_00:00:00.0000, which is loop # 1 out of 3
configflags%julyr, %julday, %gmt: 2008 134 0.0000000E+00
d01 2008-05-13_00:00:00 Timing for input 0 s.
d01 2008-05-13_00:00:00 flag_soil_layers read from met_em file is 1
i,j = 77 10
p column = 100000.0 97500.00 95000.00 92500.00
90000.00 85000.00 80000.00 75000.00 70000.00
65000.00 60000.00 55000.00 50000.00 45000.00
40000.00 35000.00 30000.00 25000.00 20000.00
15000.00 10000.00 7000.000 5000.000 3000.000
2000.000 1000.000
z column = 8.453903 310.0116 534.8116 662.3279
999.4542 1363.149 2005.928 2551.715 2914.663
3744.849 4099.316 5101.997 5458.782 6682.821
7062.484 8583.468 9040.818 10225.21 11644.86
13516.71 16140.16 18428.48 20586.80 23859.78
26482.62 31062.11
model topo = 1.092892
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: module_initialize_real.b LINE: 4812
probs with sfc p computation
-------------------------------------------
と、怒られる。
NetCDF：
-- FORTRANでのNetCDF読み書き (2010/06/24 更新)
FORTRANでNetCDFファイルの変数を書き換えたりするばあいは、
status = nf_open(ifilename, nf_write, ncid)
として、’読み書きOK’にしておく。
書き換える変数は、読み込まなくていい。
WPS：
-- 出力形式について -- (2010/06/24 更新)
WPSでの出力形式は、namelist.wps内の
"IO_FORM_GEOGRID" , "IO_FORM_GEOGRID"の番号で指定する。
出力可能形式は、
1; binary
2; NetCDF (default)
3; GRIB1
の 3つになっている。
ungribでの出力形式は、WPSか、WRFSIか、MM5になっている。
デフォルトでは、WPSなので、こちらはいじらない。
WPSでの形式については、ARWUSERSNOTE(v3.1) 3-30を参照。