爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3299|回复: 5

[求助] 如何用unformatted格式生成格点背景场

[复制链接]

新浪微博达人勋

发表于 2018-7-7 23:58:16 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 欢乐树的朋友们 于 2018-7-8 13:40 编辑

这是用binary格式生成格点背景场的fortran程序:
program gen_grid
    implicit none
    integer :: i,j
    integer,parameter :: nx=144
    integer,parameter :: ny=73
    real :: grid
    grid=1.0
    open(10,file='grid.grd',form='binary')
    do j=1,ny
        do i=1,nx
            write(10)grid
        end do
    end do
    close(10)
end program gen_grid

但是我现在用的fortran不支持binary格式,如果用unformatted格式能生成用于grads做图的背景场吗,具体程序怎么写呢?
试了一天了也没弄出来
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-7-8 11:02:53 | 显示全部楼层
open(10,file='grid.grd',form='unformatted',access='direct',recl=nx*ny*4)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-8 12:45:04 | 显示全部楼层
yuchen 发表于 2018-7-8 11:02
open(10,file='grid.grd',form='unformatted',access='direct',recl=nx*ny*4)

可以了,太谢谢你了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-8 12:45:07 | 显示全部楼层
yuchen 发表于 2018-7-8 11:02
open(10,file='grid.grd',form='unformatted',access='direct',recl=nx*ny*4)

可以了,太谢谢你了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-8 13:11:38 | 显示全部楼层
yuchen 发表于 2018-7-8 11:02
open(10,file='grid.grd',form='unformatted',access='direct',recl=nx*ny*4)

还想请教一下,如果是地面台站的资料生成grd,格式应该是怎样的呢
open(11,file='air0703.grd',status='replace',form='binary')
        tim=0.0
        nlev = 1
        flag = 1
        do i = 1,num1
                write(11)stat(i),lat(i),lon(i),tim,nlev,flag,tmp(i)
        end do
        nlev = 0
        write(11)stat(num1),lat(num1),lon(num1),tim,nlev,flag
        close(11)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-7-8 14:42:11 | 显示全部楼层
不用谢,台站试试这样:open(11,file='air0703.grd',FORM='UNFORMATTED', access='sequential')
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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