爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 21998|回复: 16

[脚本编辑] 边界层高度grads脚本求。急急急

[复制链接]

新浪微博达人勋

发表于 2017-10-12 06:33:33 | 显示全部楼层 |阅读模式

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

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

x
需要画这样的边界层高度图 请大神帮忙指点!
图片1.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-10-12 06:34:28 | 显示全部楼层
ncl也可以{:eb302:}{:eb302:}{:eb302:}
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-10-12 07:02:48 | 显示全部楼层
主要问题是 不知道如何画出这种类型的线 还有就是不知道如何将这几条线画在同一张图上{:eb302:}{:eb302:}
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-10-12 07:49:35 | 显示全部楼层
本帖最后由 tangjie199365 于 2017-10-15 06:34 编辑

进展 可以画出这样的图了 线画不出来
这是我的脚本'reinit'
'open /nuist/scratch/tj199365/wrfout/20091109/BL/test02z.ctl'
'set grads off'
'set grid off'
'set lon 119.2'
'set lat 30'

'set t 1 13'

'set z 1 10'


'set font 4'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'
'set csmooth on'
'set ylevs 0.1 0.3 0.4 0.6 0.9 1.2 1.7 2.2'


'set ccolor 2'
'set cstyle 1'
'set cmark 2'
*'set digsiz 0.12'
'set cthick 6'

'd PBLH'



'cbar_line -x 6.1 -y 7.7 -c 2 -m 2  -t "PBLH" '

'q dims'
rec=sublin(result,5)
recs=subwrd(rec,6)
'draw title PBLH'


'printim /nuist/scratch/tj199365/wrfout/20091109/BL/PBLH/pblhd2-'i'.png white'


pblhd2-i.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-10-12 08:27:06 | 显示全部楼层
把你的描述文件也放上来
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-10-14 15:53:55 | 显示全部楼层
river 发表于 2017-10-12 08:27
把你的描述文件也放上来

dset ^test02z.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.7.1 MODEL
pdef  261 279 lcc  29.761  119.552  131.000  140.000  60.00000  30.00000  119.80000   2000.000   2000.000
xdef  639 linear  116.66232   0.00900901
ydef  579 linear   27.14345   0.00900901
zdef   29 levels  
   0.04968
   0.09935
   0.19540
   0.31781
   0.47218
   0.66928
   0.91171
   1.27063
   1.73506
   2.22069
   2.72955
   3.47525
   4.44399
   5.41083
   6.37549
   7.33724
   8.29562
   9.25248
  10.21460
  11.18506
  12.16278
  13.15274
  14.14656
  15.13084
  16.10859
  17.08459
  18.06335
  19.05052
  20.04932
