WRF-Chem ./mozbc加边界场时出错

Ronin.

大家好，有个问题想要请教大家，就是在运行mozbc加边界条件时出现报错缺少变量。但是却没有显示缺少的变量是什么，出现一段空格，不知道是否有人可以帮忙解答这是我的inp
&control

do_bc     = .true.
do_ic     = .true.
domain    = 1
dir_wrf              = '/home/test/wrfchem4/WRF/test/em_real/'
dir_moz = '/mnt/hgfs/WRF-Chem-Tutorial/Data/'
fn_moz  = 'camchem-20140613-20140618.nc'
moz_var_suffix = ''
def_missing_var = .true.
spc_map = 'o3 -> O3', 'n2o -> N2O', 'no -> NO',
         'no2 -> NO2', 'nh3 -> NH3', 'hno3 -> HNO3', 'hno4 -> HO2NO2',
         'n2o5 -> N2O5', 'h2o2 -> H2O2',
         'ch4 -> CH4', 'co -> CO', 'ch3ooh -> CH3OOH',
         'hcho -> CH2O', 'ch3oh -> CH3OH', 'c2h4 -> C2H4',
         'ald -> CH3CHO', 'acet -> CH3COCH3', 'mgly -> CH3COCHO',
         'pan -> PAN', 'mpan -> MPAN', 'macr -> MACR',
         'mvk -> MVK', 'c2h6 -> C2H6', 'c3h6 -> C3H6', 'c3h8 -> C3H8',
         'c2h5oh -> C2H5OH', 'c10h16 -> MTERP',
         'isopr -> ISOP', 'acetol -> HYAC', 'mek -> MEK',
         'bigene -> BIGENE', 'bigalk -> BIGALK',
         'tol -> TOLUENE', 'benzene -> BENZENE', 'xylenes -> XYLENES',
         'cres -> CRESOL', 'dms -> DMS', 'so2 -> SO2',

          'oc_a01->0.0093*pom_a1+0.7510*soa1_a2+0.0093*soa1_a1+0.7510*soa2_a2+0.0093*soa2_a1+0.7510*soa3_a2+0.0093*soa3_a1+0.7510*soa4_a2+0.0093*soa4_a1+0.7510*soa5_a2+0.0093*soa5_a1+0.1123*pom_a1+0.2376*soa1_a2+0.1123*soa1_a1+0.2376*soa2_a2+0.1123*soa2_a1+0.2376*soa3_a2+0.1123*soa3_a1+0.2376*soa4_a2+0.1123*soa4_a1+0.2376*soa5_a2+ 0.1123*soa5_a1;1.e9',
         'oc_a02->0.3835*pom_a1+0.0113*soa1_a2+0.3835*soa1_a1+0.0133*soa2_a2+0.3835*soa2_a1+0.0113*soa3_a2+0.3838*soa3_a1+0.0113*soa4_a2+0.3838*soa4_a1+0.0113*soa5_a2+ 0.3838*soa5_a1+0.3783*pom_a1+0.0001*soa1_a2+0.3783*soa1_a1+0.0001*soa2_a2+0.3783*soa2_a1+0.0001*soa3_a2+0.3783*soa3_a1+0.0001*soa4_a2+0.3783*soa4_a1+0.0001*soa5_a2+ 0.3783*soa5_a1;1.e9',
         'oc_a03->0.1077*pom_a1+0.0000*soa1_a2+0.1077*soa1_a1+0.0000*soa2_a2+0.1077*soa2_a1+0.0000*soa3_a2+0.1077*soa3_a1+0.0000*soa4_a2+0.1077*soa4_a1+0.0000*soa5_a2+ 0.1077*soa5_a1+0.0087*pom_a1+0.0000*soa1_a2+0.0087*soa1_a1+0.0000*soa2_a2+0.0087*soa2_a1+0.0000*soa3_a2+0.0087*soa3_a1+0.0000*soa4_a2+0.0087*soa4_a1+0.0000*soa5_a2+ 0.0087*soa5_a1;1.e9',
         'oc_a04->0.0002*pom_a1+0.0000*soa1_a2+0.0002*soa1_a1+0.0000*soa2_a2+0.0002*soa2_a1+0.0000*soa3_a2+0.0002*soa3_a1+0.0000*soa4_a2+0.0002*soa4_a1+0.0000*soa5_a2+ 0.0002*soa5_a1+0.0000*pom_a1+0.0000*soa1_a2+0.0000*soa1_a1+0.0000*soa2_a2+0.0000*soa2_a1+0.0000*soa3_a2+0.0000*soa3_a1+0.0000*soa4_a2+0.0000*soa4_a1+0.0000*soa5_a2+ 0.0000*soa5_a1;1.e9',
         'bc_a01->0.1216*bc_a1+0.1216*bc_a4;1.e9',
         'bc_a02->0.7618*bc_a1+0.7618*bc_a4;1.e9',
         'bc_a03->0.1164*bc_a1+0.1164*bc_a4;1.e9',
         'bc_a04->0.0002*bc_a1+0.0002*bc_a4;1.e9',
         'so4_a01->0.9886*so4_a2+0.1216*so4_a1+0.0000*so4_a3;1.e9',
         'so4_a02->0.0114*so4_a2+0.7618*so4_a1+0.0002*so4_a3;1.e9',
         'so4_a03->0.0000*so4_a2+0.1164*so4_a1+0.0995*so4_a3;1.e9',
         'so4_a04->0.0000*so4_a2+0.0002*so4_a1+0.9003*so4_a3;1.e9',
         'nh4_a01->0.1410*so4_a2+0.0033*so4_a1+0.0000*so4_a3+0.0446*so4_a2+0.0017*so4_a1+0.0000*so4_a3;1.e9',
         'nh4_a02->0.0021*so4_a2+0.0210*so4_a1+0.0000*so4_a3+0.0000*so4_a2+0.0720*so4_a1+0.0000*so4_a3;1.e9',
         'nh4_a03->0.0000*so4_a2+0.0202*so4_a1+0.0011*so4_a3+0.0000*so4_a2+0.0001*so4_a1+0.0175*so4_a3;1.e9',
         'nh4_a04->0.0000*so4_a2+0.0000*so4_a1+0.0755*so4_a3+0.0000*so4_a2+0.0000*so4_a1+0.0935*so4_a3;1.e9',
         'no3_a01->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
         'no3_a02->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
         'no3_a03->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
         'no3_a04->0.0000*so4_a2+0.0000*so4_a1+0.0000*so4_a3;1.e9',
         'na_a01->0.2954*ncl_a2+0.0037*ncl_a1+0.0000*ncl_a3+0.0935*ncl_a2+0.0442*ncl_a1+0.0000*ncl_a3;1.e9',
         'na_a02->0.0045*ncl_a2+0.1509*ncl_a1+0.0000*ncl_a3+0.0000*ncl_a2+0.1488*ncl_a1+0.0000*ncl_a3;1.e9',
         'na_a03->0.0000*ncl_a2+0.0424*ncl_a1+0.0024*ncl_a3+0.0000*ncl_a2+0.0034*ncl_a1+0.0367*ncl_a3;1.e9',
         'na_a04->0.0000*ncl_a2+0.0000*ncl_a1+0.1582*ncl_a3+0.0000*ncl_a2+0.0000*ncl_a1+0.1960*ncl_a3;1.e9',
         'cl_a01->0.4555*ncl_a2+0.0056*ncl_a1+0.0000*ncl_a3+0.1441*ncl_a2+0.0681*ncl_a1+0.0000*ncl_a3;1.e9',
         'cl_a02->0.0068*ncl_a2+0.2326*ncl_a1+0.0000*ncl_a3+0.0000*ncl_a2+0.2295*ncl_a1+0.0000*ncl_a3;1.e9',
         'cl_a03->0.0000*ncl_a2+0.0654*ncl_a1+0.0037*ncl_a3+0.0000*ncl_a2+0.0055*ncl_a1+0.0567*ncl_a3;1.e9',
         'cl_a04->0.0000*ncl_a2+0.0001*ncl_a1+0.2439*ncl_a3+0.0000*ncl_a2+0.0000*ncl_a1+0.3023*ncl_a3;1.e9',
         'oin_a01->0.7510*dst_a2+0.0093*dst_a1+0.0000*dst_a3+0.2376*dst_a2+0.1123*dst_a1+0.0000*dst_a3;1.e9',
         'oin_a02->0.0113*dst_a2+0.3835*dst_a1+0.0000*dst_a3+0.0001*dst_a2+0.3783*dst_a1+0.0002*dst_a3;1.e9',
         'oin_a03->0.0000*dst_a2+0.1077*dst_a1+0.0061*dst_a3+0.0000*dst_a2+0.0087*dst_a1+0.0934*dst_a3;1.e9',
         'oin_a04->0.0000*dst_a2+0.0002*dst_a1+0.4020*dst_a3+0.0000*dst_a2+0.0000*dst_a1+0.4983*dst_a3;1.e9',
         'num_a01->0.9502*num_a2+0.2509*num_a1+0.0000*num_a3+0.0494*num_a2+0.4626*num_a1+0.0000*num_a3;1.0',
         'num_a02->0.0004*num_a2+0.2470*num_a1+0.0007*num_a3+0.0000*num_a2+0.0377*num_a1+0.0232*num_a3;1.0',
         'num_a03->0.0000*num_a2+0.0016*num_a1+0.1886*num_a3+0.0000*num_a2+0.0000*num_a1+0.4372*num_a3;1.0',
         'num_a04->0.0000*num_a2+0.0000*num_a1+0.2935*num_a3+0.0000*num_a2+0.0000*num_a1+0.0566*num_a3;1.0',
