爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 308|回复: 0

[图形美化] 求助GrADs格点绘图时次问题

[复制链接]
回帖奖励 1 金钱 回复本帖可获得 1 金钱奖励! 每人限 1 次(中奖概率 30%)

新浪微博达人勋

发表于 2024-6-9 14:07:32 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 时光经线 于 2024-6-10 14:45 编辑

请问为什么画出来的两个时次的图都是一样的?

                               
登录/注册后可看大图


at.jan.ctl文件:
dset D:\college\Fortran\daimh\p15\p8\at.jan.grd
dtype  station
stnmap  D:\college\Fortran\daimh\p15\p8\t1601.map
undef  -2.56E33
title 160 station temperature
tdef   61 linear  jan1951 1yr
vars  1
t 0 99 temperature
endvars


                               
登录/注册后可看大图


格点文件:
dset D:\college\Fortran\daimh\p15\p8\Grid.grd
undef -999.0
title Sample GRID Data
xdef 71 linear 70 1
ydef 41 linear 15 1
zdef 1 linear 500 1
tdef 1 linear jan1951 1yr
vars  1
g  0  99  grid data prepared for oacres function
endvars


                               
登录/注册后可看大图


gs文件:
'reinit'
*打开格点数据描述文件grid161.ctl
'open D:\college\Fortran\daimh\p15\p8\Grid161.ctl'
*打开站点数据描述文件
'open D:\college\Fortran\daimh\p15\p8\at.jan.ctl'
*调用颜色定义gs
*'define_colors'
'define_colors.gs'
*设置经纬度范围
'set lat 15 55'
'set lon 70 135'
*设置时间范围
'set t 1'
*插值
'define a1=oacres(g(t=1),t.2)'
'define b1=maskout(a1,g(t=1)-0.5)'
*绘制距平图
'enable print D:\college\Fortran\daimh\p15\p8\atem.gmf'
'set vpage 0 11 0 8.5'
'set parea 1 10 1 8'
'set grads off'
'set grid off'
'set mpdset cnhimap'
*设置时次为1973年一月
'set t 23'
*设置出图类型为填色图
'set gxout shaded'
*设置特殊等值线及颜色
'set clevs -4 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0 0.5 1.0 1.5 2.0 2.5 3 3.5 4'
'set ccols 49 48 47 46 45 44 43 42 41 61 62 63 64 65 66 67 68 69'
*显示插值后的距平b1
'd b1'
'cbarn'
*设置出图类型为等值线图
'set gxout contour'
*设置等值线间隔为0.5
'set cint 0.5'
*显示插值后的距平b1
'd b1'
*写标题
'draw title Jan 1973'
'run D:\college\Fortran\daimh\p15\p8\test.gs'
'run D:\college\Fortran\daimh\p15\p8\southsea.gs'
'print'
'c'
'set vpage 0 11 0 8.5'
'set parea 1 10 1 8'
'set grads off'
'set grid off'
*设置时次为2008年一月
'set t 58'
*设置出图类型为填色图
'set mpdset cnhimap'
'set gxout shaded'
*设置特殊等值线及颜色
'set clevs -4 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0 0.5 1.0 1.5 2.0 2.5 3 3.5 4'
'set ccols 49 48 47 46 45 44 43 42 41 61 62 63 64 65 66 67 68 69'
*显示插值后的距平b1
'd b1'
'cbarn'
*设置出图类型为等值线图
'set gxout contour'
*设置等值线间隔为0.5
'set cint 0.5'
*显示插值后的距平b1
'd b1'
*写标题
'draw title Jan 2008'
'run D:\college\Fortran\daimh\p15\p8\test.gs'
'run D:\college\Fortran\daimh\p15\p8\southsea.gs'
'print'
'c'
'disable print'
;


                               
登录/注册后可看大图


图:
QQ图片20240609140839.png
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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