请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5481|回复: 6

[图形美化] 请问我的ctl文件哪里出了问题

[复制链接]

新浪微博达人勋

发表于 2019-4-22 21:15:23 | 显示全部楼层 |阅读模式

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

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

x
问题描述:用gs文件提取nc文件,然后为提取出的grd文件写描述文件ctl,这个ctl文件是不是还要拖到grads里open一下?但是不知道哪里出了错?求大神帮忙看看
这是源文件的ctl:
dset E:\aa\uwnd.sig995.1987.nc
title mean daily NMC reanalysis
undef -9.96921e+36
dtype netcdf
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 linear 0 1
tdef 365 linear 00Z01JAN1987 1440mn
vars 1
uwnd=>uwnd  0  t,y,x  mean Daily u-wind
endvars
这是我写的gs:
'reinit'
'sdfopen E:\aa\uwnd.sig995.1987.nc'
'set gxout fwrite'
'set fwrite E:\aa\uwnd1987.grd'
'set x 105 120 '
'set y 5 15'
'set z 1'
'set t 365 '
'd uwnd'
'disable fwrite'
'reinit'
;
这是我为这个文件写的ctl:
dset E:\aa\uwnd1987.grd
title daily ltm uwnd.sig995 from the NCEP Reanalysis
undef -9.96921e+36
xdef 7 linear 105.0 2.5
ydef 5 linear 5.0 2.5
zdef 1 linearclose 0 1
tdef 365 linear 00Z01JAN1987 1440m
vars 1
uwnd=>uwnd  0  t,y,x  mean Daily u-wind
endvars
这是报错:
ga-> open E:\aa\uwnd1987.ctl
Scanning description file:  E:\aa\uwnd1987.ctl
Syntax Error:  Invalid Date/Time offset.
  Expecting yr/mo/dy/hr/mn, found m
Open Error:  Invalid time increment in tdef
use 1 for single time data
  --> The invalid description file record is:
  --> tdef 365 linear 00Z01JAN1987 1440m
  The data file was not opened.

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2019-4-22 21:23:18 | 显示全部楼层
看到错误了,但是我运行fortran的时候显示触发一个断点
program console10
    implicit none
    integer,parameter::nx=7,ny=5,nt=365,nxx=26,nyy=21,ntt=73
    real uwnd1987(nx,ny,nt),vwnd1987(nx,ny,nt),olr1(nxx,nyy,nt),olr2(nt),olr3(ntt),iabn(nt),s1(nt),s2(ntt)
    integer i,j,k,t
    open(1,file='e:\aa\uwnd1987.grd',form='binary')
    do k=1,ntt
        read(1)((uwnd1987(i,j,k),j=1,ny),i=1,nx)
    enddo
    close(1)
    print*,uwnd1987
    end
我想看一下grd文件有没有争取提取。请问为提取出的grd文件所写的ctl文件怎么使用啊,是把ctl拖进grads里就可以了吗
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-4-22 21:55:45 | 显示全部楼层
你只提取了一个时次,却告诉他有365个时次
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-4-22 23:16:25 | 显示全部楼层
zdef 1 linearclose 0 1
tdef 365 linear 00Z01JAN1987 1440m

标红的地方自己再看看
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-4-23 00:00:11 | 显示全部楼层
river 发表于 2019-4-22 23:16
zdef 1 linearclose 0 1
tdef 365 linear 00Z01JAN1987 1440m

哈哈我看到了,谢谢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-4-23 00:19:12 | 显示全部楼层
本帖最后由 hgsndg 于 2019-4-23 00:25 编辑
river 发表于 2019-4-22 23:16
zdef 1 linearclose 0 1
tdef 365 linear 00Z01JAN1987 1440m

但是我又有新问题了,这是原nc文件的ctl文件:
set E:\aa\olr.day.mean.nc
title HIRS outgoing longwave radiation, regridded to 90S to 90N
undef -9.99e+08
dtype netcdf
xdef 360 linear 0.5 1
ydef 181 linear -90 1
zdef 1 linear 0 1
tdef 12419 linear 00Z01JAN1979 1440mn
vars 1
olr=>olr  0  t,y,x  outgoing longwave radiation
endvars
这是我写的gs文件
'reinit'
'sdfopen E:\aa\olr.day.mean.nc'
'set gxout fwrite'
'set fwrite E:\aa\olr1987.grd'
'set x 105.5 120.5'
'set y 5 15'
'set z 1'
i=2922
while(i<=3287)
'set t 'i''
'd olr'
i=i+1
endwhile
'disable fwrite'
'reinit'
;
这是我写的ctl文件
dset E:\aa\olr1987.grd
title  HIRS outgoing longwave radiation, regridded to 90S to 90N
undef -9.99e+08
xdef 16 linear 0.5 1
ydef 11 linear 0 1
zdef 1 linear 0 1
tdef 365 linear 00Z01JAN1987 1440mn
vars 1
olr=>olr  0  t,y,x  outgoing longwave radiation
endvars
这是fortran 文件
program console10
    implicit none
    integer,parameter::nx=7,ny=5,nt=365,nxx=16,nyy=11,ntt=73
    real uwnd1987(nx,ny,nt),vwnd1987(nx,ny,nt),olr1(nxx,nyy,nt),olr2(nt),olr3(ntt),iabn(nt),s1(nt),s2(ntt)
    integer i,j,k,t
    open(1,file='e:\aa\olrl987.grd',form='binary')
    do k=1,nt
        read(1)((olr1(i,j,k),j=1,nyy),i=1,nxx)
    enddo
    close(1)
    print*,olr1
    end
运行时fortran会显示触发一个断点,这是怎么回事你
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-4-25 21:36:49 | 显示全部楼层
hgsndg 发表于 2019-4-23 00:19
但是我又有新问题了,这是原nc文件的ctl文件:
set E:\aa\olr.day.mean.nc
title HIRS outgoing longwa ...

你自己的循环难道不是有问题? i j 哪个是内循环,哪个是外循环,你自己看看
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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