站点数据转格点数据，画降水距平百分率

蜻蜓队长

小女子要做一个暖年发展阶段降水距平的合成分析。先已用Fortran提取了我所要的七年的降水距平百分率。写出其相应的的.ctl文件。
执行.gs文件报错太多。还请大神帮忙
rainhechengtiqu1.ctl文件如下：
DSET   d:\rain\rainhechengtiqu1.grd
DTYPE  station   
STNMAP d:\rain\cor.map
UNDEF  -999.0  
TITLE  Station Data Sample  
TDEF   7 linear 01jan1951 1yr
VARS 1  
p    0  99  rainfall anomaly percentage of 160station
endvars


站点转格点的Fortran程序：
Program  main
           parameter(nx=71, ny=41)
           real     lat(ny), lon(nx)
           real     s(nx,ny)
           open(1, file= 'd:/rain/grid.grd', form='binary')
           lat(1)=15.0
           lon(1)=70.0
           do j=1,ny-1
                lat(j+1)=lat(j)+1.0
           enddo
           do i=1,nx-1
                lon(i+1)=lon(i)+1.0
           enddo

            do i=1,nx
                 do j=1,ny
                       s(i,j)= 1
                 enddo
            enddo
            write(1)    s
            end


grid.ctl文件
dset          d:\rain\grid.grd
undef        -999.0
title            Grid data
xdef          71          linear     70  1
ydef          41          linear     15   1
zdef           1           linear     1000    1
tdef           7         linear      01jan1951     1yr
vars 1
g       0        99         grid data
endvars



gs脚本
'open d:\rain\grid.ctl'
'open d:\rain\rainhechengtiqu1.ctl'
'enable print d:\rain\eprainanomper1.gmf'
'set grads off'
'set lon 73 135'
'set lat 15 55'
'set mpdset cn  cnriver'
'define a=oacres(p.2,g,1.5)'
'define a1=maskout(a,g-0.5)'
'define aa=smth9(a1)'
'set xlops 1 10 0.18'
'set ylops 1 10 0.18'
'set gxout shaded'
'set ccolour rainbow'
'aave=ave(aa,t=1,t=7,1)'
'd aave'
'cbarn'
'set gxout contour'
'd aave'
'pull dummy'
'd:\rain\southsea.gs'
'print'
'disable print'
'reinit'
  ;

还望好心人帮我找出错误所在。

river

请在论坛grads板块搜索 站点资料作图 然后看看帖子，一步一步来

蜻蜓队长

river 发表于 2016-4-20 17:16
请在论坛grads板块搜索 站点资料作图 然后看看帖子，一步一步来

谢谢关心，我已经解决啦！

蜻蜓队长

river 发表于 2016-4-20 17:16
请在论坛grads板块搜索 站点资料作图 然后看看帖子，一步一步来

谢谢关心，我已经解决啦！

蜻蜓队长

river 发表于 2016-4-20 17:16
请在论坛grads板块搜索 站点资料作图 然后看看帖子，一步一步来

谢谢关心，我已经解决啦！

river

蜻蜓队长 发表于 2016-4-21 14:33
谢谢关心，我已经解决啦！

解决就好，下次就可以给有同样问题的同学解答问题了

dyhtt29

{:lxm_24:}{:lxm_24:}

小叶子

站点转格点程序是单独的吗

蜻蜓队长

小叶子 发表于 2017-5-11 16:40
站点转格点程序是单独的吗

一共需要两个程序吧