急~~文件名怎么循环，怎么能读数据不控制数据格式，直接用空格来控制

wyw · 发表于 2012-8-8 22:44:58

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

我正在做一个micaps第一类数据转第三类数据的程序，大块基本完成。现在要做实时调用数据的部分，每三小时调用文件一次，03-06调用02时的文件，06-09调用05时的文件，以此类推。
我现在只是先用一个文件做的，文件名：2012073014，不知道怎么调用不同文件名的数据（2012080802，2012080805。。。。2012083002。。。。。），想请教一下，多谢~·
程序1如下：
!=======================================================================
PROGRAM MAIN1
IMPLICIT NONE
!-----------------------------------------------------------------------
  integer::YEAR,MONTH,DAY,TIME !年月日录入
  !INTEGER,PARAMETER::JS=10000
  INTEGER::NUMBER,I
  INTEGER,ALLOCATABLE::NUM(:),HD(:) ！站点数，海拔
  REAL,ALLOCATABLE::LON(:),LAT(:),T(:),TD(:),TC(:),P(:),Q(:) ！经度、纬度、温度、露点、存放不规则的变量数组
  INTEGER TEMP

!-----------------------------------------------------------------------
  OPEN(12,FILE='D:\wx\12073014.000') !打开原始数据文件

  OPEN (13,FILE='D:\wx\SC.TXT',STATUS='UNKNOWN',form='formatted') !用来存放第三类数据的文件


  READ (12,*)

  READ (12,'(I2,3i3,I7)') YEAR,MONTH,DAY,TIME,NUMBER !CS 站点数
  PRINT*,NUMBER
  WRITE (13,'(A,I2,A,I2,A,I2,A,I2,A)') "diamond 3",YEAR,"年",MONTH,"月",DAY,"日" ,TIME,"时温度露点差" !题头
  WRITE (13,'(4I3,A)') YEAR,MONTH,DAY,TIME," -3"
  WRITE (13,'(A)') "    0"
  WRITE (13,'(A)') "    1    25    0"
  WRITE (13,'(A,I10)') " 1",NUMBER !题头五行
!-----------------------------------------------------------------------
  ALLOCATE(NUM(NUMBER))
  ALLOCATE(HD(NUMBER))
  ALLOCATE(LON(NUMBER))
  ALLOCATE(LAT(NUMBER))
  ALLOCATE(T(NUMBER))
  ALLOCATE(TD(NUMBER))
  ALLOCATE(TC(NUMBER))
  ALLOCATE(P(NUMBER))
  ALLOCATE(Q(NUMBER)) !给数组分配内存空间
!-----------------------------------------------------------------------------------------------------------------------------

  DO I=1,NUMBER
READ (12,'(I5,2F8.2,I5,8I5)') NUM(I),LON(I),LAT(I),HD(I),TEMP,TEMP,TEMP,TEMP,TEMP,TEMP,TEMP,TEMP
READ (12,'(F7.1,3I5,F7.1,F7.1,I5,F7.1,6I5)') Q(I),TEMP,TEMP,TEMP,TD(I),P(I),TEMP,T(I),TEMP,TEMP,TEMP,TEMP,TEMP,TEMP
  END DO !读入原始资料数据，两行是一个站点资料
  DO I=1,NUMBER
IF(T(I)==999.9000.OR.TD(I)==999.9000)THEN
TC(I)=9999
  ELSE

TC(I)=T(I)-TD(I)
END IF  !判断温度露点是否有缺省值
END DO
!-----------------------------------------------------------------------
DO I=1,NUMBER
WRITE (13,'(I5,2F8.2,I5,F7.1)') NUM(I),LON(I),LAT(I),HD(I),TC(I) !将数据读入sc.txt

END DO

CLOSE(12)
CLOSE(13)
END PROGRAM
!-----------------------------------------------------------------------
请高手指教，我是菜鸟，现学现写的程序，希望回答能写细点，能给段通用程序最好啦~~

我还有做另外一个程序，要求和上面一样，但是每列的数据格式不一样，不是列数不同，是同一列每行的数据格式都不一样，比如9999，可能突然变成0.1，下面又出现103.223，若控制输出按较大的f8.1就会读错数据，这种情况能利用空格来读入数据吗？可以的话能写的详细点吗？

以上两个数据文件挂到附件了，“.000”是程序1调用的文件
                           “.aws.txt”是格式乱的那个文件（后缀应该就是点.txt）

传不上去，我摘一部分放下面

wyw · 发表于 2012-8-8 22:52:34

2012073014.txt内容如下

