立即注册 登录
气象家园 返回首页

meteowang的个人空间 https://bbs.06climate.com/?131856 [收藏] [复制] [分享] [RSS]

日志

./mozbc < MZ-MOSAIC-4BINS.inp > mozbc.outmz 报错问题

已有 44 次阅读2022-11-15 21:38

使用MOZART-MOSAIC-4bin.inp需要将chem_opt设置为201
报错为:
chk_moz_vars: could not find OC1 in ./jjj0000.nc
NetCDF: Variable not found

首先ncdump -c jjj0000.nc > variable.txt,查看jjj0000.nc文件中有无OC1变量;并用同样的方法查看wrfinput中的变量,
发现两个文件中都没有OC1变量。

使用的inp文件如下:
spc_map =
'o3 -> O3',
'no -> NO',
'no2 -> NO2',
'no3 -> NO3',
'nh3 -> NH3',
'hno3 -> HNO3',
'hno4 -> HO2NO2',
'n2o5 -> N2O5',
'ho -> OH',
'ho2 -> HO2',
'h2o2 -> H2O2',
'ch4 -> CH4',
'co -> CO',
'ch3o2 -> CH3O2',
'ch3ooh -> CH3OOH',
'hcho -> CH2O',
'ch3oh -> CH3OH',
'c2h4 -> C2H4',
'ald -> CH3CHO',
'ch3cooh -> CH3COOH',
'acet -> CH3COCH3',
'mgly -> CH3COCHO',
'paa -> CH3COOOH',
'pan -> PAN',
'mpan -> MPAN',
'macr -> MACR',
'mvk -> MVK',
'c2h6 -> C2H6',
'c3h6 -> C3H6',
'c3h8 -> C3H8',
'c2h5oh -> C2H5OH',
'onit -> ONIT',
'onitr -> ONITR',
'isopr -> ISOP',
'isopn -> ISOPNO3',
'acetol -> HYAC',
'glyald -> GLYALD',
'hydrald -> HYDRALD',
'mek -> MEK',
'bigene -> BIGENE',
'open -> BIGALD',
'bigalk -> BIGALK',
'tol -> TOLUENE',
'cres -> CRESOL', 'dms -> DMS', 'so2 -> SO2', 'sulf -> SO4',
'oc_a01->9.60024e-2*OC1+9.60024e-2*OC2+1.152026*SOA;1.e9'
'oc_a02->.254356*OC1+.254356*OC2+3.05226*SOA;1.e9',
'oc_a03->6.00723e-2*OC1+6.00723e-2*OC2+.720865*SOA;1.e9',
'oc_a04->.00104*OC1+.00104*OC2+.012481*SOA;1.e9',
'bc_a01->6.5325E-01*CB1+6.5325E-01*CB2;0.413793103e9',
'bc_a02->2.5458E-01*CB1+2.5458E-01*CB2;0.413793103e9',
'bc_a03->4.0521E-03*CB1+4.0521E-03*CB2;0.413793103e9',
'bc_a04->1.6794E-06*CB1+1.6794E-06*CB2;0.413793103e9',
'so4_a01->2.4433e-2*SO4;3.31423e9',
'so4_a02->.47111*SO4;3.31423e9',
'so4_a03->.47857*SO4;3.31423e9',
'so4_a04->2.5727e-2*SO4;3.31423e9',
'no3_a01->2.4433e-2*NH4NO3;2.76186e9',
'no3_a02->.47111*NH4NO3;2.76186e9',
'no3_a03->.47857*NH4NO3;2.76186e9',
'no3_a04->2.5727e-2*NH4NO3;2.76186e9',
'nh4_a01->2.4433e-2*NH4;6.22876e8',
'nh4_a02->.47111*NH4;6.22876e8',
'nh4_a03->.47857*NH4;6.22876e8',
'nh4_a04->2.5727e-2*NH4;6.22876e8',
'na_a01->.2*SA1;1.e9',
'na_a02->.4*SA1;1.e9',
'na_a03->.2*SA1+.8*SA2;1.e9',
'na_a04->.8*SA3;1.e9',
'cl_a01->.3*SA1;1.e9',
'cl_a02->.6*SA1;1.e9',
'cl_a03->.3*SA1+1.2*SA2;1.e9',
'cl_a04->1.2*SA3;1.e9',
'oin_a01->0.5*[DUST1];1.e9',
'oin_a02->0.5*[DUST1]+0.24*[DUST2];1.e9',
'oin_a03->0.76*[DUST2]+1.0*[DUST3];1.e9',
'oin_a04->1.0*[DUST4];1.e9',
'num_a01-
>1.73e17*OC1+1.73e17*OC2+1.73e17*SOA+5.64e17*CB1+5.64e17*CB2+7.67e16*SO4+6.39
e16*NH4NO3+1.44e16*NH4+3.22e17*SA1+4.83e17*SA1+0.5*3.93e17*[DUST1];1.0',
'num_a02-
>1.71e16*OC1+1.71e16*OC2+1.71e16*SOA+9.91e15*CB1+9.91e15*CB2+6.06e16*SO4+5.05
e16*NH4NO3+1.14e16*NH4+2.90e16*SA1+4.36e16*SA1+0.5*1.17e16*[DUST1]+0.24*1.17e16
*[DUST2];1.0',
'num_a03-
>1.13e14*OC1+1.13e14*OC2+1.13e14*SOA+1.80e12*CB1+1.80e12*CB2+2.44e15*SO4+2.03
e15*NH4NO3+4.69e14*NH4+4.12e14*SA1+1.64e15*SA2+0.76*9.55e13*[DUST2]+9.55e13*[D
UST3];1.0',
'num_a04-
>4.09e10*OC1+4.09e10*OC2+4.09e10*SOA+3.42e07*CB1+3.42e07*CB2+3.25e12*SO4+2.71
e12*NH4NO3+6.12e11*NH4+7.60e13*SA3+1.49e12*[DUST4];1.0',
/

