D02运行icon报错

360768761

大家好，最近我在运行CMAQ5.3.2的D02层icon和bcon脚本的时候出现了如下报错，但相对而言D01运行的时候没有出现错误。具体的意思就是：层定义和CTM垂直网格类型相同，但分辨率不同。使用VGLVS进行垂直插值(如下所示)。但是对应的垂直插值并没有给出。

     The Layer Defn and CTM vertical grid types are the same, but the resolution is different.
     Vertical interpolation using VGLVS (listed below).


forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source            
ICON_v532.exe      00000000004C60E3  Unknown               Unknown  Unknown
libpthread-2.17.s  00002AC5C177C630  Unknown               Unknown  Unknown
libc-2.17.so       00002AC5C4C1C2D0  Unknown               Unknown  Unknown
ICON_v532.exe      000000000052ECF0  Unknown               Unknown  Unknown
ICON_v532.exe      0000000000502361  Unknown               Unknown  Unknown
ICON_v532.exe      00000000005005A0  Unknown               Unknown  Unknown
ICON_v532.exe      000000000041B7EC  Unknown               Unknown  Unknown
ICON_v532.exe      0000000000417B07  Unknown               Unknown  Unknown
ICON_v532.exe      00000000004148E8  Unknown               Unknown  Unknown
ICON_v532.exe      000000000040D226  Unknown               Unknown  Unknown
ICON_v532.exe      000000000040C3E2  Unknown               Unknown  Unknown
libc-2.17.so       00002AC5C4AEF555  __libc_start_main     Unknown  Unknown
ICON_v532.exe      000000000040C2E9  Unknown               Unknown  Unknown
0.010u 0.024s 0:00.04 75.0%        0+0k 0+32io 0pf+0w
exit ( )

weilin421

360768761 发表于 2023-8-18 18:25
是的，我检查了这个mcip生成的METCRO3D文件，他里边的数量如右边终端里显示，而那个cctm运行完的output文 ...

哦我知道了 cctm有一个CGRID所有层都输出的参数，你得改下，要不cctm只输出一层，但是输出所有层很大。

weilin421

d02的icon不是也是用空数据跑么，你内层模拟周期的开头有初始场浓度么。

360768761

weilin421 发表于 2023-8-18 12:46
d02的icon不是也是用空数据跑么，你内层模拟周期的开头有初始场浓度么。

您好，D01层跑的时候，我用清单分配工具制作了五个排放源的数据（emis），用这个emis跑的cctm，但是对于D02而言，第一层生成这个output即可作为D02的输入数据了嘛不是。那您说的内层模拟周期的初始场浓度就对应我这个D01层跑的这个output数据呀，是有数的。

360768761

360768761 发表于 2023-8-18 15:38
您好，D01层跑的时候，我用清单分配工具制作了五个排放源的数据（emis），用这个emis跑的cctm，但是对于D ...

有数据，只是数据不全，没有海洋的、生物的、沙尘的数据，但人为源数据是有的，所以按理说应该不是空数据。

weilin421

bcon可以用cctm d01的结果截取出每个小时的边界条件。
icon初始场d01不是也用的本底值算的么，就是d02的一个初始场和直接初始化一个背景初始场结果差不太多吧。所以你的问题是bcon切cctm d01结果有问题么。

weilin421

bcon我记得是需要编译两次，一个是用背景值初始化最外层场用的不需要输入外层数据，
一个是从外层场里切内层的边界场需要有外层的cctm结果，两种需要修改bcon build的脚本参数编译两次，不过看你这个报错大概是编译过了第二种程序了吧，你mcip输出的结果是wrf 1 2两层的输出结果，mcip好像可以压缩层数，基本上应该和你wrf的垂直分层一致吧，所以你1 2两层的wrf垂直分层不一样吗？不一样的话好像mcip转换格式的时候也可以调整，可以看下你wrf原始气象数据的垂直分层情况和mcip不同层转换数据时候的配置具体是怎么设置的。

360768761

weilin421 发表于 2023-8-18 16:54
bcon我记得是需要编译两次，一个是用背景值初始化最外层场用的不需要输入外层数据，
一个是从外层场里切内 ...

您所提到的bcon需要编译两次，这个我并没有编译两次，就直接用了同一个BLD文件中的exe来运行（不知道这个是不是必要的）。关于这个WRF的分层，应该是在namelist.input脚本中设置，当时我就用的默认的45层，没有更改垂直层数，但是相对而言我看了看这个mcip和cctm第一层运行完所生成的文件，里边有个参数是不一样的：NLAYS，在cctm中是1层，在mcip中是44层，不知道这有没有关系？

360768761

360768761 发表于 2023-8-18 17:09
您所提到的bcon需要编译两次，这个我并没有编译两次，就直接用了同一个BLD文件中的exe来运行（不知道这个 ...

然后我就在cctm运行脚本中找到NLAYS所在的选项，将其中的1更改为all layers,但是这并不能解决问题

360768761

weilin421 发表于 2023-8-18 16:49
bcon可以用cctm d01的结果截取出每个小时的边界条件。
icon初始场d01不是也用的本底值算的么，就是d02的一 ...

对，目前来看是第一次cctm运行的结果不能支撑第二层的icon和bcon正常运行