Changeset 97 for trunk/src/paper01/fig4/fig4_met_var_correction_scatter.pro

Timestamp:

08/30/11 14:38:29 (13 years ago)

Author:

pinsard

Message:

suppress blank lines trailing blank

File:

• trunk/src/paper01/fig4/fig4_met_var_correction_scatter.pro (modified) (8 diffs)

trunk/src/paper01/fig4/fig4_met_var_correction_scatter.pro

@@ -84 +84 @@
 reinitplt, /z,/invert
 key_portrait = 1
-
+;
 openps, FILENAME = project_od_env+'fig4_met_var_correction_scatter.ps'
 fi=project_id_env+'sst_correction_final.txt'
@@ -93 +93 @@
 std=reform(ff(2,*))
 mean=reform(ff(4,*))
-
+;
 splot, mean, bias, small=[2,3,1],/noer, xtitle='Mean ERAI SST', ytitle='ERAI SST bias', $
        psym=2, charsize=1.,  $
@@ -100 +100 @@
 xx=total(std)/n_elements(std)
 yy=total(bias)/n_elements(bias)
-
+;
 xyouts, 22, 0.6, 'bias=-0.05', charsize=1.1
 xyouts, 22, 0.5, 'std ratio=0.989', charsize=1.1
-
-print, xx, yy
-ab=linfit(mean, bias,yfit=yfit)
-a=float(ab(0)) 
-b=float(ab(1))
-oplot, mean, yfit, thick=3, color=50
-print, a,b
-
-;; significant test
-significance_test_99,mean,bias
-
+;
+print, xx, yy
+ab=linfit(mean, bias,yfit=yfit)
+a=float(ab(0))
+b=float(ab(1))
+oplot, mean, yfit, thick=3, color=50
+print, a,b
+;
+; significant test
+significance_test_99,mean,bias
+;
 fi=project_id_env+'t2m_correction_final.txt'
 res=read_ascii(fi,data_start=0)
@@ -121 +121 @@
 std=reform(ff(2,*))
 mean=reform(ff(4,*))
-
-xx=total(std)/n_elements(std)
-yy=total(bias)/n_elements(bias)
-print, xx, yy
-
+;
+xx=total(std)/n_elements(std)
+yy=total(bias)/n_elements(bias)
+print, xx, yy
+;
 splot, mean, bias, small=[2,3,2],/noer, xtitle='Mean ERAI T2m', ytitle='ERAI T2m bias', $
        psym=2, charsize=1.,  $
        xrange=[20,31], yrange=[-1,0.6], xmin=1, ymin=1
 oplot, [20,32], [0,0]
-
+;
 xyouts, 22, 0.5, 'bias=-0.296', charsize=1.1
 xyouts, 22, 0.4, 'std ratio=0.92', charsize=1.1
-
-ab=linfit(mean, bias,yfit=yfit)
-a=float(ab(0)) 
-b=float(ab(1))
-oplot, mean, yfit, thick=3, color=50
-print, a,b
-
-;; significant test
-significance_test_99,mean,bias
-
+;
+ab=linfit(mean, bias,yfit=yfit)
+a=float(ab(0))
+b=float(ab(1))
+oplot, mean, yfit, thick=3, color=50
+print, a,b
+;
+; significant test
+significance_test_99,mean,bias
+;
 fi=project_id_env+'q2m_correction_final.txt'
 res=read_ascii(fi,data_start=0)
@@ -150 +150 @@
 std=reform(ff(2,*))
 mean=reform(ff(4,*))
-
-xx=total(std)/n_elements(std)
-yy=total(bias)/n_elements(bias)
-print, xx, yy
-
+;
+xx=total(std)/n_elements(std)
+yy=total(bias)/n_elements(bias)
+print, xx, yy
+;
 splot, mean, bias, small=[2,3,3],/noer, xtitle='Mean ERAI Q2m', ytitle='ERAI Q2m bias', $
        psym=2, charsize=1.,  $
@@ -161 +161 @@
 xyouts, 14, 0.6, 'bias=-0.79', charsize=1.1
 xyouts, 14, 0.5, 'std ratio=0.92', charsize=1.1
-
-ab=linfit(mean, bias,yfit=yfit)
-a=float(ab(0)) 
-b=float(ab(1))
-oplot, mean, yfit, thick=3, color=50
-print, a,b
-
-;; significant test
+;
+ab=linfit(mean, bias,yfit=yfit)
+a=float(ab(0))
+b=float(ab(1))
+oplot, mean, yfit, thick=3, color=50
+print, a,b
+;
+; significant test
 significance_test_99,mean,bias
 fi=project_id_env+'ws_correction_final.txt'
@@ -177 +177 @@
 std=reform(ff(2,*))
 mean=reform(ff(4,*))
-
-xx=total(std)/n_elements(std)
-yy=total(bias)/n_elements(bias)
-print, xx, yy
-
-
+;
+xx=total(std)/n_elements(std)
+yy=total(bias)/n_elements(bias)
+print, xx, yy
+;
+;
 splot, mean, bias, small=[2,3,4],/noer, xtitle='Mean ERAI WS', ytitle='ERAI WS bias', $
         psym=2, charsize=1.,  $
@@ -189 +189 @@
 xyouts, 3, 0.7, 'bias=-0.28', charsize=1.1
 xyouts, 3, 0.6, 'std ratio=0.90', charsize=1.1
-
-ab=linfit(mean, bias,yfit=yfit)
-a=float(ab(0)) 
-b=float(ab(1))
-oplot, mean, yfit, thick=3, color=50
-
-print, a,b
-;; significant test
-significance_test_99,mean,bias
-
+;
+ab=linfit(mean, bias,yfit=yfit)
+a=float(ab(0))
+b=float(ab(1))
+oplot, mean, yfit, thick=3, color=50
+;
+print, a,b
+; significant test
+significance_test_99,mean,bias
+;
 closeps
-
-
+;
 end