Changeset 97 for trunk/src/paper01/fig10/fig10_swr_correction.pro
- Timestamp:
- 08/30/11 14:38:29 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/paper01/fig10/fig10_swr_correction.pro
r94 r97 83 83 reinitplt, /z,/invert 84 84 key_portrait = 1 85 85 ; 86 86 openps, FILENAME = project_od_env+'fig10_swr_correction.ps' 87 87 close,/all 88 88 ; 89 89 fi=project_id_env+'swr_isccp_correction.txt' 90 90 res=read_ascii(fi,data_start=0) … … 94 94 std=reform(ff(2,*)) 95 95 mean=reform(ff(4,*)) 96 97 96 ; 97 ; 98 98 splot, mean, bias, small=[2,3,1], xtitle='Mean ERA-I SWR', ytitle='ERA-I SWR bias', $ 99 99 title='', psym=2, charsize=1.2, $ 100 100 xrange=[170,290], yrange=[-10,35], xmin=1, ymin=1 101 101 ; 102 102 ab=linfit(mean, bias,yfit=yfit) 103 a=float(ab(0)) 103 a=float(ab(0)) 104 104 b=float(ab(1)) 105 105 oplot, mean, yfit, thick=3 106 106 print, a,b 107 ; ;significant test107 ; significant test 108 108 significance_test_99,mean,bias 109 109 ; 110 110 fi=project_id_env+'swr_isccp_correction.txt' 111 111 res=read_ascii(fi,data_start=0) … … 115 115 std=reform(ff(2,*)) 116 116 mean=reform(ff(4,*)) 117 117 ; 118 118 splot, mean, std, small=[2,3,3],/noer, xtitle='Mean ERA-I SWR', ytitle='ERA-I std ratio', $ 119 119 title='', psym=2, charsize=1.2, $ 120 120 xrange=[170,290], yrange=[0.4,1.4], xmin=1, ymin=1 121 121 ; 122 122 ab=linfit(mean, std,yfit=yfit) 123 a=float(ab(0)) 123 a=float(ab(0)) 124 124 b=float(ab(1)) 125 125 oplot, mean, yfit, thick=3 126 127 ; ;significant test126 ; 127 ; significant test 128 128 significance_test_99,mean,std 129 129 ; 130 130 closeps 131 132 131 ; 133 132 end
Note: See TracChangeset
for help on using the changeset viewer.