fortran语言将micaps第11类数据转换grads数据格式

whw334455

好不容易才把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)

2014-3-2 22:23 上传

点击文件名下载附件
micaps原始数据

whw334455

请各位多多指教

showme

你的意思是哪里出问题了？

lqouc

最好贴上来效果图哈

laudtim

谢谢楼主分享，来看看怎么样

mofangbao

楼主如果是提问请说明你的问题 要不然还以为你是来分享代码的 你是昨天想QQ找我问的吧？

whw334455

之前我发了一帖子求助，后来自己摸索，搞通了。现在这帖子是分享源代码程序。

whw334455

H:\micaps\micaps11tograds\uv.gif

aa22556688

楼主真是好人。

whw334455

aa22556688 发表于 2014-3-5 18:23

                               
登录/注册后可看大图

楼主真是好人。

谢谢   希望能帮到你