插值警告

wuyejincao

插值函数gs运行后提示如下，图可以出，但因为有如下警告，所以我担心图出的不对。请求指教警告是怎么回事？？？

做个霸气的木头

红色警告没事，是因为你有一句数据气候态的命令，modify那个，是正常提示。
红色上面的警告得看你的gs和ctl大家才能帮你

wuyejincao

做个霸气的木头 发表于 2012-10-10 22:35

                               
登录/注册后可看大图

红色警告没事，是因为你有一句数据气候态的命令，modify那个，是正常提示。
红色上面的警告得看你的gs和ct ...

226sm.grd (1.33 MB, 下载次数: 14)

2012-10-11 10:18 上传

点击文件名下载附件 grid.grd (2.13 MB, 下载次数: 12)

2012-10-11 10:18 上传

点击文件名下载附件 226sm.ctl (155 Bytes, 下载次数: 5)

2012-10-11 10:18 上传

点击文件名下载附件 grid.ctl (194 Bytes, 下载次数: 7)

2012-10-11 10:18 上传

点击文件名下载附件 obs.gs (5.07 KB, 下载次数: 12)

2012-10-11 10:17 上传

点击文件名下载附件
求鉴定~

做个霸气的木头

你单独画站点数据会发现，一个数都木有，说明是你生成226sm.grd的时候就出错了，你好好检查你的fortran吧~

wuyejincao

做个霸气的木头 发表于 2012-10-11 12:31

                               
登录/注册后可看大图

你单独画站点数据会发现，一个数都木有，说明是你生成226sm.grd的时候就出错了，你好好检查你的fortran吧~

我是我用22sm.ctl画的第一时次的站点图，你看应该没问题啊~

river

wuyejincao 发表于 2012-10-11 10:19

                               
登录/注册后可看大图

求鉴定~

ctl和gs贴代码的方式拿出来呗，我们又不是版主，没有那么多金钱下下来看

wuyejincao

226sm.ctl
dset D:\model\OBS\226sm.grd
dtype station
stnmap D:\model\OBS\sm.map
undef -9999.0
title sm
tdef 192 linear jan1993 1mo
vars 1
sm 0 99 sm
endvars


grid.ctl
dset D:\model\OBS\grid.grd
undef -9.99e33
title Sample Grid Data
xdef 71 linear 70 1
ydef 41 linear 15 1
zdef 1 linear 1000 1
tdef 192 linear JAN1993 1mo
vars 1
g 0 99 grid data
endvars


obs.gs
'reinit'
'open D:\model\OBS\grid.ctl'
'open D:\model\OBS\226sm.ctl'
'set mpdset hires'
'set lat 15 55'
'set lon 70 140'
'set t 1 192'
'define a=oacres(g,sm.2,1.5)'
*'define a1=maskout(a,g-0.5)'
'define aa=smth9(a)'
'set lat 15 55'
'set lon 70 140'
*'set clab forced'
*'set ccolor rainbow'
'set grads off'
'set grid off'

******开始画图*****
'D:\source\grads\define_colors.gs'
'enable print D:\model\OBS\OBS-climate-10cm(9308)test.gmf'
'set lat 15 55'
'set lon 70 140'
'set lev 1'

****标准差***
'define mean=ave(aa*1.4,t=1,t=192)'
'define sd=sqrt((sum(pow(aa*1.4-mean,2),t=1,t=192))/191)'
'set vpage 0 11 0 8.5';
'set parea 1.0 10 2 7';
'set xlopts 1 6 0.18';
'set ylopts 1 6 0.18'
'set grads off';'set grid off'
'set lon 70 140';'set lat 15 55';'set t 1'
'set xlint 10 '
'set ylint 5';
'set mpdset cnmap1'

'set gxout shaded'
'set cint 0.02'
'set clevs  0.04 0.06 0.08 0.1'
'set ccols 31 33 35 37 39 '
'd smth9(sd)'
'cbarn 0.85 0 5.50 1.05 0.65'
*'D:\source\grads\cbar_matlab.gs 1 1 0'';*****#调用色标脚本，1 1 1分别表示：隔1个颜色显示一次标注和竖线；色标尺寸；色标水平还是竖直****
'run D:\pro-software\OpenGrADS\cnbasemap.gs'

'set gxout contour'
'set cint 0.02'
'd smth9(sd)'
'run D:\pro-software\OpenGrADS\cnbasemap.gs'
'draw title Standard Deviation of OBS(10cm)'
'run D:\pro-software\OpenGrADS\scs_map.gs'
'print'
'c'

*****气候态******
'set vpage 0 11 0 8.5';
'set parea 1.0 10 2 7';
'set xlopts 1 6 0.18';
'set ylopts 1 6 0.18'
'set grads off';'set grid off'
'set lon 70 140';'set lat 15 55'
'set xlint 10 '
'set ylint 5';
*'set mpdset cnriver'
'set mpdset cnmap1'

'set gxout shaded'
'set cint 0.4';'set cmin 0'
'set clevs 2.8 3.2 3.6 4.0 4.4'
'set ccols 21 22 23 24 25 26'
'd smth9(ave(aa*1.4,t=1,t=192)*10)'
'cbarn 0.85 0 5.50 1.05 0.65';
'run D:\pro-software\OpenGrADS\cnbasemap.gs'

