grads高空站点资料画等值线

易小凯 · 发表于 2012-3-10 19:45:47

gs文件
'reinit'
'open D:\grid.ctl'
'open D:\cor.ctl'

'enable print D:\cor.gmf'
i=1
while(i<=16)
'set grads off'
'set lon 70 140'
'set lat 15 55'
'set z 'i
'set mpdset hires'
'define a=oacres(g1,cor1.2)'
'define a1=maskout(a,g1-0.5)'
'define aa=smth9(a1)'
'set gxout contour'
'd aa'
'draw title 'i
'print'
'c'
i=i+1
endwhile

'disable print'
'reinit'
;

mofangbao · 发表于 2012-3-10 23:09:16

易小凯发表于 2012-3-10 19:45

登录/注册后可看大图

gs文件
'reinit'
'open D:\grid.ctl'

fortran里面输出的部分捏

易小凯 · 发表于 2012-3-11 00:37:37

mofangbao 发表于 2012-3-10 23:09

登录/注册后可看大图

fortran里面输出的部分捏

station文件输出
tim=0.0
nlev=LEV
nflag=0
DO i=1,NN
WRITE(21) station(i).id,station(i).lat,station(i).lon,tim,nlev,nflag
DO k=1,LEV
WRITE(21) PRS(k),(station(i).cor(k,j),j=1,6)
ENDDO
ENDDO
nlev=0
WRITE(21) station(i-1).id,station(i-1).lat,station(i-1).lon,tim,nlev,nflag
CLOSE(21)

易小凯 · 发表于 2012-3-11 00:38:11

grid文件输出
g=1.0
OPEN(22,FILE='grid.grd',FORM='BINARY')
DO n=1,6
DO k=1,LEV
WRITE(22) ((g,i=1,NX),j=1,NY)
ENDDO
ENDDO

miracle · 发表于 2012-9-24 22:56:44

能看一下你将高空多层多变量资料处理成二进制文件的程序吗，我的总是不对。

whmooncool · 发表于 2013-1-18 16:38:53

可以麻烦楼主能整理下吗？我也遇到了类似的问题，高层的有问题。能否麻烦楼主吧读取站点资料的fortran程序全部贴出来？谢谢

月光竹影 · 发表于 2013-1-18 16:58:24

程序设计交流的好地方

river · 发表于 2013-1-18 18:25:19

希望楼主可以将完整过程写出来，加上测试数据就更完美了。很多人都可以借鉴了就

20009 · 发表于 2014-7-22 23:56:25

好东西。楼主辛苦了

qiuba · 发表于 2018-5-26 23:00:52

谢谢分享

		自动登录	找回密码
密码			立即注册

[分享资料] grads高空站点资料画等值线

评分