爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 6077|回复: 3

查看wrf输出结果求助

[复制链接]
发表于 2015-5-4 13:15:14 | 显示全部楼层 |阅读模式

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

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

x
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!将arwpost输出结果(.dat文件)转化为txt文件的fortran程序!!!!!!!!!!!!!!
program main
implicit none
integer,parameter:: nx=250,ny=207,nz=19,nt=1
real,allocatable :: U(:,:)
integer :: i,j,t,jrec

allocate(U(nx,ny))

jrec = 0
open(101,file='D:/d01.dat',form='unformatted',access='direct',recl=nx*ny)
open(102,file='d01.txt')

do t=1,nt
jrec=jrec+1
  read(101,rec=jrec)((U(i,j),i=1,nx),j=1,ny)  
          write(102,*)((u(i,j),i=1,nx),j=1,ny)
enddo
close(101)
close(102)
deallocate(U)
end

!!!!!!!!!!转化后地面水平风分量如图!!!!!!!!!!!!!!!!!
请问,数值明显不对吧,为什么会这样呢?多谢。



捕获.PNG
密码修改失败请联系微信:mofangbao
发表于 2015-5-4 15:41:20 | 显示全部楼层
是用错了big_endian / little_endian的问题吧
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2015-5-4 16:31:35 | 显示全部楼层
letsdo8 发表于 2015-5-4 15:41
是用错了big_endian / little_endian的问题吧

能具体说说是什么,怎么修改吗,我不太懂这个。。
密码修改失败请联系微信:mofangbao
发表于 2015-5-4 19:10:17 | 显示全部楼层
open(101,file='D:/d01.dat',form='unformatted',access='direct',recl=nx*ny)

在上面语句末尾加上“convert=big_endian”看看
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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