注:在运行中如果出现找不到变量的情况,可以直接删掉这一变量。
但是OC变量是与模拟有关的,所以这里我没有删掉,选择了其他的计算方法进行替代,修改后的inp文件为:
spc_map =
'o3 -> O3',
'no -> NO',
'no2 -> NO2',
'no3 -> NO3',
'nh3 -> NH3',
'hno3 -> HNO3',
'hno4 -> HO2NO2',
'n2o5 -> N2O5',
'ho -> OH',
'ho2 -> HO2',
'h2o2 -> H2O2',
'ch4 -> CH4',
'co -> CO',
! 'ch3o2 -> CH3O2',
! 'ch3ooh -> CH3OOH',
'hcho -> CH2O',
'ch3oh -> CH3OH',
'c2h4 -> C2H4',
'ald -> CH3CHO',
'ch3cooh -> CH3COOH',
'acet -> CH3COCH3',
'mgly -> CH3COCHO',
! 'paa -> CH3COOOH',
'pan -> PAN',
'mpan -> MPAN',
'macr -> MACR',
'mvk -> MVK',
'c2h6 -> C2H6',
'c3h6 -> C3H6',
'c3h8 -> C3H8',
'c2h5oh -> C2H5OH',
! 'onit -> ONIT',
'onitr -> ONITR',
'isopr -> ISOP',
! 'isopn -> ISOPNO3',
'acetol -> HYAC',
! 'glyald -> GLYALD',
! 'hydrald -> HYDRALD',
'mek -> MEK',
'bigene -> BIGENE',
'open -> BIGALD',
'bigalk -> BIGALK',
'tol -> TOLUENE',
'cres -> CRESOL',
'dms -> DMS',
'so2 -> SO2',
! 'sulf -> SO4',
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
'oc_a01->0.1216*pom_a1+0.9886*soa1_a2+0.1216*soa1_a1+0.9886*soa2_a2+0.1216*soa2_a1+0.9886*soa3_a2+0.12
16*soa3_a1+0.9886*soa4_a2+0.1216*soa4_a1+0.9886*soa5_a2+0.1216*soa5_a1;1.e9',
'oc_a02->0.7618*pom_a1+0.0114*soa1_a2+0.7618*soa1_a1+0.0114*soa2_a2+0.7618*soa2_a1+0.0114*soa3_a2+0.76
18*soa3_a1+0.0114*soa4_a2+0.7618*soa4_a1+0.0114*soa5_a2+0.7618*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.11
64*soa3_a1+0.0000*soa4_a2+0.1164*soa4_a1+0.0000*soa5_a2+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.00
02*soa3_a1+0.0000*soa4_a2+0.0002*soa4_a1+0.0000*soa5_a2+0.0002*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.1856*so4_a2+0.0050*so4_a1+0.0000*so4_a3;1.e9',
'nh4_a02->0.0021*so4_a2+0.0930*so4_a1+0.0000*so4_a3;1.e9',
'nh4_a03->0.0000*so4_a2+0.0203*so4_a1+0.0186*so4_a3;1.e9',
'nh4_a04->0.0000*so4_a2+0.0000*so4_a1+0.1690*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.3889*ncl_a2+0.0479*ncl_a1+0.0000*ncl_a3;1.e9',
'na_a02->0.0045*ncl_a2+0.2997*ncl_a1+0.0000*ncl_a3;1.e9',
'na_a03->0.0000*ncl_a2+0.0458*ncl_a1+0.0391*ncl_a3;1.e9',
'na_a04->0.0000*ncl_a2+0.0000*ncl_a1+0.3542*ncl_a3;1.e9',
'cl_a01->0.5996*ncl_a2+0.0737*ncl_a1+0.0000*ncl_a3;1.e9',
'cl_a02->0.0068*ncl_a2+0.4621*ncl_a1+0.0000*ncl_a3;1.e9',
'cl_a03->0.0000*ncl_a2+0.0709*ncl_a1+0.0604*ncl_a3;1.e9',
'cl_a04->0.0000*ncl_a2+0.0001*ncl_a1+0.5462*ncl_a3;1.e9',
'oin_a01->0.9886*dst_a2+0.1216*dst_a1+0.0000*dst_a3;1.e9',
'oin_a02->0.0114*dst_a2+0.7618*dst_a1+0.0002*dst_a3;1.e9',
'oin_a03->0.0000*dst_a2+0.1164*dst_a1+0.0995*dst_a3;1.e9',
'oin_a04->0.0000*dst_a2+0.0002*dst_a1+0.9003*dst_a3;1.e9',
'num_a01->0.9996*num_a2+0.6135*num_a1+0.0000*num_a3;1.0',
'num_a02->0.0004*num_a2+0.2847*num_a1+0.0239*num_a3;1.0',
'num_a03->0.0000*num_a2+0.0016*num_a1+0.6258*num_a3;1.0',
'num_a04->0.0000*num_a2+0.0000*num_a1+0.3501*num_a3;1.0',
/

