- 积分
- 1631
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-11-26
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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) |
|