Program structure (development version)

(not complete yet)

prodimo

    (Initialization)
    |
    +------ cpu_time
    |
    +------ init_nature
    |
    +------ read_parameter
    |
    +------ init_pah
    |                |
    |                +-------- pah_ipot
    |                |
    |                +-------- read_pah_resonances
    |                |
    |                +-------- setup_pah_lambda_grid
    |                |
    |                +-------- pah_cross_sections
    |
    +------ init_elements
    |
    +------ init_species
    |
    +------ init_reactions
    |                |
    |                +-------- read_umist_reactions
    |                |
    |                +-------- read_osu_reactions.f
    |                |
    |                +-------- read_KIDA_reactions.f
    |                |
    |                +-------- read_other_reactions3/read_other_reactions
    |
    +------ init_photoreactions
    |                |
    |                +-------- read_photoprocesses_data
    |
    +------ init_heatcool
    |                |
    |                +-------- read_nlte_lambda/read_nlte_woitke_ion
    |                |
    |                +-------- coll_rate_electrons
    |                |
    |                +-------- cobands
    |                |
    |                +-------- read_ntle_h2
    |                |
    |                |
    |                +-------- read_nlte_hd
    |                |
    |                +-------- stateq_it
    |                |
    |                +-------- write_nlte_lambda
    |                |
    |                +-------- read_hitran2009.f
    |                |
    |                +-------- read_nlte_water_hitran.f
    |                |
    |                +-------- read_jpl_format.f
    |
    +------ init_grid
    |
    +------ init_star_and_uv
    |
    +------ stellar_band_means
    |
    +------ stellar_band_adjust(.true.)
    |
    +------ init_xray_cs
    |
    +------ integrate_uv_cross_sections
    |
    +------ init_dust
    |                |
    |                +------ init_dust_opac (dust_opac.f)
    |                                  |
    |                                  +------ effMedium
    |                                  |                |
    |                                  |                +------ eff_func
    |                                  |                |
    |                                  |                +------ gauss
    |                                  |                |
    |                                  |                +------ eff_pullback2
    |                                  |                |
    |                                  |                +------ e2m
    |                                  |                |
    |                                  |                +------ m2e
    |                                  |
    |                                  +------ shexqnn2
    |                                  |
    |                                  +------ dmilay
    |                                  |
    |                                  +------ dust_opac
    |
    |
    +------ init_dust_photoelectric
    |
    +------ init_line_transfer
    |
    +------ init_disk
    |                |
    |                +------ init_fixed_surface_density
    |                |
    |                +------ disk_structure
    |                |                |
    |                |                +------ fixed_disk_structure
    |                |                |               |
    |                |                |               +------ disk_structure
    |                |                |                                 |
    |                |                |                                 +------ fixed_surface_density
    |                |                |                                 |
    |                |                |                                 +------ coldens
    |                |                |                                 |
    |                |                |                                 +------ disk_structure
    |                |                |
    |                |                |
    |                |                +------ coldens
    |                |                |
    |                |                +------ disk_structure
    |                |
    |                |
    |                +------ dust_structure
    |                                  |
    |                                  +------ dust_settled_structure
    |
    |
    |
    +------ init_concentrations
    |
    +------ init_from_mcfost (read_MCFOST.f90)
    |                                  |
    |                                  +------  read_mcfost
    |                                  |                 |
    |                                  |                 +------ cfitsio routines
    |                                  |
    |                                  +------ init_grid
    |                                  |
    |                                  +------ init_disk
    |                                  |
    |                                  +------ init_star_and_uv
    |                                  |
    |                                  +------ stellar_band_means
    |                                  |
    |                                  +------ integrate_uv_cross_sections
    |                                  |
    |                                  +------ incident_intensities
    |                                  |
    |                                  +------ stellar_band_adjust
    |
    +------ init_restart (restart.f)
    |                                  |
    |                                  +------ disk_structure
    |                                  |
    |                                  +------ dust_structure
    |
    +------ time_dependent_mc
    |                                  |
    |                                  +------ rate_coeffs
    |                                  |
    |                                  +------ advance_chemistry
    |


    (time dependent chemistry model, no convergence loop)
    +------ time_dependent_disk
    |                                  |
    |                                  +------ incident_intensities
    |                                  |
    |                                  +------ cpu_time
    |                                  |
    |                                  +------ tau_interpol
    |                                  |
    |                                  +------ irback
    |                                  |
    |                                  +------ heatcoolgas_it
    |                                  |
    |                                  +------ update_optical_depths
    |                                  |
    |                                  +------ update_fline
    |                                  |
    |                                  +------ ouput
    |                                                    |
    |                                                    +------ tau_interpol
    |                                                    |
    |                                                    +------ output_fline_estimate (line_flux_estimate.f)
    |                                                    |
    |                                                    +------ write_mcfost

    (disk structure iteration loop)
    |
    +------ cpu_time(t1)
    |
    +------ radiative_transfer
    |                |
    |                +------ optical_depths
    |                |
    |                +------ formal_solution
    |                |                 |
    |                |                 +------ incident_intensities
    |                |                 |
    |                |                 +------ estimate_mean_it
    |                |                 |
    |                |                 +------ output (benchmark_rt)
    |                |                 |
    |                |                 +------ cpu_time
    |                |                 |
    |                |                 +------ integrate_ray
    |                |                 |
    |                |                 +------ adjust_theta
    |                |                 |
    |                |                 +------ extrap
    |                |
    |                +------ estimate_mean_it
    |                |                   |
    |                |                   +------ incident_intensities
    |                |                                     |
    |                |                                     +------ stellar_band_means
    |                |                                     |
    |                |                                     +------ i_to_chi
    |                |                                     |
    |                |                                     +------ j_to_chi
    |                |
    |                +------ calculate_sed_mono
    |                |                   |
    |                |                   +------ kappa_interpol
    |                |                   |
    |                |                   +------ integrate_ray_sed
    |                |
    |                +------ calculate_sed2_mono
    |                |                   |
    |                |                   +------ cpu_time
    |                |                   |
    |                |                   +------ whereami
    |                |                   |
    |                |                   +------ integrate_ray_sed
    |                |
    |                +------ calculate_sed
    |                |                   |
    |                |                   +------ incident_intensities
    |                |                   |                 |
    |                |                   |                 +------ stellar_band_means
    |                |                   |                 |
    |                |                   |                 +------ i_to_chi
    |                |                   |                 |
    |                |                   |                 +------ j_to_chi
    |                |                   |
    |                |                   +------ integrate_ray
    |                |
    |                +------ calculate_sed2
    |                |                   |
    |                |                   +------ cpu_time
    |                |                   |
    |                |                   +------ wherami
    |                |                   |
    |                |                   +------ integrate_ray
    |
    +------ chemistry_and_energy_balance
    |                |
    |                +------ init_fline_estimate (line_flux_estimate.f)
    |                |
    |                +------ incident_intensities
    |                |                   |
    |                |                   +------ stellar_band_means
    |                |                   |
    |                |                   +------ i_to_chi
    |                |                   |
    |                |                   +------ j_to_chi
    |                |
    |                +------ cpu_time
    |                |
    |                +------ guess_tgas
    |                |
    |                +------ tau_interpol
    |                |
    |                +------ xtau
    |                |
    |                +------ xhx
    |                |
    |                +------ xion
    |                |
    |                +------ xdepos
    |                |
    |                +------ xyield
    |                |
    |                +------ irback
    |                |
    |                +------ point_chemistry_and_balance: for each grid point, use the Riddler's method to find the gas temperature
    |                |                   |
    |                |                   +------ determine_chemistry
    |                |                                      |
    |                |                                      +------ rate_coeffs
    |                |                                      |                 |
    |                |                                      |                 +------ ice_thermal_desorb
    |                |                                      |                 |
    |                |                                      |                 +------ ice_cr_desorb
    |                |                                      |                 |
    |                |                                      |                 +------ ice_photo_desorb_yield
    |                |                                      |                 |
    |                |                                      |                 +------ get_photorate
    |                |                                      |                 |
    |                |                                      |                 +------compute_photorate
    |                |                                      |                 |
    |                |                                      |                 +------ compute_pah_photoionization(pah_physics.f)
    |                |                                      |                 |                 |
    |                |                                      |                 |                 +------ pah_yield
    |                |                                      |                 |                 |
    |                |                                      |                 |                 +------ jnu_from_rt
    |                |                                      |                 |
    |                |                                      |                 +------ codis
    |                |                                      |                 |
    |                |                                      |                 +------ effh2 (rate_coeffs.f)
    |                |                                      |
    |                |                                      +------ guess_chemistry
    |                |                                      |
    |                |                                      +------ solve_chemistry
    |                |                                      |                 |
    |                |                                      |                 +------ chemfunc
    |                |                                      |                 |            |
    |                |                                      |                 |            +------ xyield (if xrays are on)
    |                |                                      |                 |            |
    |                |                                      |                 |            +------ xhx (if xrays are on)
    |                |                                      |                 |                     |
    |                |                                      |                 |                     +------ xcross
    |                |                                      |                 |
    |                |                                      |                 |
    |                |                                      |                 +------ chemjacobi
    |                |                                      |                 |            |
    |                |                                      |                 |            +------ xpar (if xrays are on)
    |                |                                      |                 |
    |                |                                      |                 |
    |                |                                      |                 +------ prepare_iteration
    |                |                                      |                 |
    |                |                                      |                 +------ element_conserve
    |                |                                      |                 |
    |                |                                      |                 +------ sgeir (linpack, slatec_routines.f)
    |                |                                      |                 |
    |                |                                      |                 +------ pullback (in solve_chemistry.f)
    |                |                                      |                 |
    |                |                                      |                 +------ chemanalyse
    |                |                                      |                 |
    |                |                                      |                 +------ chemqual
    |                |                                      |                 |
    |                |                                      |                 +------ is_nan
    |                |                                      |
    |                |                                      +------ advance_chemistry
    |                |                                      |                 |
    |                |                                      |                 +------ chem_rhs
    |                |                                      |                 |
    |                |                                      |                 +------ cpu_time
    |                |                                      |                 |
    |                |                                      |                 +------ limex
    |                |                                      |                 |
    |                |                                      |                 +------ chem_fcn
    |                |                                      |                 |
    |                |                                      |                 +------ chem_jac
    |                |                                      |                 |
    |                |                                      |                 +------ chemjacobi
    |                |                                      |
    |                |                                      +------ heatcool_gas_it
    |                |                                                         |
    |                |                                                         +------ compute_pah_photoionization
    |                |                                                         |
    |                |                                                         +------ stateq_it
    |                |                                                         |                 |
    |                |                                                         |                 +------ h2_ortho_para
    |                |                                                         |                 |
    |                |                                                         |                 +------ (ludcmp/lubskb/mprove)
    |                |                                                         |                 |
    |                |                                                         |                 +------ sgeir
    |                |                                                         |                 |
    |                |                                                         |                 +------ nlevel_heatcool_it
    |                |                                                         |                                   |
    |                |                                                         |                                   +------ escape
    |                |                                                         |                                   |
    |                |                                                         |                                   +------ pump
    |                |                                                         |                                   |
    |                |                                                         |                                   +------ bplanck
    |                |                                                         |
    |                |                                                         +------ nlevel_heatcool_it
    |                |                                                         |                 |
    |                |                                                         |                 +------ escape
    |                |                                                         |                 |
    |                |                                                         |                 +------ pump
    |                |                                                         |                 |
    |                |                                                         |                 +------ bplanck
    |                |                                                         |
    |                |                                                         +------ qfreifrei
    |                |                                                         |
    |                |                                                         +------ heatcool_bf
    |                |
    |                +------ update_optical_depths
    |                |                |
    |                |                +------ tauline_it
    |                |                                  |
    |                |                                  +------ h2_ortho_para
    |                |                                  |
    |                |                                  +------ stateq_it
    |                |                                                         |
    |                |                                                         +------ (see above)
    |                |
    |                +------ update_fline (line_flux_estimate.f)
    |                |                |
    |                |                +------ h2_ortho_apra
    |                |                |
    |                |                +------ tau_interpol
    |                |
    |                +------cpu_time
    |
    +------ line_radiative_transfer
    |                |
    |                +------cpu_time
    |                |
    |                +------line_precalculate
    |                |
    |                +------whereami
    |                |
    |                +------integrate_line
    |
    +------ soundspeed
    |
    +------ disk_structure
    |                |
    |                +------fixed_disk_structure
    |                |
    |                +------coldens
    |                |
    |                +------disk_rhs
    |                |
    |                +------limex (disk_fcn, disk_jacobi)

    |
    +------ dust_structure
    |                |
    |                +------dust_settled_structure (dust_settle.f)
    |                |
    |                +------dust_opac
    |
    +------ output
    |                |
    |                +------ tau_interpol
    |                |
    |                +------ output_fline_estimate (line_flux_estimate.f)
    |                |
    |                +------ write_mcfost
    |
    +------ cpu_time(t2)