求micaps第三类降水数据转GrADs的FORTRAN程序

小飘要努力

向各位高手们请教：跪求将micaps第三类降水数据转GrADs的FORTRAN程序。

御用闲人

diamond 3 10年04月10日08时24小时降水量
10 04 10 08   -2
      8
      1      5     10     25     50    100    150    200
      1     25     39
  122.7   37.3  122.5   36.7  120.8   36.2  119.3   35.0  120.4   34.3
  121.1   32.7  122.3   30.9  122.5   30.1  121.5   28.0  119.8   25.0
  116.1   22.6  111.0   21.0  111.0   18.6  109.6   18.1  108.3   18.5
  108.9   19.9  109.3   21.0  106.4   21.6  104.1   22.2  101.6   20.8
   98.7   22.0   97.1   24.1   97.7   27.0   91.7   26.6   85.7   27.7
   78.9   30.6   73.3   37.1   74.1   40.7   87.6   49.8   98.5   43.9
  110.4   44.3  122.0   54.7  135.7   48.3  131.5   42.2  125.2   39.4
  121.3   38.5  121.2   40.3  118.5   38.6  120.8   38.1
    1    732                                                                                   ************read（1，*）
50349  124.33   51.70  496    .01
59265  111.30   23.48  120     .2
59647  109.10   21.03   55    .01
52633   98.41   38.81 3362     .2
57922  107.55   25.83  971    1.8
57647  109.18   28.15  407    .01
57710  106.41   27.81  879     .2
57816  106.72   26.58 1074    .01
52707   93.51   36.66 2772    .01
54436  118.95   40.40  228    2.0
59954  110.03   18.50   12     .4
58600  115.36   28.86   80     .2
57894  114.16   26.58  838     .2
57598  114.58   29.03  147    .01
.
.
.
.
.
_________________________________________________________以上是 资料样本
以下是 fortran源码，是mofan版主的，在这里感谢版主指导的grads站点作图
Program sta2grd
Implicit none
!这里是程序的变量声明

!变量生命结束

!程序开始
   
  Character*8 stid
  Real lon,lat,var,rain,tim
   integer nlev,flag

   tim =0.0
        nlev=1
        flag=1
  Open(1,file='11031408.000',status='old')
  open(2,file='sta.grd',status='replace',form='binary')
    （这里有一段read（1,*），被我干掉了，多少个，取决于资料第一行到上面打星的行数）
    10 Read(1,*,end=100)stid,lon,lat,var,rain
     write(2)stid,lat,lon,tim,nlev,flag,rain
!         print*,stid,lon,lat,var,rain
          
         ! pause          
      Goto 10
              100 continue
                  Close(1)
                  nlev=0
                  write(2)stid,lat,lon,tim,nlev,flag
                  close(2)
!程序结束
      End

突突

帖子里以前有人发过的，你搜就有！祝好运！

小飘要努力

御用闲人 发表于 2012-4-7 13:26

                               
登录/注册后可看大图

diamond 3 10年04月10日08时24小时降水量
10 04 10 08   -2
      8

谢谢啊！我已经执行处来了

小飘要努力

突突 发表于 2012-4-7 15:04

                               
登录/注册后可看大图

帖子里以前有人发过的，你搜就有！祝好运！

谢谢

dhdh111222

十分感谢程序提供，给予我莫大的帮助

gaofengjlsy

确实很有用.

y6310985

御用闲人 发表于 2012-4-7 13:26

                               
登录/注册后可看大图

diamond 3 10年04月10日08时24小时降水量
10 04 10 08   -2
      8

开心 正不知道错哪儿呢 看到你这句： （这里有一段read（1,*），被我干掉了，多少个，取决于资料第一行到上面打星的行数）
懂了。。 谢谢~！~！

怿燃

稍作修改后运行正常，谢了

wylbdxs

不错