WRF 3DVAR边界更新

flying鹰

吴阿好大人 发表于 2016-9-23 20:04
那么请问是不是直接先更改运行时间，再real，生成新的初边条件再去更新下边界？

我前边写的很清楚了，你./real.exe之后不就重新生成wrfinput文件了么？！那你做同化干啥？！

flying鹰

吴阿好大人 发表于 2016-9-23 20:04
那么请问是不是直接先更改运行时间，再real，生成新的初边条件再去更新下边界？

我这个帖子写的也很清楚，请仔细阅读

吴阿好大人

flying鹰 发表于 2016-9-24 10:04
我这个帖子写的也很清楚，请仔细阅读

您好，我仔细看了你写的步骤，在第三步没有提到'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 .

所以我想问一下是不是用wrf_3dvar_input_d01_06来驱动6时到12时的real.exe，然后利用此时生成的wrfinput_d01 来为wrf_3dvar_input_d01_06更新下边界，然后同化，然后更新侧边界，然后在进行模拟？
自学wrfda感觉非常吃力，多次打扰您非常抱歉！谢谢！

zys1022

在学习中，谢谢楼主

wuyejincao

lee309 发表于 2015-1-9 18:14
循环同化终于成功了，但是这个结果让人失望，还不如同化之前的结果准确，这还有什么意义啊！
不知道大家都 ...

您好，我也想问跟您一样的问题，侧边界和下边界可以一起同化么？

dongqiru77

谢谢楼主，我用官网的例子很快就运行出来了，但是原理都没搞清楚，很感谢博主自己理解！

m337872010

楼主分享，您写的东西十分受用！我还想请教个问题
1、我用的观测数据是Micap2转化的Little_R格式，有高空和地面的，那么在WRFDA/var/obsproc/里面也会对应生成两个obs_gts_....3DVAR文件；那么在./da_wrfvar.exe里面是不是也要用两个.3DVAR运行两次
2、无嵌套的是不是不需要da_update_ba.exe这一步了

m337872010

楼主还在吗。。还是想问一下，我的实测数据是高空和地面的。。那么在obsproc里面能生成一obs_gts_.3DVAR文件吗？生成两个的话是链接不同的ob.ascii吗？还是需要分两次wrfdar？

flying鹰

m337872010 发表于 2017-11-11 11:57
楼主还在吗。。还是想问一下，我的实测数据是高空和地面的。。那么在obsproc里面能生成一obs_gts_.3DVAR文 ...

探空和地面常规观测先转换为成little-r格式，然后就可以用obsproc处理，出来的文件是obs_gts...3DAVR。地面探空可以同时同化，但是涉及到局地化半径不同，那就可以分两次同化

m337872010

flying鹰 发表于 2017-12-9 12:22
探空和地面常规观测先转换为成little-r格式，然后就可以用obsproc处理，出来的文件是obs_gts...3DAVR。地 ...

谢谢楼主！地面和探空的已经转成little_r了，是Micaps的第一类和第二类的数据，然后在namelist.obsproc里面如何可以同时通化这两个呢？
&record1
obs_gts_filename = 'surface_obs_r_2008-06-18_00 , obs_06_18_00',
obs_err_filename = 'obserr.txt',
gts_from_mmm_archive = .true.,
/

&record2
time_window_min  = '2008-06-17_23:00:00',
time_analysis    = '2008-06-18_00:00:00',
time_window_max  = '2008-06-18_03:00:00',
/

&record3
max_number_of_obs        = 400000,
fatal_if_exceed_max_obs  = .TRUE.,
/