爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索

WRF/DART同化流程

[复制链接]

新浪微博达人勋

发表于 2019-4-29 15:30:23 | 显示全部楼层
今朝几朵白云生 发表于 2019-4-28 22:23
我看不出你出现错误的原因。如果你已经按照WRF/DART tutorial的说明步骤,并且将input。nml拷贝到了$DART ...

好的 我试试 谢谢啦!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-4-30 14:06:26 | 显示全部楼层
嗯嗯嗯嗯,有点详细,有点明白,有点好用,谢谢!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-7-30 22:06:40 | 显示全部楼层
今朝几朵白云生 发表于 2019-4-28 22:23
我看不出你出现错误的原因。如果你已经按照WRF/DART tutorial的说明步骤,并且将input.nml拷贝到了$DART/ ...

您好~能加下您的qq或者微信吗?有一两个小问题想向您
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-7-31 12:30:07 | 显示全部楼层
您好~我有两个问题,第一个是在第二步中,运行init_ensemble_var.csh产生集合样本时,会调用first_advance.csh这个脚本,其中
set datea     = ${1}
set emember   = ${2}
set paramfile = ${3}
这个emember是什么呢?是集合个数吗

第二个问题,在rundir文件中,new_advance_model.csh脚本里 第一句set process = $1 这个应该写什么呢 process代表的什么呢?


实在是感谢!困扰好些天了,感谢感谢!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-7-31 21:17:38 | 显示全部楼层
感谢楼主,非常有用。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-8-6 10:23:54 | 显示全部楼层
halakaka 发表于 2019-7-31 12:30
您好~我有两个问题,第一个是在第二步中,运行init_ensemble_var.csh产生集合样本时,会调用first_advance. ...

您好,init_ensemble_var.csh调用first_advance.csh是通过这一条命令实现的:${SHELL_SCRIPTS_DIR}/first_advance.csh $initial_date $n ${SHELL_SCRIPTS_DIR}/$paramfile
由此可知init_ensemble_var.csh中的变量emember就是n,并且n在循环体里面表示第n个集合样本。所以emember和集合个数有一定区别。
第二个问题:
[     scripts]$ grep new_advance_model.csh *.csh
assim_advance.csh:${RUN_DIR}/new_advance_model.csh ${emember} $domains filter_control${icnum}
first_advance.csh:${RUN_DIR}/new_advance_model.csh ${emember} 1 filter_control${icnum}
param.csh:   # IMPORTANT - SHOULD MATCH those set in new_advance_model.csh  #
可以看到调用new_advance_model.csh的脚本有两个assim_advance.csh和first_advance.csh,并且第一个参数都是emember,所以process就是这两个脚本中的emember,即第n个集合样本。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-8-8 09:42:54 | 显示全部楼层
今朝几朵白云生 发表于 2019-8-6 10:23
您好,init_ensemble_var.csh调用first_advance.csh是通过这一条命令实现的:${SHELL_SCRIPTS_DIR}/first ...

十分感谢您耐心的回复~

对于第一个问题,变量emember就是第n个集合样本,那么是不是得自己做个循环:
set emember   = 1
while ( $emember <= $NUM_ENS )
@ emember++

那对于第二个问题,process就是这两个脚本中的emember,即第n个集合样本。这个意思是不是process=$emember

在进行前两步骤,执行gen_retro_icbc.csh、init_ensemble_var.csh时,其他的主要脚本里path需要改动吗(如driver.csh、assim_advance.csh等等)因为现在还没有生成prior_d01.${icnum}文件,所以很困惑是不是这些没有调用到的脚本也需要先更改了。
在【scripts】grep prior_d01  *csh 后,并没有看到哪一步是生成 prior_d01.0001等等先验集合文件的,所以好困惑~
assim_advance.csh:  echo prior_d01.${icnum}     >> ${RUN_DIR}/filter_control${icnum}
driver.csh:   ${MOVE} ${RUN_DIR}/prior_d01.${ensstring}       ${OUTPUT_DIR}/${datea}/PRIORS/.
first_advance.csh:echo prior_d01.${icnum}     >> ${RUN_DIR}/filter_control${icnum}
first_advance.csh:mv $RUN_DIR/prior_d01.${icnum} ${OUTPUT_DIR}/${datea}/PRIORS/prior_d01.${icnum}

十分感谢您~!!!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-8-8 12:32:42 | 显示全部楼层
本帖最后由 今朝几朵白云生 于 2019-8-8 12:35 编辑
halakaka 发表于 2019-8-8 09:42
十分感谢您耐心的回复~

对于第一个问题,变量emember就是第n个集合样本,那么是不是得自己做个循环:
...

是要生成每个集合样本的扰动文件,但这个循环已经在脚本init_ensemble_var.csh里写好,不需要自己写。
我记得path这些变量都是记录在param.csh文件里的,主要的脚本如driver.csh都是从param.csh里读取变量信息。
每次同化之前都要先生成集合扰动样本prior_d01.000*,但是在哪一步生成我记不大清楚了。你可以在WRF/DART教程https://www.image.ucar.edu/wrfdart/tutorial/里找到一些线索。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-8-8 12:55:17 | 显示全部楼层
今朝几朵白云生 发表于 2019-8-8 12:32
是要生成每个集合样本的扰动文件,但这个循环已经在脚本init_ensemble_var.csh里写好,不需要自己写。
...

好的好的 我再仔细看看脚本 十分感谢您啊!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-8-8 13:34:50 | 显示全部楼层
同围观dart~
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表