- 积分
- 92
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2024-1-22
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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'
;
图:
|
-
|