Ignore:
Timestamp:
08/07/08 11:32:02 (16 years ago)
Author:
pinsard
Message:

improvements of headers (alignments of IDL prompt in examples)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SRC/Interpolation/spl_keep_mean.pro

    r327 r371  
    4949;    12 monthly values of precipitations into daily values: 
    5050; 
    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])) 
     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])) 
    6868; 
    6969; @history 
Note: See TracChangeset for help on using the changeset viewer.