爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 16695|回复: 1

[脚本编辑] (已解决)请教关于WRF模拟生成的数据使用grads绘图问题

[复制链接]

新浪微博达人勋

发表于 2020-12-28 22:37:21 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 hbl1021 于 2020-12-29 21:33 编辑

    最近开始研究WRF数值模拟,摸索着把模拟数据生成了,一个ctl文件一个dat文件,但是使用grads绘图时出现了问题:使用的是NCEP1°*1°再分析资料,里面的变量如hgtprs、tmpprs等等在WRF模拟后名字都变了,找不到应该d哪个变量了,而且我试着display一下ctl里面其他的变量,都没有出图成功,请问从哪里能找到WRF模拟后用grads绘图的相关解读或者教程呢?    本人没有编程基础,使用grads也只是会基础绘图,不会fortran,对WRF模拟也理解不够深入,如果是基础知识不足的话还得麻烦大神们耐心解答,感谢!    问题已经解决,是跑WRF模拟时中心点没设置好,用ARWpost转换时变量和高度层没有设置,现在已经能够正常出图。

ctl文件:
dset D:\grads\fnl2020\test\test02.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V4.0 MODEL
pdef   73  60 lcc  40.010  117.880   37.000   30.500  60.00000  30.00000  -98.00000  30000.000  30000.000
xdef  219 linear  101.33707   0.13513513
ydef  189 linear   27.00921   0.13513513
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    2 linear 06Z26JUL2020     360MN      
VARS  149
XLAT           1  0  LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG          1  0  LONGITUDE, WEST IS NEGATIVE (degree_east)
LU_INDEX       1  0  LAND USE CATEGORY (-)
VAR_SSO        1  0  variance of subgrid-scale orography (m2)
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)
THM           19  0  either 1) pert moist pot temp=(1+Rv/Rd Qv)*(theta)-T0, or 2) pert dry pot temp=t (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)
QNICE         19  0  Ice Number concentration (  kg-1)
QNRAIN        19  0  Rain Number concentration (  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)
COSZEN         1  0  COS of SOLAR ZENITH ANGLE (dimensionless)
LAI            1  0  LEAF AREA INDEX (m-2/m-2)
DTAUX3D       19  0  LOCAL U GWDO STRESS (m s-1)
DTAUY3D       19  0  LOCAL V GWDO STRESS (m s-1)
DUSFCG         1  0  COLUMN-INTEGRATED U GWDO STRESS (Pa m s-1)
DVSFCG         1  0  COLUMN-INTEGRATED V GWDO STRESS (Pa m s-1)
VAR            1  0  OROGRAPHIC VARIANCE (-)
CON            1  0  OROGRAPHIC CONVEXITY (-)
OA1            1  0  OROGRAPHIC DIRECTION ASYMMETRY FUNCTION (-)
OA2            1  0  OROGRAPHIC DIRECTION ASYMMETRY FUNCTION (-)
OA3            1  0  OROGRAPHIC DIRECTION ASYMMETRY FUNCTION (-)
OA4            1  0  OROGRAPHIC DIRECTION ASYMMETRY FUNCTION (-)
OL1            1  0  OROGRAPHIC DIRECTION ASYMMETRY FUNCTION (-)
OL2            1  0  OROGRAPHIC DIRECTION ASYMMETRY FUNCTION (-)
OL3            1  0  OROGRAPHIC DIRECTION ASYMMETRY FUNCTION (-)
OL4            1  0  OROGRAPHIC DIRECTION ASYMMETRY FUNCTION (-)
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)
ACSWUPT        1  0  ACCUMULATED UPWELLING SHORTWAVE FLUX AT TOP (J m-2)
ACSWUPTC       1  0  ACCUMULATED UPWELLING CLEAR SKY SHORTWAVE FLUX AT TOP (J m-2)
ACSWDNT        1  0  ACCUMULATED DOWNWELLING SHORTWAVE FLUX AT TOP (J m-2)
ACSWDNTC       1  0  ACCUMULATED DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT TOP (J m-2)
ACSWUPB        1  0  ACCUMULATED UPWELLING SHORTWAVE FLUX AT BOTTOM (J m-2)
ACSWUPBC       1  0  ACCUMULATED UPWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM (J m-2)
ACSWDNB        1  0  ACCUMULATED DOWNWELLING SHORTWAVE FLUX AT BOTTOM (J m-2)
ACSWDNBC       1  0  ACCUMULATED DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM (J m-2)
ACLWUPT        1  0  ACCUMULATED UPWELLING LONGWAVE FLUX AT TOP (J m-2)
ACLWUPTC       1  0  ACCUMULATED UPWELLING CLEAR SKY LONGWAVE FLUX AT TOP (J m-2)
ACLWDNT        1  0  ACCUMULATED DOWNWELLING LONGWAVE FLUX AT TOP (J m-2)
ACLWDNTC       1  0  ACCUMULATED DOWNWELLING CLEAR SKY LONGWAVE FLUX AT TOP (J m-2)
ACLWUPB        1  0  ACCUMULATED UPWELLING LONGWAVE FLUX AT BOTTOM (J m-2)
ACLWUPBC       1  0  ACCUMULATED UPWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM (J m-2)
ACLWDNB        1  0  ACCUMULATED DOWNWELLING LONGWAVE FLUX AT BOTTOM (J m-2)
ACLWDNBC       1  0  ACCUMULATED DOWNWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM (J m-2)
SWUPT          1  0  INSTANTANEOUS UPWELLING SHORTWAVE FLUX AT TOP (W m-2)
SWUPTC         1  0  INSTANTANEOUS UPWELLING CLEAR SKY SHORTWAVE FLUX AT TOP (W m-2)
SWDNT          1  0  INSTANTANEOUS DOWNWELLING SHORTWAVE FLUX AT TOP (W m-2)
SWDNTC         1  0  INSTANTANEOUS DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT TOP (W m-2)
SWUPB          1  0  INSTANTANEOUS UPWELLING SHORTWAVE FLUX AT BOTTOM (W m-2)
SWUPBC         1  0  INSTANTANEOUS UPWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM (W m-2)
SWDNB          1  0  INSTANTANEOUS DOWNWELLING SHORTWAVE FLUX AT BOTTOM (W m-2)
SWDNBC         1  0  INSTANTANEOUS DOWNWELLING CLEAR SKY SHORTWAVE FLUX AT BOTTOM (W m-2)
LWUPT          1  0  INSTANTANEOUS UPWELLING LONGWAVE FLUX AT TOP (W m-2)
LWUPTC         1  0  INSTANTANEOUS UPWELLING CLEAR SKY LONGWAVE FLUX AT TOP (W m-2)
LWDNT          1  0  INSTANTANEOUS DOWNWELLING LONGWAVE FLUX AT TOP (W m-2)
LWDNTC         1  0  INSTANTANEOUS DOWNWELLING CLEAR SKY LONGWAVE FLUX AT TOP (W m-2)
LWUPB          1  0  INSTANTANEOUS UPWELLING LONGWAVE FLUX AT BOTTOM (W m-2)
LWUPBC         1  0  INSTANTANEOUS UPWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM (W m-2)
LWDNB          1  0  INSTANTANEOUS DOWNWELLING LONGWAVE FLUX AT BOTTOM (W m-2)
LWDNBC         1  0  INSTANTANEOUS DOWNWELLING CLEAR SKY LONGWAVE FLUX AT BOTTOM (W m-2)
OLR            1  0  TOA OUTGOING LONG WAVE (W m-2)
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 (-)
PCB            1  0  base state dry air mass in column (Pa)
PC             1  0  perturbation dry air mass in column (Pa)
LANDMASK       1  0  LAND MASK (1 FOR LAND, 0 FOR WATER) (-)
LAKEMASK       1  0  LAKE MASK (1 FOR LAKE, 0 FOR NON-LAKE) (-)
SST            1  0  SEA SURFACE TEMPERATURE (K)
SST_INPUT      1  0  SEA SURFACE TEMPERATURE FROM WRFLOWINPUT FILE (K)
pressure      19  0  Model pressure (hPa)
height        19  0  Model height (km)
tk            19  0  Temperature (K)
tc            19  0  Temperature (C)
ENDVARS



我的gs文件:
'reinit'
'open d:\grads\fnl2020\test\test02.ctl'


'set t 1'
'd hgt'



无标题.png




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

新浪微博达人勋

发表于 2021-8-26 21:51:02 | 显示全部楼层
感谢分享!!
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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