'set gxout contour'
'set cint 0.4'
'set cmin 0'
'd smth9(ave(aa*1.4,t=1,t=192)*10)'
'run D:\pro-software\OpenGrADS\cnbasemap.gs'
'draw title Monthly mean of OBS(10cm)'
'run D:\pro-software\OpenGrADS\scs_map.gs'
'print'
'c'
**************季节平均****
'set t 1 12'
'define avemon=ave(aa*1.4,t+0,t=192,12)'
'modify avemon seasonal'
'set t 1'
'define s1=(avemon(t=3)+avemon(t=4)+avemon(t=5))/3.0 '
'define s2=(avemon(t=6)+avemon(t=7)+avemon(t=8))/3.0 '
'define s3=(avemon(t=9)+avemon(t=10)+avemon(t=11))/3.0 '
'define s4=(avemon(t=1)+avemon(t=2)+avemon(t=12))/3.0 '
***********
sub.1='OBS Spring(10cm)';
sub.2='OBS Summer(10cm)';
sub.3='OBS Autumn(10cm)';
sub.4='OBS Winter(10cm)';
****  pic 1******
'set vpage 0 11 0 8.5';
'set parea 1.0 10 2 7';
'set xlopts 1 6 0.18';
'set ylopts 1 6 0.18'

'set grads off';'set grid off'
'set lon 70 140';'set lat 15 55';'set t 1'
'set xlint 10 '
'set ylint 5';
'set mpdset cnmap1'
'set gxout shaded'
'set cint 0.6';'set cmin 0'
'set clevs 2.4 3 3.6 4.2'
'set ccols 21 22 23 24 25'
'd smth9(s1*10)'
'cbarn 0.85 0 5.50 1.05 0.65';
'run D:\pro-software\OpenGrADS\cnbasemap.gs'

'set gxout contour'
'set cint 0.6';'set cmin 0'
'd smth9(s1*10)'
'run D:\pro-software\OpenGrADS\cnbasemap.gs'
'draw title 'sub.1
'run D:\pro-software\OpenGrADS\scs_map.gs'
'print'
'c'
****  pic 2******
'set vpage 0 11 0 8.5';
'set parea 1.0 10 2 7';
'set xlopts 1 6 0.18';
'set ylopts 1 6 0.18'
********
'set grads off';'set grid off'
'set lon 70 140';'set lat 15 55';'set t 1'
'set xlint 10';'set ylint 5';
'set mpdset cnmap1'
'set gxout shaded'
'set cint 0.6';'set cmin 0'
'set clevs 1.8 2.4 3 3.6 4.2'
'set ccols 21 22 23 24 25 26'
'd smth9(s2*10)'
'cbarn 0.85 0 5.50 1.05 0.65';
'run D:\pro-software\OpenGrADS\cnbasemap.gs'

'set gxout contour'
'set cint 0.6';'set cmin 0'
'd smth9(s2*10)'
'run D:\pro-software\OpenGrADS\cnbasemap.gs'
'draw title 'sub.2
'run D:\pro-software\OpenGrADS\scs_map.gs'
'print'
'c'
****  pic 3******
'set vpage 0 11 0 8.5';
'set parea 1.0 10 2 7';
'set xlopts 1 6 0.18';
'set ylopts 1 6 0.18'
********
'set grads off';'set grid off'
'set lon 70 140';'set lat 15 55';'set t 1'
'set xlint 10';'set ylint 5';
'set mpdset cnmap1'
'set gxout shaded'
'set cint 0.6';'set cmin 0'
'set clevs 2.4 3 3.6 4.2'
'set ccols 21 22 23 24 25'
'd smth9(s3*10)'
'cbarn 0.85 0 5.50 1.05 0.65';
'run D:\pro-software\OpenGrADS\cnbasemap.gs'
'set gxout contour'
'set cint 0.6';'set cmin 0'
'd smth9(s3*10)'
'run D:\pro-software\OpenGrADS\cnbasemap.gs'
'draw title 'sub.3
'run D:\pro-software\OpenGrADS\scs_map.gs'
'print'
'c'
****  pic 4******
'set vpage 0 11 0 8.5';
'set parea 1.0 10 2 7';
'set xlopts 1 6 0.18';
'set ylopts 1 6 0.18'
********
'set grads off';'set grid off'
'set lon 70 140';'set lat 15 55';'set t 1'
'set xlint 10';'set ylint 5';
'set mpdset cnmap1'
'set gxout shaded'
'set cint 0.4';'set cmin 0'
'set clevs 3.2 3.6 4 4.4 4.8'
'set ccols 21 22 23 24 25 26'
'd smth9(s4*10)'
'cbarn 0.85 0 5.50 1.05 0.65';
'run D:\pro-software\OpenGrADS\cnbasemap.gs'
'set gxout contour'
'set cint 0.4';'set cmin 0'
'd smth9(s4*10)'
'run D:\pro-software\OpenGrADS\cnbasemap.gs'
'draw title 'sub.4
'run D:\pro-software\OpenGrADS\scs_map.gs'
'print'
'c'
'disable print'
;

做个霸气的木头

wuyejincao 发表于 2012-10-11 14:09

                               
登录/注册后可看大图

我是我用22sm.ctl画的第一时次的站点图，你看应该没问题啊~

你缩到中国区看看，都是零

做个霸气的木头

river 发表于 2012-10-11 15:50

                               
登录/注册后可看大图

ctl和gs贴代码的方式拿出来呗，我们又不是版主，没有那么多金钱下下来看

钱是好东西啊，我也舍不得哇~

wuyejincao

做个霸气的木头 发表于 2012-10-11 16:48

                               
登录/注册后可看大图

钱是好东西啊，我也舍不得哇~

这个是零点几几（小数点后面好几位），因此没显示出来~