FLEXPART模型安装和模拟（二）

对方的凤飞飞

好伤心呀，没人看我的帖子呀，更没有人和我交流，大家都这么高冷吗？之前某信大的师兄问我要800元教我简单的模拟，我真的没钱呀，我只能自己摸索呀，发了很多邮件没人理我，算了，还是继续摸索吧，这次的主题是运行和模拟，下面是我的COMAND文件*********************************************************************************                                                                              *
*      Input file for the Lagrangian particle dispersion model FLEXPART        *
*                           Please select your options                         *
*                                                                              *
********************************************************************************

1. __                3X, I2
   -1      
   LDIRECT           1 FOR FORWARD SIMULATION, -1 FOR BACKWARD SIMULATION

2. ________ ______   3X, I8, 1X, I6
   20140901 000000
   YYYYMMDD HHMISS   BEGINNING DATE OF SIMULATION

3. ________ ______   3X, I8, 1X, I6
   20140905 180000
   YYYYMMDD HHMISS   ENDING DATE OF SIMULATION

4. _____             3X, I5
    3600
   SSSSS             OUTPUT EVERY SSSSS SECONDS

5. _____             3X, I5
    3600
   SSSSS             TIME AVERAGE OF OUTPUT (IN SSSSS SECONDS)

6. _____             3X, I5
     900
   SSSSS             SAMPLING RATE OF OUTPUT (IN SSSSS SECONDS)

7. _________         3X, I9
   999999999
   SSSSSSSSS         TIME CONSTANT FOR PARTICLE SPLITTING (IN SECONDS)

8. _____             3X, I5
     900
   SSSSS             SYNCHRONISATION INTERVAL OF FLEXPART (IN SECONDS)

9.  ---.--           4X, F6.4
      3.0
    CTL              FACTOR, BY WHICH TIME STEP MUST BE SMALLER THAN TL

10. ---              4X, I3
      4
    IFINE            DECREASE OF TIME STEP FOR VERTICAL MOTION BY FACTOR IFINE

11. -                4X, I1
    5  
    IOUT             1 CONCENTRATION (RESIDENCE TIME FOR BACKWARD RUNS) OUTPUT, 2 MIXING RATIO OUTPUT, 3 BOTH,4 PLUME TRAJECT., 5=1+4

12. -                4X, I1
    2  
    IPOUT            PARTICLE DUMP: 0 NO, 1 EVERY OUTPUT INTERVAL, 2 ONLY AT END

13. _                4X, I1
    1
    LSUBGRID         SUBGRID TERRAIN EFFECT PARAMETERIZATION: 1 YES, 0 NO

14. _                4X, I1
    1
    LCONVECTION      CONVECTION: 1 YES, 0 NO

15. _                4X, I1
    1
    LAGESPECTRA      AGE SPECTRA: 1 YES, 0 NO

16. _                4X, I1
    0
    IPIN             CONTINUE SIMULATION WITH DUMPED PARTICLE DATA: 1 YES, 0 NO

17. _               
    1                4X,I1
    IOFR             IOUTPUTFOREACHREL CREATE AN OUPUT FILE FOR EACH RELEASE LOCATION: 1 YES, 0 NO

18. _                4X, I1
    1
    IFLUX            CALCULATE FLUXES: 1 YES, 0 NO

19. _                4X, I1
    0
    MDOMAINFILL      DOMAIN-FILLING TRAJECTORY OPTION: 1 YES, 0 NO, 2 STRAT, 3 TRACER

20. _                4X, I1
    1
    IND_SOURCE       1=MASS UNIT , 2=MASS MIXING RATIO UNIT

21. _                4X, I1
    1
    IND_RECEPTOR     1=MASS UNIT , 2=MASS MIXING RATIO UNIT

22. _                4X, I1
    0
    MQUASILAG        QUASILAGRANGIAN MODE TO TRACK INDIVIDUAL PARTICLES: 1 YES, 0 NO

23. _                4X, I1
    0
    NESTED_OUTPUT    SHALL NESTED OUTPUT BE USED? 1 YES, 0 NO

24. _                4X, I1
    1
    LINIT_COND       INITIAL COND. FOR BW RUNS: 0=NO,1=MASS UNIT,2=MASS MIXING RATIO UNIT
我的是后向模拟，11是模拟输出文件，如果要做温度和比湿，12是要选1的，19中的DOMAIN-FILLING不懂呀，大家可以看陈斌老师的文章告诉我怎么回事，我是真不懂气象方面的东西。
********************************************************************************
*                                                                              *
*      Input file for the Lagrangian particle dispersion model FLEXPART        *
*                       Please specify your output grid                        *
*                                                                              *
********************************************************************************

1.  ------.----       4X,F11.4
        45.0000       GEOGRAFICAL LONGITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
    OUTLONLEFT        (left boundary of the first grid cell - not its centre)

2.  ------.----       4X,F11.4
       -20.0000       GEOGRAFICAL LATITUDE OF LOWER LEFT CORNER OF OUTPUT GRID
    OUTLATLOWER       (lower boundary of the first grid cell - not its centre)

3.  -----             4X,I5
      270             NUMBER OF GRID POINTS IN X DIRECTION (= No. of cells + 1)
    NUMXGRID

4.  -----             4X,I5
      140             NUMBER OF GRID POINTS IN Y DIRECTION (= No. of cells + 1)
    NUMYGRID