diamond 1 12年07月30日14时地面填图
12 07 30 14 4507
48820  105.80 21.02 6 2 7 50 1 1  -18 6 2
9999 38 6 9999 24.5 9999 3 29.6 24 9999 1 2 9999 9999
48823  106.15 20.43 3 32 7  160 4 7  -16 0 0
9999 32 7  300 25.6 9999 0 30.4 9999 9999 1 2 9999 9999
48870  109.22 13.77 6 2 8  290 1 31  -12 6 2
9999 38 8 9999 22.5 9999 21 29.5 9999 9999 1 2 9999 9999
48887  108.10 10.93 8 4 7  270 6 66  -15 0 0
9999 38 7  600 9999 20.0 0 30.2 9999 9999 1 2 9999 9999
48845  105.67 18.70 6 2 7  180 3 8  -12 0 0
   0 38 7  600 22.2 9999 0 32.1 9999 9999 1 2 9999 9999
48840  105.78 19.80 7 4 8 0 0  991  -24 0 0
9999 38 8  300 25.0 9999 0 30.4 9999 9999 1 2 9999 9999
48846  105.90 18.35 3 32 7  230 2 10  -10 0 0
   0 35 7 1500 24.5 9999 0 29.4 9999 9999 1 2 9999 9999
48855  108.18 16.03 7 4 7  320 2 20  -19 0 0
9999 38 7  600 23.0 20.0 0 33.2 9999 9999 1 2 9999 9999
48852  107.68 16.40 17 32 7  320 3 20  -14 0 0
   0 38 4  600 21.8 9999 0 31.8 23 18 1 2 9999 9999
48863  108.78 15.13 9 32 7  250 2 24  -14 0 0
9999 39 7  600 24.2 9999 0 32.2 9999 9999 1 2 9999 9999
48810  105.83 22.15  176 32 3 0 0  998  -24 0 0
9999 32 3  300 25.5 20.0 0 31.8 20 10 1 2 9999 9999
48831  105.83 21.60 32 32 4  230 2  997  -19 0 0
9999 32 4  300 23.8 20.0 0 31.5 20 10 1 2 9999 9999
54096  131.15 44.38  498 4 6  250 2  989 -6 0 0
   0 38 6  600 18.5 20.0 0 25.8 20 10 1 2 9999 9999
53764  111.10 37.50  951 16 8 70 2 21 4 9 6
20.0 39 6  600 20.1 20.0 95 20.8 27 9999 1 2 9999 9999
50973  130.55 45.75  217 32 8  320 3  997 5 0 0
   0 35 8 1500 17.3 15.0 0 20.3 9999 9999 1 2 9999 9999
54259  124.95 42.10  235 8 7  250 2 44 -6 6 0
0.4 35 7 1000 18.3 30.0 0 23.6 20 10 1 2 9999 9999
56106 93.78 31.88 4024 4 0 90 3 9999  -21 0 0
   0 30 0 9999 7.3 30.0 0 21.9 20 10 1 2 1 -2
55572 89.10 29.68 4001 32 6  250 2 9999  -12 0 0
   0 34 6 1500 9.2 30.0 0 18.4 20 10 1 2 0 0
56228 96.91 30.05 3261 32 3 70 2 9999  -30 0 0
   0 31 2 1000 3.7 30.0 0 27.2 20 12 1 2 1 -2
54326  119.70 41.38  423 16 7  340 1 54 -8 0 0
   0 31 4 1500 16.1 30.0 0 27.1 23 12 1 2 9999 9999
55690 91.95 27.98 4281 16 7  220 7 9999 -7 0 0
   0 32 6  600 7.5 30.0 0 13.1 24 12 1 2 2 -1
55696 92.47 28.42 3861 4 7  270 1 9999  -21 0 0
   0 32 3 1000 4.4 30.0 0 19.5 20 12 1 2 3 -1
56034 97.13 33.80 4418 16 7  320 4 9999  -15 0 0
   0 39 7  600 6.3 30.0 0 14.5 9999 9999 1 2 -2 -1
54195  129.76 43.33  241 16 7  220 2  998 -3 0 0
   0 35 7 1000 19.0 15.0 0 26.1 9999 9999 1 2 9999 9999
54266  125.63 42.53  341 32 7  270 4 45 4 0 0
   0 32 7  600 15.8 15.0 0 22.0 20 10 1 2 9999 9999
55472 88.63 30.95 4671 2 5  340 1 9999  -10 0 0
   0 34 5 1500 5.8 30.0 0 14.7 20 10 1 2 -1 0
55294 91.10 32.35 4801 4 7  220 3 9999 -8 0 0
   0 39 7 1000 4.9 30.0 0 11.5 9999 9999 1 2 -2 -1
55228 80.08 32.50 4279 2 5 70 3 9999  -11 0 0
   0 39 5 1500 1.1 30.0 0 22.0 20 10 1 2 1 0
