多时次站点资料作图

猫猫monty

我的fortran转换grads可识别grd文件如下：
Program sta2grd
Implicit none
!这里是程序的变量声明
Character*8 stid(2307)
Real lon(2307),lat(2307),var,tem(90,2307),tim
integer nlev,flag,i,j,k
Open(1,file='E:\data\ncep_winter_ok\std_2016\std_ok\ano\ano_tem.txt',status='old')
do k=1,2307
read(1,*) stid(k),lon(k),lat(k),(tem(i,k),i=1,90)
print*,stid(1),stid(5),stid(10)
enddo
close(1)
open(2,file='E:\data\ncep_winter_ok\std_2016\std_ok\ano\ano_ok.tem.grd',status='replace',form='binary')
Do i=1,90
tim=0.0
nlev=1
flag=1
do j=1,2307
Write(2)stid(j),lat(j),lon(j),tim,nlev,flag,tem(i,j)
enddo
print*,stid(2),lat(2),lon(2),tim,nlev,flag,tem(i,2)
nlev=0
Write(2)stid(j),lat(j),lon(j),tim,nlev,flag
Enddo
Close(2)
END
为grd文件配上ctl如下：
*生成map文件！
DSET E:/data/ncep_winter_ok/std_2016/std_ok/ano/ano_ok.tem.grd
DTYPE station
STNMAP E:/data/ncep_winter_ok/std_2016/std_ok/ano/ano_ok.tem.map
UNDEF -32766
TITLE TEM STATION Data Sample
*这里的时间可以根据上面示例文件中的第一行的时间填写
TDEF 90 linear 01Dec2015 1dy
VARS 1
tem 0 99 tem2016
ENDVARS
但是在grads里面进行简单设置：
如图，但是就是没有站点信息，检查了好久都没有看出毛病，求大神指导！！！

男紫汗

ctl缺乏经纬度

猫猫monty

男紫汗 发表于 2016-10-24 07:58
ctl缺乏经纬度

站点资料的ctl不需要经纬度啊，经纬度在grid的ctl里面

云中月影

男紫汗

猫猫monty 发表于 2016-10-24 08:04
站点资料的ctl不需要经纬度啊，经纬度在grid的ctl里面

write(2)id(j-1),lat(j-1),lon(j-1),tim,nlev,nflag

猫猫monty

男紫汗 发表于 2016-10-24 08:52
write(2)id(j-1),lat(j-1),lon(j-1),tim,nlev,nflag

还是不对。。。

猫猫monty

DSET E:\data\ncep_winter_ok\std_2016\std_ok\ano\grid0.25.grd
TITLE Grid Data Sample
UNDEF -9.99E33
XDEF 101 LINEAR 100 0.25
YDEF 61 LINEAR 20.0 0.25
ZDEF 1 LEVELS 1000
TDEF 1 LINEAR 01jan1980 1mo
VARS 1
g 0 99 Grid Data
ENDVARS


这个是网格的ctl，求大神指导！

一水天下

http://bbs.06climate.com/forum.php?mod=viewthread&tid=27685

一水天下

一水天下 发表于 2016-10-24 09:26
http://bbs.06climate.com/forum.php?mod=viewthread&tid=27685

楼主试试这个呢

男紫汗

两个ctl的时间要一样