Test simulations using Turning Bands method¶

import gstlearn as gl
import numpy as np           
import gstlearn.plot as gp
import matplotlib.pyplot as plt

nx = [100,100]
dx = [1,1]
grid = gl.DbGrid.create(nx=nx, dx=dx)
varioparam = gl.VarioParam.createMultipleFromGrid(grid, 50)

def represent(type, range=30, param=1):
    model = gl.Model.createFromParam(type=type, range=range, param=param)
    err = gl.simtub(None, grid, model, nbtuba=500)
    vario = gl.Vario.computeFromDb(varioparam, grid)
    
    # Plot
    fig = plt.figure(figsize=(10,5))
    ax1 = fig.add_subplot(1,2,1)
    ax1 = gp.raster(grid,name="Simu")
    ax2 = fig.add_subplot(1,2,2)
    ax2 = gp.varmod(vario, model)
    fig.decoration(title=type.getDescr())
    
    grid.deleteColumn("Simu")