/
然后下面是报错的内容

qqqin

可能是oc_a01那里有问题吧。可以先把oc_a01后面的全部删掉，试下能不能成功

Ronin.

qqqin 发表于 2022-5-28 16:38
可能是oc_a01那里有问题吧。可以先把oc_a01后面的全部删掉，试下能不能成功

删掉了还是不行，重新报错mapper: mstring = oc_a01       
                                      mapper: oc_a01              is invalid map

qqqin

Ronin. 发表于 2022-5-30 09:14
删掉了还是不行，重新报错mapper: mstring = oc_a01       
                                      mapper: o ...

spc_map = 'o3 -> O3', 'n2o -> N2O', 'no -> NO',
         'no2 -> NO2', 'nh3 -> NH3', 'hno3 -> HNO3', 'hno4 -> HO2NO2',
         'n2o5 -> N2O5', 'h2o2 -> H2O2',
         'ch4 -> CH4', 'co -> CO', 'ch3ooh -> CH3OOH',
         'hcho -> CH2O', 'ch3oh -> CH3OH', 'c2h4 -> C2H4',
         'ald -> CH3CHO', 'acet -> CH3COCH3', 'mgly -> CH3COCHO',
         'pan -> PAN', 'mpan -> MPAN', 'macr -> MACR',
         'mvk -> MVK', 'c2h6 -> C2H6', 'c3h6 -> C3H6', 'c3h8 -> C3H8',
         'c2h5oh -> C2H5OH', 'c10h16 -> MTERP',
         'isopr -> ISOP', 'acetol -> HYAC', 'mek -> MEK',
         'bigene -> BIGENE', 'bigalk -> BIGALK',
         'tol -> TOLUENE', 'benzene -> BENZENE', 'xylenes -> XYLENES',
         'cres -> CRESOL', 'dms -> DMS', 'so2 -> SO2'
