爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 16677|回复: 10

释放物为什么没怎么离开释放点。

[复制链接]

新浪微博达人勋

发表于 2017-9-7 08:18:00 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
想模拟下释放物的扩散趋势,输入文件如下  
用quicklook转成图片如附件
和用历史数据的结果很不一样,
历史结果中,释放物很快飘移出源释放地,
可是在这个模拟中,感觉源一直在释放。  
能帮忙找下是哪儿的原因吗?  
多谢。







=====================FORMER PATHNAMES FILE===================
/home/wrf/Work/ATM/FLEXPART-WRF/Src_flexwrf_v3.3.2/20170903/
/home/wrf/Data/WRF/Wrfdata/
/home/wrf/Data/WRF/Wrfdata/AVAILABLE_20170903
=============================================================
=====================FORMER COMMAND FILE=====================
    1                LDIRECT:          1 for forward simulation, -1 for backward simulation
    20170903 000000  YYYYMMDD HHMISS   beginning date of simulation
    20170910 000000  YYYYMMDD HHMISS   ending date of simulation
    10800             SSSSS  (int)      output every SSSSS seconds
    3600             SSSSS  (int)      time average of output (in SSSSS seconds)
    180              SSSSS  (int)      sampling rate of output (in SSSSS seconds)
    999999999        SSSSS  (int)      time constant for particle splitting (in seconds)
    180              SSSSS  (int)      synchronisation interval of flexpart (in seconds)
    10.              CTL    (real)     factor by which time step must be smaller than tl
    10               IFINE  (int)      decrease of time step for vertical motion by factor ifine
    1                IOUT              1 concentration, 2 mixing ratio, 3 both, 4 plume traject, 5=1+4
    0                IPOUT             particle dump: 0 no, 1 every output interval, 2 only at end
    0                LSUBGRID          subgrid terrain effect parameterization: 1 yes, 0 no
    0                LCONVECTION       convection: 3 yes, 0 no
    3600.            DT_CONV  (real)   time interval to call convection, seconds
    1                LAGESPECTRA       age spectra: 1 yes, 0 no
    0                IPIN              continue simulation with dumped particle data: 1 yes, 0 no
    0                IFLUX             calculate fluxes: 1 yes, 0 no
    0                IOUTPUTFOREACHREL CREATE AN OUPUT FILE FOR EACH RELEASE LOCATION: 1 YES, 0 NO
    0                MDOMAINFILL       domain-filling trajectory option: 1 yes, 0 no, 2 strat. o3 tracer
    1                IND_SOURCE        1=mass unit , 2=mass mixing ratio unit
    1                IND_RECEPTOR      1=mass unit , 2=mass mixing ratio unit
    1                NESTED_OUTPUT     shall nested output be used? 1 yes, 0 no
    0                LINIT_COND   INITIAL COND. FOR BW RUNS: 0=NO,1=MASS UNIT,2=MASS MIXING RATIO UNIT
    1                TURB_OPTION       0=no turbulence; 1=diagnosed as in flexpart_ecmwf; 2 and 3=from tke.
    1                LU_OPTION         0=old landuse (IGBP.dat); 1=landuse from WRF
    0                CBL SCHEME        0=no, 1=yes. works if TURB_OPTION=1
    0                SFC_OPTION        0=default computation of u*, hflux, pblh, 1=from wrf
    0                WIND_OPTION       0=snapshot winds, 1=mean winds,2=snapshot eta-dot,-1=w based on divergence
    1                TIME_OPTION       1=correction of time validity for time-average wind,  0=no need
    1                OUTGRID_COORD     0=wrf grid(meters), 1=regular lat/lon grid
    1                RELEASE_COORD     0=wrf grid(meters), 1=regular lat/lon grid
    0                IOUTTYPE          0=default binary, 1=ascii (for particle dump only),2=netcdf
    3                NCTIMEREC (int)   Time frames per output file, only used for netcdf
    0                VERBOSE           VERBOSE MODE,0=minimum, 100=maximum
=====================FORMER AGECLASESS FILE==================
    1                NAGECLASS        number of age classes
    999999             SSSSSS  (int)    age class in SSSSS seconds