参考的原inp文件如下:
MOZART-MOSAIC-4bin.inp
&control do_bc = .true. do_ic = .true. domain = 1 dir_wrf = '/media/small/ybjin/WRFV3.9.1.1/Build_WRF/WRFV3/run/' dir_moz = '/data/ybjin/mozbc/hdata/WAMCC/' fn_moz = 'h0003.nc' moz_var_suffix = '' def_missing_var = .true. spc_map ='o3 -> O3', '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', 'isopr -> ISOP', 'acetol -> HYAC', 'mek -> MEK', 'bigene -> BIGENE', 'bigalk -> BIGALK', 'tol -> TOLUENE', 'benzene -> BENZENE', 'cres -> CRESOL', 'dms -> DMS', 'so2 -> SO2', 'oc_a01->0.1216*pom_a1+0.9886*soa1_a2+0.1216*soa1_a1+0.9886*soa2_a2+0.1216*soa2_a1+0.9886*soa3_a2+0.1216*soa3_a1+0.9886*soa4_a2+0.1216*soa4_a1+0.9886*soa5_a2+0.1216*soa5_a1;1.e9', 'oc_a02->0.7618*pom_a1+0.0114*soa1_a2+0.7618*soa1_a1+0.0114*soa2_a2+0.7618*soa2_a1+0.0114*soa3_a2+0.7618*soa3_a1+0.0114*soa4_a2+0.7618*soa4_a1+0.0114*soa5_a2+0.7618*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.0000*soa4_a2+0.1164*soa4_a1+0.0000*soa5_a2+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.0000*soa4_a2+0.0002*soa4_a1+0.0000*soa5_a2+0.0002*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.1856*so4_a2+0.0050*so4_a1+0.0000*so4_a3;1.e9', 'nh4_a02->0.0021*so4_a2+0.0930*so4_a1+0.0000*so4_a3;1.e9', 'nh4_a03->0.0000*so4_a2+0.0203*so4_a1+0.0186*so4_a3;1.e9', 'nh4_a04->0.0000*so4_a2+0.0000*so4_a1+0.1690*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.3889*ncl_a2+0.0479*ncl_a1+0.0000*ncl_a3;1.e9', 'na_a02->0.0045*ncl_a2+0.2997*ncl_a1+0.0000*ncl_a3;1.e9', 'na_a03->0.0000*ncl_a2+0.0458*ncl_a1+0.0391*ncl_a3;1.e9', 'na_a04->0.0000*ncl_a2+0.0000*ncl_a1+0.3542*ncl_a3;1.e9', 'cl_a01->0.5996*ncl_a2+0.0737*ncl_a1+0.0000*ncl_a3;1.e9', 'cl_a02->0.0068*ncl_a2+0.4621*ncl_a1+0.0000*ncl_a3;1.e9', 'cl_a03->0.0000*ncl_a2+0.0709*ncl_a1+0.0604*ncl_a3;1.e9', 'cl_a04->0.0000*ncl_a2+0.0001*ncl_a1+0.5462*ncl_a3;1.e9', 'oin_a01->0.9886*dst_a2+0.1216*dst_a1+0.0000*dst_a3;1.e9', 'oin_a02->0.0114*dst_a2+0.7618*dst_a1+0.0002*dst_a3;1.e9', 'oin_a03->0.0000*dst_a2+0.1164*dst_a1+0.0995*dst_a3;1.e9', 'oin_a04->0.0000*dst_a2+0.0002*dst_a1+0.9003*dst_a3;1.e9', 'num_a01->0.9996*num_a2+0.6135*num_a1+0.0000*num_a3;1.0', 'num_a02->0.0004*num_a2+0.2847*num_a1+0.0239*num_a3;1.0', 'num_a03->0.0000*num_a2+0.0016*num_a1+0.6258*num_a3;1.0', 'num_a04->0.0000*num_a2+0.0000*num_a1+0.3501*num_a3;1.0', /

评论 (0 个评论)

facelist doodle 涂鸦板

您需要登录后才可以评论 登录 | 立即注册

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

返回顶部