※ WRF のインストールについては、WRFの導入へ。
新規プロジェクトのケース名を "em_9bar" とする。
em_9bar : wrf @ echo '--------------------------------------' ( cd main ; $(MAKE) MODULE_DIRS="$(ALL_MODULES)" SOLVER=em IDEAL_CASE=9bar em_ideal ) ( cd test/$@ ; /bin/rm -f wrf.exe ; ln -s ../../main/wrf.exe . ) ( cd test/$@ ; /bin/rm -f ideal.exe ; ln -s ../../main/ideal.exe . ) ( cd test/$@ ; /bin/rm -f README.namelist ; ln -s ../../run/README.namelist . ) ( cd run ; /bin/rm -f ideal.exe ; ln -s ../main/ideal.exe . ) ( cd run ; if test -f namelist.input ; then \ /bin/cp -f namelist.input namelist.input.backup ;\ fi ; \ /bin/rm -f namelist.input ; ln -s ../test/$@/namelist.input . ) ( cd run ; /bin/rm -f input_sounding ; ln -s ../test/$@/input_sounding . ) @echo "build started: $(START_OF_COMPILE)" @echo "build completed:" `date`いくつかのテストケースにあるような
( cd test/$@ ; /bin/rm -f gribmap.txt ; ln -s ../../run/gribmap.txt . ) ( cd test/$@ ; /bin/rm -f grib2map.tbl ; ln -s ../../run/grib2map.tbl . )という記述は、おそらくいらない。
module_initialize_9bar.o : \ ../frame/module_domain.o \ ../frame/module_configure.o \ ../frame/module_state_description.o \ ../frame/module_timing.o \ ../frame/module_dm.o \ ../share/module_io_domain.o \ ../share/module_model_constants.o \ ../share/module_bc.o \ module_init_utilities.o
dyn_em/module_initialize_*.F で設定する。
grid%ht(i,j)を変更する。関連する配列については、次のようにすると良さそうである。
grid%phb(i,1,j) = g*grid%ht(i,j) grid%php(i,1,j) = 0. grid%ph0(i,1,j) = grid%phb(i,1,j)このようにすると、以下の処理で上空の geopotential height も設定してくれるようである。
-> csh run_me_first.csh -> ./ideal.exeを実行した後、wrfinput_d01 (netCDFファイル) を用いて、変数 hgt を表示させれば値がわかる。
dyn_em/module_initialize_*.F で設定する。
xland(i,j)これが 1 ならば陸、2 ならば海(水の上)を表す。
dyn_em/module_initialize_*.F で設定する。
nameist.input で設定する。
periodic_x = .true., symmetric_xs = .false., symmetric_xe = .false., open_xs = .false., open_xe = .false., periodic_y = .false., symmetric_ys = .false., symmetric_ye = .false., open_ys = .true., open_ye = .true.,それぞれ、x,y 方向について、周期境界条件か、対称な境界条件か、オープンな境界条件か。
nameist.input で設定する。
mp_physics = 1
non_hydrostatic = .ture.
ra_lw_physics = 1
ra_sw_physics = 1
slope_rad = 1, top_shad = 1,斜面と地形による日射量の変化は考慮しないならば、0 にするのではなく、書かないのが正解のようだ。
sf_surface_physics = 1,※ 5層で温度を計算
sf_sfclay_physics = 5,※ 5 : MYNN
bl_pbl_physics = 5,※ 5 : MYNN 2.5
cu_physics = 14
sst_skin = 1長期間積分に便利な、SST の更新