删掉之后是这样？


另外可以ncdump看一下fn_moz和wrfinput或者wrfbdy文件里的变量，做下对比

Ronin.

qqqin 发表于 2022-5-30 09:19
spc_map = 'o3 -> O3', 'n2o -> N2O', 'no -> NO',
         'no2 -> NO2', 'nh3 -> NH3', 'hno3 -> HNO ...

解决了，把oc_a01和oc_ao2的参数改了下就跑通了

好小猴

您好，在吗，请问可以加下您的好友吗

为人恋城

Ronin. 发表于 2022-5-31 10:04
解决了，把oc_a01和oc_ao2的参数改了下就跑通了

请问具体是怎么改的oc_a01和oc_ao2的参数呀？

Ronin.

为人恋城 发表于 2023-3-17 13:28
请问具体是怎么改的oc_a01和oc_ao2的参数呀？

我现在的ao1-ao4是这样的
'oc_a01->0.1216*pom_a1+0.9886*soa1_a2+0.122*soa1_a1+0.98*soa2_a2+0.1216*soa2_a1+0.98*soa3_a2+0.1216*soa3_a1+0.98*soa4_a2+0.1216*soa4_a1+0.98*soa5_a2+0.1216*soa5_a1;1.e9',
         'oc_a02->0.7618*pom_a1+0.0114*soa1_a2+0.762*soa1_a1+0.0114*soa2_a2+0.76*soa2_a1+0.0114*soa3_a2+0.76*soa3_a1+0.0114*soa4_a2+0.76*soa4_a1+0.0114*soa5_a2+0.76*soa5_a1;1.e9',
         'oc_a03->0.1164*pom_a1+0.0000*soa1_a2+0.1164*soa1_a1+0.0000*soa2_a2+0.1164*soa2_a1+0.0000*soa3_a2+0.1164*soa3_a1+0.1164*soa4_a1+0.1164*soa5_a1;1.e9',
         'oc_a04->0.0002*pom_a1+0.0000*soa1_a2+0.0002*soa1_a1+0.0000*soa2_a2+0.0002*soa2_a1+0.0000*soa3_a2+0.0002*soa3_a1+0.0002*soa4_a1+0.0002*soa5_a1;1.e9',

为人恋城

Ronin. 发表于 2023-3-21 09:45
我现在的ao1-ao4是这样的
'oc_a01->0.1216*pom_a1+0.9886*soa1_a2+0.122*soa1_a1+0.98*soa2_a2+0.1216*s ...

非常感谢！

叶舟

oc_a0*所在的几行太长了，超出了Fortran的字符上限。