wrf-chem转出变量数值异常的问题

SG哭晕在厕所

astiny 发表于 2016-12-5 19:06
目测anthro_emis做的源吧233
这个只能用来给mozcart的那四个模式做源。或者用来做逐时次的，然后再在它 ...

我用的是2006年INTEXB的排放源目录

astiny

SG哭晕在厕所 发表于 2016-12-5 19:14
我用的是2006年INTEXB的排放源目录

prep_chem_src是emission guide 里面介绍的用来做排放源文件wrfchemi_xxZ_d0x的。。。
然后一开始提到的那个anthro_emiss是逐时次生成wrfchemi*文件的。。。
然后wrfout文件输出频次的话，分时次和输一起都是可以的。。。

SG哭晕在厕所

astiny 发表于 2016-12-5 19:22
prep_chem_src是emission guide 里面介绍的用来做排放源文件wrfchemi_xxZ_d0x的。。。
然后一开始提到的 ...

也就是说prep_chem_src是另一种生成wrfchemi的方法？之后namamelist.input的修改方法都不一样吧，那请问大神。。。有没有教程什么的能学习学习。。。因为我这种排放源还是学姐传下来的。。。

astiny

SG哭晕在厕所 发表于 2016-12-5 19:31
也就是说prep_chem_src是另一种生成wrfchemi的方法？之后namamelist.input的修改方法都不一样吧，那请问 ...

官网的emission guide啊。。。

astiny

总之，你想办法把排放源文件的e_pm25改成mol/km2/h就对了。

SG哭晕在厕所

astiny 发表于 2016-12-5 20:12
总之，你想办法把排放源文件的e_pm25改成mol/km2/h就对了。

还有一个问题。。我发现排放源文件里没有e_pm25。。。只有e_pm25i和e_pm25j，而wrfout里的PM2_5_DRY是怎么出来的。。。

夜未央

astiny 发表于 2016-12-5 19:06
目测anthro_emis做的源吧233
这个只能用来给mozcart的那四个模式做源。或者用来做逐时次的，然后再在它 ...

有chem的地方就会看到这个萌萌熊，我是来看萌萌熊的！

astiny

SG哭晕在厕所 发表于 2016-12-5 20:24
还有一个问题。。我发现排放源文件里没有e_pm25。。。只有e_pm25i和e_pm25j，而wrfout里的PM2_5_DRY是怎 ...

有的chem_opt读的是i,j,有的是没有ij的。具体可以在registry.chem里面看。i和j代表的好像是不同模态。一个是核模一个是积聚模好像。如果一定需要e_pm25的话可以把两个加起来（大概。我猜的。）

astiny

夜未央 发表于 2016-12-5 20:58
有chem的地方就会看到这个萌萌熊，我是来看萌萌熊的！

那是因为pangzi_xu大神被永封了233

say苹果有虫

astiny 发表于 2016-12-5 19:06
目测anthro_emis做的源吧233
这个只能用来给mozcart的那四个模式做源。或者用来做逐时次的，然后再在它 ...

大神，我用prep_chem_src做源，正常生成了*.bin文件，正确链接并改了名字，但是就在convert_emiss.exe的时候，出现了error：
   PGFI0-F-217/unformatted read/unit=19/attempt to read  past end of file.
        File name = wrf_gocart_backg    unformatted, swquential access  record =2
   In source file convert_emiss.f90, at line number 393
在convert_emiss.f90中：
382    IF (wrf_dm_on_monitor()) THEN
383      OPEN(19,FILE='wrf_gocart_backg',FORM='UNFORMATTED')
384    endif
385
386        IF (wrf_dm_on_monitor()) THEN
387            read(19)tmp2(ids:ide  ,jds:jde  )
388        ENDIF
389        CALL wrf_dm_bcast_bytes ( tmp2 , size ( tmp2 ) * rnum8 )
390       grid%dms_0(ips:ipe,jps:jpe) = tmp2(ips:ipe,jps:jpe)
391
392        IF (wrf_dm_on_monitor()) THEN
393            read(19)tmp3(ids:ide  ,jds:jde  ,1:3)
394        ENDIF
395        CALL wrf_dm_bcast_bytes ( tmp3 , size ( tmp3 ) * rnum8 )

求大神指点，这是哪里出了问题？我用的是PREP_CHEM_SRC_1.5, NETCDF-4.1.3, PGI-10.6, CentOS-6.8.