initmetagauss {rainfarmr}

Generate the spectral amplitudes for a metagaussian field

Description

A Fourier spectrum with prescribed logarithimic slope and zero phases is constructed.

Usage

initmetagauss(sx, ns)

Arguments

sx

spectral slope for the output field. The convention is that this is the slope of k*|A(k)|^2.

ns

size of the output field.

Value

Output complex field in Fourier space with specified spectral slope, with dimensions c(ns, ns).

Author(s)

Jost von Hardenberg, [e-mail]

Examples

# Make a synthetic rainfall field with prescribed logarithmic spectral slope
f = initmetagauss(1.7, 64)
r = metagauss(f)
# Check slope of the resulting field
fx <- fft2d(r)
fitslope(fx)
# 1.640373