D02运行icon报错

weilin421

cctm的层数设置就是mcip 3d的输出文件里的分层结构，如果不一致，检查mcip就行了，看mcip为什么不一致，是参数设的有问题还是wrfout输入给mcip的结果就不一样，我记得mcip有layer collapse，30层可以压成15层，好像也可以设eta值那么弄，垂直层分辨率不能不一样，比如同样的高度，第一层30层，第二层90层这种那肯定不行。

weilin421

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

这个少1层好像是正常的吧，wrf有些是格点数据，cctm里面转换成网格数据本身就会少1，就像5个点中间只能有四个边这个意思。
mcip本身好像也会切掉最外围一部分格子，他自己的说法是最外围一圈气象场数据不太可信就给切了。这个厚度在mcip好像可以设置nthick那块。
你看看你mcip输出的数据是不是垂直层数就不一样。

weilin421

360768761 发表于 2023-8-18 17:13
然后我就在cctm运行脚本中找到NLAYS所在的选项，将其中的1更改为all layers,但是这并不能解决问题

setenv NLAYS_PHOTDIAG 你说的是这个么，似乎和这个没关系
这个只是数据计算过程中的一个方案，你那个报错是输入数据本身有问题。两层垂直结构不一致。
#> Number of layers for PHOTDIAG2 and PHOTDIAG3 from
#>     Layer 1 to NLAYS_PHOTDIAG  [ default: all layers ]

360768761

weilin421 发表于 2023-8-18 17:50
cctm的层数设置就是mcip 3d的输出文件里的分层结构，如果不一致，检查mcip就行了，看mcip为什么不一致，是 ...

是的，我检查了这个mcip生成的METCRO3D文件，他里边的数量如右边终端里显示，而那个cctm运行完的output文件里的参数如左边终端所示，我对比着就是那个NLAYS不一致，所以我去修改了那个cctm运行脚本，您的意思是cctm如此生成的话是和mcip第一层的运行有关系吗。我第一次跑到这有很多不懂的，望您见谅。

360768761

weilin421 发表于 2023-8-18 17:57
setenv NLAYS_PHOTDIAG 你说的是这个么，似乎和这个没关系
这个只是数据计算过程中的一个方案，你那个报 ...

是的，是这个，我按照注释内容改了后发现了新的错误，应该是内存不支持。

weilin421

setenv CONC_BLEV_ELEV " 1 1"  # CONC file layer range; comment to write all layers to CONC
这个参数，不做嵌套的就可以不设成ALL
类似的还有ACONC的小时均值数据，根据需要调整输出数据的空间物种范围控制文件大小。
BCON和ICON用的是CONC的，人看结果一般用ACONC的。

360768761

weilin421 发表于 2023-8-18 18:43
setenv CONC_BLEV_ELEV " 1 1"  # CONC file layer range; comment to write all layers to CONC
这个参数 ...

因为那个icon和bcon用到的是CTM_CONC所以这里就对应先把这个地方的1 1 改成ALL，如果后续的用到别的参数再改别的嘛，all的意思是对应生成所有垂直层数据嘛。

360768761

360768761 发表于 2023-8-18 18:55
因为那个icon和bcon用到的是CTM_CONC所以这里就对应先把这个地方的1 1 改成ALL，如果后续的用到别的参数 ...

对应修改的话就是嵌套几层然后改成1 X嘛，比如说我有三层嵌套域的话就是设置为1 3，然后跑第二层的时候改成2 3，第三层即为3 3.还是说应该是直接和上边对齐改为ALL呀

weilin421

360768761 发表于 2023-8-18 18:55
因为那个icon和bcon用到的是CTM_CONC所以这里就对应先把这个地方的1 1 改成ALL，如果后续的用到别的参数 ...

好像不是改all，是把这句注了，设all未必能被他识别意思。
如果是最内层的数据，这个压成1层就行了，反正也不看conc的，conc是时刻数据就是整点那一个时刻，aconc是时均数据。0代表0点到0:59吧好像，我记得aconc一天好像是0到23 24个时间刻度。

weilin421

cctm一次就跑一层的，和嵌套几层没关系。
他那个1 1指的是从1到也就是最下面一层，全部范围就是mcip里所有层。
如果是最内层不需要再给别人提供边界条件了，那conc就可以设1 1。