请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 2827|回复: 2

降水资料粗网格化

[复制链接]

新浪微博达人勋

发表于 2015-10-28 16:53:41 | 显示全部楼层 |阅读模式
GrADS
系统平台:
问题截图: -
问题概况: 分辨率为1*1的CRU降水资料,需要将其做成4*4网格的中国降水,每个网格的值是原来1*1网格的和,而且网格之间不要连着。用grads要怎么做呀?或者用Fortran怎么做?大家可以帮我看看吗?多谢!
我看过提问的智慧: 看过
自己思考时长(天): 2

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 绯宜物语 于 2015-10-29 15:31 编辑

'open f:\pr.ctl'
'set gxout fwrite'
'set fwrite f:\pr.grd'
n = 1
while(n<=105)   
j = 19   
while(j<=51)        
i = 73        
while(i<=133)         
'set t 'n''         
'set x 73 137'         
'set y 19 55'         
'define pr=asumg(pra,x='i',x='i+4',y='j',y='j+4')'         
'd pr'         
i=i+4        
endwhile   
j=j+4     
endwhile
n=n+1
endwhile   
'disable fwrite'
'reinit';
写了一个这样的gs文件,但是输出的都是缺省值。
大家可以帮我看看吗?麻烦给我提一点建议。多谢!

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-10-28 21:43:35 | 显示全部楼层
henhaodetiezi
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2015-10-30 14:15:47 | 显示全部楼层
本帖最后由 绯宜物语 于 2015-10-30 15:30 编辑

program cmip5
  integer*4,parameter:: x=360,y=181,t=105,m=17,n=10
  integer*4 :: ix,iy,im,in,it,irec
  real*4,parameter :: undef = -9.99E+8
  real*4 pra(x,y,t),pr(m,n,t)
  
  open(1,file='f:\pr\h\cru\prcsu.grd',form='binary')
  do it=1,t
      do iy=1,y
    do ix=1,x
     read(1) pra(ix,iy,it)
    enddo
   enddo
  enddo
  close(1)
  do it=1,t
     do im=1,m
     do in=1,n
                 ix=4*im+69
                 iy=4*in+15
                 pr(im,in,it)=pra(ix,iy,it)+pra(ix,iy+1,it)+pra(ix,iy+2,it)+pra(ix,iy+3,it)+pra(ix+1,iy,it)+pra(ix+1,iy+1,it)+pra(ix+1,iy+2,it)+pra(ix+1,iy+3,it)+pra(ix+2,iy,it)+pra(ix+2,iy+1,it)+pra(ix+2,iy+2,it)+pra(ix+2,iy+3,it)+pra(ix+3,iy,it)+pra(ix+3,iy+1,it)+pra(ix+3,iy+2,it)+pra(ix+3,iy+3,it)
        enddo            
   enddo
  enddo
  
        open(3,file='f:\prmm.grd',form='unformatted',access='direct',recl=m*n*4)
        irec = 0
        do it = 1,t
           irec = irec + 1
           write(3,rec=irec) ((pr(im,in,it),im=1,m),in=1,n)
        enddo
  close(3)
end
写了Fortran程序做,但是输出的值基本都是0.
有人可以帮我看看哪有问题吗?谢谢。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表