Fortran读取micaps风场数据的问题

袜子女巫

想读micaps十一类风场的EC预报场资料，写成二进制文件再用grads画出，数据能够顺利读出，但是画出的图中却没有v风量，想求助下问什么会这样，因为用同样的程序读高空实况风场就能顺利读出，二者唯一的差别是数据存放的方式，实况的是从北极往南放，ec是从赤道往北放。
程序：program main
implicit none
integer i,j
integer,parameter :: ii=145
integer,parameter :: jj=37
character ch*100

real :: u(ii,jj),v(ii,jj)
open(10,file='d:\micaps\11\17010108.024')
read(10,*) ch
print*,ch

!先放U分量，数据按先纬向后经向放（若为直角坐标网格数据，则先X方向，后Y方向），均为浮点数。所有格点的U分量放完后再放V分量，也是按先纬向后经向放。


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='d:\micaps\11\17010108.024.grd',form='unformatted',access='direct',recl=145*37*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)
end
再贴一个数据的图片
file:///C:\Users\Administrator\AppData\Roaming\Tencent\Users\382871867\QQ\WinTemp\RichOle\9@{HVAS9]3GM(6EDBW1Y$(7.png

xdzhao2013yk

谢谢，打扰了

tulalang

你仔细看看你发的micaps数据的文件头的说明，你想要读的这个数据排列是由0-90°N，意思就是和grads的排列顺序是一样的，你改改的fortran写出的顺序就好

袜子女巫

tulalang 发表于 2017-12-26 23:28
你仔细看看你发的micaps数据的文件头的说明，你想要读的这个数据排列是由0-90°N，意思就是和grads的排列顺 ...

恩，解决了，现在只要解决把-180-180改为0-360的存储方式就好了

单纯xiao可爱

袜子女巫 发表于 2017-12-27 09:55
恩，解决了，现在只要解决把-180-180改为0-360的存储方式就好了

前辈，想问下你怎么改的呀。我弄一天了老是在v那出现问题