=====================FORMER OUTGRID FILE=====================
   110.00            OUTLONLEFT      geograhical longitude of lower left corner of output grid
    30.0              OUTLATLOWER     geographical latitude of lower left corner of output grid
    400               NUMXGRID        number of grid points in x direction (= # of cells )
    200               NUMYGRID        number of grid points in y direction (= # of cells )
    0                OUTGRIDDEF      outgrid defined 0=using grid distance, 1=upperright corner coordinate
    0.1          DXOUTLON        grid distance in x direction or upper right corner of output grid
    0.1           DYOUTLON        grid distance in y direction or upper right corner of output grid
    7                NUMZGRID        number of vertical levels
    100.0            LEVEL           height of level (upper boundary)
    400.0            LEVEL           height of level (upper boundary)
   1000.0            LEVEL           height of level (upper boundary)
   2000.0            LEVEL           height of level (upper boundary)
   5000.0            LEVEL           height of level (upper boundary)
  10000.0            LEVEL           height of level (upper boundary)
  20000.0            LEVEL           height of level (upper boundary)
================OUTGRID_NEST==========================
   120.00             OUTLONLEFT      geograhical longitude of lower left corner of output grid
    35.00             OUTLATLOWER     geographical latitude of lower left corner of output grid
    400               NUMXGRID        number of grid points in x direction (= # of cells )
    200               NUMYGRID        number of grid points in y direction (= # of cells )
    0                OUTGRIDDEF      outgrid defined 0=using grid distance, 1=upperright corner coordinate
    0.05            DXOUTLON        grid distance in x direction or upper right corner of output grid
    0.05            DYOUTLON        grid distance in y direction or upper right corner of output grid
=====================FORMER RECEPTOR FILE====================
    6                NUMRECEPTOR     number of receptors
    RN20
    116.2
    39.8
    RN37
    127.9
    26.5
    RN38
    139.0
    36.3
    RN58
    131.9
    43.7
    QH
    130.5
    42.5
    CB
    128.12
    42.4
=====================FORMER SPECIES FILE=====================
     1               NUMTABLE        number of variable properties. The following lines are fixed format
XXXX|NAME    |decaytime |wetscava  |wetsb|drydif|dryhenry|drya|partrho  |parmean|partsig|dryvelo|weight |
    Xe-133     198720.0   -9.9E-09         -9.9                 -9.9E09                   -9.99   133.00
=====================FORMER RELEEASES FILE===================
   1                NSPEC           total number of species emitted
   0                EMITVAR         1 for emission variation
   1                LINK            index of species in file SPECIES
1                 NUMPOINT        number of releases
20170903 033000   ID1, IT1        beginning date and time of release
20170903 043000   ID2, IT2        ending date and time of release
129.1100         XPOINT1 (real)  longitude [deg] of lower left corner
  41.3500         YPOINT1 (real)  latitude [deg] of lower left corner
129.1100         XPOINT2 (real)  longitude [deg] of upper right corner
  41.3500         YPOINT2 (real)  latitude [DEG] of upper right corner
        1         KINDZ  (int)  1 for m above ground, 2 for m above sea level, 3 pressure
   100.000        ZPOINT1 (real)  lower z-level
   200.000        ZPOINT2 (real)  upper z-level
  100000          NPART (int)     total number of particles to be released
1.000E+15         XMASS (real)    total mass emitted
box1              NAME OF RELEASE LOCATION


mother_0-4_anim_mother.gif
mother_0-4_anim_mother.gif
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-9-7 09:16:02 | 显示全部楼层
是不是WIND_OPTION选择0=snapshot winds的原因,我觉得如果WIND_OPTION选择1的话应该会有移动。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-9-7 19:20:24 | 显示全部楼层
刘勇的下铺 发表于 2017-9-7 09:16
是不是WIND_OPTION选择0=snapshot winds的原因,我觉得如果WIND_OPTION选择1的话应该会有移动。

确实,谢谢!  
这个在物理上的解释哪儿能查到呢?  
多谢!  
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-9-8 09:18:01 | 显示全部楼层
这个物理解释我也不是很清楚。不过我也遇到过这个问题,我也是吧0改成1就好了。我后来看那个模式介绍的文章,觉得好像就是瞬时风和平均风的区别,不过为什么选0的时候会不扩散我也不知道其中的原因。还有楼主是不是在微信群里发了这个问题哈,这个设置发到微信群里看起来格式错乱,根本看不出来哪里有问题。还是这样发个帖子或者截图一眼就能看出来。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-9-8 11:33:57 | 显示全部楼层
刘勇的下铺 发表于 2017-9-8 09:18
这个物理解释我也不是很清楚。不过我也遇到过这个问题,我也是吧0改成1就好了。我后来看那个模式介绍的文章 ...

好的,谢谢!  
多交流!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-9-8 11:44:26 | 显示全部楼层
0改成1
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2018-3-24 21:53:41 | 显示全部楼层

怎么理解呢?0和1的差别
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2019-2-27 17:03:23 | 显示全部楼层
想请教一下flexpart的结果nc文件怎么处理成上面那种效果图
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-5-25 16:30:56 | 显示全部楼层
请教楼主,FORMER OUTGRID FILE那部分和OUTGRID_NEST那部分设置的区别是什么呢?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-1-10 21:24:00 | 显示全部楼层
xiaoguoguo327 发表于 2020-5-25 16:30
请教楼主,FORMER OUTGRID FILE那部分和OUTGRID_NEST那部分设置的区别是什么呢?

一个是大网格,一个是小网格。  
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

快速回复 返回顶部 返回列表