time-step设置造成后处理不成功的问题

Mid_Farmer

在./arwpost.exe时，出现了如此帖一样的问题（http://bbs.06climate.com/forum.php?mod=viewthread&tid=5878），该帖的回复是说time-step设为162，导致在后处理过程中不能得到相应的整点输出（3600/162=22.2222）。
我的time-step也是设成的162，因为guide建议为dx的6倍，我的dx=27000，且history_interval=180,120,30，但我在后处理时遇到了三个情况：
情况1：处理区域d03时，后处理interval_seconds=1800，能够得出ctl和dat，且出图正常。
情况2：处理区域d02时，后处理interval_seconds=1800，就会得到和链接帖一样的报错。
情况3：处理区域d02时，后处理interval_seconds=3600，能够得出ctl和dat，可是出图从t=10到49的每个时刻的结果都是一模一样的。
想问几个问题：
问题1：如果time-step=162导致后处理不能整点输出，那么为什么情况1能够成功问题？
问题2：情况3为什么不报错，但结果错误？
问题3：是不是time-step最好设成180这种能被3600这种整除的数，那么dx的6倍这种建议是否需要遵守？若遵守，则dx的设置是不是局限性变大了？
万分感谢！！！

lsd1984

6倍不是必须，这是一般性的建议。时间步长在6倍时候计算不容易发散，但是计算结果不保证正确。
一般需要根据实测结果对时间步长进行微调。或者改小一系列的时间步长，做收敛性分析。
我的计算中一般3倍附近效果好点。时间步长尽量取整，这样方便后处理操作。

Mid_Farmer

lsd1984 发表于 2014-12-13 17:26
6倍不是必须，这是一般性的建议。时间步长在6倍时候计算不容易发散，但是计算结果不保证正确。
一般需要根 ...

谢谢你啦 后面几天我把时间步长改几个方案试试 看看都有什么区别
那我这次遇到到的情况1、2应该怎么解释呢？

lsd1984

wrf在输出的时候貌似会就近选择相应的结果输出（忘记了）。但是你的后处理程序可能不是这么来的。我用ncl，是指定时间间隔输出处理结果的，从未有问题。

Mid_Farmer

lsd1984 发表于 2014-12-14 10:31
wrf在输出的时候貌似会就近选择相应的结果输出（忘记了）。但是你的后处理程序可能不是这么来的。我用ncl， ...

知道了 非常感谢你！

Mid_Farmer

lsd1984 发表于 2014-12-14 10:31
wrf在输出的时候貌似会就近选择相应的结果输出（忘记了）。但是你的后处理程序可能不是这么来的。我用ncl， ...

你好 我今天遇到了一个另外的问题 我dx=30000，step_time=180 但wrf.exe在运行5分钟后会自动断掉 在rsl.out.0000里没有发现error 不过存在5 poinits exceeded cfl=2 in domain 2 at time 4.20000这样的句子 官网上只回答"when this happens,often reduce time_step can help".我把time_step改为90后，成功运行了
请问time_step在计算过程中起个什么作用，对结果有多大影响？提前感谢你，给你添麻烦了···

lsd1984

Mid_Farmer 发表于 2014-12-15 20:53
你好 我今天遇到了一个另外的问题 我dx=30000，step_time=180 但wrf.exe在运行5分钟后会自动断掉 在rsl.o ...

在计算过程中，时间步长太大的时候，在变化较大的时候容易发生数值发散。常规的方法就是减小时间步长。你如果有充足时间，可以做一点简单的收敛性分析，比如时间步长从2倍到6倍算一下，输出一些典型的点的气压、风速、气温之类的数据查看收敛性，这样可以心中有数一点。另外莫名其妙的死掉的时候，请看自己的wps和wrf里边的namelist关于网格的部分是否一致，有时候没有对应改过来，也容易发散。

Mid_Farmer

lsd1984 发表于 2014-12-16 11:59
在计算过程中，时间步长太大的时候，在变化较大的时候容易发生数值发散。常规的方法就是减小时间步长。你 ...

谢谢你了 非常感谢！！！

yujie

感谢楼主的帖子，知道不少东西，还想问一下，time step可不可以取小数啊？比如40.5之类的？

fiona0923

帮到我了！！感谢楼主