提取GrADS插值结果（oacres/gr2stn/grads2ascii）

mofangbao

burf 发表于 2013-5-3 10:30

                               
登录/注册后可看大图

对不起啊楼主，仔细看了一遍，发现还是不能理解站点背景数据是什么？我想要的就是把格点差值到站点啊，站点 ...

看我楼上给你的回复

mofangbao

传说中的谁 发表于 2013-5-3 10:36

                               
登录/注册后可看大图

就是站点文件，哪个经纬度有站点，然后随便给那个站点赋值（不能没有，但用不着）

你竟然不点回复给他。。。做好事要留名的哈

burf

传说中的谁 发表于 2013-5-3 10:36

                               
登录/注册后可看大图

就是站点文件，哪个经纬度有站点，然后随便给那个站点赋值（不能没有，但用不着）

恩 谢谢！但是刚刚实践了下，给每个站点赋值为1  结果用fnl的格点降水数据插值出来得到的txt里不是1就是0，这个明显不对啊。。。

burf

mofangbao 发表于 2013-5-3 10:43

                               
登录/注册后可看大图

看我楼上给你的回复

囧 多谢提醒。。。 但是就是我回复的那样，如果我给每个站点赋值为1  最后结果不是1就是0啊。。。

burf

mofangbao 发表于 2013-5-3 10:43

                               
登录/注册后可看大图

看我楼上给你的回复

Fortran代码是这样的。。。
program text10
   implicit none
   real r10(160)
   integer i,k,j
   r10=1
   call stntogrd(r10)
   


   end

   subroutine stntogrd(x)
   real lat(160),lon(160),x(160)
   character*8 stid(160)
   integer i,j,k
   !从文件‘china.dat’中读入160站的经纬度
   open(2,file='d:\rain\result\china.dat')
   do k=1,160
      read(2,*)lat(k),lon(k)
!     close(2)
   end do


   !给出各站的站号
   do i=1,160
      stid(i)=char(i)

    end do
    !按顺序向二进制文件中写入记录
        open(3,file='d:\rain\gr2st.grd',form='binary')      
         
               tim=0.0
               nlev=1
               nflag=1
       
              do i=1,160
               write(3)stid(i),lat(i),lon(i),tim,nlev,nflag,x(i)
              end do
       
       
           nlev=0
               write(3)stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag
           
      
         return
        end

burf

传说中的谁 发表于 2013-5-3 10:36

                               
登录/注册后可看大图

就是站点文件，哪个经纬度有站点，然后随便给那个站点赋值（不能没有，但用不着）

Fortran代码是这样的
program text10
   implicit none
   real r10(160)
   integer i,k,j
   r10=1
   call stntogrd(r10)
   


   end

   subroutine stntogrd(x)
   real lat(160),lon(160),x(160)
   character*8 stid(160)
   integer i,j,k
   !从文件‘china.dat’中读入160站的经纬度
   open(2,file='d:\rain\result\china.dat')
   do k=1,160
      read(2,*)lat(k),lon(k)
!     close(2)
   end do


   !给出各站的站号
   do i=1,160
      stid(i)=char(i)

    end do
    !按顺序向二进制文件中写入记录
        open(3,file='d:\rain\gr2st.grd',form='binary')      
         
               tim=0.0
               nlev=1
               nflag=1
       
              do i=1,160
               write(3)stid(i),lat(i),lon(i),tim,nlev,nflag,x(i)
              end do
       
       
           nlev=0
               write(3)stid(i-1),lat(i-1),lon(i-1),tim,nlev,nflag
           
      
         return
        end

mofangbao

burf 发表于 2013-5-3 10:51

                               
登录/注册后可看大图

Fortran代码是这样的
program text10
   implicit none

不要给别人贴这么长的程序，没空帮你看！
请自己检查问题

burf

mofangbao 发表于 2013-5-3 10:54

                               
登录/注册后可看大图

不要给别人贴这么长的程序，没空帮你看！
请自己检查问题

==好。。。明白了。。。非常抱歉。。

burf

传说中的谁 发表于 2013-5-3 10:36

                               
登录/注册后可看大图

就是站点文件，哪个经纬度有站点，然后随便给那个站点赋值（不能没有，但用不着）

发现问题所在了，虽然还没有解决，但应该是数据方面的问题，还是谢谢了。。。

2010fuqiang

好好好~~分享了~~