Changeset 94 for trunk/src/paper01/fig7/q2m_validation_scatter_2000_2009_basin.pro

Timestamp:

08/25/11 13:02:20 (13 years ago)

Author:

pinsard

Message:

get rid of multistatement lines

File:

• trunk/src/paper01/fig7/q2m_validation_scatter_2000_2009_basin.pro (modified) (12 diffs)

trunk/src/paper01/fig7/q2m_validation_scatter_2000_2009_basin.pro

@@ -179 +179 @@
 ;; reading data from mooring
 
-    site=sitelist(n) & csite=site
+    site=sitelist(n) 
+    csite=site
     print, csite
     x=x_site_location(site)
     y=y_site_location(site)
     if (y ge 0. and y le 30.) then y=y+360.
-    dx=0.5 & dy=0.5 & box=[y-dy, y+dy, x-dx, x+dx]
+    dx=0.5 
+    dy=0.5 
+    box=[y-dy, y+dy, x-dx, x+dx]
     read_variables_v2, csite,date1,date2,nsmooth, $
          at, sw,rh,sst,wu,wv,ws,lh
@@ -220 +223 @@
      ncep1=tropflux
 
-    ind=where(finite(q2m)) & q2m=q2m(ind) & uncr_q2m=uncr(ind)  & corr_q2m=corr(ind)
-    oafl=oafl(ind) & ncep=ncep(ind) & ncep1=ncep1(ind)
+    ind=where(finite(q2m)) 
+    q2m=q2m(ind) 
+    uncr_q2m=uncr(ind)  
+    corr_q2m=corr(ind)
+    oafl=oafl(ind) 
+    ncep=ncep(ind) 
+    ncep1=ncep1(ind)
     mean_tao=total(q2m,/nan)/n_elements(ind)
 
@@ -247 +255 @@
 lat=reform(ff(0,*))
 lon=reform(ff(1,*))
-cor_era=reform(ff(2,*))   & cor_erai=total(cor_era)/n_elements(cor_era)
-bias_era=reform(ff(3,*))  & bias_erai=total(bias_era)/n_elements(bias_era)
-std_era=reform(ff(4,*))   & std_erai=total(std_era)/n_elements(std_era)
-rmsd_era=reform(ff(5,*))  & rmsd_erai=total(rmsd_era)/n_elements(rmsd_era)
-mean_tao=reform(ff(6,*))  & mean_erai=bias_era+mean_tao
+cor_era=reform(ff(2,*))   
+cor_erai=total(cor_era)/n_elements(cor_era)
+bias_era=reform(ff(3,*))  
+bias_erai=total(bias_era)/n_elements(bias_era)
+std_era=reform(ff(4,*))   
+std_erai=total(std_era)/n_elements(std_era)
+rmsd_era=reform(ff(5,*))  
+rmsd_erai=total(rmsd_era)/n_elements(rmsd_era)
+mean_tao=reform(ff(6,*))  
+mean_erai=bias_era+mean_tao
 
 print, ''
@@ -265 +278 @@
 oplot, [11,21], [11,21]
 ab=linfit(mean_tao, mean_erai,yfit=yfit)
-a=float(ab(0)) & b=float(ab(1))
+a=float(ab(0)) 
+b=float(ab(1))
 oplot, mean_tao, yfit, color=250, thick=2
 
@@ -273 +287 @@
 lat=reform(ff(0,*))
 lon=reform(ff(1,*))
-cor_tro=reform(ff(2,*))   & cor_trop=total(cor_tro)/n_elements(cor_tro)
-bias_tro=reform(ff(3,*))  & bias_trop=total(bias_tro)/n_elements(bias_tro)
-std_tro=reform(ff(4,*))   & std_trop=total(std_tro)/n_elements(std_tro)
-rmsd_tro=reform(ff(5,*))  & rmsd_trop=total(rmsd_tro)/n_elements(rmsd_tro)
-mean_tao=reform(ff(6,*))  & mean_trop=bias_tro+mean_tao
+cor_tro=reform(ff(2,*))   
+cor_trop=total(cor_tro)/n_elements(cor_tro)
+bias_tro=reform(ff(3,*))  
+bias_trop=total(bias_tro)/n_elements(bias_tro)
+std_tro=reform(ff(4,*))   
+std_trop=total(std_tro)/n_elements(std_tro)
+rmsd_tro=reform(ff(5,*))  
+rmsd_trop=total(rmsd_tro)/n_elements(rmsd_tro)
+mean_tao=reform(ff(6,*))  
+mean_trop=bias_tro+mean_tao
 
 print, ''
@@ -292 +311 @@
 
 ab=linfit(mean_tao, mean_trop,yfit=yfit)
-a=float(ab(0)) & b=float(ab(1))
+a=float(ab(0)) 
+b=float(ab(1))
 oplot, mean_tao, yfit, color=250, thick=2
 
