Input file - Geodels/geodels-libs GitHub Wiki
About | Input file | Docker Installation | Examples
← Previous topic: About | Next topic: Docker Installation →
The file below provide the list of all the parameters that can be defined in sedPro
:
###############################################################
# Required parameters section #
##############################################################
# sedPro 2019
#
TITLE
# One-line title / comment on experiment (required)
#
Enter experiment name here
#
# Include general comments about what modules are turned on and
# what is different in this simulation than previous iterations here
#
#----------------------------------------------------------
TIME
# Time parameters (required)
#
# Simulation start time [years] End time [years]
0.0 10.0
# Use negative value to represent time before present.
#eg a simulation from 21 Ma to 19.6 Ma should be written as -21000000.0 -19600000.0
# Display interval [years] Flow sampling interval [years]
10000.0 2000.0
# Display interval decides how often the results files are updated
# Flow sampling interval decides how often fluid elemnets are released from the source
#------------------------------------------------------------------------
GRID
# Grid size definitions and geometry
#
# Grid spacing [m] Number of rows Number of columns
5000.0 60 60
# Lower left (SW) corner coordinates Base level elevation
464500.0 7930000.0 -180.0
# Topography grid file name
<file_name>.top
#
#---------------------------------------------------------
SEDIMENTS
# Sediment Parameters (required)
#
# Line 1: Diameter of each grain size [mm]
# Line 2: Density of each grain size [kg/m3]
# Line 3: 1 - suspension (normal type), 0 - bed load (new)
#
# pebble 4-64 mm, granule 2-4 mm, vcse 1-2 mm, cse 0.5-1 mm, med 0.25-0.5 mm,
# fn 0.125-0.25, vfn 0.062-0.125 mm, slt 0.0039-0.062 mm, clay < 0.0039 mm
#
# Coarse Medium Fine Finest
0.28 0.15 0.03 0.0003
2650.00 2600.00 2600.00 2550.00
1 1 1 1
#---------------------------------------------------------
#CONTINUE
# Continue a previous experiment (optional)
# Give the start time of the initial simulation and the hardness of the previous deposit (optional can leave it blank)
-140000.0 1.0
# Name of the inital simulation without .sif extention
# if a fine grid is present, the name of the previous initial simulation
samplerun1
samplerun_fine
#----------------------------------------------------------
#INTERNAL GRID
# experiment name for the internal grid.
<experiment_name>
# x,y location of the grid (SW corner)
26000 22000
# grid spacing(m), number of fine row, number of columns
250 41 41
#-----------------------------------------------------------
#INTERNAL TOPOGRAPHY GRID
# topography grid file name. (optional)
# If you have the internal grid turned on, you can specify the initial surface
# if you don't do this, one will be generated from the coarse grid
<high_res_topography_name>.top
#
#----------------------------------------------------------
#DEPOSIT
# Sediment deposited prior to this run. Do not use (optional)
# CONTINUE and DEPOSIT at the same time. Sediment initailly deposited will
# be recorded anyway
# <deposit_file>.spc
#
# Deposit file name
<initial_deposit_file_name>.spc
#-----------------------------------------------------------
#SEA LEVEL
# Define sea level curve (optional)
# if no file defined then a constant sea level will be the base level in GRID section
# Sea level curve file
<sea-level-file-name>.sl
#
#---------------------------------------------------------
#SOURCES
# Guide to flow regime input:
# A. V (m/s)
# From Hjulstrom(1939) we get:
# Taken from Krumbein &Sloss 1963) p.203
#
# The depositional flow threshold seems to be
#
# V (m/s) < 4.49e-07 + 0.01d(mm)+ 9.92e-06d^2 -8.983e-07d^3
#
# Erosional threshold for bed sediments below and including 0.1mm grain size
#
# V (m/s) > 0.0062145 * d(mm)^-0.5
#
# For grain sizes over and including 0.1mm
#
# V(m/s)= 0.020-0.00116*d +0.001d^2 -1.03e-05d^3 (R^2=1)
#
# Another rule of thumb from Ackers(1964) quoted by Schumm(1987)p.173
# only for discharges below 0.2 m3/s
# is V(m/s) = 1.92Q(m3/s)^0.15
# or Q(m3/s) = 0.0129*V(m/s)^6.67
#
# velocity as a function of discharge (Tetlaff and Harbaugh, 1989, p.58)
# for a slope of 0.01 (0.6deg) v = 0.33842Q^0.36937
# for a slope of 0.001 (0.06deg) v = 0.22543Q^0.32191
#
# B. Q (m3/s)
# From Allen 1970 p.128 we get the relationship
# Q(m3/s) = 5.17e-05*Slope^-2.27
# for the flow boundary between braided and meandering streams (for
# a given slope a higher discharge than Q gives braided streams, lower gives
# meandering streams). This can be used to give a rough ball-park check
# on discharge for a given channellized slope.
#
# C. concentration 'c' (kg/m3)
# Typical concentration values for natural river systems are:
# 0.16 kg/m3 (Colorado river - higher discharge)
# 1.6 kg/m3 (Niobrara river - lower discharge )
#
# Schumm used concentrations of around 0.44 kg/m3 with flume
# slopes of 0.008(0.45deg) and discharge rates of 0.0057 m3/s
# in river channel experiments.
# For fan-delta experiments he used 114 kg/m3 with a median grain size of
# 0.005 and0.5mm
# For turbidity currents volumic concentration ranges between a=0.1% and b=10%
# which gives a concentration between a*sed_density/(1-a) and b*sed_density/(1-b)
# with sed_density=2650 kg/m3 it gives 2.65<c<295
# Definition of sources that are constant throughout the experiment (required)
# One line per source, entries are:
# Unique Source identification number (matches with the number in the SOURCES section)
# Source location (x,y) [m]
# Velocity at source (vx,vy) [m/sec]
# Discharge rate (Q) [m3/sec]
# Sediment concentration (c) [kg/m3]
# Sediment composition (C coarse, [%]
# M medium,
# F fine,
# FF finest)
#
# Source ID# t1 t2 x y vx vy Q c %C %M %F %FF
#SE (A) source
1 -21000000.0 -19600000.0 751498 7947823 -0.2 0.2 100.0 0.2 30.0 25.0 45.0 0.0
#SE (A2)source
7 -20000000.0 -19800000.00 751356 7958405 -0.2 0.2 100.0 0.2 30.0 25.0 45.0 0.0
#South (B) source
10 -21000000.0 -19600000.0 586581 7938297 -0.15 0.05 100.0 0.2 15.0 15.0 30.0 40.0
#North E (C)source
9 -20000000.0 -19900000.0 741377 8071659 -0.05 -0.05 40.0 0.2 5.0 15.0 30.0 50.0
*
######################################################################
# Additional module parameters (modules will not run if not included)
##################################################################
#----------------------------------------------------------
SOURCE HEIGHTS
# list ends in '*'
# Unique Source identification number (matches with the number in the SOURCES section)
# flow types 0=normal 1=turbidite 2=debris flow
# for debris flows the concentration must be high enough, otherwise will
# transition to a turbidite (default transition at 60 kg/m3)
#
# Source ID# flow height(m) flow type
1 1.0 0
10 0.5 1
*
#--------------------------------------------------------
#FLUVIAL
# turn this on when you have flow in a well defined channel, which is at a higher than grid resolution (optional)
#-----------------------------------------------------------
WELL LOCATIONS
# locations in meters from sedPro-grid origin (draws them in sedview)
# use this to identify any point of interest on the output visualisation
# X Y well name
751498 7947823 source_1
586581 7938297 source_2
741377 8071659 source_3
481255 8101925 east mermaid_1
658611 8093097 wamac_1
653706 8109969 lacepede_1
* end of well data
#--------------------------------------------------------------
#LAKES
# If you use LAKES you need to define the ocean as a 'lake' in this section
# and you don't need to use the SEA LEVEL module
# For each lake provide the following information:
# Lake ID SW corner coords. rows cols lake level sal/fresh water
# saline water is 0 and fresh is 1
# you need to verify that the lake or sea level you give here is equal to
# the one in the water level curve file for the simulation starting time
1 499500.0 5000000.0 20 30 18.0 1
*
# Define water level variation for each lake
# Water level curve files must be ordered like lake definition
leman_lake_curve.dat
*
#--------------------------------------------------------------
#TECTONICS
# Define tectonic movement (optional)
#
# Tectonic movement file name
<tectonic_movement_file>.tec
#--------------------------------------------------------------
#INTERNAL GRID TECTONICS
# Define tectonic movement of the internal (optional)
#
# Tectonic movement file name
<tectonic_movement_file>.tec
#
#-----------------------------------------------------------
#SLOPE FAILURE
# Parameters for calculating slope failure (optional)
#
# Slope angles should be greater than that set in the SLOPE section.
# Typical turbidite volume concentrations : Bagnold maximum of C=0.09
# and those referred to in Piper and Savoye (1993) for the Var fans of
# 0.015 to 0.12 refer to volume concentration. Thus 0.1 volume concentration
# for a sediment of density 2650 kg/m3 , c = 265 kg/m3
# and for C = 0.05, c=132 kg/m3
#
# Maximum Subaerial and Marine slope before failure can occur (dz/dx) (dz/dx)
# Minimum and maximum height allowed to fail(m) (bigger heights will be split into multiple fluid elements)
# Concentrations slumped material after release (kg/m^3)
# Concentration of the material where it transitions from a debris flow to a tubidite (kg/m^3)
0.002 0.05
0.5 5.0
100.0
60.0
#---------------------------------------------------------
#CONTOUR CURRENT
# current calling interval [a]
100000.0
# starting time, ending time [a]
# maximum and minimum elevation where current occurs [m]
# flow rotation clockwise=1, counter-clockwise=-1,decide-by-code=0
# maximum near bottom flow strength m/s
# the latitude at Lower left (SW) corner, + North, - South [degree]
# sediment supply condition, supply 1, No-supply=0
# current supplied source sediment combination sum=1
#
# t1 t2 maxd mind clockwise velocity (m/s) lat supply %C %M %F %FF %rudist %pelagic %org1 %org2
-21000000.0 -19600000.0 -500 -600 1 0.2 -40 1 10.0 80.0 10 0 0 0 0 0
*
#---------------------------------------------------------
#CIRCULATION (CYCLIC)
# This module allows the input of external current data. the format for the
#input file is
# header line
# x velocities -one for each grid point in standard matrix format
# header line
# y velocity -one for each grid point in standard matrix format
# header line
# th value -total sheilds parameter (stress is on the sea bed)
# (th lifts sediment -velocities move it)
# Frequency factor controls Circulation calling interval.
# The bigger the factor the longer the interval
# files will cycle until the end of the simulation
# astrix ends input list
# duration [years], frequency factor, input_file_name
20.0 10.0 Current_u_v_th.dat
30.0 10.0 Current_u_v_th_2.dat
*
#---------------------------------------------------------
#CIRCULATION (LINEAR)
# This module allows the input of external current data. the format for the
#input file is
# header line
# x velocities -one for each grid point in standard matrix format
# header line
# y velocity -one for each grid point in standard matrix format
# header line
# th value -total Sheilds parameter (stress is on the sea bed)
# (th lifts sediment -velocities move it)
# Frequency factor controls Circulation calling interval.
# The bigger the factor the longer the interval
# astrix ends input list
# start_time, end_time, frequency factor, input_file_name
-21000000.0 -20000000.0 10.0 Current_u_v_th.dat
-20000000.0 -19600000.0 10.0 Current_u_v_th_2.dat
*
#-----------------------------------------------------------
#STORM
# Provide frequency and individual storm parameters (optional)
#
# Angle of mean storm direction [deg]
# Deviation of the storm direction [deg]
# Significant storm waveheight [m]
# Duration of storm [h]
# Mean storm return time [y] (put 0.0 if exact storm records known)
240.0 0.0 2.0 5.0 0.0
#if the Mean storm return time is set as 0.0, then continue read
#exact time when storm happened
# The four columns in storm data file contain
# Time[year] Incidence angle[deg] Height[m] Duration[hour]
WIND_WAVEhourly_92_97.stm
#----------------------------------------------------------
#STORM REFRACTION
# Calculate STORMs refraction,diffraction and breaking over the coastal region.
# (optional)
# Dip direction of the gereral shoreline [0.0 ~ 360.0][deg]
270.0
#----------------------------------------------------------
#WAVES
# This module must be turned on to run any form of waves
# In addition, ONE of the FOUR other wave module must be turned on to run waves
# The other four modules are in order of complexity
# WAVE (RATE)
# WAVE (HEIGHT)
# WAVE (GRID-CYCLIC)
# WAVE (GRID-LINEAR)
# in addition
# WAVE REFRACTION may be turned on with WAVE (HEIGHT) or WAVE (RATE)
#
# Sediment supply conditions at 4 Boundaries
# The Non-supply condition represents the rocky cliff shoreline
# where waves can not transport sediment into the simulation area.
# The normal condition represents normal sandy or mudy beach, sediment
# can be transported into the simulation area.
# South East North West 1==Normal, 0=No supply
1 1 1 1
#
# Berm height [m]
# Minimum mobile bed thickness [m]
# Residual thickness for deficit of erosion [m]
# Amplitude of local sea level fluctuation [m]
# Depth of mobile bed [m]
# Wind Shadow [m]
0.2 0.01 0.0001 0.5 0.02 0.0
# Note on Wind Shadow: this is the distance of the wind shadow created by land.
# Entering 0.0 will result in an infinitely long shadow
# Any distance less than or equal to 1 grid cell will result in all sea areas being exposed to waves
#
#----------------------------------------------------------
#WAVE (RATE)
# Note: WAVES module must also be turned on
# Note: All angles are specified by
# N (0 or 360)
# |
# |
# (270) W-----------E (90)
# |
# |
# S (180)
#
# Input lines of (input start at end of sequence)
# Start_time[yrs], End time[yrs], Angle of incidence[deg], Transport rate[m3/sec], Wave base[m]
-1000 -900 110 0.0000002 5.0
-900 -500 50 0.0000002 4.0
-100 0.0 210 0.0000002 6.0
*
#----------------------------------------------------------
#WAVE (HEIGHT)
# Note: WAVES module must also be turned on
# Enter a file name with the wave data
# The file contains
# Time[a], Angle of incidence[deg], Significant wave height(m),
# Note: All angles are specified by
# N (0 or 360)
# |
# |
# (270) W-----------E (90)
# |
# |
# S (180)
# The following is a example of the format
# 1990.00 180.0 2.0
# 1990.25 270.0 2.0
# 1990.30 360.0 2.0
WIND_WAVEhourly_92_97.wav
#----------------------------------------------------------
# WAVE (GRID-CYCLIC)
# Note: WAVES module must also be turned on
# Enter a list of time durations[year] and a filename containing wave data
# The list terminates in a *
# After all the files have been read, the list will repeat
# Each file contains the
# Significant wave height[m] and Angle of incidence[deg] in block format e.g.
# j=1 j=number of columns
# Hs(1,1) Hs(1,2)..........Hs(1,numcol)
# ... ... ......... ...
# ... ... ......... ...
# ... ... ......... ...
# Hs(numrow,1) ............Hs(numrow,numcol)
# Ang(1,1) Ang(1,2)..........Ang(1,numcol)
# ... ... ......... ...
# ... ... ......... ...
# ... ... ......... ...
# Ang(numrow,1) ............Ang(numrow,numcol)
# Note: All angles are specified by
# N (0 or 360)
# |
# |
# (270) W-----------E (90)
# |
# |
# S (180)
#
# 0.125 WIND_WAVEdec_jan_feb.wvs
# 0.25 WIND_WAVEmar_apr.may.wvs
# 0.25 WIND_WAVEjun_jul_aug.wvs
# 0.25 WIND_WAVEsep_oct_nov.wvs
# 0.125 WIND_WAVEdec_jan_feb.wvs
# *
#----------------------------------------------------------
#WAVE (GRID-LINEAR)
# Note: WAVES module must also be turned on
# Enter a list of start time[year] and end time[year] and a filename containing wave data
# the list terminates in a *
# After all the files have been read, the list will repeat
# Each file contains the
# Significant wave height[m] and Angle of incidence[deg] in block format e.g.
# j=1 j=number of columns
# Hs(1,1) Hs(1,2)..........Hs(1,numcol)
# ... ... ......... ...
# ... ... ......... ...
# ... ... ......... ...
# Hs(numrow,1) ............Hs(numrow,numcol)
# Ang(1,1) Ang(1,2)..........Ang(1,numcol)
# ... ... ......... ...
# ... ... ......... ...
# ... ... ......... ...
# Ang(numrow,1) ............Ang(numrow,numcol)
#Note: All angles are specified by
# N (0 or 360)
# |
# |
# (270) W-----------E (90)
# |
# |
# S (180)
#
# 2005.000 2005.125 WIND_WAVEdec_jan_feb.wvs
# 2005.125 2005.375 WIND_WAVEmar_apr_may.wvs
# 2005.375 2005.625 WIND_WAVEjun_jul_aug.wvs
# 2005.625 2005.875 WIND_WAVEsep_oct_nov.wvs
# 2005.875 2006.125 WIND_WAVEdec_jan_feb.wvs
# 2006.125 2006.375 WIND_WAVEmar_apr_may.wvs
# 2006.375 2006.625 WIND_WAVEjun_jul_aug.wvs
# 2006.625 2006.875 WIND_WAVEsep_oct_nov.wvs
# 2006.875 2007.125 WIND_WAVEdec_jan_feb.wvs
# 2007.125 2007.375 WIND_WAVEmar_apr_may.wvs
# 2007.375 2007.625 WIND_WAVEjun_jul_aug.wvs
# 2007.625 2007.875 WIND_WAVEsep_oct_nov.wvs
# 2006.875 2008.000 WIND_WAVEdec_jan_feb.wvs
# *
#-----------------------------------------------------------
#WAVE REFRACTION
# Calculate wave,refraction, diffraction and breaking over the coastal region
# (optional)
# Dip direction of the gereral shoreline [0.0 ~ 360.0][deg]
# 270.0
#--------------------------------------------------------------
#BASEMENT STRUCTURE
#name of the basement structure file (dep file)
xxx.DEP
#--------------------------------------------------------------
#LANDSCAPE
# Define rainfall cyclicity (optional)
#
# Pluviometry file name pluviometry must be in (m/y)
rain.plv
# Number of cells per side that defined a moving square frame which is used
# to average the precipitation over a bigger surface to enhance CPU time
1
# Elevation under which rain cannot be found
900
#
#--------------------------------------------------------------
#COMPACTION
# Enable compaction module; (optional)
# Depth of post depo burial by sediment (m) and depth of post depo burial by water (m)
1000 500
#--------------------------------------------------------------
#ISOSTASY
# Calculate isostatic subsidence (optional)
#
# Calling interval [y]
1000
# enter changing crustal properties, finish list with a *
# Start time [y], End time[y], Mantle density [kg/m3] Flexural rigidity [Nm2]
0.0 100000.0 3500 1.0E23
200000.0 300000.0 3500 1.0E23
*
# Time periods not coveered above will result in zero isostaic movement.
# entering a single time period with start and end time of 0, will result in the same parameters being used for the entire simulation
#
# NOTE:if continuing a simulation
# list all previous intervals when isostasy was run, otherwise the continue may be incorrect.
#-----------------------------------------------------------
#CARBONATES AND ORGANICS
# Parameters for calculating carbonate development (optional)
# (It is strongly recomended that waves are turned on with this module)
# calling interval (years)
100
# carb/org grain size diameter (one for each grain) after comminution
0.3 0.3 0.1 0.1
# carb grain density (one for each grain)
2050.00 2050.00 1800.0 1800.0
# temperature and salinity file (this is a list of year, sea surface temperature[degrees C] and salinity [ppm])
# to let sedPro interpolate between unknown values set them to -99 (warning: first value must always be defined)
# eg 2005 27.0 30000
# 2005.2 27.2 -99
# 2005.7 30.0 30250 etc.
# This line is ignored if TEMPERATURE GRID is used
#
temp_static.tmprt
#temp_1Ma.tmprt
#temp_sine_1Ma.tmprt
#
#----------------------------------------------------------
#EXPERT CARBONATES AND ORGANICS
#
# Use this if the default settings of the carbonates need modifying.
#
# the format for defining membership functions is
# function_name codevariable xpoint1 ypoint1 xpoint2 ypoint2 xpoint3 ypoint3 etc
#
# function_name - user defined variable name
# code_variable - must be one of the key words
# OUTPUT VARIABLES
# growth - defines the function as a output membership function controlling carb/org growth rate [m/yr]
# erosion - defines the function as a output membership function controlling carbonate growth rate [m/yr]
# porosity - defines the function as a output membership function controlling the porosity [0 - 1]
# hardness - defines the function as a output membership function controlling the porosity [ >1 ]
# INPUT VARIABLES
# depth - water depth [m]
# current - current velocity in the area - you must have circulation turned on)[m/s]
# shore - distance in metres to the shoreline[m]
# rivdist - distance to the nearest fluid element or river source [m]
# sedrate - sediment rate [m/year]
# temp - surface temperature[degrees celcius]
# expo - The value is calculated by ray tracing from eight different angles around the boundary
# of the simulation area. Where the ray is stopped when it encounters exposed land.
# Values vary between 0 and 1, with 0 being a completely land locked body of water (or land itself)
# through to 1 which would represent water completely exposed from all angles.
# expotime ñ time that the layer has been undisturbed on the surface [years]
# carbdist - distance to nearest surface exposed accumulation of carbonate1 [m]
# orgdist - distance to nearest surface exposed accumulation of organic1 or organic2 [metres]
# salinity - degree of salinity in the system (defined in a SALINITY section) [ppm]
# tidspd - maximum velocity of the tide (defined in a file in the TIDAL SPEED section) [m/s]
# valley - the normalise number of surrounding points uphill from the reference point
# (calclated from all grid points within two grid cells)
# -1 = all points downhill, 1 = all points uphill [-1,1]
# gradient - the average gradient of the point (calclated from all grid points within two grid cells)
# negative value indicates that most surrounding points are downhill from the point)
# grain1 - the percentage of the coarsest grain size in the top sediment layer (either basement or deposited material) [0-100]
# grain2 - the percentage of the second grain size in the top sediment layer (either basement or deposited material) [0-100]
# grain3 - the percentage of the third grain size in the top sediment layer (either basement or deposited material) [0-100]
# grain4 - the percentage of the finest grain size in the top sediment layer (either basement or deposited material) [0-100]
# carb1 - the percentage of the first carbonate type in the top sediment layer (either basement or deposited material) [0-100]
# carb2 - the percentage of the second carbonate type in the top sediment layer (either basement or deposited material) [0-100]
# org1 - the percentage of the first organic type in the top sediment layer (either basement or deposited material) [0-100]
# org2 - the percentage of the second organic type in the top sediment layer (either basement or deposited material) [0-100]
# burial - depth to which the material has been buried by subsequent deposits [meters] (affects material subsurface)
# age - time since the layer was deposited [years] (affects material subsurface)
# time - current time [years]
#
# xpoint ypoint - coordinates describing the shape of the membership function. ypoint must be between 0 and 1.
# see below for examples
#
# membership functions
# highexposure expo 0.0 0.0 1.0 1.0
# littlesed sedrate 0.0005 1.0 0.001 0.0
# vlittlesed sedrate 0.00005 1.0 0.0001 0.0
# shallow depth 0.0 0.0 0.0 1.0 15.0 1.0 25.0 0.0
# abovecarbcomp depth 0.0 0.0 0.0 1.0 2000.0 1.0 5000.0 0.0
# medtodepwat depth 20.0 0.0 50.0 1.0 65.0 1.0 105.0 0.0
# awayfromriver rivdist 10000.0 0.0 25000.0 1.0
# warmtohot temp 18.0 0.0 20.0 1.0 28.0 1.0 33.0 0.0
# temptowarm temp 12.0 0.0 14.0 1.0 20.0 1.0 22.0 0.0
# fastcurrents current 0.5 0.0 1.0 1.0
# normsalin salinity 0.0 0.0 3.3 1.0 3.5 1.0 3.6 0.0
# farfromshore shore 500.0 0.0 1000.0 1.0
# closetoreef carbdist 0.0 1.0 10.0 0.0
# overburden burial 0.0 0.0 0.1 1.0
# residence age 0 0.0 10 1.0
#
# slow_growth growth 0.0 0.0 5.0e-4 1.0
# v_slow_growth growth 0.0 0.0 4.0e-5 1.0
# quick_erosion erosion 0.0 0.0 4.0e-5 1.0
# slowdiagen porosity 0.2 1.0 1.0 0.0
# fastdiagen porosity 0.1 1.0 0.9 0.0
# hardening hardness 10 0.0 1000 1.0
# *
#
# the format for defining fuzzy rules is
# car/org_type = output_function_name = input_function_name1 and input_function_name2 and etc
# carb/org_type - must be one of "carbonate1", "carbonate2","organic1" or "organic2"
# output_function_name - must be one of the output functions defined earlier (with code_variable "growth")
# input_function_name - must be one of the functions defined earlier (without the code_variable "growth")
# rules
# carbonate1 = slow_growth = temptowarm and shallow
# carbonate1 = slowdiagen = longexposure
# carbonate1 = fastdiagen = overburden
#
# carbonate2 = v_slow_growth = warmtohot and shallow and closetoreef
# carbonate2 = slowdiagen = longexposure
# carbonate2 = fastdiagen = overburden
#
# siliciclast = fastdiagen = overburden
#
# organic2 = slow_growth = warmtohot and awayfromriver and vlittlesed and highexposure and shallow
# organic1 = slow_growth = temptowarm and awayfromriver and shallow
# *
#
#
# membership functions
shallow depth 0.0 0.0 0.0 1.0 15.0 1.0 25.0 0.0
medtodepwat depth 20.0 0.0 50.0 1.0 65.0 1.0 105.0 0.0
abovesealevel depth 0.0 1.0 0.0 0.0
warmtohot temp 18.0 0.0 20.0 1.0 28.0 1.0 33.0 0.0
temptowarm temp 12.0 0.0 14.0 1.0 20.0 1.0 22.0 0.0
longexposure expotime 1000.0 0.0 10000 1.0
closetoreef carbdist 0.0 1.0 10.0 0.0
awayfromriver rivdist 10000.0 0.0 25000.0 1.0
vlittlesed sedrate 0.00005 1.0 0.0001 0.0
highexposure expo 0.0 0.0 1.0 1.0
littleoverburden burial 5.0 0.0 20.0 1.0
overburden burial 20.0 0.0 120.0 1.0
oldage age 30 0.0 100 1.0
#
slow_growth growth 0.0 0.0 4.0e-4 1.0
v_slow_growth growth 0.0 0.0 2.0e-5 1.0
vv_slow_growth growth 0.0 0.0 1.0e-6 1.0
quickerosion erosion 0.0 0.0 3.0e-5 1.0
fastdiagen porosity 0.1 1.0 0.5 0.0
slowdiagen porosity 0.2 1.0 0.9 0.0
hardening hardness 10 0.0 1000 1.0
*
# rules
carbonate1 = slow_growth = temptowarm and shallow
carbonate1 = quickerosion = longexposure and abovesealevel
carbonate1 = slowdiagen = longexposure
carbonate1 = fastdiagen = littleoverburden and oldage
carbonate1 = hardening = longexposure
#
carbonate2 = v_slow_growth = warmtohot and medtodepwat and closetoreef
carbonate2 = quickerosion = longexposure and abovesealevel
carbonate2 = fastdiagen = littleoverburden and oldage
#
organic1 = vv_slow_growth =warmtohot and awayfromriver and vlittlesed and highexposure and shallow
organic1 = hardening = longexposure
#
siliciclast = fastdiagen = overburden and oldage
*
#-----------------------------------------------------------
#TIDAL SPEED
# this module is used in conjunction with organics and carbonate growth. (when tidspd is specified)
# the file contains the typical file speeds for each grid location
filename.tid
#-----------------------------------------------------------
#SALINITY
# this module is used in conjunction with organics and carbonate growth. (when salinity is specified)
# the file contains the typical slainity values for each grid location
filename.sal
#-----------------------------------------------------------
#TEMPERATURE GRID
# this module is used in conjunction with organics and carbonate growth. (when temperature is specified)
# the file contains the surface temperature values for each grid location
# when this file module is used, the temperature file specified in ORGANICS or CARBONATES is not used
filename.tpt
#-----------------------------------------------------------
#AEOLIAN
# Are there any sand sources ? (O/1 : no/yes)
1
# Define the sediment sources
# Definition of source
# One line per event, entries are:
# Event start and end time [sec]
# Source location (x,y) [m]
# Sediment height (H) [m]
# Sediment composition (C) [%]
#
# nb t1 t2 x y H %C %M %F %FF
1 0 10 -200 -199 2 0.0 5.0 15.0 80.0
*
# Define:
# - display wind time step,
# - transport parameter aliasing,
# - aeolian parametric sampling interval.
0.05 1 0.001
# Define a global wind regime for the entire area
#
# file containing time and wind velocity values for the region
#
# note that from 10 June 2009 the .wnd file has a preservation factor
# that is the fraction of a dune at any grid node that cannot be eroded
# by wind - this miics the effect of water-table preservation
# in wet aeolian systems
# 0.0 means all the dune is moved
# 1.0 all is preserved,
# 0.2 means 20% of the height of the dune at that grid node
# is preserved and the rest is moved.
#
# t1 t2 Vx(m/s) Vy(m/s) max_gust_velocity(m/s) preservation_factor(%) water_table_elevation(m)
windfile.wnd
#---------------------------------------------------------
#EXPERT AEOLIAN
# Is there some sources ? (O/1 : no/yes)
1
# Define the sediment sources
# Definition of source
# One line per event, entries are:
# Event start and end time [sec]
# Source location (x,y) [m]
# Sediment height (H) [m]
# Sediment composition (C) [%]
#
# nb t1 t2 x y H %C %M %F %FF
1 0 10 -200 -199 2 0.0 5.0 15.0 80.0
*
# Define:
# - display wind time step,
# - transport parameter aliasing,
# - aeolian parametric sampling interval.
0.05 1 0.001
# Wind measurement height Zu [m]
10
# Aerodynamic roughness length of the surface zo [m]
0.01
# Zero displacement height D [m]
0.03
# Vegetation cover growth function
# Vegetation cover updated time [year]
0.1
# number of points defining the curve
5
# define the curve for the considered vegetation type
# each point defined defined sedimentation balance and according annual concentration growth
-1.5 -1.0 0.1 0.0 0.6 0.4 2.0 0.0 2.2 -1.0
# Critical vegetation density over wich sand movement decreases sharply and maximum efficiency of suppression parameter
0.75 2.8
# input of initial vegetation cover (percent of cell area)
vegetation.dat
# This module allows the input of external wind data. the format for the
# input file is
# header line
# x velocities -one for each grid point in standard matrix format
# header line
# y velocity -one for each grid point in standard matrix format
# header line
# Frequency factor controls Aeolian calling interval.
# The bigger the factor the longer the interval
# astrix ends input list
# start_time, end_time, frequency factor, preservation_factor(%), water_table_elevation(m), input_file_name
-21000000.0 -20000000.0 10.0 0.0 -1000 wind_1.dat
-20000000.0 -19600000.0 10.0 0.2 140 wind_2.dat
*
##################################################################
# Optional parameters section. Defaults will be used if not entered
###############################################################
#--------------------------------------------------------
#ACCURACY FACTOR
# accuracy determines to what percentage accuracy the solution is solved to.
# importantly a small accuracy does not necessarily lead to longer runtimes
# as the solution stays more stable recommended=0.001 maximum=0.01
0.001
#------------------------------------------------------------------------
#BOUNDARY
#specify the boundary types for the coarse grid. (optional)
# 0 = regular boundary 1= insert a wall to prevent sediment and fluid loss
# without this section the default is regular boundaries
# North South East West
0 0 1 0
#--------------------------------------------------------
#DENSITIES
# Water densities (optional)
# fresh water 1000.0 kg/m3, sea water with 30000ppm NaCl 1027 kg/m3
# density of fluid entering the system [kg/m3] sea density [kg/m3]
1015.0 1027.0
#----------------------------------------------------------
#PARAMETRIC SAMPLING INTERVAL
# (optional)
#
# This should be seen in relation to temporal resolution - ie. think aliasing
# Sampling interval for sea-level, rainfall and tectonics [y]
1000.0
#----------------------------------------------------------
#SLOPE ANGLES
# Define maximum slope per grain size (optional)
#
#tg(1.0)=0.017; tg(2.0)=0.035
#
# a rough guide to appropriate slopes derived from deltas around the world:
# slope(deg)=0.4+0.5log(sand/shale ratio)
#
# Example slopes:
# a. based on sand/shale ratio for ratios 0.2-0.8 - use Sabesi equation :
# slope (gradient) = tan(0.4 + 0.52 log(sand/shale ratio))
# b. based on median grain size -
# from Dean(1983) quoted in Olsen 1990 p.38 (NTH Diplom)
# h(x) = Ax^m
# where :
# A =0.24 + 0.254*log(d50mm)
# m= 2/3
# Modified by Griffiths (1993) to agree more closely with
# Short(1979)'s observations as below.
# (Short noted that 0.5mm+ gave 5-6deg, 0.25-0.5mm gave between 1.5-5.5,
# and <0.25 mm gave less than 1.5deg)
# Therefore to give a slope in degrees from a grain size
# between 0.1 and 50 mm use
# Slope(deg) = atan(((0.28 + 0.23*log(d mm))*100^(2/3))/100)
#or
# Slope(gradient) = ((0.28+ 0.23*log(d mm))*100^(2/3))/100
#
# Max. slope of four grain sizes below sea level (dz/dx)
# Max. slope of four grain sizes above sea level (dz/dx)
# Max. slope of four reworked grain sizes below sea level (dz/dx)
0.007 0.0055 0.005 0.003
0.0001 0.0001 0.0001 0.0001
0.005 0.004 0.003 0.0005
# max slope carb/org grains below sea level (dz/dx) 4 grains
# max slope carb/org grains above sea level (dz/dx) 4 grains
0.005 0.004 0.005 0.004
0.0001 0.0001 0.0001 0.0001
# Minimum slope (dz/dx)
# Slope module calling interval [years]
0.0001
1000.0
#-----------------------------------------------------------
#SLOPE PARAMETERS
# Parameters for calculating equilibrium slope (optional)
#
# Number of diffusion cycles during one pass
# Maximum number of iterations per diffusion cycle
# Diffusion residual [m]
4 1000 0.0001
#--------------------------------------------------------
#SEDIMENT TRANSPORT PARAMETERS
# Limiting factors for sediment transport (optional)
#
# Sedimentation time step factor (sedimentation time step =
# flow time step * sedimentation time step factor)
# Maximum depth of fluid elements [m]
# Minimum velocity of fluid elements [m/s]
# Minimum ratio of sediment load fluid element to average sediment
# load at source [kg/m3]
# Basement hardness factor
1.0e4 1.0 0.05 0.10 1.0E2
#
#----------------------------------------------------------
#MANNING
# Manning's coefficients (optional)
#
# Coefficient for open-channel flow
# Coefficient for hyperpycnal flow
# Coefficient for hypopycnal flow
# Coefficient for debris flow
#default: 0.020 0.010 0.070 0.080
0.030 0.030 0.070 0.070
#--------------------------------------------------------------
#POROSITY TABLE
# Porosity function (optional)
#
# Number of entries in effective pressure look-up table
6
# Effective pressure look-up table [MPa]
0.0 10.0 20.0 30.0 40.0 50.0
# Number of entries in fine-to-coarse-ratio look-up table
12
# Fine-to-coarse-ratio look-up table
0.0 0.05 0.10 0.15 0.20 0.25 0.30 0.40 0.50 0.65 0.85 1.0
# Porosity look-up table
# must have the size of (pressures*ratios)
# rows: constant fine-to-coarse ratio
# columns: constant effective pressure
0.39 0.36 0.35 0.34 0.33 0.32
0.36 0.33 0.31 0.30 0.295 0.29
0.34 0.31 0.28 0.27 0.26 0.25
0.335 0.28 0.23 0.22 0.21 0.20
0.35 0.25 0.21 0.19 0.18 0.17
0.36 0.26 0.20 0.18 0.16 0.14
0.39 0.27 0.21 0.17 0.15 0.13
0.43 0.31 0.24 0.18 0.16 0.12
0.48 0.34 0.27 0.22 0.18 0.14
0.53 0.39 0.31 0.24 0.21 0.17
0.58 0.44 0.35 0.27 0.22 0.19
0.61 0.47 0.37 0.29 0.24 0.20
*
# Porosity table one row for each carbonate and organic
.70 .65 .50 .49 .48 .47
.30 .30 .30 .30 .30 .30
.60 .50 .40 .30 .20 .20
.20 .16 .13 .10 .05 .05
#
# Linear weighting coefficients for 4 grain sizes:
# r=Sum(h*w)/Sum(h)
# r fine-to-coarse ratio
# h thickness of individual grain size
# w linear weighting coefficient of individual grain size
# 4 coefficients(one for each siliciclastic grain):
0.0 0.0 1.0 1.0
#
#------------------------------------------------------------------------
#VISUALISATION
# Give the frequency at which you want to ouput the results : (optional)
# for example: -> 1 every time intervals
# -> 10 each 10s time interval
10
# Is this an Underworld simulation? (0-no ; 1-yes)
0
# *** End of example input file ***