登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
求助:wrf资料画流线图,出图线条特别密集,已经调节strmden至10,还是没什么效果,调了半天也没发现解决办法,是数据本身的精度问题么?
附ctl,gs和图:
dset D:\531\wrf03d2014053100.dat
options byteswapped
undef 1.e30
title OUTPUT FROM WRF V3.4 MODEL
pdef 288 288 lcc 32.236 117.107 144.500 144.500 60.00000 30.00000 117.30000 3000.000 3000.000
xdef 744 linear 112.06458 0.01351351
ydef 606 linear 28.11144 0.01351351
zdef 17 levels
1000.00000
975.00000
950.00000
925.00000
900.00000
850.00000
800.00000
750.00000
700.00000
650.00000
600.00000
550.00000
500.00000
400.00000
300.00000
200.00000
100.00000
tdef 73 linear 00Z31May2013 60MN
VARS 26
U 17 0 x-wind component (m s-1)
V 17 0 y-wind component (m s-1)
W 17 0 z-wind component (m s-1)
Q2 1 0 QV at 2 M (kg kg-1)
T2 1 0 TEMP at 2 M (K)
QVAPOR 17 0 Water vapor mixing ratio (kg kg-1)
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)
XLAT 1 0 LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG 1 0 LONGITUDE, WEST IS NEGATIVE (degree_east)
height 17 0 Model height (km)
tc 17 0 Temperature (C)
theta 17 0 Potential Temperature (K)
td 17 0 Dewpoint Temperature (C)
rh 17 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 (%)
u10m 1 0 Rotated wind component (m s-1)
v10m 1 0 Rotated wind component (m s-1)
slp 1 0 Sea Levelp Pressure (hPa)
dbz 17 0 Reflectivity (-)
cape 17 0 CAPE (J/kg)
ENDVARS
GS文件部分内容如下:
'reinit'
'open D:\531\wrfoutd03.ctl'
'enable print D:\531\1\stream\850stream.gmf'
'set xlopts 1 7 0.17'
'set ylopts 1 7 0.17'
'set parea 1 10 1 8'
'set mpdset cnworld cnriver'
'set grid off'
'set grads off'
tt=1
while(tt<=2)
'set t ' tt
'set lev 850'
'set lat 28.1114 36.2871'
'set lon 112.065 122.105'
'set gxout stream'
'set cthick 1'
'set poli on'
'set strmden 1'
'd u;v'
'q dim'
a=sublin(result,5)
a=subwrd(a,6)
'draw title stream ' a
'print'
|