height_above_vent = [
30000, 20000, 10000, 5000, 2000, 1000
] # m (if point, if suzuki or if plume where Height_or_MFR = Height)
A = 4 # (suzuki only)
L = 5 # (suzuki only)
height_or_MFR = 'MFR' # plume only
MFR_minimum = 1e7 # kg/s (plume only)
MFR_maximum = 1e9 # kg/s (plume only)
exit_velocity = 100 # m/s (plume only)
exit_temperature = 1073 # K (plume only)
exit_volatile_fraction = 0 # % (plume only)
# Fall3D (Volcanological input file)
terminal_velocity_model = 'ganser' # Possibilites are ARASTOOPOR/GANSER/WILSON/DELLINO
vertical_turbulence_model = 'similarity' # Possibilites are CONSTANT/SIMILARITY
horizontal_turbulence_model = 'rams' # Possbilities are CONSTANT/RAMS
vertical_diffusion_coefficient = 100 # m2/s
horizontal_diffusion_coefficient = 1000 # m2/s
value_of_CS = 0.1 # RAMS only
# Contouring: True, False, number or list of numbers
thickness_contours = True
load_contours = True
thickness_units = 'cm' # mm/cm/m
# Run model using specified parameters
if __name__ == '__main__':
from aim import run_scenario
run_scenario(__file__)
Map_deposit_thickness = 'Yes' # Yes/No
Map_thickness_units = 'cm' # Possibilities (mm, cm, m)
Map_thickness_compaction_factor = 0.7 # Degree of compaction (i.e. 0.7)
Thickness_contours = [0.1, 1, 5, 10, 50, 100, 500] # List contour intervals (PS format only)
Map_total_concentration = 'No' # Yes/No
Map_total_concentration_z_cuts = [1000, 2000] # List height in meters of each z-cut required (i.e. 1000 2000)
Total_concentration_contours = [1e-5, 1e-4] # List contour intervals (PS format only)
Map_z_cummulative_concentration = 'No' # Yes/No
Cummulative_concentration_contours = [0.01, 0.1, 1, 10] # List contour intervals (PS format only)
Map_Z_maximum_concentration = 'No' # Yes/No
Maximum_concentration_contours = [1e-4, 1e-3] # List contour intervals (Ps format only)
Topography_grid = 'merapi_topography.txt' # Specify ASCII topography grid to use.
# If empty, AIM will look for a topography grid named
# <scenario_name>_topography.txt
# Run model using specified parameters
if __name__ == '__main__':
from aim import run_scenario
run_scenario(__file__, dircomment=Eruption_comment)
exit_volatile_fraction = 0 # % (plume only)
# Fall3D
terminal_velocity_model = 'ganser' # Possibilites are ARASTOOPOR/GANSER/WILSON/DELLINO
vertical_turbulence_model = 'similarity' # Possibilites are CONSTANT/SIMILARITY
horizontal_turbulence_model = 'rams' # Possbilities are CONSTANT/RAMS
vertical_diffusion_coefficient = 100 # m2/s
horizontal_diffusion_coefficient = 1000 # m2/s
value_of_CS = 0.1 # RAMS only
# Contouring: True, False, number or list of numbers
thickness_contours = True
load_contours = True
thickness_units = 'cm' # mm/cm/m
# Run model using specified parameters
if __name__ == '__main__':
from aim import run_scenario
run_scenario(__file__)
Topography_grid = 'merapi_topography.txt' # Specify ASCII topography grid to use.
# If empty, AIM will look for a topography grid named
# <scenario_name>.top (surfer GRD format)
# Contouring:
# False: Disabled
# True: Provide a fixed number of contours covering entire range
# Number: Fixed (vertical) interval between contours
# List of numbers: Exact contour levels
Thickness_contours = [1, 2, 5, 50, 100] # True, False, number or list of numbers
Thickness_units = 'cm' # mm/cm/m
Load_contours = 2000 # True, False, number or list of numbers
# Run model using specified parameters
if __name__ == '__main__':
import time, random
from aim import run_scenario
run_scenario(__file__,
timestamp_output=True,
dircomment=Eruption_comment)
# Fall3D (Volcanological input file)
terminal_velocity_model = 'ganser' # Possibilites are ARASTOOPOR/GANSER/WILSON/DELLINO
vertical_turbulence_model = 'similarity' # Possibilites are CONSTANT/SIMILARITY
horizontal_turbulence_model = 'rams' # Possbilities are CONSTANT/RAMS
vertical_diffusion_coefficient = 100 # m2/s
horizontal_diffusion_coefficient = 1000 # m2/s
value_of_CS = 0.1 # RAMS only
# Contouring: True, False, number or list of numbers
thickness_contours = [1, 2, 5, 10, 20, 50, 100] # Options: 'True', 'False' or [1, 2] (a list of contour values)
load_contours = [0.1, 90, 150, 300] # Options: 'True', 'False' or [1, 2] (a list of contour values)
thickness_units = 'cm' # mm/cm/m
# Run model using specified parameters
if __name__ == '__main__':
import time, random
from aim import run_scenario
run_scenario(__file__,
timestamp_output=True)
MFR_maximum = 1e9 # kg/s (plume only)
exit_velocity = 100 # m/s (plume only)
exit_temperature = 1073 # K (plume only)
exit_volatile_fraction = 0 # % (plume only)
# Fall3D (Volcanological input file)
terminal_velocity_model = 'ganser' # Possibilites are ARASTOOPOR/GANSER/WILSON/DELLINO
vertical_turbulence_model = 'similarity' # Possibilites are CONSTANT/SIMILARITY
horizontal_turbulence_model = 'rams' # Possbilities are CONSTANT/RAMS
vertical_diffusion_coefficient = 100 # m2/s
horizontal_diffusion_coefficient = 1000 # m2/s
value_of_CS = 0.1 # RAMS only
# Contouring: True, False, number or list of numbers
thickness_contours = [
1, 2, 5, 10, 20, 50, 100
] # Options: 'True', 'False' or [1, 2] (a list of contour values)
load_contours = [
0.1, 90, 150, 300
] # Options: 'True', 'False' or [1, 2] (a list of contour values)
thickness_units = 'cm' # mm/cm/m
# Run model using specified parameters
if __name__ == '__main__':
import time, random
from aim import run_scenario
run_scenario(__file__, timestamp_output=True)
Value_of_CS = 0.1 # RAMS only
# Output (Volcanological input file)
Postprocess_time_interval = 1 # Hours
Postprocess_3D_variables = 'No' # Yes/No
Postprocess_classes = 'No' # Yes/No
Track_points = 'No' # Yes/No
Topography_grid = 'merapi_topography.txt' # Specify ASCII topography grid to use.
# If empty, AIM will look for a topography grid named
# <scenario_name>.top (surfer GRD format)
# Contouring:
# False: Disabled
# True: Provide a fixed number of contours covering entire range
# Number: Fixed (vertical) interval between contours
# List of numbers: Exact contour levels
Thickness_contours = [1, 2, 5, 50,
100] # True, False, number or list of numbers
Thickness_units = 'cm' # mm/cm/m
Load_contours = 2000 # True, False, number or list of numbers
# Run model using specified parameters
if __name__ == '__main__':
import time, random
from aim import run_scenario
run_scenario(__file__, timestamp_output=True, dircomment=Eruption_comment)