topspl_keep_mean
Return value
y2: the meean value between two consecutive values of x2. This
array has one element less than y2. y2 has double precision.
Parameters
x
in
required
yin
in
required
x2
in
required
Keywords
YP0
YPN_1
GE0
Examples
y2 = spl_keep_mean(x, y, x2)
12 monthly values of precipitations into daily values:
yr1 = 1990
yr2 = 1992
nyr = yr2-yr1+1
n1 = 12*nyr+1
x = julday(1+findgen(n1), replicate(1, n1) $
, replicate(yr1, n1), fltarr(n1))
n2 = 365*nyr + total(leapyr(yr1+indgen(nyr))) + 1
x2 = julday(replicate(1, n2), 1+findgen(n2) $
, replicate(yr1, n2), fltarr(n2))
y = abs(randomn(0, n1-1))
y2 = spl_keep_mean(x, y, x2, /ge0)
print, min(x, max = ma), ma
print, min(x2, max = ma), ma
print, vairdate([min(x, max = ma), ma])
print, total(y*(x[1:n1-1]-x[0:n1-2]))
print, total(y2*(x2[1:n2-1]-x2[0:n2-2]))
Version history
History
Sebastien Masson (smasson@lodyc.jussieu.fr): May 2005
Known issues
Restrictions
It might be possible that y2 has very small negative values
(amplitude smaller than 1.e-6)...