Run a modal code simulation.

Parameters

[in,out] file Spec file

See Input File Format: Spec File Format for the input file text format.

A modal simulation spec file has the parameters:

run_type (string): must be modal
output_prefix (string): prefix of the output filenames — the filenames will be of the form PREFIX_SSSSSSSS.nc where SSSSSSSS is the eight-digit output index (starting at 1 and incremented each time the state is output)
t_max (real, unit s): total simulation time
del_t (real, unit s): timestep size
t_output (real, unit s): the interval on which to output data to disk (see Output File Format)
t_progress (real, unit s): the interval on which to write summary information to the screen while running
input_format_radius_bin_grid
do_camp_chem (logical): must be no (not currently supported)
do_tchem (logical): must be no (not currently supported)
gas_data (string): name of file from which to read the gas material data — the file format should be Input File Format: Gas Material Data
aerosol_data (string): name of file from which to read the aerosol material data (single species only) — the file format should be Input File Format: Aerosol Material Data
do_fractal (logical): whether to use fractal treatment (must be no; not currently supported)
aerosol_init (string): filename for the initial lognormal aerosol distribution — the file format should be Input File Format: Aerosol Distribution
Input File Format: Scenario
Input File Format: Environment State
do_coagulation (logical): must be no (not currently supported)
do_condensation (logical): must be no (not currently supported)
do_mosaic (logical): must be no (not currently supported)
do_optical (logical): must be no (not currently supported)
do_nucleation (logical): must be no (not currently supported)
do_immersion_freezing (logical): must be no (not currently supported)
do_parallel (logical): must be no (not currently supported)

Note: Modal runs currently support only the drydep and none loss functions. Aerosol emissions and background dilution are not supported; their rates must be zero in the scenario files. The options do_camp_chem, do_tchem, do_coagulation, do_condensation, do_mosaic, do_optical, do_nucleation, do_immersion_freezing, and do_parallel are included for future extensibility but must currently be set to no — setting any of them to yes will produce a fatal error.

Example:

 run_type modal           # modal run
 output_prefix out/modal  # prefix of output files

 t_max 28800    # total simulation time (s)
 del_t 60       # timestep (s)
 t_output 3600  # output interval (0 disables) (s)
 t_progress 0   # progress printing interval (0 disables) (s)

 n_bin 1000    # number of bins (for processing purposes)
 d_min 4e-8    # minimum diameter (m)
 d_max 2.5e-3  # maximum diameter (m)

 do_camp_chem no  # whether to use CAMP for chemistry (yes/no)
 do_tchem no      # whether to use TChem for chemistry (yes/no)

 gas_data gas_data.dat            # file containing gas data
 aerosol_data aero_data.dat       # file containing aerosol data
 do_fractal no                    # use fractal treatment (yes/no)
 aerosol_init aero_init_dist.dat  # aerosol initial condition file

 temp_profile temp.dat            # temperature profile file
 pressure_profile pres.dat        # pressure profile file
 height_profile height.dat        # height profile file
 gas_emissions gas_emit.dat       # gas emissions file
 gas_background gas_back.dat      # background gas concentrations file
 aero_emissions aero_emit.dat     # aerosol emissions file
 aero_background aero_back.dat    # aerosol background file
 loss_function drydep             # loss function specification
 drydep_params drydep_params.dat  # dry deposition parameters

 rel_humidity 0.95  # initial relative humidity (1)
 latitude 0         # latitude (degrees, -90 to 90)
 longitude 0        # longitude (degrees, -180 to 180)
 altitude 0         # altitude (m)
 start_time 21600   # start time (s since 00:00 UTC)
 start_day 200      # start day of year (UTC)

 do_coagulation no         # whether to do coagulation (yes/no)
 do_condensation no        # whether to do condensation (yes/no)
 do_mosaic no              # whether to do MOSAIC (yes/no)
 do_optical no             # whether to compute optical props (yes/no)
 do_nucleation no          # whether to do nucleation (yes/no)
 do_immersion_freezing no  # whether to do freezing (yes/no)
 do_parallel no            # whether to run in parallel (yes/no)