MICAPS第四类格点数据转换为grads格点数据求助

hanqianli

想把MICAPS欧洲中心细网格第四类格点数据转换为grads格点数据（多时次，多层次，多变量），写了一个Fortran程序，第一个变量转换为grads数据文件是能够正确转换的，但第二个和第三个变量在数据转换过程中，高度层次有问题，求论坛里的各位大侠指教！！！
program main
integer,parameter::x=361,y=281,z=7,g=25
real u(x,y,z,g),v(x,y,z,g),t(x,y,z,g)
character(8)::day
character(4)::cen(z)
character(4)::hour(g)
character(6)::typ(4)
data cen/'1000','925','850','700','500','200','100'/
data hour/'.000','.006','.012','.018','.024','.030','.036','.042','.048','.054','.060','.066','.072','.078','.084','.090','.096','.102','.108',&
'.114','.120','.126','.132','.138','.144'/
data typ/'U','V','T','R'/
print*,'输入日期（格式:yymmddhh,15012820）'
read(*,*)day
! print*,day
print*,'输入文件数（1-25）'
read(*,*)gg
! print*,g

  do k=1,z
   do l=1,gg
    open(10,file='H:\ecmwf_thin\'//trim(typ(1))//'\'//trim(cen(k))//'\'//day//hour(l),status='old')
    read(10,*)
    read(10,*)
    read(10,*)
    do i=1,y
    read(10,*) (u(j,i,k,l),j=1,x)
    enddo
    close(10)
    open(10,file='H:\ecmwf_thin\'//trim(typ(2))//'\'//trim(cen(k))//'\'//day//hour(l),status='old')
    read(10,*)
    read(10,*)
    read(10,*)
    do i=1,y
    read(10,*) (v(j,i,k,l),j=1,x)
    enddo
    close(10)
    open(10,file='H:\ecmwf_thin\'//trim(typ(3))//'\'//trim(cen(k))//'\'//day//hour(l),status='old')
    read(10,*)
    read(10,*)
    read(10,*)
    do i=1,y
    read(10,*) (t(j,i,k,l),j=1,x)
    enddo
    close(10)
   enddo
  enddo
  open(12,file='D:\ectin\d\'//day//'.grd',status='unknown',form='binary')
  do l=1,gg
   do k=1,z
    do i=y,1,-1
    write(12) (u(j,i,k,l),j=1,x)
    enddo
    do i=y,1,-1
    write(12) (v(j,i,k,l),j=1,x)
    enddo
    do i=y,1,-1
    write(12) (t(j,i,k,l),j=1,x)
    enddo
   enddo
  enddo
  close(12)
end

hanqianli

DSET  d:\ectin\d\15020820.grd
TITLE Grid Data Sample
UNDEF 9999
XDEF 361 LINEAR 60  0.25
YDEF 281 LINEAR -10  0.25
ZDEF  7 LEVELS 1000 925 850 700 500 200 100  
TDEF  25 LINEAR 00Z04feb2015 6hr  
VARS  3
u     7  99  Grid Data
v     7  99  Grid Data
t     7  99  Grid Data
ENDVARS

hanqianli

'reinit'
'open d:/ectin/d/h.ctl'
i=1
while(i<=25)
'set cint 5'
'set z 1'
'set t 'i
'd t'
'printim d:/ectin/d/t/'i'.gif  x1000 y800 white'
'c'
i=i+1
endwhile
;

mofangbao

遇到问题先搜索：http://bbs.06climate.com/forum.php?mod=viewthread&tid=1598

hanqianli

感谢2楼回复，我仔细看了你的帖子，好像层次问题和时次问题没有解决啊，我想把多时次，多层次，多变量数据统一放在一个文件里以方便grads绘图。

瓶子

请问楼主解决了吗？是fortran脚本有问题？还是ctl有问题呀？

边缘人糊涂乐

厉害厉害{:eb502:}{:eb502:}{:eb502:}{:eb502:}