Changeset 97 for trunk/src/paper01/fig6/statistics.pro
- Timestamp:
- 08/30/11 14:38:29 (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/paper01/fig6/statistics.pro
r94 r97 60 60 cor, bias, std, rmsd 61 61 @cm_project 62 ; ;mean1= mean of mooring variable std1=std of mooring variable63 ; ;mean2= mean of erai variable std2=std of erai variable64 x=moor 62 ; mean1= mean of mooring variable std1=std of mooring variable 63 ; mean2= mean of erai variable std2=std of erai variable 64 x=moor 65 65 y=erai 66 ind1=where(finite(x,/nan)) 66 ind1=where(finite(x,/nan)) 67 67 ind2=where(finite(y,/nan)) 68 68 ind=inter(ind1,ind2) 69 69 if (ind ne -1 ) then begin 70 x(ind)=!Values.f_nan 70 x(ind)=!Values.f_nan 71 71 y(ind)=!Values.f_nan 72 72 endif 73 73 ; 74 74 ind=where(finite(x)) 75 x=x(ind) 75 x=x(ind) 76 76 y=y(ind) 77 77 ; 78 78 nn=n_elements(ind) 79 79 ; 80 80 cor=correlate(x,y) 81 81 xxx=moment(x) 82 mean1=xxx(0) 82 mean1=xxx(0) 83 83 std1=sqrt(xxx(1)) 84 84 xxx=moment(y) 85 mean2=xxx(0) 85 mean2=xxx(0) 86 86 std2=sqrt(xxx(1)) 87 bias=mean2-mean1 87 bias=mean2-mean1 88 88 std=std2/std1 89 89 rmsd=sqrt(total((x-y)*(x-y))/n_elements(x)) 90 ; 90 91 end
Note: See TracChangeset
for help on using the changeset viewer.