请问我的ctl文件哪里出了问题

hgsndg

问题描述：用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.

hgsndg

看到错误了，但是我运行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里就可以了吗

尽头的尽头

你只提取了一个时次，却告诉他有365个时次

river

zdef 1 linearclose 0 1
tdef 365 linear 00Z01JAN1987 1440m

标红的地方自己再看看

hgsndg

river 发表于 2019-4-22 23:16
zdef 1 linearclose 0 1
tdef 365 linear 00Z01JAN1987 1440m

哈哈我看到了，谢谢

hgsndg

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会显示触发一个断点，这是怎么回事你

river

hgsndg 发表于 2019-4-23 00:19
但是我又有新问题了，这是原nc文件的ctl文件：
set E:\aa\olr.day.mean.nc
title HIRS outgoing longwa ...

你自己的循环难道不是有问题？ i j 哪个是内循环，哪个是外循环，你自己看看