WRF 3DVAR边界更新

flying鹰

到WRF-3DVAR这一部分，我自己认为看手册更新边界有点难理解，论坛里又没有，我自己的理解，分享给大家～如有认为错误的请批评指正～～～
！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！
手册说明：
cold-start run                      ——只有侧边界（lateral boundary）需要更新,无需更新下边界（low boundary）
warm-start（cycling） run  ——侧边界（lateral boundary）、下边界（low boundary）都要更新
                                             ——但是下边界（low boundary）要在运行WRFDA前更新
嵌套时：dommain_id=1  侧边界（lateral boundary）、下边界（low boundary）都要更新
              dommain_id>1  只有下边界（low boundary）需要更新
下边几句是WRF官网论坛里的（网址：http://forum.wrfforum.com/viewtopic.php?f=28&t=371）
&Use wrfbdy_d01 from WPS/real.exe.
&da_update_bc.exe only processes the first two time periods of data in wrfbdy_d01.
&New boundary tendency of time1 and time2 is recalculated
&using the WRF-Var updated analysis valid at time1 and the existing time2 value in wrfbdy_d01.
以同化一重区域，00：00和 06：00为例
（1）～～～～～～～～～更新侧边界～～～～～～～～～～～～～～～～～～～～～～～～～～～～～
00：00是cold-start run
cd $WRFDA_DIR/var/test/update_bc /00  
cp $WRFDA_DIR/var/test/update_bc/parame.in .                                                           
cp  /home/lishyi/WRFV3run/wrfbdy_d01 .  （ from WPS and WRF real ( make a copy of wrfbdy_d01, as the wrf_bdy_file will be over-written by da_update_bc.exe)
cp  /home/lishyi/WRFV3/run/wrfinput_d01 . （ from WPS and WRF real）  
cp  /home/lishyi/WRFDA/workdir/00/wrfvar_output .（wrfvar_output是同化后输出文件）
ln -sf $WRFDA_DIR/var/da/da_update_bc.exe .

     parame.in文件设置
    domain_id= 1  
    cycling = .false .                                                                        
    update_lateral_bdy = .true.                                          
    update_low_bdy     = .false.   
                                                         
然后运行 ./da_update_bc.exe
此步更新侧边界 (wrf_bdy_file will be overwritten by da_update_bc).
（2）～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～
将同化输出的wrfvar_output及更新的边界拷贝到……WRFV3/run/目录下
cp $WRFDA_DIR/var/test/update_bc/00/wrfbdy_d01 .
cp  /home/lishyi/WRFDA/workdir/00/wrfvar_output wrfinput_d01
在namelist.input中我设置积分时间是00：00到06：00
同时有（可参看官网手册）
write_input = .true.
inputout_interval = 360
input_outname = 'wrf_3dvar_input_d<domain>_<date>'【我设置的是‘wrf_3dvar_input_d01_06’】
inputout_begin_h = 6
inputout_end_h = 6
运行过后输出文件有wrf_3dvar_input_d<domain>_<date>【我的文件名'wrf_3dvar_input_d01_06】
（3）～～～～～～更新下边界～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～
06：00是warm-start run
在namelist.input里设置06：00到12：00，real.exe后生成wrfbdy文件和wrfinput文件
cd $WRFDA_DIR/var/test/update_bc /06
cp $WRFDA_DIR/var/test/update_bc/parame.in .
cp  /home/lishyi/WRFV3/run/wrfbdy_d01 . （ from WPS and WRF real）
cp  /home/lishyi/WRFV3/run/wrfinput_d01 .（ from WPS and WRF real）
cp  /home/lishyi/WRFV3/run/wrfvar_input_d01_06  wrfvar_output  
ln -sf $WRFDA_DIR/var/da/da_update_bc.exe .

parame.in文件设置
cycling = .true .                                                                        
update_lateral_bdy = .false.                                                               
update_low_bdy     = .true.                                                        
然后运行 ./da_update_bc.exe
此步更新下边界 (da_file will be overwritten by da_update_bc).
～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～～
然后进入/home/lishyi/WRFDA/workdir/06
将更新下边界后的wrfvar_output链接为fg进行同化
之后就是重复（1）与（2）中的步骤更新侧边界与进行积分

天边的那边

在学习中，谢谢楼主

freekiller

——但是下边界（low boundary）要在运行WRFDA前更新
运行wrfda之前用啥去更新下边界？ 再琢磨一下吧

flying鹰

更新下边界手册上有说到，至于原因……我就说不明白了

flying鹰

freekiller 发表于 2014-6-20 16:34
——但是下边界（low boundary）要在运行WRFDA前更新
运行wrfda之前用啥去更新下边界？ 再琢磨一下吧

更新下边界手册上有说到，至于原因……我就说不明白了

flying鹰

freekiller 发表于 2014-6-20 16:34
——但是下边界（low boundary）要在运行WRFDA前更新
运行wrfda之前用啥去更新下边界？ 再琢磨一下吧

更新下边界手册上有说到，至于原因……我就说不明白了

燕雀要变鸿鹄

楼主请问一下，你的parame.in文件中有update_lateral_bdy = .true. 吗？还是你自己加上的呢？我的文件里有update_lsm=.false.这一设置，这是什么意思呢？谢谢啦~

flying鹰

燕雀要变鸿鹄 发表于 2014-7-9 16:20
楼主请问一下，你的parame.in文件中有update_lateral_bdy = .true. 吗？还是你自己加上的呢？我的文件里有u ...

update_lateral_bdy = .true.是我自己加上去的

update_lsm=.true.: the LSM predicted variables: TSLB, SMOIS, SNOW, SH2O, RHOSN, CANWAT, SNOWH
                                  will be updated based on wrf_input file
false表示不更新
上述是$WRFDA/var/README.basics的解释
我运行的时候设置的是update_lsm=.true.

燕雀要变鸿鹄

flying鹰 发表于 2014-7-9 19:37
update_lateral_bdy = .true.是我自己加上去的

update_lsm=.true.: the LSM predicted variables: TSL ...

谢谢你~我还想问一下如果不加上update_lateral_bdy = .true.，是不是边界条件不会更新呢？谢啦

燕雀要变鸿鹄

flying鹰 发表于 2014-7-9 19:37
update_lateral_bdy = .true.是我自己加上去的

update_lsm=.true.: the LSM predicted variables: TSL ...

还有就是想问一下楼主是不是把low_bdy_only=.false.去掉了，然后换成update_lateral_bdy和update_low_bdy了呢？谢谢楼主