有人用过NASA MERRA的数据么？不知道有没有人画出过它的剖面图。。。。。

晚晴山

因为气象家园最近添加不了图片和附件了
所以我把这种资料的nc文件的ctl给手动打出来了。。。。。。。。。我不知道这个问题有没有人遇到过，也可能没有答案。只是请有过类似经历或者懂的大神指点一下，之前我一直以为是因为多nc合并之后出的错，没想到单个文件也画不出垂直剖面图，二维的时间高度剖面图也作不出来，其他图正常。@兰溪之水

dset G:/MERRA/MERRA2014/MERRA300.prod.assim.inst6_3d_ana_Np.20140620.SUB.nc
title
undef 1e+15
xdef 540 levels -180 -179.333 -178.667 -178 -177.333 -176.667 -176 -175.333
-174.667 -174 -173.333 -172.667 -172 -171.333 -170.667 -170 -169.333 -168.667
-168 -167.333 -166.667 -166 -165.333 -164.667 -164 -163.333 -162.667 -162
-161.333 -160.667 -160 -159.333 -158.667 -158 -157.333 -156.667 -156 -155.333
-154.667 -154 -153.333 -152.667 -151.333 -150.667 -150 -149.333 -148.667
-148 -147.333 -146.667 -146 -145.333 -144.667 -144 -143.333 -142.667 -142
-141.333 -140.667 -140 -139.333 -138.667 -138 -137.333 -136.667 -136 -135.333
-134.667 -134 -133.333 -126.667 -126 -125.333 -124.667 -124 -123.333 -122.667 -122
-121.333 -120.667 -120 -119.333 -118.667 -118 -117.333 -116.667 -116 -115.333
-114.667 -114 -113.333 -112.667 -112 -111.333 -110.667 -110 -109.333 -108.667
-108 -107.333 -106.667 -106 -105.333 -104.667 -104 -103.333 -102.667 -102
-101.333 -100.667 -100 -99.3333 -98.6667 -98 -97.3333 -96.6667 -96 -95.3333
-94.6667 -94 -93.3333 -92.6667 -92 -91.3333 -90.6667 -90 -89.3333 -88.6667
-88 -87.3333 -86.6667 -86 -85.3333 -84.6667 -84 -83.3333 -82.6667 -82
-81.3333 -80.6667 -80 -79.3333 -78.6667 -78 -77.3333 -76.6667 -76 -75.3333
-74.6667 -74 -73.3333 -72.6667 -72 -71.3333 -70.6667 -70 -69.3333 -68.6667
-68 -67.3333 -66.6667 -66 -65.3333 -64.6667 -64 -63.3333 -62.6667 -62
-61.3333 -60.6667 -60 -59.3333 -58.6667 -58 -57.3333 -56.6667 -56 -55.3333
-54.6667 -54 -53.3333 -52.6667 -52 -51.3333 -50.6667 -50 -49.3333 -48.6667
-48 -47.3333 -46.6667 -46 -45.3333 -44.6667 -44 -43.3333 -42.6667 -42
-41.3333 -40.6667 -40 -39.3333 -38.6667 -38 -37.3333 -36.6667 -36 -35.3333
-34.6667 -34 -33.3333 -32.6667 -32 -31.3333 -30.6667 -30 -29.3333 -28.6667
-28 -27.3333 -26.6667 -26 -25.3333 -24.6667 -24 -23.3333 -22.6667 -22
-21.3333 -20.6667 -20 -19.3333 -18.6667 -18 -17.3333 -16.6667 -16 -15.3333
-14.6667 -14 -13.3333 -12.6667 -12 -11.3333 -10.6667 -10 -9.33333 -8.66667
-8 -7.33333 -6.66667 -6 -5.33333 -4.66667 -4 -3.33333 -2.66667 -2
-1.33333 -0.666672 0 0.666672 1.33333 2 2.66667 3.33333 4 4.66667
5.33333 6 6.66667 7.33333 8 8.66667 9.33333 10 10.6667 11.3333
12 12.6667 13.3333 14 14.6667 15.3333 16 16.6667 17.3333 18
18.6667 19.3333 20 20.6667 21.3333 28 28.6667 29.3333 30 30.6667 31.3333
32 32.6667 33.3333 40 40.6667 41.3333 42 42.6667 43.3333 44 44.6667
45.3333 46 46.6667 47.3333 48 48.6667 49.3333 50 50.6667 51.3333
52 52.6667 53.3333 54 54.6667 55.3333 56 56.6667 57.3333 58
58.6667 59.3333 60 60.6667 61.3333 62 62.6667 63.3333 64 64.6667
65.3333 66 66.6667 67.3333 68 68.6667 69.3333 70 70.6667 71.3333
72 72.6667 73.3333 74 74.6667 75.3333 76 76.6667 77.3333 78
78.6667 79.3333 80 80.6667 81.3333 82 82.6667 83.3333 84 84.6667
85.3333 86 86.6667 87.3333 88 88.6667 89.3333 90 90.6667 91.3333
92 92.6667 93.3333 100 100.667 101.333 102 102.667 103.333 104 104.667
105.333 106 106.667 107.333 108 108.667 109.333 110 110.667 111.333
112 112.667 113.333 114 114.667 115.333 116 116.667 117.333 118
118.667 119.333 120 120.667 121.333 122 122.667 123.333 124 124.667
125.333 126 126.667 127.333 128 128.667 129.333 130 130.667 131.333
132 132.667 133.333 134 134.667 135.333 136 136.667 137.333 138
138.667 139.333 140 140.667 141.333 142 142.667 143.333 144 144.667
145.333 146 146.667 147.333 148 148.667 149.333 150 150.667 151.333
152 152.667 153.333 154 154.667 155.333 156 156.667 157.333 158
158.667 159.333 160 160.667 161.333 162 162.667 163.333 164 164.667
165.333 166 166.667 167.333 168 168.667 169.333 170 170.667 171.333
172 172.667 173.333 174 174.667 175.333 176 176.667 177.333 178
178.667 179.333
ydef 361 linear -90 0.5
zdef 18 levels 1000 925 900 850 800 750 700 650 600 550 500 400 350 300 250 200 150 100
tdef 4 linear 00Z20JUN2014 360mn
vars 4
t=>t 18  t,z,y,x  Air temperature
u=>u 18  t,z,y,x  Eastward wind component
v=>v 18  t,z,y,x  Northward wind component
qv=>qv 18  t,z,y,x  Specific humidity
endvars

