从WRF中删除旋转帧的技术

Question

我是WRF V3.5的用户，主要用于气象模拟。
有时，我会运行模型很长一段时间．

我的尝试

我把长周期(e。G 1年)分成几个短的。

例如:

## WPS名称列表。wps start_date = '2014-01-01_06:00:00' end_date = '2014-01-31_23:00:00'在。csh文件中设置WRF运行周期。Set year = "2014" foreach strtime (010106010512 010500011012 011000011512 011500012012 012000012512 012500013012 013000020100) Set smon = ' echo ${strtime}|cut -c1-2 ' Set sday = ' echo ${strtime}|cut -c3-4 ' Set SHR = ' echo ${strtime}|cut -c5-6 ' Set emon = ' echo ${strtime}|cut -c9-10 ' Set ehr = ' echo ${strtime}|cut -c11-12 ' &time_control run_days = 0, run_hours = 132, run_minutes = 0, run_seconds = 0, start_year = ${year}，${year}，Start_month = ${smon}，${smon}， start_day = ${sday}，${sday}， start_hour = ${shr}，${shr}， start_minute = 00,00, start_second = 00,00, end_year = ${year}，${year}， end_month = ${emon}，${emon}， end_day = ${eday}，${ehr}， end_hour = ${ehr}， end_minute = 00,00, end_second = 00, interval_seconds = 21600

在这短短的时间里，我将前12小时作为旋转的时间(除了开始的时间)。

这是我从网上剪下来的时序建模设置的示意图。

http://i12.tietuku.cn/afd906ddd823ddd6.png

结果

它生成8个.nc文件，每个文件包含132帧数据。
现在，我使用netCDF运营商(NCO)来处理输出文件。

1. 使用ncks要分别删除每个文件的前12小时帧。

2. 使用ncrcat将所有8个.nc文件连接在一起，得到一个文件，代表整个模拟周期。

问题

我曾经用过这种方法。但是当文件数量很大时，它实际上是一种转储方法。

有人提供您处理旋转数据删除的方法吗?

我的目标

• 该模型只能生成一个包含所有仿真周期的文件。
• 开机时间可以被上一个文件的帧自动覆盖。

更新:希望你分享你的方法!

Answer 1

有一个特点- d选项，用于调用ncks多余的。但是，从NCO版本4.2.1(2012年8月)开始就可以使用了。

ncrcat -d时间，13，，132,120 wrfout*摘要。Nc # -d dim，[min][，[max][，[stride][，[subcycle]]]] .

这一特性深深隐藏在在线手册中。你可以在这里找到它:nco.sourceforge.net/nco.html #次循环．

解释:

• 13:我们考虑的第一步
• [空]:直到最后
• 132:我们的循环长度为132步
• 120:我们取每个周期的前120个值(从13号开始)

如果我们写0132年- d, 13日，等于0132年- d, 13日,1，我们只能得到从第13个开始的第132个值。

我不处理WRF输出，但我们使用化学传输模型CMAQ工作。我们得到每小时解析浓度输出文件的每一天。这些文件有25个时间步长(当前的0点钟和第二天)，我们需要删除一个)点钟时间步长并追加文件。为此，我使用ncrcat -d TSTEP,0, 25,24 [infiles] [outfile]．

最直接的方法是为此目的编写一个Fortran程序。但是，我也更喜欢军士。也许你可以试试cdo。code.zmaw.de /项目/ cdo)作为替代。然而，我使用cdo的解决方案比NCO解决方案更笨重。

Answer 2

这是我目前的方法。

1. 所有wrfout-files都位于同一路径上。

2. 使用pynco在python循环中拆分旋转数据帧

   import glob nco = nco () for fn in glob.iglob('*.nc'): filename, extname = os.path.splitext(fn) output_fn = filename[11:22] + str(int(filename[22:24]) + 12) + extname nco。nks (input=fn, output=output_fn, dimension='Time,13 ')

输入如下:2014 - 04 - 04 - _12.nc(内装120帧)
输出如下所示:wrfout_2014-04-04_00:00:00.nc(内装133帧)