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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3335|回复: 1

[分享资料] 求助!多层次站点资料的读取

[复制链接]
回帖奖励 1 金钱 回复本帖可获得 1 金钱奖励! 每人限 2 次(中奖概率 50%)

新浪微博达人勋

发表于 2013-1-18 17:07:57 | 显示全部楼层 |阅读模式

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

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

x
看了一些关于站点资料的读取的帖子,但是多个层次的读取方法基本没有涉及到,自己试着按照家园整理的《站点作图详解》写了程序,但是仍然有问题,求指点!!!
station资料的读取:
program station
      parameter (mt=427,nz=4)    ! mt 站点数量 ,nz层次
      character stid(mt)*8
      integer flag,nlev
      real lat(mt),lon(mt),tim,high(mt),number(mt),lev(nz)
     real height(mt,nz),temp(mt,nz),ttd(mt,nz),ud(mt,nz),us(mt,nz)


      character stationfile(4)*8
   data stationfile/'850','700','500','200'/
  
   do k=1,nz

   open(k,file=''//trim(stationfile(k))//'_20.000')
   read(k,*)
   read(k,*)
   do i=1,mt
   
      read(k,*)stid(i),lon(i),lat(i),high(i),number(i),height(i,k),temp(i,k),ttd(i,k),ud(i,k),us(i,k)

   end do   

end do  !结束k的循环

      tim=0.0
      flag=1
      nlev=4   
   open(40,file='20.grd',form='binary')  

   do i=1,mt
  write(40)stid(i),lat(i),lon(i),tim,nlev,flag,high(i)
  
  do k=1,nz
  lev(k)=k
  write(40)lev(k),height(i,k),temp(i,k)
  write(*,*)lev(k),height(i,k)  
  end do

end do
  
nlev=0
    write(40)stid(mt-1),lat(mt-1),lon(mt-1),tim,nlev,flag

close(40)
end
显示运行无误
ctl描述文件如下:
dset f:\20.grd
dtype station
stnmap f:\20.map
undef 9999.0
title 500hpa 20

zdef 4 levels 850 700 500 200
tdef 1  linear  29jul2012 1DY
vars 3

high   0    99   high
h      4    99      height
t      4    99      temp
endvars

随后在grads中输入!stnmap -i f:/20.ctl,运行出map文件;
最后用grads画出站点资料,只有high显示有数据,而变量h和t显示空白。
请问问题出在哪?请大家指点下,谢谢~~



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

新浪微博达人勋

发表于 2013-11-7 17:16:02 | 显示全部楼层
write 语句没有将3个变量都输出,只输出了2个变量
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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