关于prep-chem-src的用法

央哥儿

从本科起就在家园学到很多东西，去年开始做模式，一开始简直晕头转向，在家园里获得了很多帮助，今天答应别人写prep-chem-src的用法，为了回馈家园的帮助，也为了帮助更多学习模式的同学。虽然我也经验不够充分，写出来给大家参考，一起讨论。prep-chem-src是一个清单前处理模型，用在WRF/Chem, convert_emiss.exe之前。
安装方法 http://bbs.06climate.com/forum.php?mod=viewthread&tid=51340 这里讲得十分清楚。
主要用法是：
bin文件夹下：prep_chem_sources.inp （作用类似于namelist）和 prep_chem_sources_RADM_WRF_FIM_.exe （这里我主要使用RADM作为化学机制）
prep_chem_sources.inp：
RP_INPUT
!#################################################################
!  CCATT-BRAMS/MCGA-CPTEC/WRF-Chem/FIM-Chem emission models CPTEC/INPE
!  version 1.5: Mar 2015
!  contact: gmai@cptec.inpe.br   - http://meioambiente.cptec.inpe.br
!#################################################################


!---------------- grid_type of the grid output
   grid_type= 'lambert',
   rams_anal_prefix = '../ANL/OPQUE',
!---------------- date of emission
    ihour=00,
    iday=26,
    imon=09,
    iyear=2016,

!---------------- select the sources datasets to be used（这里有不同的排放清单，大部分从prep-chem研发团队的网站上可以下载，下载后放入特定路径中，然后将路径分写入下面xxx_dir=''）
   use_retro=1,  ! 1 = yes, 0 = not
   retro_data_dir='/pathway/RETRO/anthro',

   use_edgar =3,  ! 0 - not,
                  ! 1 - Version 3,
                  ! 2 - Version 4 for some species
                  ! 3 - Version HTAP

   edgar_data_dir='/pathway/EDGAR-HTAP',

   use_gocart=1,
   gocart_data_dir='/pathway/emissions',

   use_streets =1,
   streets_data_dir='pathway/STREETS',

   use_seac4rs =1,
   seac4rs_data_dir='/pathway/SEAC4RS',


   use_fwbawb =1,

fwbawb_data_dir ='/pathway/Emissions_Yevich_Logan',

   use_bioge =2, ! 1 - geia, 2 - megan
   ! ######
   ! # BIOGENIC = 0
   bioge_data_dir ='/pathway/MEGAN/2000',
   ! # MEGAN = 1
   ! ######
   !bioge_data_dir='/home/Yang/wrf_chem/emissions/source_data/MEGAN/2000',
   ! ######

   use_gfedv2=1,
   gfedv2_data_dir='/home/Yang/wrf_chem/emissions/source_data/GFEDv2-8days',

   use_bbem=0,
   use_bbem_plumerise=0,

!--------------------------------------------------------------------------------------------------

!---------------- if  the merging of gfedv2 with bbem is desired (=1, yes, 0 = no)
   merge_GFEDv2_bbem =0,

!---------------- Fire product for BBBEM/BBBEM-plumerise emission models
   bbem_wfabba_data_dir   ='./datain/Emission_data/fires_data/GOES/f',
   bbem_modis_data_dir    ='./datain/Emission_data/fires_data/MODIS/Fires',
   bbem_inpe_data_dir     ='./datain/Emission_data/fires_data/DSA/Focos',
   bbem_extra_data_dir    ='NONE',

!---------------- veg type data set (dir + prefix)
   veg_type_data_dir      ='./datain/surface_data/GL_IGBP_MODIS_INPE/MODIS',


!---------------- vcf type data set (dir + prefix)
  use_vcf = 0,
  vcf_type_data_dir      ='./datain/surface_data/VCF/data_out/2005/VCF',
!---------------- olson data set (dir + prefix)
  olson_data_dir      ='./datain/Emission_data/OLSON2/OLSON',

!---------------- carbon density data set (dir + prefix)

   carbon_density_data_dir='./datain/surface_data/GL_OGE_INPE/OGE',

   fuel_data_dir      ='./datain/Emission_data/Carbon_density_Saatchi/amazon_biomass_final.gra',


!---------------- gocart background
   use_gocart_bg=1,
   gocart_bg_data_dir='/home/Yang/wrf_chem/emissions/source_data/GOCART',

!---------------- volcanoes emissions
   use_volcanoes =0,
   volcano_index =1143, !REDOUBT

   use_these_values='NONE',
! define a text file for using external values for INJ_HEIGHT, DURATION,
! MASS ASH (units are meters - seconds - kilograms) and the format for
   begin_eruption='201303280000',  !begin time UTC of eruption YYYYMMDDhhmm

!---------------- degassing volcanoes emissions
   use_degass_volcanoes =0,
   degass_volc_data_dir ='./datain/Emission_data/VOLC_SO2',

!---------------- user specific  emissions directory
!---------------- Update for South America megacities
   user_data_dir='NONE',


!--------------------------------------------------------------------------------------------------
   pond=1,   ! mad/mfa  0 -> molar mass weighted
             !          1 -> Reactivity weighted