5.  ------.---        4X,F10.3
         0.500        GRID DISTANCE IN X DIRECTION
    DXOUTLON

6.  ------.---        4X,F10.3
         0.500        GRID DISTANCE IN Y DIRECTION
    DYOUTLAT

12. -----.-           4X, F7.1
     1000.0
    LEVEL 8           HEIGHT OF LEVEL (UPPER BOUNDARY)

Outgrid文件是设置经纬网网格，0.5乘以0.5，模拟的高度，这里的是是一层的高度，还是从地面到1000m积分，求会的告知。
*************************************************************************
*                                                                       *
*                                                                       *
*                                                                       *
*   Input file for the Lagrangian particle dispersion model FLEXPART    *
*                        Please select your options                     *
*                                                                       *
*                                                                       *
*                                                                       *
*************************************************************************
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  1      
___                        i3    Total number of species emitted

24
___                        i3    Index of species in file SPECIES

=========================================================================
20140901 000000
________ ______            i8,1x,i6 Beginning date and time of release

20140905 180000
________ ______            i8,1x,i6 Ending date and time of release

118.95
____.____                  f9.4  Longitude [DEG] of lower left corner

  32.12
____.____                  f9.4  Latitude [DEG] of lower left corner

118.95
____.____                  f9.4  Longitude [DEG] of upper right corner

  32.12
____.____                  f9.4  Latitude [DEG] of upper right corner

       2
_________                  i9    1 for m above ground, 2 for m above sea level,3 for pressure in hPa

  800.0
_____.___                  f10.3 Lower z-level (in m agl or m asl)

1000.0
_____.___                  f10.3 Upper z-level (in m agl or m asl)

    20000
_________                  i9    Total number of particles to be released

12.0000E00
_.____E__                  e9.4  Total mass emitted

RELEASE_TEST1
________________________________________   character*40 comment
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

设置粒子发射时间高度范围和粒子总数，粒子发射的经纬度，我设置的大概是南京的范围。接着看，往下看！
1 添加pathname文件：设置资料输入输出数据的路径以及数据描述文件的路径
/home/four/flexpart_9.0.3/options/
/mnt/hgfs/model1/output/
/mnt/hgfs/model/
/home/four/lc/AVAILABLE
2 创建output文件夹：输出数据的位置
这就不说了吧
3 创建input文件夹：输入数据的位置
同上
4 添加AVAILABLE文件：输入数据的描述文件
DATE     TIME         FILENAME     SPECIFICATIONS
YYYYMMDD HHMISS
________ ______      __________      _________
20140901 000000      fnl_20140901_00_00.grib2      ON DISC
20140901 060000      fnl_20140901_06_00.grib2      ON DISC
20140901 120000      fnl_20140901_12_00.grib2      ON DISC
20140901 180000      fnl_20140901_18_00.grib2      ON DISC
20140902 000000      fnl_20140902_00_00.grib2      ON DISC
20140902 060000      fnl_20140902_06_00.grib2      ON DISC
20140902 120000      fnl_20140902_12_00.grib2      ON DISC
20140902 180000      fnl_20140902_18_00.grib2      ON DISC
20140903 000000      fnl_20140903_00_00.grib2      ON DISC
20140903 060000      fnl_20140903_06_00.grib2      ON DISC
20140903 120000      fnl_20140903_12_00.grib2      ON DISC
20140903 180000      fnl_20140903_18_00.grib2      ON DISC
20140904 000000      fnl_20140904_00_00.grib2      ON DISC
20140904 060000      fnl_20140904_06_00.grib2      ON DISC
20140904 120000      fnl_20140904_12_00.grib2      ON DISC
20140904 180000      fnl_20140904_18_00.grib2      ON DISC
20140905 000000      fnl_20140905_00_00.grib2      ON DISC
20140905 060000      fnl_20140905_06_00.grib2      ON DISC
20140905 120000      fnl_20140905_12_00.grib2      ON DISC
20140905 180000      fnl_20140905_18_00.grib2      ON DISC

5 修改option文件夹内内容
上面说了
6 ./FLEXPART_GFS_GFORTRAN
结果有恭喜你字样！
7.结果展示
数据处理我用到了matlab,python中的quicklook，但是这个quicklook是个傻瓜软件，还有一个python包在linux中一直安装不了！F2py安装不了！！！！！！！！！！！！！！！！！！下面是我的图的展示！希望大家和我交流怎么处理！

菜鸟的梦想

可以请教一下如何分析结果吗？
如果方便的话，可以有偿请教

对方的凤飞飞

菜鸟的梦想 发表于 2018-11-20 21:15
可以请教一下如何分析结果吗？
如果方便的话，可以有偿请教

你加我扣扣984008820

haozl07

关注，感兴趣，曾经运行到生成数据文件那一步，卡在了可视化显示那一步！

bobbycao2008

很想学FLEXPART，请问楼主多长时间能掌握基本操作，简单运行出结果

对方的凤飞飞

交流一下还是很快的

imastrid

想问一下楼主用过flexpartwrf吗~

imastrid

楼主，你反向模拟的数据有试过全是0的吗~

beijing135830

您的数据范围是 20140901 为什么输出结果是2013年的呢

Grace123

楼主你太搞笑了，话说现在找你交流学习一下可以吗？