爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3466|回复: 0

[脚本编辑] 用FORTRAN存储GrADS数据dat

[复制链接]

新浪微博达人勋

发表于 2016-5-8 10:15:02 | 显示全部楼层 |阅读模式

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

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

x
各位大神,你们好!我把原始的dat数据经过处理后,得到了一个三维的数组PRES(nx,ny,nt)和一个二维数组NET(nz,nt),其中nx是经向方向,有240个格点,ny是纬向方向,有202个格点,nz是垂直方向,有19层,nt是时间,有33个时次。
下面是我用FORTRAN程序把PRES和NET数据写进二进制文件用来GrADS出图

OPEN (11, FILE="G:\F\OUTPUT.DAT", FORM="UNFORMATTED", ACCESS="DIRECT", RECL=NX*NY)
JREC = 0
DO IT = 1, NT
    JREC = JREC + 1
    WRITE (11, REC = JREC) (NET(IZ, IT), IZ = 1, NZ)
    JREC = JREC + 1
    WRITE (11, REC = JREC) ((PRES(IX, IY, IT), IX = 1, NX), IY = 1, NY)
END DO

程序运行没有报错,也生成了dat文件,相应的写了ctl文件

DSET G:\F\OUTPUT.DAT
UNDEF 1.e30
TITLE NET WATER VAPOR FLUX
XDEF 240 LINEAR 92.13150 0.27505
YDEF 202 LINEAR 1.543289 0.274933
ZDEF  19 LEVELS 1000950 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100
TDEF  33 LINEAR 12Z23AUG2012 180MN
VARS 2
NET      19  99   NET WATER VAPOR
PRES      0  99   SURFACE PRESSURE
ENDVARS

问题是用GrADS作图时,有一些时次或层次是有图片生成,可是有一些时次或层次又没有图。
请问我的FORTRAN程序和ctl文件哪里出错了?
谢谢~
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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