请教GRADS里的PDEF插值函数

alexclimate

请教大家个问题，有1*1度的网格数据，想插值成2.5*2.5的，听说GRADS里通过PDEF可以设置，但是没弄过，怎么设置呀，谢谢！

CTL如下

dset D:\h500mean01to0608.grd
undef 9.999e+20
title hgt-JJA
xdef 360 linear 0 1
ydef 181 linear -90 1
zdef 1 linear 0 1
tdef 1 linear 00Z01JAN1982 1mn
vars 1
x1  0  t,y,x  Geopotential Height
endvars

维克多·蓝晶

RE的使用

一、使用方法
1．re(expr,nx,'linear',lon,dlon,ny,'linear',lat,dlat, ['ig',nyig],['ba'|'bl'|'bs'|'vt',vtmax,vtmin|'ma',min]

2．re(expr,nx,'linear',lon,dlon,ny,'gaus',gstart,njog, ['ig',nyig],['ba'|'bl'|'bs'|'vt',vtmax,vtmin|'ma',min]

Erpr：变量
Nx：x方向上或经向格点数
‘linear’和’gaus’表示目标变量存放的属性为“线性”和“高斯”格点
Lon：起始经度
Dlon：经度间隔
Ny：y方向上或纬向格点数
Lat：起始纬度
Dlat：纬度间隔
Gstart：第一个高斯格点数，即高斯格点起始纬(经)度
Njog：输出格点的全球高斯格点纬度数，一般与ny相同
['ig',nyig]：输入的格点为高斯格点，nyig表示高斯纬度数，必须≥8且为4的倍数，这个参数是为了在高斯网格间进行更为精确的边界计算。
['ba'|'bl'|'bs'|'vt',vtmax,vtmin|'ma',min]：为选择插值方案，ba区域平均，bl双线性插值，bs三阶Bessel插值
二、举例说明
1．re(expr,nx,'linear',lon,dlon,ny,'linear',lat,dlat, ['ig',nyig],['ba'|'bl'|'bs'|'vt',vtmax,vtmin|'ma',min]
1.1 原始数据ctl文件
dset g:\yang\sst.dat
undef -999.0
title sst
xdef 360 linear -180.0  1.0
ydef 180 linear  -90.  1.0
zdef 1 levels 1000
tdef 12 linear 00z01jan1990  1mon
vars 1
sst 0 99  globle ocean sst
endvars
1.2 转化为2°x 2°格点gs文件
'reinit'
'open G:\ sst.ctl'
'set gxout fwrite'
'set fwrite G:\ sst2.grd'

'set t 1'
'set lev 1000'
'set x 1 360'
'set y 1 180'
'd re(sst,180,linear,-180,2,91,linear,-90,2,bl)'
'disable fwrite'
;
1.3输出文件的ctl文件
dset g:\sst2.grd
undef -999.0
title  sst
xdef 180 linear -180  2.0
ydef 91 linear  -90.  2.0
zdef 1 levels 1000
tdef 1 linear 00z01jan1960  1mon
vars 1
sst 0 -999  globle ocean sst
endvars
2．re(expr,nx,'linear',lon,dlon,ny,'gaus',gstart,njog, ['ig',nyig],['ba'|'bl'|'bs'|'vt',vtmax,vtmin|'ma',min]
2.1 读取gs文件，x方向为线性格点，y方向为高斯格点
'open sst.ctl'
'set gxout fwrite'
'set fwrite sst3.dat'
'set t 1'
'set z 1'
'd re(sst,256,'linear',0,1.40625,128,'gaus',1,128)'
'disable fwrite'
2.3输出文件的ctl文件
dset  e:\sst3.dat
undef  -9.99e8
title  sst
xdef 256  linear 0 1.40625
ydef 128  levels -88.9277        -87.5387        -86.1415        -84.7424        -83.3426        -81.9425
        -80.5421        -79.1417        -77.7412        -76.3406        -74.94        -73.5394        -72.1387
        -70.7381        -69.3374        -67.9367        -66.536        -65.1353        -63.7345        -62.3338
        -60.9331        -59.5323        -58.1316        -56.7309        -55.3301        -53.9294        -52.5286
        -51.1279        -49.7271        -48.3264        -46.9256        -45.5248        -44.1241        -42.7233
        -41.3226        -39.9218        -38.5211        -37.1203        -35.7195        -34.3188        -32.918
        -31.5172        -30.1165        -28.7157        -27.315        -25.9142        -24.5134        -23.1127
        -21.7119        -20.3111        -18.9104        -17.5096        -16.1088        -14.7081        -13.3073
        -11.9065        -10.5058        -9.105        -7.7042        -6.3035        -4.9027        -3.5019
        -2.1012        -0.7004        0.7004        2.1012        3.5019        4.9027        6.3035        7.7042        9.105        10.5058
        11.9065        13.3073        14.7081        16.1088        17.5096        18.9104
        20.3111        21.7119        23.1127        24.5134        25.9142        27.315        28.7157
        30.1165        31.5172        32.918        34.3188        35.7195        37.1203        38.5211
        39.9218        41.3226        42.7233        44.1241        45.5248        46.9256        48.3264
        49.7271        51.1279        52.5286        53.9294        55.3301        56.7309        58.1316
        59.5323        60.9331        62.3338        63.7345        65.1353        66.536        67.9367
        69.3374        70.7381        72.1387        73.5394        74.94        76.3406        77.7412
        79.1417        80.5421        81.9425        83.3426        84.7424
        86.1415        87.5387        88.9277
zdef  1  levels   1000
tdef  1  linear  0z01Jan1961  1mo
vars   1
sst   0   99   sst
endvars
注：
1.在输出时可不对经纬度进行设置，若需设置时则尽量设置为x,y，设置为lat,lon容易发生漂移
2.减少使用define，尽量在d之后直接re()。使用define也容易产生漂移现象
3.在linear和gaus尽量加单引号‘’

alexclimate

非常感谢！

sa13226139

谢楼主无私分享！

liujxt

学习，学习

梦纯

维克多·蓝晶 发表于 2015-1-8 15:49
RE的使用

一、使用方法

楼主，这个re命令是哪个版本grads带的呀？我用的grads2.0.1没有哎。多谢了哈

梦纯

alexclimate 发表于 2015-1-8 16:09
非常感谢！

请问re命令是哪个版本的grads中的呀？我用的grads2.0.1，没有这个命令呀，多谢了哈

alexclimate

2.0.a9.oga.1

小胡00

又学习到了

野田妹酱

真的太感谢了