<<prev file | next file >>    view single page | view frames    summary: fields | routine    details: routine

Interpolation/

extrapsmooth.pro

similar to extrapolate but could to the job in a better way because the extrapolated values are smoothed... takes more time than extrapolate. extrapolate data where mskin is equal 0 by filling step by step the coastline points with the mean value of the 8 neighbourgs.

topextrapsmooth Interpolation

result = extrapsmooth(in, mskin, X_PERIODIC=X_PERIODIC, MINVAL=scalar, MAXVAL=scalar, GE0=scalar 0 or 1, _EXTRA=_EXTRA)

Return value

the extrapolated array with no more masked values

Parameters

in        in required type: 2d array

data to be extrapolate

mskin        in required type: 2d array or -1

a 2D array, the land-sea mask of the output data (1 on ocean, 0 on land) put -1 if input data are not masked

Keywords

X_PERIODIC       

MINVAL        type: scalar default: not used

to specify a minimum value to the extrapolated values

MAXVAL        type: scalar default: not used

to specify a maximum value to the extrapolated values

GE0        type: scalar 0 or 1 default: 0

put 1 to force the extrapolated values to be larger than 0, same as using minval=0.

_EXTRA       

to be able to call extrapsmooth with _extra keyword

Examples

IDL> a=extrapsmooth(dist(jpi,jpj),tmask[*,*,0],/x_periodic) IDL> tvplus, a compare to extrapolate result: IDL> b=extrapolate(dist(jpi,jpj),tmask[*,*,0],/x_periodic) IDL> tvplus, b, window = 1

Version history

Version

$Id: extrapsmooth.pro 271 2007-08-30 12:44:23Z smasson $

History

January 2007: Sebastien Masson (smasson@lodyc.jussieu.fr)

Known issues

Restrictions

You cannot specify the number of iterations done in the extrapolation process

Statistics

McCabe cyclic 14
McCabe essential 1
McCabe modular design 1
Produced by IDLdoc 2.0.