lqouc

贴一下报错和你用的命令吧，光靠一个ctl啥也猜不出的。

晚晴山

lqouc 发表于 2015-8-6 14:56
贴一下报错和你用的命令吧，光靠一个ctl啥也猜不出的。

你怎么换头像了。。。。。。。。

晚晴山

lqouc 发表于 2015-8-6 14:56
贴一下报错和你用的命令吧，光靠一个ctl啥也猜不出的。

不报错 就是用MI将几个时段的nc合并了，然后用GrADS画时间高度剖面，出来的图类似密集恐惧症的。。。。总之一看就是错的剖面。气象家园不知怎么上传不了图片了。。所以不好描述。最新进展表明有可能是合并后的文件过于巨大（4G）GrADS可能搞不定了。。。。。。。。。

lqouc

晚晴山 发表于 2015-8-6 16:03
不报错 就是用MI将几个时段的nc合并了，然后用GrADS画时间高度剖面，出来的图类似密集恐惧症的。。。。总 ...

。。。。。好吧，是有可能因为文件过大。
其实附件的事，只有发帖时不好用，回帖时是可以的。

兰溪之水

你这啥也没有让人家怎么猜。。

晚晴山

lqouc 发表于 2015-8-6 16:06
。。。。。好吧，是有可能因为文件过大。
其实附件的事，只有发帖时不好用，回帖时是可以的。

我已经用CDO解决了 谢谢大神~

晚晴山

兰溪之水 发表于 2015-8-6 17:18
你这啥也没有让人家怎么猜。。

哈哈 我已经用CDO解决了 结果是MI的合并有问题。

jr921102

手打了这个ctl。。。真是辛苦了