从全国降水的站点资料中提取某个区域的值

polestar7

求各位大神帮帮忙吧，我这里有一个降水资料（txt)，全国525个台站，有台站编号以及台站经纬度，降水资料从1959年1月到2009年12月，资料是月平均的，现在想截取lat(21-30)，lon(97-107)这个区间的台站降水资料，关于截取部分不太明白，txt文件格式及我的编程如下，请各位帮忙看一下吧。

全国降水资料

program rain
   implicit none
   integer:: iyr(52),imo(12)
    integer:: i,j,k
   real :: lon(525),lat(525),rr(525,12,52),
   integer,allocatable::a(:,:,:)，rm(:,:,:)
   character*5 stn(525)
   character*55 filein,fileout
   filein="D:\data to stud\rainfall\station\1\rain_quanguo_mon_grads.txt"
   fileout="D:\data to stud\rainfall\station\1\rain_djf_525stn.txt"
   open(10,file=filein,status='old',form='formatted')
   open(11,file=fileout,status='old',form='formatted')
    rr=0.0
    do k=1,51
do j=1,12
do i=1,525
    read(10,400)iyr(k),imo(j),stn(i),lat(i),lon(i),rr(i,j,k)
    enddo
enddo
enddo
write(*,400)iyr(51),imo(12),stn(525),lat(525),lon(525),rr(525,12,1)
write(*,*) rr(525,1,2),rr(525,2,2)
400 format(i4,2x,i2,2x,a5,f8.2,f8.2,2x,f8.1)  !逐月数据读入格式

  
  do k=1,51
  do j=1,12
  do i=1,525
  if(lat(i)>=21.00.and.lat(i)<=30.00.and.lat(i)>=97.00.and.lat(i)<=107.00)
    write(11,400)iyr(i),imo(12),stn(i),lat(i),lon(i),r(:,12,51)
end if
enddo
   enddo
enddo
   
close(10)
close(11)
stop
end

残月

额 直接没明白你在问什么。。是结果不对？程序不对？既然你自己写的为什么不明白？

polestar7

残月 发表于 2012-3-26 14:43

                               
登录/注册后可看大图

额 直接没明白你在问什么。。是结果不对？程序不对？既然你自己写的为什么不明白？

就是程序不对，下面截取区域部分出错，不知道该怎么弄

mofangbao

残月 发表于 2012-3-26 14:43

                               
登录/注册后可看大图

额 直接没明白你在问什么。。是结果不对？程序不对？既然你自己写的为什么不明白？

真希望楼主提问前自己先弄清楚自己已经有的东西

lotina

我没有看楼主的程序，但是我知道你的意思。不用程序反而简单一点啊，直接就在TXT里把你需要的范围内的台站号复制到另外一个TXT里就好了，你的范围不大，台站不多。然后直接用程序读取新建的TXT里的台站号为一个数组。在你读取源文件的时候，IF语句中就直接用这个为条件读资料就好了。

言深深

我看了一下，基本上程序是没有什么问题的，应该能跑出来！
不知道你想知道啥？程序主要由两部分构成，400以前的是读入数据，并且写程序的人还在屏幕输出了，这让你跑的应该比较吃力；2，if语句判断经纬度在你需要提取的范围以内，是正确的！
跑不出来的话就改改文件路径什么的吧，在不对你就把错误截图发上来吧

残月

mofangbao 发表于 2012-3-26 17:56

                               
登录/注册后可看大图

真希望楼主提问前自己先弄清楚自己已经有的东西

额 这个不是在说我吧我错了，我刚开始是真不明白他在问什么

残月

言深深 发表于 2012-3-26 20:28

                               
登录/注册后可看大图

我看了一下，基本上程序是没有什么问题的，应该能跑出来！
不知道你想知道啥？程序主要由两部分构成，400以 ...

按楼主所需要的数据的话，这句  write(11,400)iyr(i),imo(12),stn(i),lat(i),lon(i),r(:,12,51)
有点不恰当吧？？再者，莫非写数据之前 fileout 已经创建了？（status='old'）
还有400之前的那个只是测试下吧，就几个数，也不算是刷屏吧？

ps：深深我上次去祝贺来着，结果没给我红包~故意的

polestar7

言深深 发表于 2012-3-26 20:28

                               
登录/注册后可看大图

我看了一下，基本上程序是没有什么问题的，应该能跑出来！
不知道你想知道啥？程序主要由两部分构成，400以 ...

非常感谢，下午接着弄了一段时间，已经出结果了

言深深

polestar7 发表于 2012-3-26 21:49

                               
登录/注册后可看大图

非常感谢，下午接着弄了一段时间，已经出结果了

好，搞定了，可以改一下你的帖子内容和标题，一面再度有别人回复