grads提取wrf数据

泡沫学天气 · 发表于 2017-7-8 09:59:10

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

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

x

grads提取wrf数据，查找了很多帖子，也有人没有解决的，用了各种方法发现都出错，提取的文件很大，可是数据都是0
求大神指点
以下是gs

'reinit'
'open F:\2017\wrfout_d02_2017-03-24.ctl'
'set gxout fwrite'
'set fwrite f:\t.dat'
'set t 4'
z=1
  while(z<=19)
'set z 'z
'd T'
z=z+1
endwhile
'disable fwrite'
return
********************************************************
以下是wrf的ctl文件
dset F:\newwork\shiQ\wrfout_d02_2017-03-24.dat
undef 9.999E+20
title  OUTPUT FROM WRF V3.6 MODEL
xdef  260 levels
111.420792
111.451950
111.483109
111.514267
111.545425
111.576584
111.607742
111.638901
111.670059
111.701218
111.732376
111.763535
111.794693
111.825851
111.857010
111.888168
111.919327
111.950485
111.981644
112.012802
112.043961
112.075119
112.106277
112.137436
112.168594
112.199753
112.230911
112.262070
112.293228
112.324387
112.355545
112.386703
112.417862
112.449020
112.480179
112.511337
112.542496
112.573654
112.604813
112.635971
112.667130
112.698288
112.729446
112.760605
112.791763
112.822922
112.854080
112.885239
112.916397
112.947556
112.978714
113.009865
113.041023
113.072182
113.103340
113.134499
113.165657
113.196815
113.227974
113.259132
113.290291
113.321449
113.352608
113.383766
113.414925
113.446083
113.477242
113.508400
113.539558
113.570717
113.601875
113.633034
113.664192
113.695351
113.726509
113.757668
113.788826
113.819984
113.851143
113.882301
113.913460
113.944618
113.975777
114.006935
114.038094
114.069252
114.100410
114.131569
114.162727
114.193886
114.225044
114.256203
114.287361
114.318520
114.349678
114.380836
114.411995
114.443153
114.474312
114.505470
114.536629
114.567787
114.598946
114.630096
114.661255
114.692413
114.723572
114.754730
114.785889
114.817047
114.848206
114.879364
114.910522
114.941681
114.972839
115.003998
115.035156
115.066315
115.097473
115.128632
115.159790
115.190948
115.222107
115.253265
115.284424
115.315582
115.346741
115.377899
115.409058
115.440216
115.471375
115.502533
115.533691
115.564850
115.596008
115.627167
115.658325
115.689484
115.720642
115.751801
115.782959
115.814117
115.845276
115.876434
115.907593
115.938751
115.969910
116.001068
116.032227
116.063385
116.094543
116.125702
116.156860
116.188019
116.219177
116.250336
116.281494
116.312645
116.343803
116.374962
116.406120
116.437279
116.468437
116.499596
116.530754
116.561913
116.593071
116.624229
116.655388
116.686546
116.717705
116.748863
116.780022
116.811180
116.842339
116.873497
116.904655
116.935814
116.966972
116.998131
117.029289
117.060448
117.091606
117.122765
117.153923
117.185081
117.216240
117.247398
117.278557
117.309715
117.340874
117.372032
117.403191
117.434349
117.465508
117.496666
117.527824
117.558983
117.590141
117.621300
117.652458
117.683617
117.714775
117.745934
117.777092
117.808250
117.839409
117.870567
117.901726
117.932877
117.964035
117.995193
118.026352
118.057510
118.088669
118.119827
118.150986
118.182144
118.213303
118.244461
118.275620
118.306778
118.337936
118.369095
118.400253
118.431412
118.462570
118.493729
118.524887
118.556046
118.587204
118.618362
118.649521
118.680679
118.711838
118.742996
118.774155
118.805313
118.836472
118.867630
118.898788
118.929947
118.961105
118.992264
119.023422
119.054581
119.085739
119.116898
119.148056
119.179214
119.210373
119.241531
119.272690
119.303848
119.335007
119.366165
119.397324
119.428482
119.459641
119.490799
ydef  240 levels
26.1235123
26.1514893
26.1794510
26.2074127
26.2353592
26.2633133
26.2912445
26.3191833
26.3470993
26.3750229
26.4029312
26.4308395
26.4587326
26.4866333
26.5145111
26.5423889
26.5702591
26.5981293
26.6259766
26.6538315
26.6816711
26.7095108
26.7373352
26.7651672
26.7929916
26.8207932
26.8485870
26.8763885
26.9041748
26.9319611
26.9597397
26.9875183
27.0152740
27.0430145
27.0707779
27.0985184
27.1262512
27.1539764
27.1816940
27.2094116
27.2371140
27.2648239
27.2925110
27.3201981
27.3478775
27.3755493
27.4032135
27.4308777
27.4585266
27.4861755
27.5138092
27.5414429
27.5690613
27.5966873
27.6242905
27.6519012
27.6794815
27.7070770
27.7346573
27.7622375
27.7898026
27.8173676
27.8449249
27.8724594
27.9000015
27.9275436
27.9550629
27.9825745
28.0101013
28.0376053
28.0651016
28.0925827
28.1200714
28.1475525
28.1750259
28.2024841
28.2299423
28.2573853
28.2848282
28.3122711
28.3396835
28.3671188
28.3945236
28.4219360
28.4493256
28.4767303
28.5041046
28.5314789
28.5588531
28.5862122
28.6135788
28.6409225
28.6682739
28.6955948
28.7229309
28.7502518
28.7775726
28.8048782
28.8321686
28.8594666
28.8867416
28.9140244
28.9413071
28.9685593
28.9958115
29.0230637
29.0503159
29.0775375
29.1047745
29.1319885
29.1592026
29.1864090
29.2136078
29.2407913
29.2679749
29.2951660
29.3223343
29.3494873
29.3766479
29.4038010
29.4309311
29.4580765
29.4851913
29.5123215
29.5394363
29.5665359
29.5936279
29.6207199
29.6478043
29.6748810
29.7019424
29.7290192
29.7560654
29.7831116
29.8101425
29.8371811
29.8642120
29.8912201
29.9182281
29.9452286
29.9722290
29.9992218
30.0261993
30.0531616
30.0801392
30.1070938
30.1340485
30.1609879
30.1879196
30.2148590
30.2417755
30.2686844
30.2955933
30.3224945
30.3493881
30.3762741
30.4031525
30.4300232
30.4568787
30.4837418
30.5105820
30.5374222
30.5642624
30.5910797
30.6178970
30.6447144
30.6715164
30.6983185
30.7251053
30.7518768
30.7786484
30.8054276
30.8321762
30.8589325
30.8856659
30.9124146
30.9391327
30.9658508
30.9925766
31.0192795
31.0459747
31.0726624
31.0993576
31.1260376
31.1527023
31.1793671
31.2060165
31.2326584
31.2592926
31.2859268
31.3125534
31.3391647
31.3657761
31.3923874
31.4189682
31.4455490
31.4721375
31.4987030
31.5252762
31.5518265
31.5783691
31.6049118
31.6314392
31.6579742
31.6844864
31.7109985
31.7375031
31.7639999
31.7904892
31.8169785
31.8434372
31.8699112
31.8963547
31.9228210
31.9492493
31.9756851
32.0021133
32.0285339
32.0549393
32.0813599
32.1077499
32.1341476
32.1605225
32.1868896
32.2132492
32.2396164
32.2659683
32.2923050
32.3186417
32.3449783
32.3712845
32.3976135
32.4239044
32.4502029
32.4765015
32.5027847
32.5290527
32.5553207
32.5815735
32.6078339
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 00Z24MAR2017    60MN
VARS  107
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          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)
QNWFA2D       1  0  Surface aerosol number conc emission (kg-1 s-1)
QNCCN       19  0  CCN Number concentration (  kg(-1))
QNCLOUD    19  0  cloud water 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)
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    19  0  Radar reflectivity (lamda = 10 cm) (dBZ)
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)
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 (-)
AFWA_CLOUD    1  0  AFWA Diagnostic: Cloud ceiling (m)
AFWA_CAPE    1  0  AFWA Diagnostic: Convective Avail Pot Energy (J kg-1)
AFWA_ZLFC    1  0  AFWA Diagnostic: Level of Free Convection (m)
AFWA_PLFC    1  0  AFWA Diagnostic: Pressure of LFC (Pa)
MIDRH_MIN    1  0  Min Mid-level relative humidity (%)
AFWA_HAIL    1  0  AFWA Diagnostic: Hail Diameter (Weibull) (mm)
AFWA_LLWS    1  0  AFWA Diagnostic: 0-2000 ft wind shear (m s-1)
AFWA_TORNA    1  0  AFWA Diagnostic: Tornado wind speed (Weibull) (m s-1)
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)
ENDVARS
@ global String comment TITLE =  OUTPUT FROM WRF V3.6 MODEL
@ global String comment START_DATE = 2017-03-24_00:00:00
@ global String comment SIMULATION_START_DATE = 2017-03-23_18:00:00
@ global String comment WEST-EAST_GRID_DIMENSION = 261
@ global String comment SOUTH-NORTH_GRID_DIMENSION = 241
@ global String comment BOTTOM-TOP_GRID_DIMENSION = 40
@ global String comment DX =    3000.00
@ global String comment DY =    3000.00
@ global String comment STOCH_FORCE_OPT =    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 = 16
@ 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 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 =       18.00
@ global String comment RADT =       3.00
@ global String comment BLDT =       0.00
@ global String comment CUDT =       5.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 GRAV_SETTLING =    0
@ global String comment DFI_OPT =    0
@ global String comment WEST-EAST_PATCH_START_UNSTAG =    1
@ global String comment WEST-EAST_PATCH_END_UNSTAG = 260
@ global String comment WEST-EAST_PATCH_START_STAG =    1
@ global String comment WEST-EAST_PATCH_END_STAG = 261
@ global String comment SOUTH-NORTH_PATCH_START_UNSTAG =    1
@ global String comment SOUTH-NORTH_PATCH_END_UNSTAG = 240
@ global String comment SOUTH-NORTH_PATCH_START_STAG =    1
@ global String comment SOUTH-NORTH_PATCH_END_STAG = 241
@ global String comment BOTTOM-TOP_PATCH_START_UNSTAG =    1
@ global String comment BOTTOM-TOP_PATCH_END_UNSTAG = 39
@ global String comment BOTTOM-TOP_PATCH_START_STAG =    1
@ global String comment BOTTOM-TOP_PATCH_END_STAG = 40
@ global String comment GRID_ID =    2
@ global String comment PARENT_ID =    1
@ global String comment I_PARENT_START = 101
@ global String comment J_PARENT_START = 116
@ global String comment PARENT_GRID_RATIO =    5
@ global String comment CEN_LAT =       29.42
@ global String comment CEN_LON =    115.46
@ global String comment TRUELAT1 =       30.00
@ global String comment TRUELAT2 =       60.00
@ global String comment MOAD_CEN_LAT =       27.50
@ global String comment STAND_LON =    116.40
@ global String comment POLE_LAT =       90.00
@ global String comment POLE_LON =       0.00
@ global String comment GMT =       0.00
@ global String comment JULYR =  2017
@ global String comment JULDAY = 83
@ global String comment MAP_PROJ =    3
@ global String comment MAP_PROJ_CHAR = Mercator
@ 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
************************************************************************

astiny · 发表于 2017-7-8 14:11:57

目测ctl里面加big edian

river · 发表于 2017-7-8 16:40:31

你自己变得描述文件呢？缺测值设置的多少？

xibao · 发表于 2017-9-6 18:30:14

有没有人会画剖面图的啊，能不能教教我

xibao · 发表于 2017-9-6 18:30:29

最近毕业论文要画剖面图啊，希望大神帮帮我

ymm9106 · 发表于 2017-10-9 13:03:41

同求解决方法大神

		自动登录	找回密码
密码			立即注册

[分享资料] grads提取wrf数据

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