- 积分
- 2796
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-5-10
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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文件哪里出错了?
谢谢~
|
|