ARWpost插值的问题（无法插值到等压面层）

youxianchan

最近刚开始学习wrf,在使用wrf后处理模式ARWpost时遇到一些问题，就是无法将变量插值到等压面层。我的namelist.ARWpost如下:
&datetime
start_date = '2018-08-05_03:00:00',
end_date   = '2018-08-05_03:00:00',
interval_seconds = 10800,
tacc = 0,
debug_level = 0,
/
&io
input_root_name = './wrfout'
output_root_name = './test'
plot = 'all'
fields = 'height,pressure,tk,tc'
mercator_defs = .False.
/
split_output = .False.
frames_per_outfile = 2

&interp
interp_method = 1
interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,
/
extrapolate = .True.


其中， interp_method = 1 这里我已经改为1，而且也设置了interp_levels。但是我的输出test.ctl如下：
dset ^test.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM *             PROGRAM:WRF/CHEM V3.8.1 MODEL
pdef  222 201 lcc  39.500  114.800  111.500  101.000  60.00000  30.00000  114.80000  12000.000  12000.000
xdef  724 linear   95.21928   0.05405406
ydef  447 linear   27.02565   0.05405406
zdef   24 linear 1 1  
tdef    1 linear 03Z05AUG2018     180MN      
VARS  234
XLAT           1  0  LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG          1  0  LONGITUDE, WEST IS NEGATIVE (degree_east)
LU_INDEX       1  0  LAND USE CATEGORY (-)
VAR_SSO        1  0  variance of subgrid-scale orography (m2)
U             37  0  x-wind component (m s-1)
V             37  0  y-wind component (m s-1)
W             24  0  z-wind component (m s-1)
PH            24  0  perturbation geopotential (m2 s-2)
PHB           24  0  base-state geopotential (m2 s-2)
T             24  0  perturbation potential temperature (theta-t0) (K)
MU             1  0  perturbation dry air mass in column (Pa)
MUB            1  0  base state dry air mass in column (Pa)
NEST_POS       1  0  - (-)
P             24  0  perturbation pressure (Pa)
PB            24  0  BASE STATE PRESSURE (Pa)
P_HYD         24  0  hydrostatic pressure (Pa)
Q2             1  0  QV at 2 M (kg kg-1)
T2             1  0  TEMP at 2 M (K)
TH2            1  0  POT TEMP at 2 M (K)
PSFC           1  0  SFC PRESSURE (Pa)
U10            1  0  U at 10 M (m s-1)
V10            1  0  V at 10 M (m s-1)
QVAPOR        24  0  Water vapor mixing ratio (kg kg-1)
QCLOUD        24  0  Cloud water mixing ratio (kg kg-1)
QRAIN         24  0  Rain water mixing ratio (kg kg-1)
QICE          24  0  Ice mixing ratio (kg kg-1)
QSNOW         24  0  Snow mixing ratio (kg kg-1)
QGRAUP        24  0  Graupel mixing ratio (kg kg-1)
QNICE         24  0  Ice Number concentration (  kg-1)
QNRAIN        24  0  Rain Number concentration (  kg(-1))
SHDMAX         1  0  ANNUAL MAX VEG FRACTION (-)
SHDMIN         1  0  ANNUAL MIN VEG FRACTION (-)
SNOALB         1  0  ANNUAL MAX SNOW ALBEDO IN FRACTION (-)
TSLB           4  0  SOIL TEMPERATURE (K)
SMOIS          4  0  SOIL MOISTURE (m3 m-3)
SH2O           4  0  SOIL LIQUID WATER (m3 m-3)
SMCREL         4  0  RELATIVE SOIL MOISTURE (-)
SEAICE         1  0  SEA ICE FLAG (-)
XICEM          1  0  SEA ICE FLAG (PREVIOUS STEP) (-)
SFROFF         1  0  SURFACE RUNOFF (mm)
UDROFF         1  0  UNDERGROUND RUNOFF (mm)
IVGTYP         1  0  DOMINANT VEGETATION CATEGORY (-)
ISLTYP         1  0  DOMINANT SOIL CATEGORY (-)
VEGFRA         1  0  VEGETATION FRACTION (-)
GRDFLX         1  0  GROUND HEAT FLUX (W m-2)
ACGRDFLX       1  0  ACCUMULATED GROUND HEAT FLUX (J m-2)
ACSNOM         1  0  ACCUMULATED MELTED SNOW (kg m-2)
SNOW           1  0  SNOW WATER EQUIVALENT (kg m-2)
SNOWH          1  0  PHYSICAL SNOW DEPTH (m)
CANWAT         1  0  CANOPY WATER (kg m-2)
SSTSK          1  0  SKIN SEA SURFACE TEMPERATURE (K)
COSZEN         1  0  COS of SOLAR ZENITH ANGLE (dimensionless)
LAI            1  0  LEAF AREA INDEX (m-2/m-2)
VAR            1  0  OROGRAPHIC VARIANCE (-)
MAPFAC_M       1  0  Map scale factor on mass grid (-)
MAPFAC_MX      1  0  Map scale factor on mass grid, x direction (-)
MAPFAC_MY      1  0  Map scale factor on mass grid, y direction (-)
MF_VX_INV      1  0  Inverse map scale factor on v-grid, x direction (-)
F              1  0  Coriolis sine latitude term (s-1)
E              1  0  Coriolis cosine latitude term (s-1)
SINALPHA       1  0  Local sine of map rotation (-)
COSALPHA       1  0  Local cosine of map rotation (-)
HGT            1  0  Terrain Height (m)
TSK            1  0  SURFACE SKIN TEMPERATURE (K)
RAINC          1  0  ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
RAINSH         1  0  ACCUMULATED SHALLOW CUMULUS PRECIPITATION (mm)
RAINNC         1  0  ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
SNOWNC         1  0  ACCUMULATED TOTAL GRID SCALE SNOW AND ICE (mm)
GRAUPELNC      1  0  ACCUMULATED TOTAL GRID SCALE GRAUPEL (mm)
HAILNC         1  0  ACCUMULATED TOTAL GRID SCALE HAIL (mm)
CLDFRA        24  0  CLOUD FRACTION (-)
SWDOWN         1  0  DOWNWARD SHORT WAVE FLUX AT GROUND SURFACE (W m-2)
GLW            1  0  DOWNWARD LONG WAVE FLUX AT GROUND SURFACE (W m-2)
SWNORM         1  0  NORMAL SHORT WAVE FLUX AT GROUND SURFACE (SLOPE-DEPENDENT) (W m-2)
OLR            1  0  TOA OUTGOING LONG WAVE (W m-2)
ALBEDO         1  0  ALBEDO (-)
CLAT           1  0  COMPUTATIONAL GRID LATITUDE, SOUTH IS NEGATIVE (degree_north)
ALBBCK         1  0  BACKGROUND ALBEDO (-)
EMISS          1  0  SURFACE EMISSIVITY (-)
NOAHRES        1  0  RESIDUAL OF THE NOAH SURFACE ENERGY BUDGET (W m{-2})
TMN            1  0  SOIL TEMPERATURE AT LOWER BOUNDARY (K)
XLAND          1  0  LAND MASK (1 FOR LAND, 2 FOR WATER) (-)
UST            1  0  U* IN SIMILARITY THEORY (m s-1)
PBLH           1  0  PBL HEIGHT (m)
HFX            1  0  UPWARD HEAT FLUX AT THE SURFACE (W m-2)
QFX            1  0  UPWARD MOISTURE FLUX AT THE SURFACE (kg m-2 s-1)
LH             1  0  LATENT HEAT FLUX AT THE SURFACE (W m-2)
ACHFX          1  0  ACCUMULATED UPWARD HEAT FLUX AT THE SURFACE (J m-2)
ACLHF          1  0  ACCUMULATED UPWARD LATENT HEAT FLUX AT THE SURFACE (J m-2)
SNOWC          1  0  FLAG INDICATING SNOW COVERAGE (1 FOR SNOW COVER) (-)
SR             1  0  fraction of frozen precipitation (-)
EROD           3  0  fraction of erodible surface in each grid cell (0-1) (none)
CLDFRA2       24  0  CLOUD FRACTION (-)
RAINPROD      24  0  TOTAL RAIN PRODUCTION RATE (s-1)
EVAPPROD      24  0  RAIN EVAPORATION RATE (s-1)
UST_T          1  0  Threshold Friction Velocity (m s-1)
ROUGH_COR      1  0  roughness elements correction (-)
SMOIS_COR      1  0  soil moisture correction (-)
DUSTLOAD_1     1  0  dust loading for size 1 (ug/m2)
DUSTLOAD_2     1  0  dust loading for size 2 (ug/m2)
DUSTLOAD_3     1  0  dust loading for size 3 (ug/m2)
DUSTLOAD_4     1  0  dust loading for size 4 (ug/m2)
DUSTLOAD_5     1  0  total dust loading (ug/m2)
DRYDEPVEL      1  0  dust dry deposition velocity (m/s)
GRASET_1       1  0  dust gravitational settling  for size 1 (ug/m2/s)
GRASET_2       1  0  dust gravitational settling  for size 2 (ug/m2/s)
GRASET_3       1  0  dust gravitational settling  for size 3 (ug/m2/s)
GRASET_4       1  0  dust gravitational settling  for size 4 (ug/m2/s)
GRASET_5       1  0  dust gravitational settling  for size 5 (ug/m2/s)
DRYDEP_1       1  0  dust dry deposition for size 1 (ug/m2/s)
DRYDEP_2       1  0  dust dry deposition for size 2 (ug/m2/s)
DRYDEP_3       1  0  dust dry deposition for size 3 (ug/m2/s)
DRYDEP_4       1  0  dust dry deposition for size 4 (ug/m2/s)
DRYDEP_5       1  0  dust dry deposition for size 5 (ug/m2/s)
EDUST1         1  0  DUST emissions bin1 (-)
EDUST2         1  0  DUST emissions bin2 (-)
EDUST3         1  0  DUST emissions bin3 (-)
EDUST4         1  0  DUST emissions bin4 (-)
EDUST5         1  0  DUST emissions bin5 (-)
EXTCOF55      24  0  Extinction coefficients for .55um (km^-1)
DRY_DEP_LE    24  0  dry deposition velocity (cm/s)
DRY_DEP_FL    24  0  dry deposition flux (mol or ug m^-2)
WET_DEP_FL    24  0  column wet scavening flux (mmol or ug m^-2)
EBIO_ISO       1  0  Actual biog emiss (mol km^-2 hr^-1)
EBIO_API       1  0  Actual biog emiss (mol km^-2 hr^-1)
LAI_VEGMAS     1  0  MODIS LAI vegetation mask for this date; 0=no dust produced (vegetation) (none)
dvel_o3        1  0  O3 deposition velocity (cm/s)
PM2_5_DRY     24  0  pm2.5 aerosol dry mass (ug m^-3)
PM10          24  0  pm10 dry mass (ug m^-3)
DMS_0          1  0  dms oceanic concentrations (nM/L)
PHOTR201      24  0  cl2 photolysis rate (min{-1})
PHOTR202      24  0  hocl photolysis rate (min{-1})
PHOTR203      24  0  fmcl photolysis rate (min{-1})
so2           24  0  SO2 mixing ratio (ppmv)
sulf          24  0  SULF mixing ratio (ppmv)
no2           24  0  NO2 mixing ratio (ppmv)
no            24  0  NO mixing ratio (ppmv)
o3            24  0  O3 mixing ratio (ppmv)
hno3          24  0  HNO3 mixing ratio (ppmv)
h2o2          24  0  H2O2 mixing ratio (ppmv)
ald           24  0  ALD mixing ratio (ppmv)
hcho          24  0  HCHO mixing ratio (ppmv)
op1           24  0  OP1 mixing ratio (ppmv)
op2           24  0  OP2 mixing ratio (ppmv)
ora1          24  0  ORA1 mixing ratio (ppmv)
ora2          24  0  ORA2 mixing ratio (ppmv)
nh3           24  0  NH3 mixing ratio (ppmv)
n2o5          24  0  N2O5 mixing ratio (ppmv)
no3           24  0  NO3 mixing ratio (ppmv)
pan           24  0  PAN mixing ratio (ppmv)
eth           24  0  ETH mixing ratio (ppmv)
co            24  0  CO mixing ratio (ppmv)
ol2           24  0  OL2 mixing ratio (ppmv)
olt           24  0  OLT mixing ratio (ppmv)
oli           24  0  OLI mixing ratio (ppmv)
tol           24  0  TOL mixing ratio (ppmv)
xyl           24  0  XYL mixing ratio (ppmv)
hono          24  0  HONO mixing ratio (ppmv)
hno4          24  0  HNO4 mixing ratio (ppmv)
ket           24  0  KET mixing ratio (ppmv)
mgly          24  0  MGLY mixing ratio (ppmv)
onit          24  0  ONIT mixing ratio (ppmv)
csl           24  0  CSL mixing ratio (ppmv)
iso           24  0  ISO mixing ratio (ppmv)
ho            24  0  HO mixing ratio (ppmv)
ho2           24  0  HO2 mixing ratio (ppmv)
hcl           24  0  HCL mixing ratio (ppmv)
ch3o2         24  0  CH3O2 mixing ratio (ppmv)
ethp          24  0  ETHP mixing ratio (ppmv)
ch3oh         24  0  CH3OH mixing ratio (ppmv)
c2h5oh        24  0  C2H5OH mixing ratio (ppmv)
par           24  0  PAR mixing ratio (ppmv)
to2           24  0  TO2 mixing ratio (ppmv)
cro           24  0  CRO mixing ratio (ppmv)
open          24  0  OPEN mixing ratio (ppmv)
op3           24  0  OP3 mixing ratio (ppmv)
c2o3          24  0  C2O3 mixing ratio (ppmv)
ro2           24  0  RO2 mixing ratio (ppmv)
ano2          24  0  ANO2 mixing ratio (ppmv)
nap           24  0  NAP mixing ratio (ppmv)
xo2           24  0  XO2 mixing ratio (ppmv)
xpar          24  0  XPAR mixing ratio (ppmv)
isoprd        24  0  ISOPRD mixing ratio (ppmv)
isopp         24  0  ISOPP mixing ratio (ppmv)
isopn         24  0  ISOPN mixing ratio (ppmv)
isopo2        24  0  ISOPO2 mixing ratio (ppmv)
so4_a01       24  0  Sulfate, aerosol bin 01 (ug/kg-dryair)
no3_a01       24  0  Nitrate, aerosol bin 01 (ug/kg-dryair)
cl_a01        24  0  Chloride, aerosol bin 01 (ug/kg-dryair)
nh4_a01       24  0  Ammonium, aerosol bin 01 (ug/kg-dryair)
na_a01        24  0  Sodium, aerosol bin 01 (ug/kg-dryair)
oin_a01       24  0  Other inorganics, aerosol bin 01 (ug/kg-dryair)
oc_a01        24  0  Organic carbon, aerosol bin 01 (ug/kg-dryair)
bc_a01        24  0  Black carbon, aerosol bin 01 (ug/kg-dryair)
hysw_a01      24  0  Hysteresis water, aerosol bin 01 (ug/kg-dryair)
water_a01     24  0  Water, aerosol bin 01 (ug/kg-dryair)
num_a01       24  0  Number, aerosol bin 01 ( /kg-dryair)
so4_a02       24  0  Sulfate, aerosol bin 02 (ug/kg-dryair)
no3_a02       24  0  Nitrate, aerosol bin 02 (ug/kg-dryair)
cl_a02        24  0  Chloride, aerosol bin 02 (ug/kg-dryair)
nh4_a02       24  0  Ammonium, aerosol bin 02 (ug/kg-dryair)
na_a02        24  0  Sodium, aerosol bin 02 (ug/kg-dryair)
oin_a02       24  0  Other inorganics, aerosol bin 02 (ug/kg-dryair)
oc_a02        24  0  Organic carbon, aerosol bin 02 (ug/kg-dryair)
bc_a02        24  0  Black carbon, aerosol bin 02 (ug/kg-dryair)
hysw_a02      24  0  Hysteresis water, aerosol bin 02 (ug/kg-dryair)
water_a02     24  0  Water, aerosol bin 02 (ug/kg-dryair)
num_a02       24  0  Number, aerosol bin 02 ( /kg-dryair)
so4_a03       24  0  Sulfate, aerosol bin 03 (ug/kg-dryair)
no3_a03       24  0  Nitrate, aerosol bin 03 (ug/kg-dryair)
cl_a03        24  0  Chloride, aerosol bin 03 (ug/kg-dryair)
nh4_a03       24  0  Ammonium, aerosol bin 03 (ug/kg-dryair)
na_a03        24  0  Sodium, aerosol bin 03 (ug/kg-dryair)
oin_a03       24  0  Other inorganics, aerosol bin 03 (ug/kg-dryair)
oc_a03        24  0  Organic carbon, aerosol bin 03 (ug/kg-dryair)
bc_a03        24  0  Black carbon, aerosol bin 03 (ug/kg-dryair)
hysw_a03      24  0  Hysteresis water, aerosol bin 03 (ug/kg-dryair)
water_a03     24  0  Water, aerosol bin 03 (ug/kg-dryair)
num_a03       24  0  Number, aerosol bin 03 ( /kg-dryair)
so4_a04       24  0  Sulfate, aerosol bin 04 (ug/kg-dryair)
no3_a04       24  0  Nitrate, aerosol bin 04 (ug/kg-dryair)
cl_a04        24  0  Chloride, aerosol bin 04 (ug/kg-dryair)
nh4_a04       24  0  Ammonium, aerosol bin 04 (ug/kg-dryair)
na_a04        24  0  Sodium, aerosol bin 04 (ug/kg-dryair)
oin_a04       24  0  Other inorganics, aerosol bin 04 (ug/kg-dryair)
oc_a04        24  0  Organic carbon, aerosol bin 04 (ug/kg-dryair)
bc_a04        24  0  Black carbon, aerosol bin 04 (ug/kg-dryair)
hysw_a04      24  0  Hysteresis water, aerosol bin 04 (ug/kg-dryair)
water_a04     24  0  Water, aerosol bin 04 (ug/kg-dryair)
num_a04       24  0  Number, aerosol bin 04 ( /kg-dryair)
LANDMASK       1  0  LAND MASK (1 FOR LAND, 0 FOR WATER) (-)
LAKEMASK       1  0  LAKE MASK (1 FOR LAND, 0 FOR WATER) (-)
SST            1  0  SEA SURFACE TEMPERATURE (K)
SST_INPUT      1  0  SEA SURFACE TEMPERATURE FROM WRFLOWINPUT FILE (K)
ENDVARS
@ global String comment TITLE =  OUTPUT FROM *             PROGRAM:WRF/CHEM V3.8.1 MODEL
以下ctl文件中的内容省略。。。
。。。。。。。
请问是否是我的namelist设置的不对或者哪里写的不对？
看到之前（很多年前）的帖子中提到ARWpost对化学模式的结果识别的不好，需要修改一些参数，不知道是不是和这个有关系，如果有请问怎么修改？

youxianchan

灰溜溜的回来自问自答。
果然是因为ARWpost对化学模式不识别。在运行的过程中有一个warning被我忽略了：
   WARNING --- I do not recognize this data.
               OUTPUT FROM *             PROGRAM:WRF/CHEM V3.8.1 MODEL
               Will make an attempt to read it.
查找了一下，发现需要修改src/input_module.f90 。在大概194行左右，增加一行内容：
IF ( INDEX(title,'OUTPUT FROM *             PROGRAM:WRF/CHEM V3.8.1 MODEL') /= 0 )   iprogram = 8
这样就可以识别这个结果了。

通过这件事我们可以得到如下结论：
1、多思考再来提问
2、不要忽视warning
红红火火恍恍惚惚哈哈哈哈哈哈

chenli

感谢楼主，遇到同样问题，已参考解决。

qqluckqq

改了以后没有WARNING了
但显示  Segmentation fault (core dumped)

z.biye

大家都能成功，为啥我改了后还是一样的警告