关于站点数据多时次处理和插值保存问题

Sylvia

lulu 发表于 2013-11-19 22:18

                               
登录/注册后可看大图

那请问这种格式的台站经纬度应该怎么读取啊？求指教

我也在处理这个数据，应该是跟160站一样的，但是麻烦点吧

tmw

正在学习处理站点数据，正在收罗各种资料，谢谢楼主提供

wpp851118

你好，我按照你给的程序生成格点数据，但是grads打不开其对应的ctl呢，老是报错怎么回事呢

Sylvia

wpp851118 发表于 2013-12-17 15:42

                               
登录/注册后可看大图

你好，我按照你给的程序生成格点数据，但是grads打不开其对应的ctl呢，老是报错怎么回事呢

你看下是不是你设置的时次不是同一个时次，最好都按照一个时次来弄，如果实在不行的话，改天你发给我，我帮你看看，但这几天估计没空，我得回去找工作咯，呵呵

wpp851118

Sylvia 发表于 2013-12-18 09:42

                               
登录/注册后可看大图

你看下是不是你设置的时次不是同一个时次，最好都按照一个时次来弄，如果实在不行的话，改天你发给我，我 ...

谢谢啦，我知道哪里错了呵呵，是没把数据弄成对应的数组，解决了呵呵

chuanhao

program ok3
implicit none
integer i,j,k,n,l,x,p,t,y,a,b,c,d
parameter(n=39,l=65,k=50)
integer pps(k,n),vpps(n,k),yu(61),stpps,s(161),nflag,nlev
character(8) sta(n),aaa(n)
real lon(n),lat(n),tim


open(13,file='f:\yuzhifenbu\gzpc\filename.txt')
open(14,file='F:\xin\svd\39zhanjingwei.txt')
open(15,file='f:\xin\svd\vji.grd',form='binary')

do a=1,n
read(13,*) AAA(a)
enddo
close(13)
do i=1,n
open(100+i,file='F:\yuzhifenbu\gzpc\'//trim(aaa(i)),form='formatted')
do j=1,50
read(100+i,*) pps(j,i)
enddo
enddo
print*,'ok'

do b=1,n
read(14,*) sta(b),lon(b),lat(b)
enddo
do y=1,50
do x=1,39
vpps(x,y)=pps(y,x)
enddo
enddo
print*,'ok'
!print*,sta,lon,lat
do y=1,50
do x=1,39
tim=0.0
nflag=1
nlev=1


write(15) sta(x),lat(x),lon(x),tim,nlev,nflag,vpps(x,y)
print*,sta(x),lat(x),lon(x),tim,nlev,nflag,vpps(x,y)
enddo

tim=0.0
nflag=1
nlev=0

write(15) sta(x-1),lat(x-1),lon(x-1),tim,nlev,nflag

enddo

end

dset    f:\xin\svd\vji.grd
dtype   station
stnmap   f:\xin\svd\rain.map
undef   -999.0
title   rain
tdef    50  linear jul1951 1yr
vars    1
p  0  99  rainfall  data
endvars

program  grid
parameter(nx=41,ny=31,nt=50)
real lat(ny),lon(nx)
real s(nx,ny,nt)
open(1,file='F:\xin\svd\kongbai\grid.grd',form='binary')
lat(1)=20.0
lon(1)=80.0
do j=1,ny-1
   lat(j+1)=lat(j)+1.0
enddo
do i=1,nx-1
   lon(i+1)=lon(i)+1.0
enddo
do k=1,50
  do i=1,nx
    do j=1,ny
      s(i,j,k)=1
enddo
  enddo
enddo
write(1) s
end

dset F:\xin\svd\kongbai\grid.grd
undef -9.99e33
title Sample Grid Data
xdef 41 linear 80 1
ydef 31 linear 20 1
zdef 1 linear 1000 1
tdef 50 linear jul1951 1yr
vars 1
g 0 99 grid data
endvars

'REINIT'
'open F:\xin\svd\kongbai\grid.ctl'
'open F:\xin\svd\vji.ctl'
'enable print F:\xin\svd\jiangshuifenbu.gmf'
i=1
while(i<=50)
'set t 'i''

'set xlopts 1 5 0.2'
'set ylopts 1 5 0.2'
'set grads off'


'SET YLINT 2'
'set lon 90 112'
'set lat 25 40'
'set mpdset hires cnworld'
'define aa=oacres(g,p.2,1.5)'
'define aa1=maskout(aa,g-0.5)'

'set grid off'
'set gxout contour'
'set clab forced'
*'set cint 3.0'
*'set cmin 28'

'd aa1'

为什么报错呢？

chuanhao

麻烦楼主帮我看看。

chuanhao

画出的图是这样的

Sylvia

chuanhao 发表于 2013-12-26 20:43

                               
登录/注册后可看大图

画出的图是这样的

我看你这个图，应该是数据没有插上去，只是画了站点。你报了什么错，也没有贴出来。

chuanhao

就是这个错误