def init(): global minthick, maxlayers, cent, threshdist, defmaxlayers global fineBottomFlag fineBottomFlag = 1 minthick = 5 defmaxlayers = 2 cent = (8700, 17093.376121) threshdist = 6000 ncd.nc_init("gmsh/islay-jura-bathy.grd")
def init(finebottom=0, layers=5): global minthick, maxlayers, cent, threshdist, defmaxlayers global fineBottomFlag print "*** Initialising netcdf bathy data" fineBottomFlag = finebottom minthick = 5 defmaxlayers = layers cent = (8700, 17093.376121) threshdist = 5000 ncd.nc_init("gmsh/islay-jura-bathy.grd")
def init(): global u_tau, nu mu = 1.8e-4 density = 1.027 nu = mu / density kappa = 0.41 ks = 3.5 * 5 y0 = ks / 30.0 # Assume flow speed 2.0 m/s, at 30m above seabed udef = 2.0 ydef = 30.0 u_tau = kappa * udef / math.log(ydef / y0) global dist_init try: dist_init except NameError: import nc_depth_io as ncd import math ncd.nc_init("gmsh/islay-jura-bathy.grd") print "*** Initialising distance to wall field" init() dist_init = 1