如何将NCEP R1的sigma层数据转换成等压面数据？

kftqui

下载了NCEP的非绝热加热率资料，它是由T62全球谱模式输出的sigma层数据，垂直方向分为18个不等距的sigma层，ctl文件如下：

1. dset d:/2013/201311.grb3d1630
   
2. index d:/2013/201311.grb3d1630.idx
   
3. undef 9.999E+20
   
4. title d:/2013/201311.grb3d1630
   
5. *  produced by grib2ctl v0.9.12.5p39c
   
6. dtype grib 255
   
7. options yrev
   
8. xdef 192 linear 0.000000 1.875
   
9. ydef 94 levels
   
10. -88.542 -86.653 -84.753 -82.851 -80.947 -79.043 -77.139 -75.235 -73.331 -71.426
    
11. -69.522 -67.617 -65.713 -63.808 -61.903 -59.999 -58.094 -56.189 -54.285 -52.380
    
12. -50.475 -48.571 -46.666 -44.761 -42.856 -40.952 -39.047 -37.142 -35.238 -33.333
    
13. -31.428 -29.523 -27.619 -25.714 -23.809 -21.904 -20.000 -18.095 -16.190 -14.286
    
14. -12.381 -10.476  -8.571  -6.667  -4.762  -2.857  -0.952   0.952   2.857   4.762
    
15.    6.667   8.571  10.476  12.381  14.286  16.190  18.095  20.000  21.904  23.809
    
16.   25.714  27.619  29.523  31.428  33.333  35.238  37.142  39.047  40.952  42.856
    
17.   44.761  46.666  48.571  50.475  52.380  54.285  56.189  58.094  59.999  61.903
    
18.   63.808  65.713  67.617  69.522  71.426  73.331  75.235  77.139  79.043  80.947
    
19.   82.851  84.753  86.653  88.542
    
20. tdef 60 linear 06Z16nov2013 6hr
    
21. *  z has 28 levels, for sig
    
22. zdef 28 levels
    
23. 9950 9821 9644 9425 9159 8838 8458 8014 7508 6943 6329 5681 5017 4357 3720 3125 2582 2101 1682 1326 1028 782 580 418 288 183 101 27
    
24. vars 11
    
25. CNVHRsig 28 242,107,0 ** (profile) Deep convective heating [K/s]
    
26. CNVMRsig 28 243,107,0 ** (profile) Deep convective moistening [kg/kg/s]
    
27. LRGHRsig 28 241,107,0 ** (profile) Large scale condensation heating [K/s]
    
28. LWHRsig 28 251,107,0 ** (profile) Longwave radiative heating [K/s]
    
29. SHAHRsig 28 244,107,0 ** (profile) Shallow convective heating [K/s]
    
30. SHAMRsig 28 245,107,0 ** (profile) Shallow convective moistening [kg/kg/s]
    
31. SWHRsig 28 250,107,0 ** (profile) Solar radiative heating [K/s]
    
32. VDFHRsig 28 246,107,0 ** (profile) Vertical diffusion heating [K/s]
    
33. VDFMRsig 28 249,107,0 ** (profile) Vertical diffusion moistening [kg/kg/s]
    
34. VDFUAsig 28 247,107,0 ** (profile) Vertical diffusion zonal accel [m/s^2]
    
35. VDFVAsig 28 248,107,0 ** (profile) Vertical diffusion meridional accel [m/s^2]
    
36. ENDVARS
    

复制代码

我需要将该资料转换成等压面资料，应该怎么做？

夏朗的芒果

如果楼主你会ncl的话，这个脚本可以帮到你，从sigma坐标插值到p坐标

kftqui

夏朗的芒果 发表于 2014-11-20 18:19
如果楼主你会ncl的话，这个脚本可以帮到你，从sigma坐标插值到p坐标

谢谢您！听说NCL功能很强大，正打算学。