54472  122.71 40.88 27 32 5  340 2 52 -6 0 0
   0 31 5 1000 19.6 20.0 0 26.2 20 10 1 2 9999 9999
54476  122.15 40.16 38 32 6  360 2 51 -1 0 0
   0 31 1 1000 19.9 20.0 0 26.8 20 12 1 2 9999 9999
54486  123.28 40.28 81 32 2 20 3 38 -6 0 0
   0 31 2  600 20.2 12.0 0 29.0 20 10 1 2 9999 9999
50968  127.97 45.22  191 16 8  220 3 12 15 8 0
3.0 35 8 1000 18.8 12.0 60 19.3 9999 9999 1 2 9999 9999
53772  112.55 37.78  779 2 8  250 1 12 -6 6 0
   .1 30 8 2500 21.8 8.0 60 25.5 22 9999 1 2 9999 9999
58608  115.55 28.06 30 16 3  140 3 9  -16 0 0
   0 31 3  600 22.0 20.0 0 35.8 20 10 1 2 9999 9999
56317 94.21 29.21 2952 32 6  360 2 9999  -26 0 0
   0 31 2 1500 10.1 30.0 0 25.0 20 12 1 2 4 -3
55597 91.68 29.03 3741 32 2  340 1 9999  -15 0 0
   0 31 2 1000 5.7 30.0 0 21.4 20 10 1 2 1 0
56202 93.28 30.66 4490 16 2  140 2 9999  -16 0 0
   0 31 2 1000 4.2 30.0 0 16.8 20 10 1 2 4 -2
50955  126.30 45.38  167 32 8  320 2 15 10 6 0
5.0 35 8 1500 18.4 14.0 60 18.7 9999 9999 1 2 9999 9999
56109 93.78 31.48 3941 32 5  110 2 9999  -22 0 0
   0 31 4 1500 5.9 30.0 0 24.6 20 12 1 2 1 -1
56223 95.83 30.75 3640 32 6 40 3 9999  -21 0 0
   0 32 6 1000 2.5 30.0 0 24.4 20 10 1 2 -1 -1
56331 97.83 29.66 3781 32 5 0 0 9999  -21 0 0
   0 32 5 1000 3.2 30.0 0 21.5 20 10 1 2 1 -2
56137 97.17 31.15 3307 1 7  220 2 9999  -26 0 0
   0 32 2 1000 4.1 30.0 0 26.3 20 12 1 2 1 -2
52895  104.68 36.56 1399 16 8 0 0 47 2 6 6
2.0 35 6 1500 15.7 30.0 21 19.6 27 9999 1 2  -10 5
56093  104.01 34.43 2316 16 8 40 3 9999 3 6 0
5.0 35 8  600 10.0 15.0 0 16.4 9999 9999 1 2  -10 2
54808  115.58 36.03 47 16 8  290 3 9 -7 0 0
.01 30 8 2500 24.1 12.0 60 29.7 27 9999 1 2 9999 9999
54823  116.98 36.68 58 2 6  110 3 11  -15 0 0
   0 33 2  600 24.4 25.0 0 32.9 20 12 1 2 9999 9999
52868  101.43 36.03 2238 16 8  290 1 9999 -1 8 6
   .6 35 8 1500 14.9 30.0 60 19.2 9999 9999 1 2  -10 4
56046 99.65 33.75 3968 2 5  320 2 9999  -17 0 0
   0 39 5  600 8.4 30.0 0 18.7 20 10 1 2 -1 -2
50745  123.92 47.38  148 2 7 40 2 48 7 6 0
1.0 35 6 1500 17.2 14.0 60 18.9 20 12 1 2 9999 9999
50953  126.77 45.75  143 1 7  360 1 7 5 0 0
   0 35 7 1500 17.5 15.0 0 21.1 9999 9999 1 2 9999 9999
50442  124.11 50.40  372 16 7 40 4 52 -2 0 0
   0 31 1 1000 10.8 20.0 0 27.6 20 12 1 2 9999 9999
50468  127.45 50.25  166 16 7 20 3 43 0 0 0
   0 30 0 9999 12.9 20.0 0 26.1 20 12 1 2 9999 9999
56116 95.60 31.42 3874 4 3 90 2 9999  -22 0 0
   0 31 2 1000 5.8 30.0 0 21.5 20 12 1 2 0 -2
53463  111.68 40.82 1065 2 8 70 2 48  -12 8 0
5.0 39 8  600 16.5 10.0 80 19.9 9999 9999 1 2 -8 1
55589 90.98 29.30 3556 32 7 20 2 9999  -14 0 0
   0 32 1 1500 7.9 30.0 0 23.0 20 18 1 2 2 -1
55248 84.41 32.15 4416 4 7  110 2 9999  -11 0 0
   .1 39 7 1500 5.0 30.0 25 17.0 9999 9999 1 2 1 0
