Fortran进行t检验后使用GrADS画打点图

前段时间写本科论文，使用到了对于海温sst差值场的显著性t检验，但苦于不会用GrADS画打点图，网上也找不到画打点的.gs代码，5月1放假自己在楼里硬憋出来一段代码，供大家参考，实测可以运行，输入数据为差值场和t检验结果，其核心是循环判断并且画点“draw sting”。

'reinit'
'open D:\lunwen\nczhuangrd\grdfenlei\sst_jianyan\chazhi12.ctl'
'open D:\lunwen\nczhuangrd\grdfenlei\sst_jianyan\t12.ctl''set grid off'
'set grads off'
'set cthick 6'
* These are the BLUE shades蓝色
'set rgb  16    0    0  255'
'set rgb  17   55   55  255'
'set rgb  18  110  110  255'
'set rgb  19  165  165  255'
'set rgb  20  220  220  255'
* These are the RED shades红色
'set rgb  21  255  220  220'
'set rgb  22  255  165  165'
'set rgb  23  255  110  110'
'set rgb  24  255   55   55'
'set rgb  25  255    0    0'
*阴影'set xlint 10'
'set xlopts 1 4 0.15'
'set ylint 5'
'set ylopts 1 4 0.15''set x 105 255'
'set y 90 145'
*'set gxout shaded'
*'set ccols 14'
*'set clevs -2.654 -1.997 1.997 2.654'
*'set ccols 18 20 0 23 25'
*'d sst.2''set gxout shaded'
'set clevs -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4'
'set ccols 16 17 18 19 20 21 22 23 24 25'
'd sst.1'
'cbar 1.0 0'y2=90
while (y2<145)
x2=105
while(x2<255)
'set x 'x2
'set y 'y2
'd sst.2'
aaa=subwrd(result,4)
if(aaa=-9.99e+08)
aaa=0
endif
if(aaa<0)
aaa=-aaa
endif
*if(aaa>2.654)
if(aaa>1.997)
'q w2xy 'x2-0.5' ' y2-90.5''
abx=subwrd(result,3)
aby=subwrd(result,6)
'set string 1 c 9'
'set strsiz 0.01'
'draw string 'abx' 'aby' `32'
endif
x2=x2+1
endwhile
y2=y2+1
endwhile'draw title difference field and t test type2-type1'
'printim D:\lunwen\nczhuangrd\grdfenlei\sst_jianyan\chazhi_and_t12.png white'
*'c'
;

图片如图，希望有所帮助。

Fortran进行t检验后使用GrADS画打点图相关推荐

Grads画等值线（一）-----心得感言
之前研究了如何用surfer画等值线,总体来说还算不错,在实际项目中得到了应用.但是在surfer等软件中使用的插值法都是纯数学差值法,在孤立点附近容易产生一些不合理区域.后来知道Grads中有在气象 ...
用grads读取txt中经纬度画站点图
用grads读取txt中经纬度,用于画站点图 txt中的站点信息格式如下: stn_num lon lat data1 data2 .... (每个变量用空格格开一行为一个记录) ********** ...
54部经典国外动画片：送给70－80年代的朋友 90年代后的一边画圈圈去41-50
导读: 动漫音乐区 TS服务另类图片汽车俱乐部 54部经典国外动画片:送给70-80年代的朋友 90年代后的一边画圈圈去41-50 41.日本动画片<天空战记> [原片名]天空戦記 ...
三个变量存在一个协整方程_三变量误差修正模型：两个解释变量，用JJ协整检验后得到以下结果，如何看这个结果？怎样做误差修正模型？...
三变量误差修正模型:两个解释变量,用JJ协整检验后得到以下结果,如何看这个结果?怎样做误差修正模型?0 vldxhg2013.04.12浏览2007次分享举报 Date: 01/13/12 Time: ...
matlab如何画极零图,用MATLAB画零极点图.ppt
用MATLAB画零极点图用MATLAB将传递函数化为零极点增益模型并绘制零极点图将传递函数化为零极点增益模型并绘制零极点图 >> num=[3 2 8]; >> den=[ ...
origin中文版散点图拟合曲线_origin中两组数已经画出散点图，如何新增类似于excel中的趋势线？而且是幂函式的，求高手指点...
origin中两组数已经画出散点图,如何新增类似于excel中的趋势线?而且是幂函式的,求高手指点以下文字资料是由(历史新知网www.lishixinzhi.com)小编为大家搜集整理后发布的内容,让 ...
用Grads画等值线（二）------.map生成
要生成.map文件,需要以下四步: (一)准备一个.ch文件,其中有需要的站点数据: (二)编写一个简单的程序把.ch中的数据写成二进制格式的.dat文件.(三)编写一个站点控制文件.ctl.(四)运 ...
grads中画站点图
1.把不规则的数据按照其经纬读写成stream型data(其fortran程序如下) parameter(n=160) character*8 stid(n),sti(n) real rlat(n), ...
斗鱼html5插件音画不同步,直播开发完成后关于音画不同步的原因及办法
如果说直播过程中出现卡顿.延时是很正常的,那么唯有音画不同步是最令人难以忍受的.对于直播开发完成之后出现音画不同步的现象,也是非常令人头疼的.我们可以从以下几个方面来简单分析一下,音画不同步现象产生的 ...

Fortran进行t检验后使用GrADS画打点图

Fortran进行t检验后使用GrADS画打点图相关推荐

最新文章

热门文章