格点提取中国区域数据

言深深

写在最后········声明：刚刚在清风的提醒下才发现，原来手里的那份数据原来也只是一个规则的区域29*17=493
这样程序可以直接简化如下：

parameter(n=73*144)
open(2,file="data.txt")
do i=1,n
read (2,*) (x(i,j),j=1,3)
enddo
close(2)

open(3,file="picc.txt")
do ii=1,n
if (x(ii,1).ge.70.and.x(ii,1).le.140.and.x(ii,2).ge.15.and.x(ii,2).le.55)then
write(3,"(2f12.1,f12.2)") (x(ii,j),j=1,3)
endif
enddo
close(3)


end

忏悔去了{:soso__4029110217240599864_2:}
！！！！！！！！！！！！！！

昨天晚上看到大家在群聊讨论如何从格点数据中提取中国区域内的数据，正巧手里有一份已经提取好的中国境内数据，于是简单写了个小程序实现从2.5*2.5全球格点数据中提取适合中国的部分用于研究，希望对大家有帮助。
中国国内一共493个格点，数据存放在 china.txt (10.59 KB, 下载次数: 361)

2011-11-28 13:59 上传

点击文件名下载附件 文件里面，全球格点数据存放在 data.txt (390.09 KB, 下载次数: 250)

2011-11-28 13:59 上传

点击文件名下载附件 里面，为了方便阅读，所有的格点资料已经转换成经纬度，大家可以看一下。用程序 World 2 China.f90 (1 KB, 下载次数: 212)

2011-11-28 14:09 上传

点击文件名下载附件 选择读取就ok啦

!---------------------------------------------------------------------------------
!        从ncep/ncar2.5*2.5格点资料提取国内数据
!---------------------------------------------------------------------------------
parameter(m=493,n=73*144)
real china(m,2),x(n,3)
!---------------------------------------------------------------------------------
open(1,file="china.txt")
do i=1,m
read(1,*) (china(i,j),j=1,2)                                                                !        读取国内格点
enddo
close(1)
!---------------------------------------------------------------------------------
open(2,file="data.txt")
do i=1,n
read (2,*) (x(i,j),j=1,3)                                                                        !        读取全球格点
enddo
close(2)
!---------------------------------------------------------------------------------
open(3,file="pic.txt")
do i=1,m
do ii=1,n
if (x(ii,1).eq.china(i,1).and.x(ii,2).eq.china(i,2))then        !        选择国内格点格点
write(3,"(2f12.1,f12.2)") (x(ii,j),j=1,3)
endif
enddo
enddo
close(3)
!---------------------------------------------------------------------------------
print*,"----------------------"
print*,"程序来自http://bbs.06climate.com"
print*,"----------------------"
end

!---------------------------------------------------------------------------------

结果放在 pic.txt (18.29 KB, 下载次数: 89)

2011-11-28 13:59 上传

点击文件名下载附件 里面，下面是全球格点数据绘图

提取的中国部分

尽头的尽头

深深又一大作，这个沙发我必须来，哈哈

topmad

我来坐深深的沙发，看了深深的东东，感觉 好多东西要学习啊
深深，你要好好教教我啊

言深深

先占一楼再说

topmad

我来坐深深的沙发，看了深深的东东，感觉 好多东西要学习啊
深深，你要好好教教我啊

言深深

尽头的尽头 发表于 2011-11-28 14:12

                               
登录/注册后可看大图

深深又一大作，这个沙发我必须来，哈哈

不是我做的沙发么？肿么又没占到

传说中的谁

死鬼，周末不配伦家逛街。原来在折腾这个啊

mofangbao

先顶一下  然后表示这个提取出来和grads中是一样的，

言深深

传说中的谁 发表于 2011-11-28 14:14

                               
登录/注册后可看大图

死鬼，周末不配伦家逛街。原来在折腾这个啊

555555555

言深深

mofangbao 发表于 2011-11-28 14:24

                               
登录/注册后可看大图

先顶一下  然后表示这个提取出来和grads中是一样的，

好吧，诚恳的接受批评······