请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3619|回复: 4

[脚本编辑] 【求助】用GrADS作u随时间的变化图

[复制链接]

新浪微博达人勋

发表于 2017-3-29 16:31:52 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 微笑的橙子 于 2017-3-30 07:51 编辑

刚学GrADS不久的小白一枚。。用的方法比较蠢。。先用fortran把数据算出来了然后再写ctl文件进GrADS里面跑,但是不知道为什么出不来图
求各位大神帮看!!!!要求是在300个时间步长的时间段内在每个时间点上,闭合圈上360个格点的速度。。
fortran程序:
program main
implicit none
  real u(360,301)
  real c,dx,dt,pi
  integer i,j,l
  l=1
  pi=3.1415927

  write(*,*) "x steps:"
  read(*,*) dx
  write(*,*) "t steps:"
  read(*,*) dt
  write(*,*) "c:"
  read(*,*) c

  do i=1,360
    u(i,1)=20.0*cos((3.0*(i-1))*pi/180)
  end do

  do j=2,360
    u(j,2)=u(j,1)+c*dt*(u(j-1,1)-u(j+1,1))/(2*dx)
  end do

  u(1,2)=c*dt*(u(360,1)-u(2,1))/(2*dx)+u(360,1)
  u(360,2)=c*dt*(u(359,1)-u(1,1))/(2*dx)+u(1,1)

  do j=3,301
    do i=2,359
      u(i,j)=c*dt*(u(i-1,j-1)-u(i+1,j-1))/dx+u(i,j-2)
    end do
    u(1,j)=2*c*dt*(u(360,j-1)-u(2,j-1))/(2*dx)+u(1,j-2)
    u(360,j)=2*c*dt*(u(359,j-1)-u(1,j-1))/(2*dx)+u(360,j-2)
  end do

  open(unit=10,file='desktop/result.grd',form='unformatted',status='new',access='direct',recl=360*4)
  do j=1,301
    write(10,rec=l) (u(i,j),i=1,360)
    l=l+1
  end do
  close(10)

stop
end



ctl文件:
dset desktop/result.grd
title homework3
undef -9999.0
xdef 360 linear 1. 1.
ydef 1 linear 0. 2.5
zdef 301 linear 1. 1.
tdef 1 linear jan1981 5mn
vars 1
u 1 99 windspeed
endvars


d了u后除了坐标轴什么都没有 求帮助求帮助!!

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-3-29 17:00:57 | 显示全部楼层
随时间变化?t只有一个时次。。你确定垂直方向有301层?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-3-30 07:45:42 | 显示全部楼层
男紫汗 发表于 2017-3-29 17:00
随时间变化?t只有一个时次。。你确定垂直方向有301层?

按照您的建议改了!
但是还是不行诶。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2017-3-30 08:48:02 | 显示全部楼层
闭合圈360个格点什么意思?环绕地球一圈?那你的横轴是时间,纵轴应该是经度了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2017-3-30 09:37:48 | 显示全部楼层
男紫汗 发表于 2017-3-30 08:48
闭合圈360个格点什么意思?环绕地球一圈?那你的横轴是时间,纵轴应该是经度了

解决啦~~
感谢你的帮助!!
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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