运行SMOKE的点源程序point时的错误

lwkraul2019

大家好，我用SMOKE3.5，能运行smk_area.csh并生成正常的结果，但在跑smk_point.csh的时候未成功，在grdmat.point.log中出现了以下的错误：
=============================================
     File "GRIDDESC" opened for input on unit:  96
     /home/cmaq246/SMOKE3.5.1/data/ge_dat/GRIDDESC-baise

     Value for IOAPI_GRIDNAME_1:  'baise'

     File "GRIDDESC" opened for input on unit:  96
     /home/cmaq246/SMOKE3.5.1/data/ge_dat/GRIDDESC-baise

     NOTE: Grid settings initialized using baise            in
          grid description file.
     NOTE: Output grid "baise" set; described as
          No description available
     Value for IOAPI_ISPH:  '12'
     INITSPHERES:  input sphere WGS 84
     Value for IOAPI_GRIDNAME_1:  'baise'

     *** ERROR ABORT in subroutine GRDMAT
     No source-cell intersections found.

========================================
根据以上错误信息寻找了源代码，出错原因是点源没有落在模拟区域之内，但我反复核查了我的区域和点源（只有1个），肯定在区域之内。我的GRIDDESC情况如下：
========================================
' '
'baise'
  2        15.000        40.000       106.170       106.170        27.500
' '
'baise'
'baise'   -177001.516   -563697.250      3000.000      3000.000 118 118   1
' '=======================================

我的点源位置是（24N，106.17E），选择的是整个模拟区域的正中央，恳求各位高手帮忙看看问题

weilin421

xenos_an 发表于 2020-6-30 14:00
问题解决了吗？我也遇到同样的问题

这个脚本里有一个参数 会把西经弄成正数，东经弄成负数，就是对经度取相反数，默认例子脚本这个参数是开启的，你得把这个参数关了，东经才会真的落在东经的位置上，就不会报这个了。
否则你写的东经经度，它都进来默认按西经计算，自然没有。

xenos_an

问题解决了吗？我也遇到同样的问题

xenos_an

weilin421 发表于 2020-6-30 14:19
这个脚本里有一个参数 会把西经弄成正数，东经弄成负数，就是对经度取相反数，默认例子脚本这个参数是开 ...

真是感谢呢，果然有这么个参数存在：
setenv WEST_HSPHERE         N   # Y converts longitudes to negative values

{:eb511:}{:eb511:}

xenos_an

weilin421 发表于 2020-6-30 14:19
这个脚本里有一个参数 会把西经弄成正数，东经弄成负数，就是对经度取相反数，默认例子脚本这个参数是开 ...

还有个问题需要请教，我使用MCIP输出的GRIDDESC输入SMOKE，SMOKE输出再输入到BCON的时候，显示 和MCIP的网格坐标不一致：

     Maximum current record number        67
          Differences found between files CTM_CONC_1       and MET_CRO_3D_CRS  :
          CTM_CONC_1 XORIG = -3.24000E+05   MET_CRO_3D_CRS XORIG = -2.80500E+05
          CTM_CONC_1 YORIG = -3.28500E+05   MET_CRO_3D_CRS YORIG = -2.85000E+05
          CTM_CONC_1 VGTOP = -9.99900E+36   MET_CRO_3D_CRS VGTOP =  5.00000E+03
          CTM_CONC_1 VGLVS( 1 ) =  0.00000E+00MET_CRO_3D_CRS VGLVS( 1 ) =  1.00000E+00
          CTM_CONC_1 VGLVS( 2 ) =  0.00000E+00MET_CRO_3D_CRS VGLVS( 2 ) =  9.98800E-01

     *** ERROR ABORT in subroutine CK_MET_FL
     MET_CRO_3D_CRS file inconsistent: stopping

0.027u 0.027s 0:00.05 80.0%     0+0k 0+0io 0pf+0w
exit ( )

一时没想明白

weilin421

你mcip 输出的griddesc那几个参数是多少。

weilin421

你这是cmaq的griddesc和mcip输出的气象文件不匹配呀。你是不是用了自己的griddesc但是气象数据还指向的例子数据。

xenos_an

weilin421 发表于 2020-7-6 18:20
你这是cmaq的griddesc和mcip输出的气象文件不匹配呀。你是不是用了自己的griddesc但是气象数据还指向的例子 ...

从头跑了一遍后OK了，感谢支持