- 积分
- 9000
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2013-3-14
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看，没有帐号？立即注册
x
本帖最后由 亽小笨蜑■ɑ 于 2014-7-15 16:55 编辑
我用wrf模拟了十天的数据，没有报错，ncdump -v Times wrfout。。。。。之后，每个时刻都有数据，经过后处理arwpost转换后，用GRADS画图，只能画出前8天的数据，后面就一直报错，感觉是GRADS的原因，一直解决不了，求高手帮忙： log中 d01 2013-11-12_00:00:00 wrf: SUCCESS COMPLETE WRF
ARWPOST.namelist
&datetime
start_date = '2013-11-02_00:00:00',
end_date = '2013-11-12_00:00:00',
interval_seconds = 3600,
tacc = 0,
debug_level = 0,
/
&io
io_form_input = 2,
input_root_name = '/home/programfile/WRF/wrfout/20131130/wrfout_d01_2013-11-02_00:00:00'
output_root_name = '/home//programfile/WRF/ARWpost/20131130'
plot = 'all_list'
fields = 'height,pressure,tk,tc,slp'
output_type = 'grads'
mercator_defs = .true.
/
split_output = .true.
frames_per_outfile = 2
output_type = 'grads'
output_type = 'v5d'
&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.
ctl文件
dset F:\20131130.dat
options byteswapped
undef 1.e30
title OUTPUT FROM WRF V3.3.1 MODEL
pdef 140 165 lcc 30.000 110.000 70.500 83.000 60.00000 30.00000 110.00000 30000.000 30000.000
xdef 468 linear 78.37687 0.13513514
ydef 353 linear 6.27682 0.13513514
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 241 linear 00Z02NOV2013 60MN
VARS 95
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)
QICE 19 0 Ice mixing ratio (kg kg-1)
QSNOW 19 0 Snow mixing ratio (kg kg-1)
QGRAUP 19 0 Graupel 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)
H_DIABATIC 19 0 MICROPHYSICS LATENT HEATING (K s-1)
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 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)
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 (-)
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)
AVG_FUEL_F 1 0 fuel remaining averaged to atmospheric grid (1)
UAH 1 0 wind at fire_wind_height (m/s)
VAH 1 0 wind at fire_wind_height (m/s)
pressure 19 0 Model pressure (hPa)
height 19 0 Model height (km)
tk 19 0 Temperature (K)
tc 19 0 Temperature (C)
slp 1 0 Sea Levelp Pressure (hPa)
ENDVARS
@ global String comment TITLE = OUTPUT FROM WRF V3.3.1 MODEL
@ global String comment SIMULATION_START_DATE = 2013-11-02_00:00:00
@ global String comment GRIDTYPE = C
@ global String comment MMINLU = USGS
@ global String comment WEST-EAST_GRID_DIMENSION = 141
@ global String comment SOUTH-NORTH_GRID_DIMENSION = 166
@ global String comment BOTTOM-TOP_GRID_DIMENSION = 28
@ global String comment MAP_PROJ = 1
@ global String comment DX = 30000.00
@ global String comment DY = 30000.00
@ global String comment CEN_LAT = 30.00
@ global String comment CEN_LON = 110.00
@ global String comment TRUELAT1 = 30.00
@ global String comment TRUELAT2 = 60.00
@ global String comment MOAD_CEN_LAT = 30.00
@ global String comment STAND_LON = 110.00
@ global String comment DIFF_OPT = 1
@ global String comment KM_OPT = 1
@ global String comment DAMP_OPT = 0
@ global String comment KHDIF = 0.00
@ global String comment KVDIF = 0.00
@ global String comment MP_PHYSICS = 2
@ 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 = 1
@ 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 = 1
@ global String comment PARENT_ID = 0
@ global String comment I_PARENT_START = 1
@ global String comment J_PARENT_START = 1
@ global String comment PARENT_GRID_RATIO = 1
@ global String comment DT = 60.00
@ global String comment ISWATER = 16
@ global String comment ISICE = 24
@ global String comment ISURBAN = 1
@ global String comment ISOILWATER = 14
GS文件
'reinit'
'open F:\20131130.ctl'
'set lon 98 126'
'set lat 21 49'
'set zlog on'
'set grads off'
'set grid off'
i=193
while(i<241)
'set t 'i''
'set map 1 1 10'
'set gxout shaded'
'set rgb 16 0 0 143'
'set rgb 18 0 0 207'
'set rgb 20 0 16 255'
'set rgb 24 0 128 255'
'set rgb 28 0 223 255'
'set rgb 32 48 255 207'
'set rgb 34 96 255 159'
'set rgb 36 143 255 112'
'set rgb 40 239 255 16'
'set rgb 42 255 239 0'
'set rgb 43 255 223 0'
'set rgb 44 255 207 0'
'set rgb 45 255 191 0'
'set rgb 46 255 175 0'
'set rgb 47 255 159 0'
'set rgb 48 255 143 0'
'set rgb 49 255 128 0'
'set rgb 50 255 112 0'
'set rgb 52 255 80 0'
'set rgb 53 255 64 0'
'set rgb 54 255 48 0'
'set rgb 55 255 32 0'
'set rgb 56 255 16 0'
'set rgb 57 255 0 0'
'set clevs 1008 1010 1012 1014 1016 1018 1020 1022 1024 1026 1028 1030'
'set ccols 16 20 24 28 32 36 40 43 46 50 54 56 57'
'd slp'
'cbarn'
'set gxout contour'
'd slp'
'draw title 2013.11.'math_int((i-1)/24)+10'.'math_mod(i-1,24)*1''
'printim F:\2013.11\2013.11.'math_int((i-1)/24)+10'.'math_mod(i-1,24)*1'.png x1000 y800 white'
'c'
i=i+1
endwhile
;
|
-
|
窥视卡
雷达卡
发表于 2014-7-15 16:32:43
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
楼主