1D interface

TODO: check if up-to-date, also mention special 1D interfaces. Also various parameters can be set via the 1D interface.

It is also possible to use ProDiMo with prescribed radial surface density profiles for the gas and the dust including the grain size distribution. One has to provide those values as a function of radius in a file.

To create the surface density input file we highly recommend the utility routines of prodimopy (see interface1D)

To use the 1D interface set the following Parameters in the Parameter.in

.true.        ! interface1D          activate the 1D interface
sdprofile.in  ! surf_density_file    the name of the file with the surface densities

The sdprofile.in defines the initial surface density profile for the dust and the gas and if requested the dust size distribution function. The vertical structure is still controlled via the settings in Parameter.in (like in a normal model). Also, dust settling etc. is still possible.

The given surface density profile also determines the total disk mass, the inner and outer radius, the numerical radial grid and if requested the grain size grid. (i.e. the user is responsible for that). As a consequence parameters like the disk mass, inner/outer radius or grain size distributions from Parameter.in are ignored (overwritten, check the output/log file of your ProDiMo run)

Example ProDiMo models with user-defined surface density profiles can be found in examples/DiskStructure of your ProDiMo installation.

References: Rab+ 2020 (paper that uses 1D input profiles).

Old DEPRECATED version

TODO: Maybe just remove it

WARNING: Should not be used anymore

There is still an older version included in ProDiMo that uses a simple input file but also has less functionality.

To use the 1D interface set the following Parameters in the Parameter.in

.true.        ! fixed_surface_density   activate the 1D interface
sdprofile.in  ! surf_density_file       the name of the file with the surface densities

The old format for the sdprofile.in has to look like this :

# number of radii
250
# Rin(au)
0.410652
# Rout(au)
605.4412
# normalization factor
-2.0
#    R(cm)      Profile    gas2dust
 6.14326648e+12  2.42569800e+02  8.90598782e+04
 6.14326798e+12  2.42569852e+02  8.90598782e+04
 6.14328144e+12  2.42569158e+02  8.90598780e+04
 6.14338017e+12  2.42564568e+02  8.90598930e+04
 6.14411620e+12  2.42532473e+02  8.90599065e+04
 6.14962589e+12  2.42291200e+02  8.90600724e+04
 6.19082065e+12  2.40501795e+02  8.90598593e+04
 6.27872436e+12  2.36765549e+02  8.90596087e+04
 6.39441139e+12  2.32014139e+02  8.90598274e+04
 6.53416871e+12  2.26509700e+02  8.90598127e+04

....