!---------------- for grid type 'll' or 'gg' only
   grid_resolucao_lon=1.0,
   grid_resolucao_lat=1.0,


   nlat=320,          ! if gg (only global grid)
   lon_beg   = -180., ! (-180.:+180.) long-begin of the output file
   lat_beg   =  -90., ! ( -90.:+90. ) lat -begin of the output file
   delta_lon =  360, ! total long extension of the domain (360 for global)
   delta_lat =  180, ! total lat  extension of the domain (180 for global)

!---------------- For regional grids (polar or lambert)

   NGRIDS   = 2,            ! Number of grids to run（domian数量）

   NNXP     = 91,79,86,46,        ! Number of x gridpoints （同namelist.wps中e_ws）
   NNYP     = 73,61,74,46,        ! Number of y gridpoints
   NXTNEST  = 0,1,1,1,          ! Grid number which is the next coarser grid （同namelist.wps中parent_id）
   DELTAX   = 45000.,
   DELTAY   = 45000.,         ! X and Y grid spacing （最外面一层domain的间距）

   ! Nest ratios between this grid and the next coarser grid.
   NSTRATX  = 1,3,3,4,           ! x-direction （domain之间间距的比例）
   NSTRATY  = 1,3,3,4,           ! y-direction

   NINEST = 1,37,0,0,        ! Grid point on the next coarser （同namelist.wps中i_parent_start）
   NJNEST = 1,23,0,0,        !  nest where the lower southwest
                             !  corner of this nest will start.
                             !  If NINEST or NJNEST = 0, use CENTLAT/LON
   POLELAT  =  28.14, !-89.99,          ! If polar, latitude/longitude of pole point (同namelist.wps中ref_lat)
   POLELON  =  112.771,         ! If lambert, lat/lon of grid origin (x=y=0.)

   STDLAT1  = 30.,           ! If polar for BRAMS, use 90.0 in STDLAT2
   STDLAT2  = 60.,         ! If lambert, standard latitudes of projection
                            !(truelat2/truelat1 from namelist.wps, STDLAT1 < STDLAT2)
                            ! If mercator STDLAT1 = 1st true latitude
   CENTLAT  = 28.14,!-89.99,  -23., 27.5,  27.5, （与上面polelat一样）
   CENTLON  = 112.771,  -46.,-80.5, -80.5,



!---------------- model output domain for each grid (only set up for rams)
   lati =  -90.,  -90.,   -90.,

   latf =  +90.,  +90.,   +90.,
   loni = -180., -180.,  -180.,
   lonf =  180.,  180.,   180.,

!---------------- project rams grid (polar sterogr) to lat/lon: 'YES' or 'NOT'
   proj_to_ll='YES',

!---------------- output file prefix (may include directory other than the current) （这里可以修改输出文件名的前缀）
   chem_out_prefix = 'WRF',
   chem_out_format = 'vfm',
!---------------- convert to WRF/CHEM (yes,not)
  special_output_to_wrf = 'YES',

$END


修改好.inp后就可以./prep_chem_sources_RADM_WRF_FIM_.exe
接着会生成文件
WRF-T-2016-11-15-000000-g3-ab.bin （人为排放）
WRF-T-2016-11-15-000000-g3-bb.bin  （biomass）
WRF-T-2016-11-15-000000-g3-gocartBG.bin （全球背景排放）
将这些文件分别链接为
emissopt3_d01
emissfire_d01
wrf_gocart_bg
然后./convert_emiss.exe
最后根据namelist.input里的定义修改清单的文件名，这里注意清单需要一层一层的转
我目前也就会这样简单的用法啦，个人感觉这个清单处理模型更偏向于全球清单，对于区域尺度的怎么处理我还没研究出来，有研究出来的小伙伴能提供帮助吗，特别是MICS-ASIA，问了几个人都说要自己写程序转，我还没研究出来这个程序怎么搞，可以给点提示吗~~~

Dahoe

楼主好样的，请问您知道怎么把meic清单应用到WRF-Chem里了吗？

guodp

你好！
我运行prep-chem-src的时候提示，WF_ABBA loohing for none2017032
checking directory-none2017032
error opening /tme/fire_filelist

xyan88

请问运行第二层的时候 只需要把namelist.inpu第二列换到第一列 其他任何地方都不需要改吗，我这样子做得到的结果数据好像不对

Kimmy2014

非常感谢楼主，根据楼主的总结终于跑通一次，鸡冻~

Kimmy2014

请问楼主可以贴一下用convert_emiss.exe时候的namelist.input设置嘛？运行的时候还是有点问题

Sun丶Yang

不知楼主现在研究了怎样样用区域尺度的源清单处理方法没，是否可以交流一下呢？

Legendary

xyan88 发表于 2019-5-16 20:49
请问运行第二层的时候 只需要把namelist.inpu第二列换到第一列 其他任何地方都不需要改吗，我这样子做得到 ...

因为生成的文件是同样的文件名，你需要改之前大domain生成文件的文件名，或者把他们转移到另外一个目录里，然后再生成新的domain02的文件（文件名其实还是对应domian01的），最后把真正的domain 01的文件再搬运回来。max_dom =  1也需要改成3以后再一起跑一次。

yangxiajie

请问楼主，这个prep_chem_src可以处理SMOKE数据吗？谢谢~

央哥儿

yangxiajie 发表于 2019-10-27 14:16
请问楼主，这个prep_chem_src可以处理SMOKE数据吗？谢谢~

没有用过smoke的数据，smoke里不是有一套自己的格点处理方式吗