请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5037|回复: 18

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

[复制链接]

新浪微博达人勋

发表于 2012-4-5 13:18:48 | 显示全部楼层 |阅读模式

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

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
向各位高手们请教:跪求将micaps第三类降水数据转GrADs的FORTRAN程序。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-4-7 13:26:45 | 显示全部楼层
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


密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-4-7 15:04:53 | 显示全部楼层
帖子里以前有人发过的,你搜就有!祝好运!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-4-7 22:58:25 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-4-7 22:59:23 | 显示全部楼层
突突 发表于 2012-4-7 15:04
帖子里以前有人发过的,你搜就有!祝好运!

谢谢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-7-26 17:46:30 | 显示全部楼层
十分感谢程序提供,给予我莫大的帮助
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-10-16 21:33:56 | 显示全部楼层
确实很有用.
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-3-18 16:40:13 | 显示全部楼层
御用闲人 发表于 2012-4-7 13:26
diamond 3 10年04月10日08时24小时降水量
10 04 10 08   -2
      8

开心 正不知道错哪儿呢 看到你这句: (这里有一段read(1,*),被我干掉了,多少个,取决于资料第一行到上面打星的行数)
懂了。。 谢谢~!~!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-4-2 17:01:00 | 显示全部楼层
稍作修改后运行正常,谢了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-5-2 16:09:43 | 显示全部楼层
不错
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表