交叉变量画图求助

infocert

各位前辈，小弟在处理mksrf_soitex.10level.c010119.nc文件过程中遇到问题，请各位伸出援助之手，非常感谢！！！

1. netcdf mksrf_soitex.10level.c010119 {
   
2. dimensions:
   
3.     lon = 4320 ;
   
4.     lat = 2160 ;
   
5.     number_of_layers = 10 ;
   
6.     number_of_mapunits = 4931 ;
   
7.     max_value_mapunit = 6998 ;
   
8. variables:
   
9.     float LON(lon) ;
   
10.         LON:long_name = "lon" ;
    
11.         LON:units = "degrees east" ;
    
12.     float LAT(lat) ;
    
13.         LAT:long_name = "lat" ;
    
14.         LAT:units = "degrees north" ;
    
15.     float LONGXY(lat, lon) ;
    
16.         LONGXY:long_name = "longitude-2d" ;
    
17.         LONGXY:units = "degrees east" ;
    
18.     float LATIXY(lat, lon) ;
    
19.         LATIXY:long_name = "latitude-2d" ;
    
20.         LATIXY:units = "degrees north" ;
    
21.     float LANDMASK(lat, lon) ;
    
22.         LANDMASK:long_name = "land mask" ;
    
23.         LANDMASK:units = "unitless" ;
    
24.     float DZSOI(number_of_layers) ;
    
25.         DZSOI:long_name = "soil layer thickness" ;
    
26.         DZSOI:units = "m" ;
    
27.     float ZSOI(number_of_layers) ;
    
28.         ZSOI:long_name = "soil layer depth" ;
    
29.         ZSOI:units = "m" ;
    
30.     float EDGEN ;
    
31.         EDGEN:long_name = "northern edge of surface grid" ;
    
32.         EDGEN:units = "degrees north" ;
    
33.     float EDGEE ;
    
34.         EDGEE:long_name = "eastern edge of surface grid" ;
    
35.         EDGEE:units = "degrees east" ;
    
36.     float EDGES ;
    
37.         EDGES:long_name = "southern edge of surface grid" ;
    
38.         EDGES:units = "degrees north" ;
    
39.     float EDGEW ;
    
40.         EDGEW:long_name = "western edge of surface grid" ;
    
41.         EDGEW:units = "degrees east" ;
    
42.     float MAPUNITS(lat, lon) ;
    
43.         MAPUNITS:long_name = "igbp soil mapunit" ;
    
44.         MAPUNITS:units = "unitless" ;
    
45.     float PCT_SAND(number_of_layers, max_value_mapunit) ;
    
46.         PCT_SAND:long_name = "percent sand" ;
    
47.         PCT_SAND:units = "unitless" ;
    
48.     float PCT_CLAY(number_of_layers, max_value_mapunit) ;
    
49.         PCT_CLAY:long_name = "percent clay" ;
    
50.         PCT_CLAY:units = "unitless" ;
    
51. 
    
52. // global attributes:
    
53.         :data_type = "igbp_soil_texture_data" ;
    
54. }
    

复制代码

小弟写的ctl如下
　

1. DEET sand_soitex.10level.c010119.nc
   
2. XDEF lon 4320 LINEAR -179.9583 0.0833
   
3. YDEF lat 2160 LINEAR -89.95834 0.08333
   
4. ZDEF number_of_layers 10 LINEAR 0 1
   

复制代码

我想显示PCT_SAND的图像，应该怎么搞？小弟新手，多谢各位前辈

infocert

怎么没前辈回复呢，写得不清楚吗？
自己顶，不要沉了

river

兰溪有一个帖子讲得比较清楚，你自己搜索一下看看吧

infocert

多谢了，那个FNL_1x1资料下载不下来，不知道为啥，可能是我的网络的问题。

infocert

不好意思，我还想问问，这个PCT_SAND变量有没有办法直接画出来？
还是必须要写程序自己转成GrADS的格式，比如改写成Lon,Lat,PCT_SAND * 10这种数据格式？