ncl读grd数据

liulinaacc

grd的ctl如下：
dset    ^wrfchem2index_20140715.grd
undef    -999.9
title    9km OUTPUT FROM WRF/CHEM V3.4 MODEL
xdef     284  linear   99.0   0.081
ydef     161  linear   24.0   0.081
zdef     1    linear   1      1
tdef     3   linear 00z16JUL2014  24hr
VARS     8
pm25     1  0   IAQI of pm2.5
pm10     1  0   IAQI of pm10
so2      1  0   IAQI of SO2
no2      1  0   IAQI of NO2
co       1  0   IAQI of CO
o3_1h    1  0   IAQI of O3-1hr
o3_8h    1  0   IAQI of O3-8hr
aqi      1  0   Air Quality Index (AQI)
ENDVARS
用ncl读取aqi第八个变量的时候应噶如何读取？
dims = (/8,3,161,284/)
a=fbindirread(f0+"wrfchem2index_20140715.grd",0,dims,"float")
我这样读对不对？画出来的图打不开是错的

心如止水的饭团

楼主，你的问题解决了吗？后面的变量都要怎么读的

liulinaacc

心如止水的饭团 发表于 2015-7-10 15:25
楼主，你的问题解决了吗？后面的变量都要怎么读的

已经解决了，
dims = (/3,8,161,284/)
a=fbindirread(f0+"wrfchem2index_20140715.grd",0,dims,"float") 这样读就行了

liulinaacc

心如止水的饭团 发表于 2015-7-10 15:25
楼主，你的问题解决了吗？后面的变量都要怎么读的

已经解决了，
dims = (/3,8,161,284/)
a=fbindirread(f0+"wrfchem2index_20140715.grd",0,dims,"float") 这样读就行了

心如止水的饭团

我现在有个问题，我的grd里面有20个层次，8个变量，但是有一两个变量只有一层，有些是有20层的
XDEF  180 LINEAR 1. 1.0
YDEF  170 LINEAR 1. 1.0
ZDEF   12 LINEAR 1. 1.0
TDEF  73 LINEAR 16:00Z21JUN2015 60MN
VARS   25
U              12    99   U             m/s
V              12    99   V             m/s
T              12    99   T             degree celsius
P              12    99   PRESS         hpa
HGT            12    99   HEIGHT        km
RAIN            1    99   RAIN          m
PBL             1    99   PBL           m
RH             12    99   RH            %
SO2            12    99   SO2           ugm-3
NO2            12    99   NO2           ugm-3
CO             12    99   CO            mgm-3
O3             12    99   O3            ugm-3
ASO4           12    99   ASO4          ugm-3
ANO3           12    99   ANO3          ugm-3
ANH4           12    99   ANH4          ugm-3
BC             12    99   BC            ugm-3
OC             12    99   OC            ugm-3
PPMF           12    99   PPMF          ugm-3
PPMC           12    99   PPMC          ugm-3
SOA            12    99   SOA           ugm-3
TPM25          12    99   TPM25         ugm-3
TPM10          12    99   TPM10         ugm-3
O3_8H          12    99   O3_8H         ugm-3
VISIB          12    99   visibility in km
AOD             1    99   aerosol Optical DEPTH
ENDVARS
这样的ctl要怎么读？

liulinaacc

心如止水的饭团 发表于 2015-7-22 15:44
我现在有个问题，我的grd里面有20个层次，8个变量，但是有一两个变量只有一层，有些是有20层的
XDEF  180  ...

这个我也没读过呢？应该各个变量是相同层次的就好读了

xiaocaoqiqiao

心如止水的饭团 发表于 2015-7-22 15:44
我现在有个问题，我的grd里面有20个层次，8个变量，但是有一两个变量只有一层，有些是有20层的
XDEF  180  ...

请问你解决了吗？我也有这样的困惑

心如止水的饭团

xiaocaoqiqiao 发表于 2016-7-12 16:23
请问你解决了吗？我也有这样的困惑

后来我直接用cdo把它转成nc格式了

zly4814624

请问你这个读取，时间是怎么做循环的。

藕茭白空心菜

心如止水的饭团 发表于 2016-7-13 14:28
后来我直接用cdo把它转成nc格式了

CDO可以把一堆grid数据（a年*12月）变成一个nc吗