爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 5969|回复: 6

[分享资料] 请教:用Grads直接读取WRF模式原始数据问题

[复制链接]
发表于 2014-4-4 13:13:49 | 显示全部楼层 |阅读模式

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

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

x
之前在学校的时候一直使用ARWpost将模式数据转为dat和ctl文件后再用grads画图
现在毕业了想重新处理一下以前的数据
发现之前用ARWpost输出数据的时间间隔比较大(1小时输出一次) 现在想用10分钟的数据(模式原始数据的时间分辨率是10s)
苦于无赖  现在没有linux环境  不知道怎么转dat 和ctl  文件

上网搜了一下  说grads也可以直接处理WRF原始数据,诸如  wrfout_d03_2013-04-17_12_00_00
现在出现了一些  问题
不管用什么变量画图  画出来都是常数0

——————————————————————————————————————
现将相关ctl文件贴上来 请教各位高手 错误出在哪里了?

CTL文件:
dset H:\20130417\wrfout_d03_2013-04-17_12_00_00
options  sequential big_endian
undef 1.e30
title  OUTPUT FROM WRF V3.4 MODEL
pdef  360 291 lcc  28.083  113.182  180.500  146.000  28.00000  28.00000  115.00000   2000.000   2000.000
xdef  852 linear  109.29938     0.01
ydef  610 linear  25.29202   0.00900901  0.01
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   7 linear 12Z17APR2013      10MN      
VARS  118
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             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)
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)
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)
TKE_PBL       19  0  TKE from PBL (m2 s-2)
EL_PBL        19  0  Length scale from PBL (m)
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 (-)
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 (-)
LANDMASK       1  0  LAND MASK (1 FOR LAND, 0 FOR WATER) (-)
SST            1  0  SEA SURFACE TEMPERATURE (K)
pressure      19  0  Model pressure (hPa)
height        19  0  Model height (km)
tk            19  0  Temperature (K)
tc            19  0  Temperature (C)
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)
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)
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           19  0  Reflectivity (-)
max_dbz        1  0  Max Reflectivity (-)
cape          19  0  CAPE (J/kg)
cin           19  0  CIN (J/kg)
mcape          1  0  MCAPE (J/kg)
mcin           1  0  MCIN (J/kg)
lcl            1  0  LCL (meters AGL)
lfc            1  0  LFC (meters AGL)
ENDVARS

——————————————————————————————————————————————————————————
接下来直接在grads窗口中输入命令,出现一下结果: 1.png
2.png


密码修改失败请联系微信:mofangbao
发表于 2014-4-4 23:06:01 | 显示全部楼层
把options那行去掉看看对不对
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-4-8 09:02:17 | 显示全部楼层
燕雀要变鸿鹄 发表于 2014-4-4 23:06
把options那行去掉看看对不对

去掉试了 还是出现相同的问题
我个人猜测是不是插值的问题
密码修改失败请联系微信:mofangbao
发表于 2014-4-8 14:46:28 | 显示全部楼层
繁星一点 发表于 2014-4-8 09:02
去掉试了 还是出现相同的问题
我个人猜测是不是插值的问题

我记得当时我也出现了这个问题,然后我就改动了options  sequential big_endian
undef 1.e30这两行,就没问题了,不知道你这个是什么情况
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-4-8 19:43:37 | 显示全部楼层
燕雀要变鸿鹄 发表于 2014-4-8 14:46
我记得当时我也出现了这个问题,然后我就改动了options  sequential big_endian
undef 1.e30这两行,就 ...

好的 谢谢 我再看看有没有其他方法
密码修改失败请联系微信:mofangbao
发表于 2015-7-11 18:29:18 | 显示全部楼层
xdef  852 linear  109.29938     0.01
ydef  610 linear  25.29202   0.00900901  0.01
不是说不是等网格点的吗    为什么还可以这样表示
密码修改失败请联系微信:mofangbao
发表于 2017-5-3 10:31:54 | 显示全部楼层
你好,请问一下,你的ctl文件里面的u和v是怎么产生的呢,为什么我的ctl里面就没有呢
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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