- 积分
- 729
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-7-13
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看，没有帐号？立即注册
x
本帖最后由 油麦面疙瘩 于 2014-4-4 14:49 编辑
今早又全部试了一遍，我发现问题所在了，就是我模式模拟出来的所有结果里面，除了rainc和rainnc在低层有值以外。其他的输出变量，主要是风变量在低层都没有值，随着高度的升高，有值的区域逐渐扩大。。。
我不知道是什么原因导致了这样。用xmanager看了一下，wrfout文件是正常的，低层也是有值的。估计就是后处理出了问题。但是按理来说，降雨量的值都有，都是对的，整个文件也应该是对的才对啊。
急求大神帮助~~~~急哭了都要。
以下是我的namelist.ARWpost文件
&datetime
start_date = '2012-06-04_12:00:00',
end_date = '2012-06-05_12:00:00',
interval_seconds = 3600,
tacc = 0,
debug_level = 0,
/
&io
io_form_input = 2,
input_root_name = '/public1/home/CAS/lzhpc36/WRFV3/run/bmj-wsm3/wrfout_d02*'
output_root_name = '/public1/home/CAS/lzhpc36/dataout/bmj-wsm3/wrfout-bmj-wsm3_d02'
plot = 'all_list'
fields = 'height,pressure,tk,geopt,theta,tc,td,td2,rh,rh2,clfr,umet,vmet,wdir,wspd,wd10,ws10'
output_type = 'grads'
mercator_defs = .true.
/
split_output = .False.
frames_per_outfile = 1
output_type = 'grads'
Noutput_type = 'v5d'
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 = 1,
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,
最近在用wrf做暴雨的数值模拟，打算画一下垂直速度，涡度，散度的垂直剖面图。先用垂直速度试了试，结果出来图上都是0.不知道怎么回事，求大神帮我看看，看问题出在哪里。
以下是我的gs文件。
'reinit'
'open f:\grib1\wrfout\bmj-ferrier\wrfout-bmj-ferrier_d02.ctl'
'enable print f:\grib1\wrfout\hb-bmj-ferrier-w02.gmf'
'set grid off'
'set grads off'
'set font 1'
'set mpdset hires'
'set xlint 1'
'set ylint 100'
'set xlopts 1 4 0.18'
'set ylopts 1 4 0.18'
'set clopts 1 -1 0.15'
'set lon 90 110'
'set lat 30 45'
'set t 6'
'set z 1 17'
'set lon 92 101'
'set lat 40'
'set clab on'
'set zlog on'
'set clskip 2'
'set cint 0.3'
'set ccolor 1'
'set cthick 6'
'set csmooth on'
'd W'
'print'
'disable print'
W用的就直接是wrfout里的输出变量，垂直风分量，应该就是垂直速度的意思吧。。。
虽然问题比较小白，但是自己也是搜了很多，也没找到类似的问题或者解决办法，只好求助了。希望大家帮帮我。
ctl文件如下：
dset f:\grib1\wrfout\bmj-ferrier\wrfout-bmj-ferrier_d02.dat
options byteswapped
undef 1.e30
title OUTPUT FROM WRF V3.2.1 MODEL
pdef 84 66 lcc 39.249 94.577 42.500 33.500 60.00000 30.00000 94.90000 18000.000 18000.000
xdef 260 linear 84.00488 0.08108108
ydef 158 linear 32.68425 0.08108108
zdef 19 levels
1000.00000
950.00000
900.00000
850.00000
800.00000
750.00000
700.00000
650.00000
600.00000
550.00000
500.00000
450.00000
400.00000
350.00000
300.00000
250.00000
200.00000
150.00000
100.00000
tdef 25 linear 12Z04JUN2012 60MN
VARS 105
LU_INDEX 1 0 LAND USE CATEGORY (-)
U 19 0 x-wind component (m s-1)
V 19 0 y-wind component (m s-1)
W 19 0 z-wind component (m s-1)
PH 19 0 perturbation geopotential (m2 s-2)
PHB 19 0 base-state geopotential (m2 s-2)
T 19 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 19 0 perturbation pressure (Pa)
PB 19 0 BASE STATE PRESSURE (Pa)
P_HYD 19 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 19 0 Water vapor mixing ratio (kg kg-1)
QCLOUD 19 0 Cloud water mixing ratio (kg kg-1)
QRAIN 19 0 Rain water mixing ratio (kg kg-1)
QSNOW 19 0 Snow mixing ratio (kg kg-1)
CWM 19 0 Total condensate mixing ratio (kg kg-1)
LANDMASK 1 0 LAND MASK (1 FOR LAND, 0 FOR WATER) (-)
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)
SNOW 1 0 SNOW WATER EQUIVALENT (kg m-2)
SNOWH 1 0 PHYSICAL SNOW DEPTH (m)
CANWAT 1 0 CANOPY WATER (kg m-2)
SST 1 0 SEA SURFACE TEMPERATURE (K)
SSTSK 1 0 SKIN SEA SURFACE TEMPERATURE (K)
LAI 1 0 Leaf area index (area/area)
F_ICE_PHY 19 0 FRACTION OF ICE (-)
F_RAIN_PHY 19 0 FRACTION OF RAIN (-)
F_RIMEF_PH 19 0 MASS RATIO OF RIMED ICE (-)
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)
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 19 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)
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 (-)
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 (-)
POTEVP 1 0 accumulated potential evaporation (W m-2)
SNOPCX 1 0 snow phase change heat flux (W m-2)
SOILTB 1 0 bottom soil temperature (K)
pressure 19 0 Model pressure (hPa)
geopt 19 0 Geopotential (m2/s2)
height 19 0 Model height (km)
tk 19 0 Temperature (K)
tc 19 0 Temperature (C)
theta 19 0 Potential Temperature (K)
td 19 0 Dewpoint Temperature (C)
td2 1 0 Dewpoint Temperature at 2m (C)
rh 19 0 Relative Humidity (%)
clflo 1 0 Low Cloud Fraction (%)
clfmi 1 0 Mid Cloud Fraction (%)
clfhi 1 0 High Cloud Fraction (%)
rh2 1 0 Relative Humidity at 2m (%)
wspd 19 0 Wind Speed (m s-1)
wdir 19 0 Wind Direction (Degrees)
ws10 1 0 Wind Speed at 10 M (m s-1)
wd10 1 0 Wind Direction at 10 M (Degrees)
umet 19 0 Rotated wind component (m s-1)
vmet 19 0 Rotated wind component (m s-1)
ENDVARS
@ global String comment TITLE = OUTPUT FROM WRF V3.2.1 MODEL
@ global String comment SIMULATION_START_DATE = 2012-06-03_18:00:00
@ global String comment GRIDTYPE = C
@ global String comment MMINLU = USGS
@ global String comment WEST-EAST_GRID_DIMENSION = 85
@ global String comment SOUTH-NORTH_GRID_DIMENSION = 67
@ global String comment BOTTOM-TOP_GRID_DIMENSION = 28
@ global String comment MAP_PROJ = 1
@ global String comment DX = 18000.00
@ global String comment DY = 18000.00
@ global String comment CEN_LAT = 39.25
@ global String comment CEN_LON = 94.58
@ global String comment TRUELAT1 = 30.00
@ global String comment TRUELAT2 = 60.00
@ global String comment MOAD_CEN_LAT = 39.50
@ global String comment STAND_LON = 94.90
@ global String comment DIFF_OPT = 1
@ global String comment KM_OPT = 4
@ global String comment DAMP_OPT = 0
@ global String comment KHDIF = 0.00
@ global String comment KVDIF = 0.00
@ global String comment MP_PHYSICS = 5
@ global String comment RA_LW_PHYSICS = 1
@ global String comment RA_SW_PHYSICS = 1
@ global String comment SF_SFCLAY_PHYSICS = 1
@ global String comment SF_SURFACE_PHYSICS = 2
@ global String comment BL_PBL_PHYSICS = 1
@ global String comment CU_PHYSICS = 2
@ global String comment SURFACE_INPUT_SOURCE = 1
@ global String comment SST_UPDATE = 0
@ global String comment GRID_FDDA = 0
@ global String comment FEEDBACK = 1
@ global String comment SMOOTH_OPTION = 0
@ global String comment W_DAMPING = 0
@ global String comment OBS_NUDGE_OPT = 0
@ global String comment GRID_ID = 2
@ global String comment PARENT_ID = 1
@ global String comment I_PARENT_START = 32
@ global String comment J_PARENT_START = 17
@ global String comment PARENT_GRID_RATIO = 3
@ global String comment DT = 100.00
@ global String comment ISWATER = 16
@ global String comment ISICE = 24
@ global String comment ISURBAN = 1
@ global String comment ISOILWATER = 14
|
|
窥视卡
雷达卡
发表于 2014-4-3 15:33:58
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
发表于 2014-4-3 15:43:41
楼主