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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 20911|回复: 32

[源代码] fortran语言将micaps第11类数据转换grads数据格式

[复制链接]

新浪微博达人勋

发表于 2014-3-2 22:27:34 | 显示全部楼层 |阅读模式

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

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

x
好不容易才把micaps第11类数据转换成grads的格式,并成功出图;1、fortran程序如下(用fortran4.0编译的);若用6.0或者6.6版本请把(ecl=361*241*4)中乘于4删除
program main
implicit none
integer i,j
integer,parameter :: ii=361
integer,parameter :: jj=241
character ch*100

real :: u(ii,jj),v(ii,jj)
open(10,file='H:\micaps\micaps11tograds\13111408.000')
read(10,*) ch
print*,ch
read(10,*) ch
print*,ch
read(10,*) ch
print*,ch

do j=1,jj
read(10,*) (u(i,j),i=1,ii)
enddo

do j=1,jj
read(10,*) (v(i,j),i=1,ii)
enddo

close(10)
!!!!注意micaps是从北向南排列数据,而grads是从南向北排列数据。
open(20,file='H:\micaps\micaps11tograds\13111408.grd',form='unformatted',access='direct',recl=361*241*4)
write(20,rec=1) ((u(i,j),i=1,ii),j=jj,1,-1)
write(20,rec=2) ((v(i,j),i=1,ii),j=jj,1,-1)
close(20)
stop
end


2、配套的ctl
dset H:\micaps\micaps11tograds\13111408.grd
undef -999
title dattest
xdef 361 linear 60.0   0.25
ydef 241 linear 0.0    0.25
zdef 1 levels 500
tdef 1 linear 08Z14nov2013 1dy
vars 2
u  0 99 uwind
v  0 99 vwind
endvars



3、配套的gs文件
'reinit'
'open h:\micaps\micaps11tograds\13111408.ctl'
'enable print H:\micaps\micaps11tograds\uv.gmf'
'set grads off'
'set vpage 0 11 0 8.5'
'set parea 0.5 10.5 1.8 7.6'
'set strsiz 0.15'
*'set lat 30 40';'set lon 80 100'
'set mpdset  lowres'
'set gxout stream'
*'set cint 4'
'd u;v'
'print'
'disable print‘






13111408.000 (1021.08 KB, 下载次数: 107)

评分

参与人数 3金钱 +30 贡献 +8 收起 理由
暴走小摔哥 + 5 很给力!
mofangbao + 15 + 5
lqouc + 10 + 3 感谢分享

查看全部评分

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

新浪微博达人勋

 楼主| 发表于 2014-3-2 22:30:42 | 显示全部楼层
请各位多多指教
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-3-2 23:20:54 | 显示全部楼层

回帖奖励 +1 金钱

你的意思是哪里出问题了?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-3-3 07:56:28 | 显示全部楼层
最好贴上来效果图哈
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-3-3 08:31:10 | 显示全部楼层
谢谢楼主分享,来看看怎么样
密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2014-3-3 09:17:23 | 显示全部楼层
楼主如果是提问请说明你的问题 要不然还以为你是来分享代码的 你是昨天想QQ找我问的吧?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-3-3 13:43:30 | 显示全部楼层
之前我发了一帖子求助,后来自己摸索,搞通了。现在这帖子是分享源代码程序。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-3-3 13:51:07 | 显示全部楼层
H:\micaps\micaps11tograds\uv.gif
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-3-5 18:23:05 | 显示全部楼层
楼主真是好人。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2014-3-5 18:47:25 | 显示全部楼层
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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