54727  117.56 36.73 69 32 8  110 3 22 -2 9 8
14.0 39 8  600 24.1 10.0 80 26.6 9999 9999 1 2 9999 9999
57741  109.18 27.71  281 32 3  180 2  995  -16 0 0
   0 32 3  600 23.1 30.0 0 35.1 20 10 1 2 9999 9999
57932  108.53 25.97  287 16 6  250 2  997  -21 0 0
.01 32 6 1000 23.9 30.0 25 33.3 24 10 1 2 9999 9999
57902  105.18 25.43 1379 4 7  140 2 10  -12 0 0
   0 38 7  600 20.2 25.0 0 25.4 9999 9999 1 2 1 -3
50877  129.58 46.30  101 16 8 0 0  999 3 0 0
   0 39 8 1000 17.0 25.0 0 20.6 9999 9999 1 2 9999 9999
50844  123.42 46.40  150 16 7 20 7 38 3 6 0
   .4 35 7 1500 16.4 12.0 60 20.0 9999 9999 1 2 9999 9999
53588  113.53 39.03 2898 4 9  360 11 9999 -8 9 6
18.0 9999 9999 9999 12.8    .0 81 13.3 9999 9999 1 2 9999 9999
53578  112.43 39.30 1092 32 8 0 0 54 -3 8 6
2.0 37 5 1000 17.2 8.0 60 19.1 22 9999 1 2 9999 9999

言深深 · 发表于 2012-8-9 07:19:18

最简单的办法是将需要的文件名写在一个txt文档里面，然后用字符数组读取，然后程序执行循环，对不同文件进行处理：
do l=1,name
open(1,file=filename(l))
??
close(1)
enddo

也可以用trim语句进行拼接，具体可以参考mofangbao的一片帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=3079

mofangbao · 发表于 2012-8-9 08:42:27

读取时候格式符改为*，就可以默认用空格来分隔的。
批处理一种方法就像深深说的写入文档，另外你的这个情况，也可以写一个专门的函数来生成文件名，用switch根据时间段选择一下就行了，仅提供思路，程序自己写吧

wyw · 发表于 2012-8-9 22:35:53

言深深发表于 2012-8-9 07:19

登录/注册后可看大图

最简单的办法是将需要的文件名写在一个txt文档里面，然后用字符数组读取，然后程序执行循环，对不同文件进行 ...

我的是要动态生成文件名，你给的清风的帖子虽然很有用，但是我这种恰好只给了思路，菜鸟死路一条啊~~
彭国伦的书里面有讲关于获取系统时间函数和文件名循环的函数吗？大概是哪个部分，我真的是边学边写的程序，书只是做参考用的。。。。。

wyw · 发表于 2012-8-9 22:43:25

mofangbao 发表于 2012-8-9 08:42

登录/注册后可看大图

读取时候格式符改为*，就可以默认用空格来分隔的。
批处理一种方法就像深深说的写入文档，另外你的这个情况 ...

非常感谢你的帮助，用空格分隔这个问题我已经解决了，但是时间这个还是不太会。。。。（菜鸟我也很痛苦）思路我都有了解，但是就是不会操作
另外就是程序中要进行时间对比，知道要用gettime、cpu_time(time)，但是没用具体的格式套用举例，比如输出就是大致“print*，A” （A，任意需输出对象），那么gettime这个是怎么使用的？彭国伦的书有写吗，

mofangbao · 发表于 2012-8-9 22:46:21

你问的问题论坛里面都有讨论过，一个是我发帖的，是一个按指定字符来分割字符串的帖子
另一个是topmad的计算fortran运算耗时的帖子里面提到的用法，你可以自己搜索查看一下

言深深 · 发表于 2012-8-10 08:33:05

wyw 发表于 2012-8-9 22:35

登录/注册后可看大图

我的是要动态生成文件名，你给的清风的帖子虽然很有用，但是我这种恰好只给了思路，菜鸟死路一条啊~~
彭 ...

最简单的办法是将需要的文件名写在一个txt文档里面，然后用字符数组读取，然后程序执行循环，对不同文件进行处理：
do l=1,name
open(1,file=filename(l))
??
close(1)
enddo
用这个方法试试看吧

Forcast · 发表于 2012-8-15 13:07:01

FORTRAN处理日期还没有弄过。

不语-承受 · 发表于 2012-8-15 16:34:07

dir */b
看看是不是你需要的格式？这就是批处理的基本命令了。
写入文件：
dir */b>a.txt
看看a.txt的内容。
接下来用上面的文件读取函数就可以了。
这是最快最简单代码都全了的。

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

[求助] 急~~文件名怎么循环，怎么能读数据不控制数据格式，直接用空格来控制

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