grd文件ctl出错求解

追风的阳光 · 发表于 2013-5-22 21:39:26

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

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

x

program main
parameter(nx=21,ny=11,nz=26,nt=13,delta_lamda=3.14/180,delta_phi=3.14/180)! nx-lon(90E-140E);ny-lat(15N-60N);nt-time
real h(nx,ny,nz,nt),u(nx,ny,nz,nt),v(nx,ny,nz,nt),vor1(nx,ny,nz,nt),vor2(nx,ny,nz,nt),div(nx,ny,nz,nt),svbit(nx,ny,nz,nt)
real f(ny)
a=6371000.0    !radiu of earth
open(21,file='d:\zssdata\u.dat',form='binary')
do k=1,nt
do s=1,nz
do j=1,ny
read(21,rec=21)(u(i,j,s,k),i=1,nx)
enddo
end do
end do
close(21)
.........
open(5,file='d:\zssdata\data.grd',form='binary')
do k=1,nt
do s=1,nz
write(5)((svbit(i,j,s,k),i=1,nx),j=1,ny)
enddo
enddo

上面是程序的读写部分，下面是输出文件的ctl
ctl:
dset d:\zssdata\data.grd
undef 9.999E+20
title svbit
options template
xdef 21 linear 110 1
ydef 11 linear 40 1
zdef 26 levels
1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10
tdef 13 linear 00:00Z11jul2011 6hr
vars 1
svbit 0 99 vorticity using hgt
endvars

为什么运行ctl显示错误呢？是哪儿有问题？

松下风 · 发表于 2013-5-22 22:11:42

最好能把grads报错的图发上来，另外检查一下时间描述是否正确。

夏朗的芒果 · 发表于 2013-5-22 22:40:41

把报错图贴出来吧

river · 发表于 2013-5-23 07:11:44

ctl是需要使用open命令打开的

蔚蓝的天 · 发表于 2013-5-27 12:46:16

我也来顶一个

一闭一睁 · 发表于 2013-6-4 09:14:40

看帖，必回！！！！！！

eeyore · 发表于 2013-6-4 09:26:27

你的Z有26层，变量svbit后面怎么是0呢，应该是写26

Madao · 发表于 2016-4-23 19:46:49

学习写ctl文件中/////

l308417179 · 发表于 2016-4-24 22:28:53

我也是正在学习中

扫眉才子 · 发表于 2016-11-21 16:20:42

在那遥远的地方发表于 2016-11-21 10:41
您好，想请教您一下，grads处理grd格式的文件对应的ctl是怎么生成的呢？

自己写啊，建议多看看手册

		自动登录	找回密码
密码			立即注册

[求助] grd文件ctl出错求解

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