爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 17876|回复: 13

【求助】wrfout使用ARWPOST转换后没有温度数据

[复制链接]

新浪微博达人勋

发表于 2013-8-30 19:55:54 | 显示全部楼层 |阅读模式

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

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

x
最近在学习wrf。在运行完之后得到了wrfout_d01,使用ARWPOST转换之后的ctl中没有显示温度数据,只有位温,求各位指教。
下面是我运行ARWPOST时的namelist(wrf是按照官网上wrf-chem练习1做的http://ruc.noaa.gov/wrf/WG11/wrf_tutorial_exercises_v35/exercise_1.html):
&datetime
start_date = '2010-07-14_00:00:00',
end_date   = '2010-07-16_00:00:00',
interval_seconds = 21600,
tacc = 0,
debug_level = 0,
/
&io
input_root_name = '/home/pub/model/wrf/WRFV3/test/em_real/wrfout_d01_2010-07-14_00:00:00'
output_root_name = './erod'
plot = 'all_list'
fields = 'height,pressure,tk,tc'
mercator_defs = .true.
/
split_output = .true.
frames_per_outfile = 2

plot = 'all'
plot = 'list'
plot = 'all_list'
! Below is a list of all available diagnostics
fields = 'height,geopt,theta,tc,tk,td,td2,rh,rh2,umet,vmet,pressure,u10m,v10m,wdir,wspd,wd10,ws10,slp,mcape,mcin,lcl,lfc,cape,cin,dbz,max_dbz,clfr'

&interp
interp_method = 0,
interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
/
extrapolate = .true.
interp_method = 0,     ! 0 is model levels, -1 is nice height levels, 1 is user specified pressure/height
interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
interp_levels = 0.25, 0.50, 0.75, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0, 20.0,

搜索了一下类似的问题,发现大家都说是修改field里面的变量就可以了,但是我的field里面是有温度的(应该是tc和tk吧)。
下面是转换完成之后的ctl变量部分:
VARS  107
LU_INDEX       1  0  LAND USE CATEGORY (-)
VAR_SSO        1  0  variance of subgrid-scale orography (m2)
LAP_HGT        1  0  Laplacian of orography (m)
U             30  0  x-wind component (m s-1)
V             30  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)
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)
LAI            1  0  Leaf area index (area/area)
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)
REFL_10CM     24  0  Radar reflectivity (lamda = 10 cm) (dBZ)
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)
XLAT           1  0  LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG          1  0  LONGITUDE, WEST IS NEGATIVE (degree_east)
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})
FLX4           1  0  sensible heat from canopy (W m{-2})
FVB            1  0  fraction of vegetation with snow below (-)
FBUR           1  0  fraction of vegetation covered by snow (-)
FGSN           1  0  fraction of ground covered by snow (-)
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)
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)
DUST_1        24  0  dust size bin 1: 0.5um effective radius (ug/kg-dryair)
DUST_2        24  0  dust size bin 2: 1.4um effective radius (ug/kg-dryair)
DUST_3        24  0  dust size bin 3: 2.4um effective radius (ug/kg-dryair)
DUST_4        24  0  dust size bin 4: 4.5um effective radius (ug/kg-dryair)
DUST_5        24  0  dust size bin 5: 8.0um effective radius (ug/kg-dryair)
LANDMASK       1  0  LAND MASK (1 FOR LAND, 0 FOR WATER) (-)
SST            1  0  SEA SURFACE TEMPERATURE (K)
ENDVARS
这107个废物里面没有温度
求各位前辈指点如何得到气温数据,还有为啥我的数据是24层的,别人都是27层啊{:soso_e141:},难道是wps时就设置的?

ps:帮忙解决问题的有奖励哈,拜托了各位。

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-8-30 20:20:39 | 显示全部楼层
tc是摄氏温度,tk是绝对温度吧。这个应该是你的namelist.input设定了25了吧?为什么一定要和别人的一样呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-8-31 13:02:30 | 显示全部楼层
nterp_method = 0,
把这一项改成1试试看
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-10-4 15:55:47 | 显示全部楼层
楼主试试直接'plot = list '看看输出的变量是不是4个?
还有mercator_defs = .true.,这个true是墨卡托投影的,是否需要修改?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-10-4 16:01:26 | 显示全部楼层

这个已经改好了,问题在于arwpost对于化学版本的wrfout不识别,需要修改其中某一模块中的参数,然后重新编译就好了。
感谢你的热心帮助。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-4-7 19:22:06 | 显示全部楼层
lqouc 发表于 2013-10-4 16:01
这个已经改好了,问题在于arwpost对于化学版本的wrfout不识别,需要修改其中某一模块中的参数,然后重新 ...

请问是修改的wrfchem的参数还是ARWpost的参数呢?我也遇到同样的问题
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-4-7 19:40:17 | 显示全部楼层
nxchenxin 发表于 2014-4-7 19:22
请问是修改的wrfchem的参数还是ARWpost的参数呢?我也遇到同样的问题

ARWpost的参数
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-4-9 23:38:07 | 显示全部楼层

请问要怎么改才会输出位温theta数据呢?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-10-17 18:40:44 | 显示全部楼层
请教版主:看到您帖子“wrfout使用ARWPOST转换后没有温度数据”,这个问题最后如何解决的?wrfchem的输出结果能否用arwpost转化为等高面或等压面?各层温度如何输出?期待您的解答。非常感谢!!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-10-20 11:22:45 | 显示全部楼层
tripitaka 发表于 2015-10-17 18:40
请教版主:看到您帖子“wrfout使用ARWPOST转换后没有温度数据”,这个问题最后如何解决的?wrfchem的输出结 ...

输出温度的方法就是在arwpost编译的时候的要修改一部分东西,但是当时是同学帮我改的,我也说不清应该怎么弄。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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