前段時間寫本科論文，使用到了對于海溫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画打点图的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。