登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本帖最后由 今朝几朵白云生 于 2019-3-11 21:55 编辑
前段时间做了一些集合同化的工作,现在整理出来给刚接触WRF/DART的朋友一些参考。如有不当或错误的地方,敬请指正。 1. DART下载与安装 官方推荐的安装方式是通过SVN来管理你的文件,但是绝大多数计算机群处于安全的考虑都不会允许普通用户通过网络和服务器交换数据。因此他们还提供下载压缩包用于本地安装的方法。DART压缩包的下载地址在 这里,我用的版本为 DART_Manhattan_r12989.tar.gz。解压之后(设置当前路径为$DART)第一步就是要配置$DART/build_templates/mkmf.template文件。具体内容是要指定netcdf库的路径和使用何种编译器,根据自己机器的实际情况进行配置即可。之后就是编译同化过程中要用到的可执行文件。进入$DART/models/wrf/work,将模板文件input.nml拷贝到当前目录,执行批量编译脚本quickbuild.csh,完成之后共有27个可执行文件。至此,DART安装完成。 2. Testdata的说明 DART提供的Testdata下载地址在 这里,解压之后(设置当前路径为$WORKDIR)要用到的子文件夹有icbc, obs_diag, rundir, scripts, output, perts, 和template。对于部分目录有必要做一些说明: icbc(initial condition and boundary condition)包含提供初边界条件的grib文件,这里面提供GFS 2017042700-2017043018一日四次每个时刻00时和06时的预报文件。 rundir 集合同化的工作目录,add_bank_perts.ncl用于添加扰动;{input/output}_list_d01.txt记录背景/分析场文件名;input.nml记录控制滤波过程的namelist;new_advance_model.csh主要完成相邻两次观测时刻的模拟过程,将WRFOUT转换为DART格式的状态向量,为下次滤波过程提供prior文件。 output 包含各个时刻的观测文件,也是将来保存同化结果的目录。 perts 扰动库文件,共100个。add_bank_perts.ncl会随机调用这里的扰动库文件来生成集合样本。 template 包含WRF、WRF-3DVAR、集合滤波的namelist 1) 启动 将$DART/models/wrf/shell_scripts 下的所有脚本复制到$WORKDIR/scripts下 将WRF/run目录下除namelist.input之外的所有文件复制到$WORKDIR/rundir/WRF_RUN/ 将WRFDA/var/run/be.dat.cv3(复制到目标路径下后记得改名为be.dat)以及da_wrfvar.exe也都复制到$WORKDIR/rundir/WRF_RUN/,如果用WRFDARANDOM生成扰动样本需要这样做,Testdata里用扰动库文件生成集合样本,可以略过这一步。 将之前编译好的DART可执行文件复制到rundir 2) 初始条件 根据scripts/param.csh脚本的描述修改变量,以适应自己机器的配置。 运行gen_retro_icbc.csh产生初边界条件 运行init_ensemble_var.csh产生集合样本,并输出为prior文件,为集合同化提供输入文件。 3) 准备观测资料 观测资料要处理成DART能够读取的格式,DART提供多种观测资料转换程序。后续有时间的话,我再把这部分内容补全。这里我们直接使用处理好的观测序列文件obs_seq.out 4) 循环同化 scripts/driver.csh是驱动循环同化的脚本,它里面调用了为滤波准备初始条件的prep_ic.csh脚本;产生滤波分析场的assimilate.csh脚本;产生集合平均的assim_advance_mem.csh脚本。注意:第一次运行集合滤波的时候,filter会读取inflation信息,一种是用fill_inflation_restart产生input_priorinf_{mean/sd}.nc,另一种是在input.nml里设置inf_initial_from_restart=false和inf_sd_initial_from_restart=false,即从input.nml读取inflation的相关信息。这里建议使用后者会方便一些,但是接的在之后的cycle里要将前两条改回true。
作为一个WRF/DART的粗略介绍,这里面没有包含资料处理和分析同化结果的内容,需要了解这部分内容的朋友可以从以上的链接中获取帮助。20190311 |