爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 13489|回复: 16

[求助] fortran读取.dat文件

[复制链接]

新浪微博达人勋

发表于 2012-12-26 15:05:49 | 显示全部楼层 |阅读模式

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

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

x
各位先贤,最近在尝试读取ARWpost处理后的数据文件,很是迷茫,希望有懂得朋友给点帮助,不胜感激。我自己尝试写的数据如下:
Program main
  Implicit None
  integer ::a1,a2,a3,a4
  integer ::i
  Open( 12 , File = 'newbeijing7.dat' , Access = 'Direct' , Form = 'Unformatted' , RecL =4 )
  do i=1,100
  Read( 12 , Rec = i) a1,a2,a3,a4
         Write( * , * )a1,a2,a3,a4        
  enddo
  Close( 12 )
End Program main
新手,写的不太专业,对Recl这个参数不知道是多少,看论坛上有填xy方向个点数乘积的,还有说填pdef后面两个数乘积的,不懂了,但感觉读出的数据都不对,看了看grads手册,不知道怎么具体的读取。读取的数据感觉总是不对。或者有做个这个的朋友能给一些代码。多谢了!
.ctl描述如下:
dset ^test1.dat
options  byteswapped
undef 1.e30
title  OUTPUT FROM WRF V3.4.1 MODEL
pdef   73  60 lcc  34.830  -81.030   37.000   30.500  60.00000  30.00000  -98.00000  30000.000  30000.000
xdef  224 linear  -94.99892   0.13513514
ydef  166 linear   23.26908   0.13513514
zdef   27 linear 1 1
tdef    1 linear 00Z01JUL2012     360MN     
VARS   91
LU_INDEX       1  0  LAND USE CATEGORY (-)
VAR_SSO        1  0  variance of subgrid-scale orography (m2)
LAP_HGT        1  0  Laplacian of orography (m)
U             27  0  x-wind component (m s-1)
V             27  0  y-wind component (m s-1)
W             27  0  z-wind component (m s-1)
PH            27  0  perturbation geopotential (m2 s-2)
PHB           27  0  base-state geopotential (m2 s-2)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2012-12-26 15:36:07 | 显示全部楼层
关于RECL的含义,可以好好看彭国伦那本书,RECL没有固定值,也没有固定的算法,完全和你自己想要一次性读取的数据量的大小有关,说白了就好比你有一瓶水,你一口喝的量就是RECL,你一次喝的多,就少喝几口(少循环几次),一次喝的多,就得多喝几口,所以RECL的值和循环的次数是成比例的,和你数组的大小要相适应。
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

 成长值: 0
发表于 2012-12-26 15:40:03 | 显示全部楼层

喝大了···
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-12-26 16:19:52 | 显示全部楼层
mofangbao 发表于 2012-12-26 15:36
关于RECL的含义,可以好好看彭国伦那本书,RECL没有固定值,也没有固定的算法,完全和你自己想要一次性读取 ...

recl 就是一次要读取得数据量是这个含义吧,看书上说这个是记录数据时每次记录占用的字节数,我个人比较同意这个说法,因为每次我设定不同的值,读出的数据是不同的,也就是把不同字节的数据归为一次记录读出的数据是不同的。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2012-12-26 21:36:00 | 显示全部楼层
L0725 发表于 2012-12-26 16:19
recl 就是一次要读取得数据量是这个含义吧,看书上说这个是记录数据时每次记录占用的字节数,我个人比较同 ...

是这个意思
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-1-8 06:26:13 | 显示全部楼层
支持支持!!!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-9-21 09:32:41 | 显示全部楼层
学习学习!!
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2016-6-27 14:06:42 | 显示全部楼层
mofangbao 发表于 2012-12-26 15:36
关于RECL的含义,可以好好看彭国伦那本书,RECL没有固定值,也没有固定的算法,完全和你自己想要一次性读取 ...

先问下 arwpost之后的数据 recl 里面用不用*4
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-10-11 15:04:24 | 显示全部楼层
问题有解决吗?我也遇到了类似的问题
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-10-25 16:22:58 | 显示全部楼层
问题有解决吗?我也遇到了类似的问题
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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