【已解决】关于scorr函数中的一个问题

开始起飞1

错误信息提示

river 发表于 2013-4-18 09:33

                               
登录/注册后可看大图

你都已经用set dfile命令了干吗还在变量后面加那个区分码，运行的时候肯定提示你a.3 a.4不是有效的变量名 ...

批评的是，我把set t 1 给加上了，这是提示的错误信息，麻烦了啊 file:///C:/Documents%20and%20Settings/Administrator/Application%20Data/Tencent/Users/617039055/QQ/WinTemp/RichOle/ZPR%7BHWDPLB5@X%7DBDR1M%7BV~K.jpg

开始起飞1

river 发表于 2013-4-18 09:33

                               
登录/注册后可看大图

你都已经用set dfile命令了干吗还在变量后面加那个区分码，运行的时候肯定提示你a.3 a.4不是有效的变量名 ...

那个错误信息提示，我粘贴了一个，又上传了一个，希望能看到，若下载要金钱的话我报销哈，呵呵

river

开始起飞1 发表于 2013-4-18 13:22

                               
登录/注册后可看大图

那个错误信息提示，我粘贴了一个，又上传了一个，希望能看到，若下载要金钱的话我报销哈，呵呵

错误提示说的很明白了，不匹配的网格，你站点资料的网格的分辨率和另外那个资料的网格是不同的吧？

开始起飞1

river 发表于 2013-4-18 18:05

                               
登录/注册后可看大图

错误提示说的很明白了，不匹配的网格，你站点资料的网格的分辨率和另外那个资料的网格是不同的吧？

我现在蒙圈了，这些资料的网格分辨率都是1.25*1的，方便加下QQ聊么？我的是617039055，有劳这位仁兄了，感激不尽啊

river

开始起飞1 发表于 2013-4-18 21:11

                               
登录/注册后可看大图

我现在蒙圈了，这些资料的网格分辨率都是1.25*1的，方便加下QQ聊么？我的是617039055，有劳这位仁兄了， ...

画站点资料的等值线图的时候不是先要生成一个背景网格吗，那个分辨率应该也弄成1.25*1的才行吧

开始起飞1

river 发表于 2013-4-18 22:03

                               
登录/注册后可看大图

画站点资料的等值线图的时候不是先要生成一个背景网格吗，那个分辨率应该也弄成1.25*1的才行吧

dset c:\GrADS19\data\grid.grd
undef -999.0
title sample grid data
xdef 71 linear 70 1.25
ydef 41 linear 15 1.0
zdef 1 linear 500 1
tdef 1 linear jun1980 1mo
vars 1
g 0 99 grid data
endvars
那个格点数据的范围是全球的，不知道这个范围的不一致对结果影不影响呢？

river

开始起飞1 发表于 2013-4-18 22:05

                               
登录/注册后可看大图

dset c:\GrADS19\data\grid.grd
undef -999.0
title sample grid data

全球的那个资料分辨率也是1.25*1的吗？可能是那个资料的分辨率问题吧，楼主把这个grid.grd的分辨率生成和那个资料一样的再试试

开始起飞1

river 发表于 2013-4-18 22:12

                               
登录/注册后可看大图

全球的那个资料分辨率也是1.25*1的吗？可能是那个资料的分辨率问题吧，楼主把这个grid.grd的分辨率生成和 ...

Longitudes:  288 bins centered on 179.375 W to 179.375 E  (1.25 degree steps)  
Latitudes :  180 bins centered on  89.5   S to  89.5   N  (1.00 degree steps)
我根据这个改了下grid.grd和ctl描述文件，ctl如下：
dset c:\GrADS19\data\grid.grd
undef -999.0
title sample grid data
xdef 288 linear -179.375 1.25
ydef 180 linear -89.5 1.0
zdef 1 linear 500 1
tdef 1 linear jan1980 1mo
vars 1
g 0 99 grid data
endvars
gs文件为：'reinit'
'open c:\GrADS19\data\1980grid.ctl'
'open c:\GrADS19\data\22.ctl'
'set grads off'
'define a1=oacres(g,x.2)'
'define aaa=smth9(a1)'
'open 1980data.ctl'
'set dfile 3'
'define aa=ave(a,t=1,t=366)'
'open 1981data.ctl'
'set dfile 4'
'define a2=ave(a,t=1,t=365)'
'set mpdset cnscs'
'set lon 73 135.5'
'set lat 17 55'
'set t 1'
'd scorr(a1,(aa+a2)/2,lon=73,lon=120.5,lat=17,lat=55)'
'enable print xiangguan.gmf'
'print'
'disable print'
;
倒没发现什么错误信息，但运行了还是什么都没有，只是有个类似于图的东西闪了一下file:///C:/Documents%20and%20Settings/Administrator/Application%20Data/Tencent/Users/617039055/QQ/WinTemp/RichOle/1L2527HT(_WQAUWH574)C08.jpg

river

开始起飞1 发表于 2013-4-18 22:37

                               
登录/注册后可看大图

Longitudes:  288 bins centered on 179.375 W to 179.375 E  (1.25 degree steps)  
Latitudes :  180  ...

我不是说要把背景格点数整的和另外那个资料一模一样，我是说分辨率一样应该就行，而且背景格点是需要重新生成的，改ctl文件好像不太行啊。你再试试吧

开始起飞1

river 发表于 2013-4-19 07:40

                               
登录/注册后可看大图

我不是说要把背景格点数整的和另外那个资料一模一样，我是说分辨率一样应该就行，而且背景格点是需要重新 ...

我知道grid.grd也是要重新生成的，这是新修改的生成网格的fortran程序，麻烦仁兄看下有啥问题么
parameter(nx=288,ny=180)
real lat(ny),lon(nx)
real s(nx,ny)
open(1,file='c:\GrADS19\data\grid.grd',form='binary')
lat(1)=-89.5
lon(1)=-179.375
do j=1,ny-1
lat(j+1)=lat(j)+1.0
enddo
do i=1,nx-1
lon(i+1)=lon(i)+1.25
enddo
do i=1,nx
do j=1,ny
s(i,j)=1
enddo
enddo
write(1)s
end