tdef   13 linear 00Z09NOV2009      1hr      
VARS  136
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)
LAP_HGT        1  0  Laplacian of orography (m)
U             29  0  x-wind component (m s-1)
V             29  0  y-wind component (m s-1)
W             29  0  z-wind component (m s-1)
PH            29  0  perturbation geopotential (m2 s-2)
PHB           29  0  base-state geopotential (m2 s-2)
T             29  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             29  0  perturbation pressure (Pa)
PB            29  0  BASE STATE PRESSURE (Pa)
P_HYD         29  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        29  0  Water vapor mixing ratio (kg kg-1)
QCLOUD        29  0  Cloud water mixing ratio (kg kg-1)
QRAIN         29  0  Rain water mixing ratio (kg kg-1)
QICE          29  0  Ice mixing ratio (kg kg-1)
QSNOW         29  0  Snow mixing ratio (kg kg-1)
QGRAUP        29  0  Graupel mixing ratio (kg kg-1)
QHAIL         29  0  Hail mixing ratio (kg kg-1)
QNCLOUD       29  0  cloud water Number concentration (  kg(-1))
QNDROP        29  0  Droplet number mixing ratio (  kg-1)
QNRAIN        29  0  Rain Number concentration (  kg(-1))
QNICE         29  0  Ice Number concentration (  kg-1)
QNSNOW        29  0  Snow Number concentration (  kg(-1))
QNGRAUPEL     29  0  Graupel Number concentration (  kg(-1))
QNHAIL        29  0  Hail Number concentration (  kg(-1))
QVGRAUPEL     29  0  Graupel Particle Volume (m(3) kg(-1))
QVHAIL        29  0  Hail Particle Volume (m(3) 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)
VAR            1  0  OROGRAPHIC VARIANCE (-)
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)
REFL_10CM     29  0  Radar reflectivity (lamda = 10 cm) (dBZ)
CLDFRA        29  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)
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 (-)
REFD_MAX       1  0  MAX DERIVED RADAR REFL (dbZ)
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      29  0  Model pressure (hPa)
geopt         29  0  Geopotential (m2/s2)
height        29  0  Model height (km)
tk            29  0  Temperature (K)
tc            29  0  Temperature (C)
theta         29  0  Potential Temperature (K)
td            29  0  Dewpoint Temperature (C)
td2            1  0  Dewpoint Temperature at 2m (C)
rh            29  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          29  0  Wind Speed (m s-1)
wdir          29  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          29  0  Rotated wind component (m s-1)
vmet          29  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           29  0  Reflectivity (-)
max_dbz        1  0  Max Reflectivity (-)
cape          29  0  CAPE (J/kg)
cin           29  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
@ global String comment TITLE =  OUTPUT FROM WRF V3.7.1 MODEL
@ global String comment START_DATE = 2009-11-09_00:00:00
@ global String comment SIMULATION_START_DATE = 2009-11-09_00:00:00
@ global String comment WEST-EAST_GRID_DIMENSION =   262
@ global String comment SOUTH-NORTH_GRID_DIMENSION =   280
@ global String comment BOTTOM-TOP_GRID_DIMENSION =    30
@ global String comment DX =      2000.00
@ global String comment DY =      2000.00
@ global String comment SKEBS_ON =     0
@ global String comment SPEC_BDY_FINAL_MU =     0
@ global String comment USE_Q_DIABATIC =     0
@ global String comment GRIDTYPE = C
@ global String comment DIFF_OPT =     1
@ global String comment KM_OPT =     4
@ global String comment DAMP_OPT =     0
@ global String comment DAMPCOEF =         0.20
@ global String comment KHDIF =         0.00
@ global String comment KVDIF =         0.00
@ global String comment MP_PHYSICS =    17
@ 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 =     0
@ global String comment SF_LAKE_PHYSICS =     0
@ global String comment SURFACE_INPUT_SOURCE =     1
@ global String comment SST_UPDATE =     0
@ global String comment GRID_FDDA =     0
@ global String comment GFDDA_INTERVAL_M =     0
@ global String comment GFDDA_END_H =     0
@ global String comment GRID_SFDDA =     0
@ global String comment SGFDDA_INTERVAL_M =     0
@ global String comment SGFDDA_END_H =     0
@ global String comment HYPSOMETRIC_OPT =     2
@ global String comment USE_THETA_M =     0
@ global String comment SF_URBAN_PHYSICS =     0
@ global String comment SHCU_PHYSICS =     0
@ global String comment MFSHCONV =     0
@ global String comment FEEDBACK =     1
@ global String comment SMOOTH_OPTION =     0
@ global String comment SWRAD_SCAT =         1.00
@ global String comment W_DAMPING =     0
@ global String comment DT =        12.00
@ global String comment RADT =         2.00
@ global String comment BLDT =         0.00
@ global String comment CUDT =         0.00
@ global String comment AER_OPT =     0
@ global String comment SWINT_OPT =     0
@ global String comment AER_TYPE =     1
@ global String comment AER_AOD550_OPT =     1
@ global String comment AER_ANGEXP_OPT =     1
@ global String comment AER_SSA_OPT =     1
@ global String comment AER_ASY_OPT =     1
@ global String comment AER_AOD550_VAL =         0.12
@ global String comment AER_ANGEXP_VAL =         1.30
@ global String comment AER_SSA_VAL =         0.00
@ global String comment AER_ASY_VAL =         0.00
@ global String comment MOIST_ADV_OPT =     1
@ global String comment SCALAR_ADV_OPT =     1
@ global String comment TKE_ADV_OPT =     1
@ global String comment DIFF_6TH_OPT =     0
@ global String comment DIFF_6TH_FACTOR =         0.12
@ global String comment OBS_NUDGE_OPT =     0
@ global String comment BUCKET_MM =        -1.00
@ global String comment BUCKET_J =        -1.00
@ global String comment PREC_ACC_DT =         0.00
@ global String comment SF_OCEAN_PHYSICS =     0
@ global String comment ISFTCFLX =     0
@ global String comment ISHALLOW =     0
@ global String comment ISFFLX =     1
@ global String comment ICLOUD =     1
@ global String comment ICLOUD_CU =     0
@ global String comment TRACER_PBLMIX =     1
@ global String comment SCALAR_PBLMIX =     0
@ global String comment YSU_TOPDOWN_PBLMIX =     0
@ global String comment GRAV_SETTLING =     0
@ global String comment DFI_OPT =     0
@ global String comment SIMULATION_INITIALIZATION_TYPE = REAL-DATA CASE
@ global String comment WEST-EAST_PATCH_START_UNSTAG =     1
@ global String comment WEST-EAST_PATCH_END_UNSTAG =   261
@ global String comment WEST-EAST_PATCH_START_STAG =     1
@ global String comment WEST-EAST_PATCH_END_STAG =   262
@ global String comment SOUTH-NORTH_PATCH_START_UNSTAG =     1
@ global String comment SOUTH-NORTH_PATCH_END_UNSTAG =   279
@ global String comment SOUTH-NORTH_PATCH_START_STAG =     1
@ global String comment SOUTH-NORTH_PATCH_END_STAG =   280
@ global String comment BOTTOM-TOP_PATCH_START_UNSTAG =     1
@ global String comment BOTTOM-TOP_PATCH_END_UNSTAG =    29
@ global String comment BOTTOM-TOP_PATCH_START_STAG =     1
@ global String comment BOTTOM-TOP_PATCH_END_STAG =    30
@ global String comment GRID_ID =     2
@ global String comment PARENT_ID =     1
@ global String comment I_PARENT_START =    28
@ global String comment J_PARENT_START =    37
@ global String comment PARENT_GRID_RATIO =     3
@ global String comment CEN_LAT =        29.76
@ global String comment CEN_LON =       119.55
@ global String comment TRUELAT1 =        30.00
@ global String comment TRUELAT2 =        60.00
@ global String comment MOAD_CEN_LAT =        29.60
@ global String comment STAND_LON =       119.80
@ global String comment POLE_LAT =        90.00
@ global String comment POLE_LON =         0.00
@ global String comment GMT =         0.00
@ global String comment JULYR =  2009
@ global String comment JULDAY =   313
@ global String comment MAP_PROJ =     1
@ global String comment MAP_PROJ_CHAR = Lambert Conformal
@ global String comment MMINLU = USGS
@ global String comment NUM_LAND_CAT =    24
@ global String comment ISWATER =    16
@ global String comment ISLAKE =    -1
@ global String comment ISICE =    24
@ global String comment ISURBAN =     1
@ global String comment ISOILWATER =    14
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-10-14 15:54:12 | 显示全部楼层
river 发表于 2017-10-12 08:27
把你的描述文件也放上来

是这个吗
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-10-14 18:30:36 | 显示全部楼层

'set ylevs 0.1 0.3 0.4 0.6 0.9 1.2 1.7 2.2' 和'set z 1 10' 去掉看看。你这个运行没有错误提示吗?有的话一起发上来
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-10-15 06:30:16 | 显示全部楼层
本帖最后由 tangjie199365 于 2017-10-15 06:32 编辑
river 发表于 2017-10-14 18:30
'set ylevs 0.1 0.3 0.4 0.6 0.9 1.2 1.7 2.2' 和'set z 1 10' 去掉看看。你这个运行没有错误提示吗?有 ...

是z 1 10 那里的问题 。因为ctl文件里面只有一层的数据 画不出来 1到10 层的。谢谢你把纵坐标设置的值也去掉之后 系统自动给画出的图纵坐标赋值 不明白这个纵坐标是什么意思?是高度的意思吗?C:\Users\Fury\Desktop\付烨2个例\pblhd2-6
pblhd2-6.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-10-15 06:32:08 | 显示全部楼层
本帖最后由 tangjie199365 于 2017-10-15 06:33 编辑

求回复求回复
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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