Changeset 371 for trunk/SRC/Interpolation/spl_keep_mean.pro
- Timestamp:
- 08/07/08 11:32:02 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/Interpolation/spl_keep_mean.pro
r327 r371 49 49 ; 12 monthly values of precipitations into daily values: 50 50 ; 51 ; IDL> yr1 = 199052 ; IDL> yr2 = 199253 ; IDL> nyr = yr2-yr1+154 ; IDL> n1 = 12*nyr+155 ; IDL> x = julday(1+findgen(n1), replicate(1, n1) $56 ; IDL> , replicate(yr1, n1), fltarr(n1))57 ; IDL> n2 = 365*nyr + total(leapyr(yr1+indgen(nyr))) + 158 ; IDL> x2 = julday(replicate(1, n2), 1+findgen(n2) $59 ; IDL> , replicate(yr1, n2), fltarr(n2))60 ; IDL> y = abs(randomn(0, n1-1))61 ; IDL> y2 = spl_keep_mean(x, y, x2, /ge0)62 63 ; IDL> print, min(x, max = ma), ma64 ; IDL> print, min(x2, max = ma), ma65 ; IDL> print, vairdate([min(x, max = ma), ma])66 ; IDL> print, total(y*(x[1:n1-1]-x[0:n1-2]))67 ; IDL> print, total(y2*(x2[1:n2-1]-x2[0:n2-2]))51 ; IDL> yr1 = 1990 52 ; IDL> yr2 = 1992 53 ; IDL> nyr = yr2-yr1+1 54 ; IDL> n1 = 12*nyr+1 55 ; IDL> x = julday(1+findgen(n1), replicate(1, n1) $ 56 ; IDL> , replicate(yr1, n1), fltarr(n1)) 57 ; IDL> n2 = 365*nyr + total(leapyr(yr1+indgen(nyr))) + 1 58 ; IDL> x2 = julday(replicate(1, n2), 1+findgen(n2) $ 59 ; IDL> , replicate(yr1, n2), fltarr(n2)) 60 ; IDL> y = abs(randomn(0, n1-1)) 61 ; IDL> y2 = spl_keep_mean(x, y, x2, /ge0) 62 ; 63 ; IDL> print, min(x, max = ma), ma 64 ; IDL> print, min(x2, max = ma), ma 65 ; IDL> print, vairdate([min(x, max = ma), ma]) 66 ; IDL> print, total(y*(x[1:n1-1]-x[0:n1-2])) 67 ; IDL> print, total(y2*(x2[1:n2-1]-x2[0:n2-2])) 68 68 ; 69 69 ; @history
Note: See TracChangeset
for help on using the changeset viewer.