一、问题产生
在做一个月的浓度预测时，由于GFS只能预报16天左右，因此需要使用CFS气象数据来驱动WRF模型，但CFS在WRF4.x版本有问题，因此重新装了WRF3.9.1版本的WRF。
而我这里的CMAQ ICON会去取前一天的CCTM CONC制作新的ICON文件，由于MCIP出来的结果不同，导致了CMAQ最终报错了：
二、解决过程
从报错图中可以看出，"The COORD.EXT and CTM vertical grid types are different"表明coord type不匹配导致了错误的产生。WRF3.9.1的VGTYP=7，而WRF4.3.1的则为-9999。（WRF3.9.1是 terrain-following dry hydrostatic pressure，而WRF4.3.1 是 hybrid sigma-pressure vertical coordinate）
最好的方法就是，不使用第一天的化学场作为初始场，直接重新生成初始场即可。但这里想找到其他的直接修改全局变量的方法来解决。
于是使用ncatted(用法：https://linux.die.net/man/1/ncatted)进行VGTYP和CONC(VGLVLS)的替换，将CWRF模拟时段的前一天的METCRO3D文件中的VGTYP改为7：
（1）修改VGTYP
ncatted -O -a VGTYP,global,m,l,7 METCRO3D_d01.2022113016.nc_wrf-431 METCRO3D_d01.2022113016.nc
（2）修改VGLVLS
ncatted -O -a VGTYP,global,m,l,7 -a VGLVLS,global,m,f,"1.0, 0.9930, 0.9830, 0.970, 0.9540, 0.9340, 0.9090, 0.880, 0.84066630, 0.80133270, 0.7619990, 0.72266530, 0.65257550, 0.58773610, 0.52781920, 0.4725140, 0.42152620, 0.37457750, 0.33140440, 0.29175790, 0.25540260, 0.22211620, 0.19168880, 0.16392220, 0.13862970, 0.11563510, 0.095250160, 0.077334810, 0.061589830, 0.047752310, 0.035591150, 0.024903280, 0.01551020, 0.0072550590, 0.0" METCRO3D_d01.2022113016.nc_wrf-431 METCRO3D_d01.2022113016.nc
修改完成后，重新运行，CMAQ就能正常运行了：
ncatted命令的用法后面熟练了再做介绍：
https://linux.die.net/man/1/ncatted
本人独自运营了微信公众号，用于分享个人学习及工作生活趣事，大家可以关注一波。（微信搜索“微思研”）
