登录后查看更多精彩内容~
您需要 登录 才可以下载或查看，没有帐号？立即注册
x
本帖最后由 Legendary 于 2020-5-1 10:17 编辑
我目前在学习使用Flexpart-wrf。咨询关于Flexpart-wrf网格设置的问题。下文46-56行来自flexwrf.input，我理解的是填写WRFout文件的参数
6 =====================FORMER COMMAND FILE===================== 7 -1 LDIRECT: 1 for forward simulation, -1 for backward simulation 8 20170602 000000 YYYYMMDD HHMISS beginning date of simulation 9 20170608 000000 YYYYMMDD HHMISS ending date of simulation 10 3600 SSSSS (int) output every SSSSS seconds 11 3600 SSSSS (int) time average of output (in SSSSS seconds) 12 180 SSSSS (int) sampling rate of output (in SSSSS seconds) 13 999999999 SSSSS (int) time constant for particle splitting (in seconds) 14 180 SSSSS (int) synchronisation interval of flexpart (in seconds) 15 10. CTL (real) factor by which time step must be smaller than tl 16 10 IFINE (int) decrease of time step for vertical motion by factor ifine 17 5 IOUT 1 concentration, 2 mixing ratio, 3 both, 4 plume traject, 5=1+4 18 1 IPOUT particle dump: 0 no, 1 every output interval, 2 only at end 19 1 LSUBGRID subgrid terrain effect parameterization: 1 yes, 0 no 20 3 LCONVECTION convection: 3 yes, 0 no 21 3600. DT_CONV (real) time interval to call convection, seconds 22 0 LAGESPECTRA age spectra: 1 yes, 0 no 23 0 IPIN continue simulation with dumped particle data: 1 yes, 0 no 24 0 IFLUX calculate fluxes: 1 yes, 0 no 25 1 IOUTPUTFOREACHREL CREATE AN OUPUT FILE FOR EACH RELEASE LOCATION: 1 YES, 0 NO 26 0 MDOMAINFILL domain-filling trajectory option: 1 yes, 0 no, 2 strat. o3 tracer 27 2 IND_SOURCE 1=mass unit , 2=mass mixing ratio unit 28 2 IND_RECEPTOR 1=mass unit , 2=mass mixing ratio unit 29 0 NESTED_OUTPUT shall nested output be used? 1 yes, 0 no 30 2 LINIT_COND INITIAL COND. FOR BW RUNS: 0=NO,1=MASS UNIT,2=MASS MIXING RATIO UNIT 31 1 TURB_OPTION 0=no turbulence; 1=diagnosed as in flexpart_ecmwf; 2 and 3=from tke. 32 1 LU_OPTION 0=old landuse (IGBP.dat); 1=landuse from WRF 33 1 CBL SCHEME 0=no, 1=yes. works if TURB_OPTION=1 34 0 SFC_OPTION 0=default computation of u*, hflux, pblh, 1=from wrf 35 0 WIND_OPTION 0=snapshot winds, 1=mean winds,2=snapshot eta-dot,-1=w based on divergence 36 0 TIME_OPTION 1=correction of time validity for time-average wind, 0=no need 37 1 OUTGRID_COORD 0=wrf grid(meters), 1=regular lat/lon grid 38 1 RELEASE_COORD 0=wrf grid(meters), 1=regular lat/lon grid 39 2 IOUTTYPE 0=default binary, 1=ascii (for particle dump only),0=netcdf 40 500 NCTIMEREC (int) Time frames per output file, only used for netcdf 41 0 VERBOSE VERBOSE MODE,0=minimum, 100=maximum :set number 19,10 Top 42 =====================FORMER AGECLASESS FILE================== 43 2 NAGECLASS number of age classes 44 7200 SSSSSS (int) age class in SSSSS seconds 45 999999 SSSSSS (int) age class in SSSSS seconds 46 =====================FORMER OUTGRID FILE===================== 47 -95.0444 OUTLONLEFT geograhical longitude of lower left corner of output grid 48 39.3059 OUTLATLOWER geographical latitude of lower left corner of output grid 49 99 NUMXGRID number of grid points in x direction (= # of cells ) 50 81 NUMYGRID number of grid points in y direction (= # of cells ) 51 0 OUTGRIDDEF outgrid defined 0=using grid distance, 1=upperright corner coordinate 52 12000 DXOUTLON grid distance in x direction or upper right corner of output grid 53 12000 DYOUTLON grid distance in y direction or upper right corner of output grid 54 2 NUMZGRID number of vertical levels 55 100.0 LEVEL height of level (upper boundary) 56 20000.0 LEVEL height of level (upper boundary) 57 =====================FORMER RECEPTOR FILE==================== 58 0 NUMRECEPTOR number of receptors 59 =====================FORMER SPECIES FILE===================== 60 2 NUMTABLE number of variable properties. The following lines are fixed format 61 XXXX|NAME |decaytime |wetscava |wetsb|drydif|dryhenry|drya|partrho |parmean|partsig|dryvelo|weight | 62 AIRTRACER -999.9 -9.9E-09 -9.9 -9.9E09 -9.99 29.00 63 Cs-137 -999.9 1.0E-04 0.80 -9.9 2.5E03 6.0E-7 3.0E-1 -9.99 -9.99 64 =====================FORMER RELEEASES FILE=================== 65 1 NSPEC total number of species emitted 66 0 EMITVAR 1 for emission variation 67 1 LINK index of species in file SPECIES 68 1 NUMPOINT number of releases 69 20170603 100000 ID1, IT1 beginning date and time of release 70 20170603 110000 ID2, IT2 ending date and time of release 71 -87.6298 XPOINT1 (real) longitude [deg] of lower left corner 72 41.8781 YPOINT1 (real) latitude [deg] of lower left corner 73 -87.6298 XPOINT2 (real) longitude [deg] of upper right corner 74 41.8781 YPOINT2 (real) latitude [DEG] of upper right corner 75 1 KINDZ (int) 1 for m above ground, 2 for m above sea level, 3 pressure 76 10 ZPOINT1 (real) lower z-level 77 2000 ZPOINT2 (real) upper z-level 78 300 NPART (int) total number of particles to be released 79 .1000E+01 XMASS (real) total mass emitted 80 release1 NAME OF RELEASE LOCATION
运行后得到一堆报错信息：
Pre-generating random numbers Calling readinput Opening 'flexwrf.input_shujian-chuan ' for reading Reading pathnames option_verbose= 0 Opening file: /qw/feng/Flexpart/available/AVAILABLE_LMOS_0602 for reading
test_xyindex_to_ll_wrf -- lgrid, rmserr (km) = 0 1.76E-01
# of vertical levels in WRF data n_bottom_top & "true" nuvz: 32 32 nwz & "augmented" nuvz: 33 33
nwzmax, nuvzmax, nzmax : 61 61 61 nconvlevmax, nconvlev : 60 31 nx, ny, nxmax, nymax : 99 81 250 190
Mother domain: east-west range: 0.0 1176000.0 Grid distance: 12000.0 south-north range: 0.0 960000.0 Grid distance: 12000.0
Reading OUTGRID
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 1.188E+06 9.720E+05 x_i, y_j = -3.822E+02 9.100E+02 Reading RECEPTORS Reading SPECIES option_verbose: 0 option_verbose: 0 Reading RELEASES
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 5.050E+02 6.393E+02 x_i, y_j = 1.368E+03 6.132E+03
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 1.705E+03 6.393E+02 x_i, y_j = -4.624E+03 -2.762E+03
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 2.905E+03 6.393E+02 x_i, y_j = 5.749E+03 -3.787E+01
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 4.105E+03 6.393E+02 x_i, y_j = -3.840E+03 4.765E+03
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 5.305E+03 6.393E+02 x_i, y_j = 1.906E+02 -5.173E+03
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 6.505E+03 6.393E+02 x_i, y_j = 3.725E+03 4.952E+03
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 7.705E+03 6.393E+02 x_i, y_j = -5.614E+03 -3.186E+02
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 8.905E+03 6.393E+02 x_i, y_j = 4.880E+03 -2.527E+03
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 1.010E+04 6.393E+02 x_i, y_j = -1.544E+03 6.060E+03
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 1.130E+04 6.393E+02 x_i, y_j = -2.389E+03 -4.630E+03
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 5.050E+02 1.839E+03 x_i, y_j = 1.774E+02 1.098E+03
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 1.705E+03 1.839E+03 x_i, y_j = -4.058E+02 2.327E+02
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 2.905E+03 1.839E+03 x_i, y_j = 6.037E+02 4.979E+02
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 4.105E+03 1.839E+03 x_i, y_j = -3.295E+02 9.653E+02
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 6.505E+03 1.839E+03 x_i, y_j = 4.068E+02 9.835E+02
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 7.705E+03 1.839E+03 x_i, y_j = -5.022E+02 4.705E+02
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 8.905E+03 1.839E+03 x_i, y_j = 5.192E+02 2.556E+02
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 1.010E+04 1.839E+03 x_i, y_j = -1.060E+02 1.091E+03
*** ll_to_xyindex_wrf error -- x_i, y_j out of bounds x_lon, y_lat = 1.130E+04 1.839E+03 x_i, y_j = -1.883E+02 5.091E+01
Program received signal SIGSEGV: Segmentation fault - invalid memory reference.
Backtrace for this error: #0 0x2BA1CE5D4697 #1 0x2BA1CE5D4CDE #2 0x2BA1CF0673AF #3 0x484A90 in outgrid_init_reg_ #4 0x445A68 in MAIN__ at flexwrf.f90:? Segmentation fault (core dumped)
下面是我WRF-Chem的网格参数（namelist.wps）:&share wrf_core = 'ARW', max_dom = 1, start_date = '2017-06-02_00:00:00','2017-06-02_00:00:00','2017-06-05_00:00:00','2017-05-31_00:00:00', end_date = '2017-06-08_00:00:00','2017-06-16_00:00:00','2017-06-11_00:00:00','2017-06-31_06:00:00', interval_seconds = 21600, io_form_geogrid = 2, /
&geogrid ! parent_id = 0, 1, 2, 3, ! parent_grid_ratio = 1, 3, 3, 3, ! i_parent_start = 1, 35, 37, 25, ! j_parent_start = 1, 23, 21, 23, ! e_we = 100, 100, 85, 100, ! e_sn = 81, 82, 127, 250, parent_id = 0, 1, 2, 3, parent_grid_ratio = 1, 3, 3, 3, i_parent_start = 1, 35, 37, 25, j_parent_start = 1, 23, 21, 23, e_we = 100, 100, 85, 100, e_sn = 82, 82, 127, 250,
geog_data_res = '30s','30s','30s','30s', dx =12000, dy = 12000, map_proj = 'lambert', ref_lat = 44, ref_lon = -88, truelat1 = 30.0, truelat2 = 60.0, stand_lon = -88.0, geog_data_path = '/qw/feng/Build_WRF/WPS_GEOG/' /
&ungrib out_format = 'WPS', prefix = 'FILE', /
我困惑的地方是，WRF是设置模拟中心点的经纬度和网格格距，所以最后区域的边框经纬度并不统一，如下图所示：地图呈现北面经线收紧南面分散的特点，那么这样的网格在flexwrf.input里应该怎么设置呢？
| 
窥视卡
雷达卡
发表于 2020-5-1 10:01:34
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
楼主