@@ -300 +320 @@
 lat=reform(ff(0,*))
 lon=reform(ff(1,*))
-cor_oaf=reform(ff(2,*))   & cor_oafl=total(cor_oaf)/n_elements(cor_oaf)
-bias_oaf=reform(ff(3,*))  & bias_oafl=total(bias_oaf)/n_elements(bias_oaf)
-std_oaf=reform(ff(4,*))   & std_oafl=total(std_oaf)/n_elements(std_oaf)
-rmsd_oaf=reform(ff(5,*))  & rmsd_oafl=total(rmsd_oaf)/n_elements(rmsd_oaf)
-mean_tao=reform(ff(6,*))  & mean_oafl=bias_oaf+mean_tao
+cor_oaf=reform(ff(2,*))   
+cor_oafl=total(cor_oaf)/n_elements(cor_oaf)
+bias_oaf=reform(ff(3,*))  
+bias_oafl=total(bias_oaf)/n_elements(bias_oaf)
+std_oaf=reform(ff(4,*))   
+std_oafl=total(std_oaf)/n_elements(std_oaf)
+rmsd_oaf=reform(ff(5,*))  
+rmsd_oafl=total(rmsd_oaf)/n_elements(rmsd_oaf)
+mean_tao=reform(ff(6,*))  
+mean_oafl=bias_oaf+mean_tao
 
 print, ''
@@ -319 +344 @@
 
 ab=linfit(mean_tao, mean_oafl,yfit=yfit)
-a=float(ab(0)) & b=float(ab(1))
+a=float(ab(0)) 
+b=float(ab(1))
 oplot, mean_tao, yfit, color=250, thick=2
 
@@ -328 +354 @@
 lat=reform(ff(0,*))
 lon=reform(ff(1,*))
-cor_nce=reform(ff(2,*))   & cor_ncep=total(cor_nce)/n_elements(cor_nce)
-bias_nce=reform(ff(3,*))  & bias_ncep=total(bias_nce)/n_elements(bias_nce)
-std_nce=reform(ff(4,*))   & std_ncep=total(std_nce)/n_elements(std_nce)
-rmsd_nce=reform(ff(5,*))  & rmsd_ncep=total(rmsd_nce)/n_elements(rmsd_nce)
-mean_tao=reform(ff(6,*))  & mean_ncep=bias_nce+mean_tao
+cor_nce=reform(ff(2,*))   
+cor_ncep=total(cor_nce)/n_elements(cor_nce)
+bias_nce=reform(ff(3,*))  
+bias_ncep=total(bias_nce)/n_elements(bias_nce)
+std_nce=reform(ff(4,*))   
+std_ncep=total(std_nce)/n_elements(std_nce)
+rmsd_nce=reform(ff(5,*))  
+rmsd_ncep=total(rmsd_nce)/n_elements(rmsd_nce)
+mean_tao=reform(ff(6,*))  
+mean_ncep=bias_nce+mean_tao
 
 print, ''
@@ -347 +378 @@
 
 ab=linfit(mean_tao, mean_ncep,yfit=yfit)
-a=float(ab(0)) & b=float(ab(1))
+a=float(ab(0)) 
+b=float(ab(1))
 oplot, mean_tao, yfit, color=250, thick=2
 
@@ -355 +387 @@
 lat=reform(ff(0,*))
 lon=reform(ff(1,*))
-cor_nce=reform(ff(2,*))   & cor_ncep=total(cor_nce)/n_elements(cor_nce)
-bias_nce=reform(ff(3,*))  & bias_ncep=total(bias_nce)/n_elements(bias_nce)
-std_nce=reform(ff(4,*))   & std_ncep=total(std_nce)/n_elements(std_nce)
-rmsd_nce=reform(ff(5,*))  & rmsd_ncep=total(rmsd_nce)/n_elements(rmsd_nce)
-mean_tao=reform(ff(6,*))  & mean_ncep=bias_nce+mean_tao
+cor_nce=reform(ff(2,*))   
+cor_ncep=total(cor_nce)/n_elements(cor_nce)
+bias_nce=reform(ff(3,*))  
+bias_ncep=total(bias_nce)/n_elements(bias_nce)
+std_nce=reform(ff(4,*))   
+std_ncep=total(std_nce)/n_elements(std_nce)
+rmsd_nce=reform(ff(5,*))  
+rmsd_ncep=total(rmsd_nce)/n_elements(rmsd_nce)
+mean_tao=reform(ff(6,*))  
+mean_ncep=bias_nce+mean_tao
 
 print, ''
@@ -374 +411 @@
 
 ab=linfit(mean_tao, mean_ncep,yfit=yfit)
-a=float(ab(0)) & b=float(ab(1))
+a=float(ab(0)) 
+b=float(ab(1))
 oplot, mean_tao, yfit, color=250, thick=2