八步完成，精确处理台站资料~到插值

dcy_23 · 发表于 2012-7-27 21:58:27

正需要,谢谢分享!

mikado · 发表于 2012-7-27 22:19:27

辛苦了，支持一个

月流火 · 发表于 2012-8-31 15:41:49

这个一定要顶一下。。。有本事再来个穷神，哼!

liuzf406 · 发表于 2012-9-7 20:56:21

版主您好, 我用您的方法来处理从国家气候中心下载的1954年1月-2012年7的月降水数据, 第3步出现了问题, 能帮我诊断一下吗?

代码如下:
program ex00
implicit none
integer nlev,nflag,i,j
real lat(160),lon(160),tim,sd(160)
character*8 stid(160)
tim=0.0
nflag=1
open(10,file='61year_160id.txt')
open(20,file='160_sd_1951_2011.txt')
open(30,file='160_sd_1951_2011.grd',form='binary')
do i=1,160
read(10,*) stid(i),lat(i),lon(i)
enddo
close(10)
do i=1,739
   nlev=1
   do j=1,160
read(20,*) sd(j)
write(30) stid(j),lat(j),lon(j),tim,nlev,nflag,sd(j)
  enddo
  nlev=0
  write(30) stid(160),lat(160),lon(160),tim,nlev,nflag
  write(*,*) "Finish writting:",i
enddo
close(20)
close(30)
end

运行结果如下:

做个霸气的木头 · 发表于 2012-9-8 08:20:51

本帖最后由做个霸气的木头于 2012-9-8 08:36 编辑

liuzf406 发表于 2012-9-7 20:56

登录/注册后可看大图

版主您好, 我用您的方法来处理从国家气候中心下载的1954年1月-2012年7的月降水数据, 第3步出现了问题, 能帮 ...

问题出现在20号文件，就是你存储降水数据的文件，读到第38个时间点，提示在读取时文件已结束。说明你数据文件数据不是按照我的方式排列的，长度不够，我的数据要求是只有一列数据，739*160行的数据，数据存放顺序文件有说明的，你的应该不是~

liuzf406 · 发表于 2012-9-8 09:56:39

做个霸气的木头发表于 2012-9-8 08:20

登录/注册后可看大图

问题出现在20号文件，就是你存储降水数据的文件，读到第38个时间点，提示在读取时文件已结束。说明你数据 ...

谢谢版主指点, 我是从气候中心下载的月降水数据, 该如何存放呢?
如1月份的为:
3 4 2 4 2 1 4 5 0 0 1 0 1 0 1 4 0 1 2 0
0 6 1 3 3 4 3 1 2 0 0 1 8 2 7 6 4 9 11 7
9 22 13 9 18 17 12 24 32 5 14 13 30 20 33 27 26 19 66 34
..........

2月份的为:
3 6 2 4 2 3 4 7 5 22 16 13 2 48 23 5 26 26 3 26
41 13 3 3 5 4 8 11 25 28 6 1 10 8 14 4 7 22 10 12
10 16 15 8 27 59 27 31 33 21 17 70 36 94 48 79  112 85 99 99
  ..........

看以前站内的帮助说是 "一般要将所要的要素，如日照时间，单个要素放到一个txt文件中，数据的顺序是：同一个时间点371个站数据，然后放下一个时间点的371个要素值，如此直到最后时刻"

那是不是这些数据应该这样放啊:
3 4 2 4 2 1 4 5 0 0 1 0 1 0 1 4 0 1 2 0
0 6 1 3 3 4 3 1 2 0 0 1 8 2 7 6 4 9 11 7
9 22 13 9 18 17 12 24 32 5 14 13 30 20 33 27 26 19 66 34
..........(1月份的)
3 6 2 4 2 3 4 7 5 22 16 13 2 48 23 5 26 26 3 26
41 13 3 3 5 4 8 11 25 28 6 1 10 8 14 4 7 22 10 12
10 16 15 8 27 59 27 31 33 21 17 70 36 94 48 79  112 85 99 99
  ..........(2月份的)
  3月份数据....

就这样排列放完12个月,放在一个文件内?

结果是数据有5912行*20列=118240.  从1月到7月, 每个月有496行*20列=9920, 从8月到12月, 每个月488行*20列=9760 (没有2012年的数据).

请指教!

做个霸气的木头 · 发表于 2012-9-8 11:08:29

liuzf406 发表于 2012-9-8 09:56

登录/注册后可看大图

谢谢版主指点, 我是从气候中心下载的月降水数据, 该如何存放呢?
如1月份的为:
3 4 2 4 ...

你这些数字把我绕晕乎了。你没看懂我的回复吧亲
我现在给你举个数据例子，看下面的截图，我程序要调用的数据文件是这样的搜狗截图_2012-09-08_10-58-43.png

文件中数据只能有一列数，放的顺序是：第一个数据 15 是1954年1月第一个站的，第二个数据23是1954年1月份第二个站的，依次类推，直到1954年1月最后一个站，接着开始二月份的，52是1954年2月第一个站的，85是1954年2月第二个站的，如此类推，直到2012年12月的最后一个站的降水。

没办法，只有数据时这样排列才能用我的程序，你得通过程序将数据改成这样的形式，或者修改读取部分的程序，不过站点数据的排列就是这样一个顺序。

liuzf406 · 发表于 2012-9-8 11:22:15

做个霸气的木头发表于 2012-9-8 11:08

登录/注册后可看大图

你这些数字把我绕晕乎了。你没看懂我的回复吧亲
我现在给你举个数据例子，看下面的截图，我程序要调用的 ...

谢谢, 这下明白你原始数据的存放格式了.
对于气候中心下载的数据, 用你原来的代码处理, 还需要做很多工作啊!

谢谢!

做个霸气的木头 · 发表于 2012-9-8 11:24:59

liuzf406 发表于 2012-9-8 11:22

登录/注册后可看大图

谢谢, 这下明白你原始数据的存放格式了.
对于气候中心下载的数据, 用你原来的代码处理, 还需要做很多工 ...

一般通用的程序都较少，主要是咱俩处理的数据不一样，但是最终形式一样，grads就能识别了~

liuzf406 · 发表于 2012-9-8 11:46:33

做个霸气的木头发表于 2012-9-8 11:24

登录/注册后可看大图

一般通用的程序都较少，主要是咱俩处理的数据不一样，但是最终形式一样，grads就能识别了~

是啊, 谢谢您的指教!

		自动登录	找回密码
密码			立即注册