grads格点资料画图（初学）

佩佩 · 发表于 2012-4-23 18:20:56

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

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

x

资料为6层151*151的资料，文本格式,只有一个变量
fortran为

parameter(nx=151,ny=151,s=6)
real,dimension(nx,ny,6,1)::dbz
integer it,k,i,j,irec
open(1,file="jg1.txt")
open(2,file="jg1.grd",form="binary",status="replace",access="direct",recl=nx*ny*4)
irec=0
do it=1,1
do k=1,6
read(1,100)((dbz(i,j,k,it),i=1,nx),j=1,ny)
irec=irec+1
write(2,rec=irec)((dbz(i,j,k,it),i=1,nx),j=1,ny)
enddo
enddo
close(1)
close(2)
100 format(151f12.6)
end

ctl为

DSET E:\grads\jg1.grd
TITLE LAGERANGE
UNDEF -999.0
XDEF 151 LINEAR -150 2
YDEF 151 LINEAR -150 2
ZDEF 6 LEVELS 2 3 4 5 6 7
TDEF 1 LINEAR 00Z02JAN1987 1DY
VARS 1
DBZ 6 99 基本反射率
ENDVARS

打开grads open e:\grads\jg1.grd
出错：missing or invalid dimension size
the invalid description file record is xdef 151 linear -150 2
the data file was not opened

我想画直角坐标系的，不想用经纬度的，不知道可不可以？连文件都打不开！

随缘 · 发表于 2012-4-23 18:46:06

XDEF 151 LINEAR -150 2
YDEF 151 LINEAR -150 2
你这个描述不知道什么意思啊？这个算下来得有151多个格点吧？你改成-150 0试试看

佩佩 · 发表于 2012-4-23 18:54:03

随缘发表于 2012-4-23 18:46

登录/注册后可看大图

XDEF 151 LINEAR -150 2
YDEF 151 LINEAR -150 2
你这个描述不知道什么意思啊？这个算下来得有151多个格点 ...

就是想要151个格点呢，改了也不行啊！求救啊！

mofangbao · 发表于 2012-4-23 18:59:52

你先不要输出到二进制，你先按照这个格式输出到文本文件看一下你的数据写出有没有什么问题，这种格点资料可以不用recl写的，直接write(1)var 这样就可以了，用recl不同的编译器可能还有区别

佩佩 · 发表于 2012-4-23 19:09:16

mofangbao 发表于 2012-4-23 18:59

登录/注册后可看大图

你先不要输出到二进制，你先按照这个格式输出到文本文件看一下你的数据写出有没有什么问题，这种格点资料可 ...

我试过了，输出文本格式的话，和原文件是一样的，我试试不用recl的、。。。。
还是不行啊，我怀疑是不是格点太多，grads不认啊？

mofangbao · 发表于 2012-4-23 19:24:27

佩佩发表于 2012-4-23 19:09

登录/注册后可看大图

我试过了，输出文本格式的话，和原文件是一样的，我试试不用recl的、。。。。
还是不行啊，我怀疑是不是 ...

你这个ydef乱写的吧？描述文件y方向最多只能有180度，但是格点个数是不限制的，你这个分辨率写的有问题！你不想用经纬网是在脚本画图的时候设置的，而不是在ctl文件中设置，具体命令你可以去查手册

随缘 · 发表于 2012-4-23 20:04:25

佩佩发表于 2012-4-23 18:54

登录/注册后可看大图

就是想要151个格点呢，改了也不行啊！求救啊！

XDEF 151 LINEAR 0 1
YDEF 151 LINEAR 0 0.5
你改成这个试试呢？

佩佩 · 发表于 2012-4-23 20:51:34

mofangbao 发表于 2012-4-23 19:24

登录/注册后可看大图

你这个ydef乱写的吧？描述文件y方向最多只能有180度，但是格点个数是不限制的，你这个分辨率写的有问题！ ...

我改了一下文件，ydef为76个点，还是不行啊

佩佩 · 发表于 2012-4-23 21:00:05

mofangbao 发表于 2012-4-23 19:24

登录/注册后可看大图

你这个ydef乱写的吧？描述文件y方向最多只能有180度，但是格点个数是不限制的，你这个分辨率写的有问题！ ...

哦，刚想明白你的意思，试试，谢啦

佩佩 · 发表于 2012-4-23 21:24:39

随缘发表于 2012-4-23 20:04

登录/注册后可看大图

XDEF 151 LINEAR 0 1
YDEF 151 LINEAR 0 0.5
你改成这个试试呢？

仍然不行啊，为什么会有这个错误呢？

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

[分享资料] grads格点资料画图（初学）

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

相关帖子