dustdensity

Module

Description

$Id$

This module takes care of everything related to dust density.

** AUTOMATIC CPARAM.INC GENERATION ************************ Declare (for generation of cparam.inc) the number of f array variables and auxiliary variables added by this module

CPARAM logical, parameter :: ldustdensity = .true.

MVAR CONTRIBUTION 1 MAUX CONTRIBUTION 0

PENCILS PROVIDED glnnd(3,ndustspec); gmi(3,ndustspec); gmd(3,ndustspec) PENCILS PROVIDED gnd(3,ndustspec); grhod(3,ndustspec) PENCILS PROVIDED ad(ndustspec); md(ndustspec); mi(ndustspec); nd(ndustspec) PENCILS PROVIDED rhod(ndustspec); rhod1(ndustspec); epsd(ndustspec) PENCILS PROVIDED udgmi(ndustspec); udgmd(ndustspec); udglnnd(ndustspec) PENCILS PROVIDED udgnd(ndustspec); glnnd2(ndustspec) PENCILS PROVIDED sdglnnd(3,ndustspec); del2nd(ndustspec); del2rhod(ndustspec) PENCILS PROVIDED del6nd(ndustspec); del2md(ndustspec) PENCILS PROVIDED del2mi(ndustspec); del6lnnd(ndustspec) PENCILS PROVIDED gndglnrho(ndustspec); glnndglnrho(ndustspec) PENCILS PROVIDED udrop(3,ndustspec); udropgnd(ndustspec) PENCILS PROVIDED fcloud; ccondens; ppwater; ppsat PENCILS PROVIDED ppsf(ndustspec); mu1; udropav(3) PENCILS PROVIDED glnrhod(3,ndustspec); PENCILS PROVIDED rhodsum; rhodsum1; grhodsum(3); glnrhodsum(3)


Quick access

Variables:

adpeak, adref_diffnd, advec_ddensity, amplnd, amplnd_rel, bordernd, coag_kernel, coeff_smooth, copy_bcs_dust_short, deltavd_const, deltavd_imposed, deriv_size, diffmd, diffmi, diffnd, diffnd_anisotropic, diffnd_exponent, diffnd_hyper3, diffnd_hyper3_mesh, diffnd_law, diffnd_ndustspec, diffnd_shock, dkern, dkern_cst, dlnad, dlnmd, dndfac_sum, dndfac_sum2, droplet_init, droplet_redistr, dt_substep, dust_coagulation, dust_condensation, dust_condensation_lmdvar, dust_condensation_nolmdvar, dustdensity_floor, dustdensity_floor_log, enum_bordernd, enum_self_collisions, epsz1_smooth, f_lucky, g_condensparam, get_ccondens, get_deltavd_turbu, get_mfluxcond, gs_condensparam, gs_condensparam0, hepsd, hnd, iadvec_ddensity, idiag_adm, idiag_epsdm, idiag_epsdmax, idiag_epsdmin, idiag_epsdrms, idiag_kk2m, idiag_kkm, idiag_mdm, idiag_mdmtot, idiag_mmxm, idiag_mmym, idiag_mmzm, idiag_nd2m, idiag_ndm, idiag_ndmax, idiag_ndmin, idiag_ndmt, idiag_ndmx, idiag_ndmxy, idiag_ndmz, idiag_rhodm, idiag_rhodmax, idiag_rhodmin, idiag_rhodmt, idiag_rhodmxy, idiag_rhodmz, idiag_rhoimt, idiag_rmom, idiag_ssrm, idiag_ssrmax, idiffd, iglobal_nd, init_distr_ki, initnd, initnd_lognormal, kern_max, kern_min, kernel_mean, kx_nd, ky_nd, kz_nd, latm_chemistry, lcalcdkern, ldeltavd_thermal, ldeltavd_turbulent, ldiffd_dusttogasratio, ldiffd_hyper3, ldiffd_hyper3_mesh, ldiffd_hyper3_polar, ldiffd_hyper3lnnd, ldiffd_shock, ldiffd_simpl_anisotropic, ldiffd_simplified, ldust_cdtc, ldustcoagulation_simplified, ldustcondensation_simplified, ldustcontinuity, ldustnucleation, ldustnulling, lkeepinitnd, lkernel_mean, lmice, lmomcons, lmomcons2, lmomcons3, lmomcons3b, lmomconsb, lno_deltavd, lnoaerosol, lnocondens_term, lpiecewise_constant_kernel, lradius_binning, lresetuniform_dustdensity, lself_collisions, lsemi_chemistry, lsubstep, ludstickmax, lupw_ndmdmi, lzero_upper_kern, mdave0, mi, momcons_sum_x, momcons_sum_y, momcons_sum_z, momcons_term_frac, nd0_luck, nd_const, nd_reuni, ndiffd_max, ndmin_for_mdvar, null_dust_vars, phase_nd, ppsf_full, r_collected, r_lucky, radius_nd, redist_mdbins, reinitialize_nd, self_collision_factor, self_collisions, set_border_dustdensity, sigmad, supsatfac, supsatratio_given, supsatratio_given0, supsatratio_omega, teta, teta1, tl0, tl01, ueta, ul0, widthnd, xblob_nd, yblob_nd, z0_smooth, z1_smooth, zblob_nd

Routines:

calc_diagnostics_dustdensity(), calc_pencils_dustdensity(), dndmd_dt(), dustdensity_after_boundary(), dustdensity_before_boundary(), get_slices_dustdensity(), impose_dustdensity_floor(), init_nd(), initialize_dustdensity(), pencil_criteria_dustdensity(), pencil_interdep_dustdensity(), read_dustdensity_init_pars(), read_dustdensity_run_pars(), register_dustdensity(), rprint_dustdensity(), write_dustdensity_init_pars(), write_dustdensity_run_pars()

Needed modules

Variables

  • dustdensity/a0 [real,private/optional/default=0.0]
  • dustdensity/a1 [real,private/optional/default=0.0]
  • dustdensity/aa [real,private/optional/default=6.6e-05]
  • dustdensity/adpeak [real,private/optional/default=0.0005]
  • dustdensity/adref_diffnd [real,private/optional/default=0.0]
  • dustdensity/advec_ddensity [character,private/optional/default='normal']
  • dustdensity/amplnd [real,private/optional/default=1.0]
  • dustdensity/amplnd_rel (ndustspec) [real,private/optional/default=0.]
  • dustdensity/beta_glnrho_scaled (*) [real,private/pointer]
  • dustdensity/bordernd [character,private/optional/default='nothing']
  • dustdensity/coag_kernel [private]
  • dustdensity/coeff_smooth (6) [real,private/optional/default=0.0]
  • dustdensity/constant_richardson [private]
  • dustdensity/copy_bcs_dust_short [private]
  • dustdensity/deltamd [real,private/pointer]
  • dustdensity/deltavd_const [real,private/optional/default=1.0]
  • dustdensity/deltavd_imposed [real,private/optional/default=0.0]
  • dustdensity/deriv_size [private]
  • dustdensity/diffmd [real,private/optional/default=0.0]
  • dustdensity/diffmi [real,private/optional/default=0.0]
  • dustdensity/diffnd [real,private/optional/default=0.0]
  • dustdensity/diffnd_anisotropic (3) [real,private/optional/default=0.0]
  • dustdensity/diffnd_exponent [real,private/optional/default=0.0]
  • dustdensity/diffnd_hyper3 [real,private/optional/default=0.0]
  • dustdensity/diffnd_hyper3_mesh [real,private/optional/default=5.0]
  • dustdensity/diffnd_law [character,private/optional/default='const']
  • dustdensity/diffnd_ndustspec (ndustspec) [real,private]
  • dustdensity/diffnd_shock [real,private/optional/default=0.0]
  • dustdensity/dkern (nx,ndustspec,ndustspec) [real,private]
  • dustdensity/dkern_cst [real,private/optional/default=0.0]
  • dustdensity/dlnad [real,private]
  • dustdensity/dlnmd [real,private]
  • dustdensity/dndfac_sum [real,private]
  • dustdensity/dndfac_sum2 [real,private]
  • dustdensity/droplet_init [private]
  • dustdensity/droplet_redistr [private]
  • dustdensity/dsize (ndustspec) [real,private]
  • dustdensity/dt_substep [real,private/optional/default=2e-07]
  • dustdensity/dust_coagulation [private]
  • dustdensity/dust_condensation [private]
  • dustdensity/dust_condensation_lmdvar [private]
  • dustdensity/dust_condensation_nolmdvar [private]
  • dustdensity/dustbin_width [real,private/pointer]
  • dustdensity/dustdensity_floor [real,private/optional/default=-1]
  • dustdensity/dustdensity_floor_log [real,private]
  • dustdensity/dwater [real,private/optional/default=0.220784]
  • dustdensity/enum_bordernd [integer,private/optional/default=0]
  • dustdensity/enum_self_collisions [integer,private/optional/default=0]
  • dustdensity/eps1 [real,private/optional/default=0.5]
  • dustdensity/eps_dtog [real,private/optional/default=0.0]
  • dustdensity/epsz1_smooth [real,private/optional/default=0.0]
  • dustdensity/f_lucky [real,private/optional/default=0.0]
  • dustdensity/g_condensparam [real,private/optional/default=0.0]
  • dustdensity/gamma [real,private]
  • dustdensity/get_ccondens [private]
  • dustdensity/get_deltavd_turbu [private]
  • dustdensity/get_mfluxcond [private]
  • dustdensity/gs_condensparam [real,private]
  • dustdensity/gs_condensparam0 [real,private]
  • dustdensity/hepsd [real,private/optional/default=1.0]
  • dustdensity/hnd [real,private/optional/default=1.0]
  • dustdensity/iadvec_ddensity [integer,private/optional/default=0]
  • dustdensity/idiag_adm [integer,private/optional/default=0]
  • dustdensity/idiag_admom (25) [integer,private/optional/default=0]
  • dustdensity/idiag_divud2m (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_epsdm (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_epsdmax (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_epsdmin (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_epsdrms (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_kk2m [integer,private/optional/default=0]
  • dustdensity/idiag_kkm [integer,private/optional/default=0]
  • dustdensity/idiag_mdm (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_mdmtot [integer,private/optional/default=0]
  • dustdensity/idiag_mmxm [integer,private/optional/default=0]
  • dustdensity/idiag_mmym [integer,private/optional/default=0]
  • dustdensity/idiag_mmzm [integer,private/optional/default=0]
  • dustdensity/idiag_nd2m (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_ndm (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_ndmax (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_ndmin (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_ndmt [integer,private/optional/default=0]
  • dustdensity/idiag_ndmx (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_ndmxy [integer,private/optional/default=0]
  • dustdensity/idiag_ndmz (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_rhodm (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_rhodmax (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_rhodmin (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_rhodmt [integer,private/optional/default=0]
  • dustdensity/idiag_rhodmxy [integer,private/optional/default=0]
  • dustdensity/idiag_rhodmz (ndustspec) [integer,private/optional/default=0]
  • dustdensity/idiag_rhoimt [integer,private/optional/default=0]
  • dustdensity/idiag_rmom (25) [integer,private/optional/default=0]
  • dustdensity/idiag_ssrm [integer,private/optional/default=0]
  • dustdensity/idiag_ssrmax [integer,private/optional/default=0]
  • dustdensity/idiffd (4) [character,private/optional/default='']
  • dustdensity/iglobal_nd [integer,private/optional/default=0]
  • dustdensity/ilnrho_cs2 [integer,private/parameter/optional/default=5]
  • dustdensity/ilnrho_ee [integer,private/parameter/optional/default=2]
  • dustdensity/ilnrho_eth [integer,private/parameter/optional/default=14]
  • dustdensity/ilnrho_lntt [integer,private/parameter/optional/default=4]
  • dustdensity/ilnrho_pp [integer,private/parameter/optional/default=3]
  • dustdensity/ilnrho_ss [integer,private/parameter/optional/default=1]
  • dustdensity/ilnrho_tt [integer,private/parameter/optional/default=9]
  • dustdensity/init_distr (mx,ndustspec) [real,private]
  • dustdensity/init_distr2 (ndustspec) [real,private]
  • dustdensity/init_distr_ki (ndustspec,ndustspec0) [real,private]
  • dustdensity/init_x1 [real,private/optional/default=0.0]
  • dustdensity/init_x2 [real,private/optional/default=0.0]
  • dustdensity/initnd (ninit) [character,private/optional/default='nothing']
  • dustdensity/initnd_lognormal [private]
  • dustdensity/ipp_cs2 [integer,private/parameter/optional/default=12]
  • dustdensity/ipp_ss [integer,private/parameter/optional/default=11]
  • dustdensity/irho_cs2 [integer,private/parameter/optional/default=6]
  • dustdensity/irho_ee [integer,private/parameter/optional/default=15]
  • dustdensity/irho_eth [integer,private/parameter/optional/default=13]
  • dustdensity/irho_lntt [integer,private/parameter/optional/default=8]
  • dustdensity/irho_pp [integer,private/parameter/optional/default=16]
  • dustdensity/irho_ss [integer,private/parameter/optional/default=7]
  • dustdensity/irho_tt [integer,private/parameter/optional/default=10]
  • dustdensity/kern_max [real,private/optional/default=0.0]
  • dustdensity/kern_min [real,private/optional/default=0.0]
  • dustdensity/kernel_mean (ndustspec,ndustspec) [real,private]
  • dustdensity/kx_nd [real,private/optional/default=1.0]
  • dustdensity/ky_nd [real,private/optional/default=1.0]
  • dustdensity/kz_nd [real,private/optional/default=1.0]
  • dustdensity/latm_chemistry [logical,private/optional/default=.false.]
  • dustdensity/lcalcdkern [logical,private/optional/default=.true.]
  • dustdensity/lcondensing_species [logical,private/optional/default=.false.]
  • dustdensity/ldeltavd_thermal [logical,private/optional/default=.false.]
  • dustdensity/ldeltavd_turbulent [logical,private/optional/default=.false.]
  • dustdensity/ldiffd_dusttogasratio [logical,private/optional/default=.false.]
  • dustdensity/ldiffd_hyper3 [logical,private/optional/default=.false.]
  • dustdensity/ldiffd_hyper3_mesh [logical,private/optional/default=.false.]
  • dustdensity/ldiffd_hyper3_polar [logical,private/optional/default=.false.]
  • dustdensity/ldiffd_hyper3lnnd [logical,private/optional/default=.false.]
  • dustdensity/ldiffd_shock [logical,private/optional/default=.false.]
  • dustdensity/ldiffd_simpl_anisotropic [logical,private/optional/default=.false.]
  • dustdensity/ldiffd_simplified [logical,private/optional/default=.false.]
  • dustdensity/ldust_cdtc [logical,private/optional/default=.false.]
  • dustdensity/ldustcoagulation_simplified [logical,private/optional/default=.false.]
  • dustdensity/ldustcondensation_simplified [logical,private/optional/default=.false.]
  • dustdensity/ldustcontinuity [logical,private/optional/default=.true.]
  • dustdensity/ldustnucleation [logical,private/optional/default=.false.]
  • dustdensity/ldustnulling [logical,private/optional/default=.false.]
  • dustdensity/lfree_molecule [logical,private/optional/default=.false.]
  • dustdensity/lkeepinitnd [logical,private/optional/default=.false.]
  • dustdensity/lkernel_mean [logical,private/optional/default=.false.]
  • dustdensity/llin_radiusbins [logical,private/pointer]
  • dustdensity/llog10_for_admom_above10 [logical,private/optional/default=.true.]
  • dustdensity/llog_massbins [logical,private/pointer]
  • dustdensity/lmice [logical,private/optional/default=.false.]
  • dustdensity/lmomcons [logical,private/optional/default=.false.]
  • dustdensity/lmomcons2 [logical,private/optional/default=.false.]
  • dustdensity/lmomcons3 [logical,private/optional/default=.false.]
  • dustdensity/lmomcons3b [logical,private/optional/default=.false.]
  • dustdensity/lmomconsb [logical,private/optional/default=.false.]
  • dustdensity/lno_deltavd [logical,private/optional/default=.false.]
  • dustdensity/lnoaerosol [logical,private/optional/default=.false.]
  • dustdensity/lnocondens_term [logical,private/optional/default=.false.]
  • dustdensity/lpiecewise_constant_kernel [logical,private/optional/default=.false.]
  • dustdensity/lradius_binning [logical,private/optional/default=.false.]
  • dustdensity/lresetuniform_dustdensity [logical,private/optional/default=.false.]
  • dustdensity/lself_collisions [logical,private/optional/default=.false.]
  • dustdensity/lsemi_chemistry [logical,private/optional/default=.false.]
  • dustdensity/lsubstep [logical,private/optional/default=.false.]
  • dustdensity/ludstickmax [logical,private/optional/default=.false.]
  • dustdensity/lupw_ndmdmi [logical,private/optional/default=.false.]
  • dustdensity/lzero_upper_kern [logical,private/optional/default=.false.]
  • dustdensity/m_w [real,private/optional/default=18.0]
  • dustdensity/mdave0 [real,private/optional/default=1.0]
  • dustdensity/mi (ndustspec) [real,private]
  • dustdensity/mmom [integer,private/parameter/optional/default=24]
  • dustdensity/momcons_sum_x [real,private]
  • dustdensity/momcons_sum_y [real,private]
  • dustdensity/momcons_sum_z [real,private]
  • dustdensity/momcons_term_frac [real,private/optional/default=1.0]
  • dustdensity/nd0_luck [real,private/optional/default=0.0]
  • dustdensity/nd_const [real,private/optional/default=1.0]
  • dustdensity/nd_reuni [real,private/optional/default=0.0]
  • dustdensity/ndiffd_max [integer,private/parameter/optional/default=4]
  • dustdensity/ndmin_for_mdvar [real,private/optional/default=0.0]
  • dustdensity/ntot [real,private/optional/default=1.0]
  • dustdensity/null_dust_vars [private]
  • dustdensity/phase_nd [real,private/optional/default=0.0]
  • dustdensity/ppsf_full (nx,ndustspec,ndustspec0) [real,private]
  • dustdensity/r_collected [real,private/optional/default=0.0]
  • dustdensity/r_lucky [real,private/optional/default=0.0]
  • dustdensity/radius (ndustspec) [real,private]
  • dustdensity/radius_nd [real,private/optional/default=1.0]
  • dustdensity/redist_mdbins [private]
  • dustdensity/reinitialize_nd [logical,private/optional/default=.false.]
  • dustdensity/rgas [real,private/optional/default=83100000.0]
  • dustdensity/rgas_unit_sys [real,private]
  • dustdensity/rho_w [real,private/optional/default=1.0]
  • dustdensity/ri0 [real,private/optional/default=1.0]
  • dustdensity/rotat_position [real,private/optional/default=0.0]
  • dustdensity/self_collision_factor [real,private/optional/default=1.0]
  • dustdensity/self_collisions [character,private/optional/default='nothing']
  • dustdensity/set_border_dustdensity [private]
  • dustdensity/sigmad [real,private/optional/default=1.0]
  • dustdensity/spot_number [integer,private/optional/default=1]
  • dustdensity/supsatfac [real,private/optional/default=1.0]
  • dustdensity/supsatratio_given [real,private/optional/default=0.0]
  • dustdensity/supsatratio_given0 [real,private/optional/default=0.0]
  • dustdensity/supsatratio_omega [real,private/optional/default=0.0]
  • dustdensity/teta [real,private/parameter/optional/default=0.0]
  • dustdensity/teta1 [real,private/parameter/optional/default=teta/(teta+tini)]
  • dustdensity/tl0 [real,private/parameter/optional/default=0.0]
  • dustdensity/tl01 [real,private/parameter/optional/default=tl0/(tl0+tini)]
  • dustdensity/tstart_droplet_coagulation [real,private/optional/default=impossible]
  • dustdensity/ueta [real,private/parameter/optional/default=0.0]
  • dustdensity/ul0 [real,private/parameter/optional/default=0.0]
  • dustdensity/widthnd [real,private/optional/default=1.0]
  • dustdensity/xblob_nd [real,private/optional/default=0.0]
  • dustdensity/yblob_nd [real,private/optional/default=0.0]
  • dustdensity/z0_smooth [real,private/optional/default=0.0]
  • dustdensity/z1_smooth [real,private/optional/default=0.0]
  • dustdensity/zblob_nd [real,private/optional/default=0.0]

Subroutines and functions

subroutine  dustdensity/register_dustdensity()

Initialise variables which should know that we solve the compressible hydro equations: ind; increase nvar accordingly.

4-jun-02/axel: adapted from hydro

Use :

farraymanager (farray_register_pde(), farray_index_append(), farray_register_auxiliary()), general (itoa()), sharedvariables (put_shared_variable())

Call to:

svn_id(), get_gamma_etc(), fatal_error(), init_nd(), warning(), set_init_parameters(), request_border_driving(), sinwave_phase(), blob(), hat3d(), initial_condition_nd(), dot2_mn(), multmv_mn(), del6(), dot_mn(), identify_bcs(), cond_spec_cond(), cond_spec_nucl(), del2fj(), dot2fj(), der6(), special_calc_dustdensity(), calc_diagnostics_dustdensity(), max_mn_name(), xysum_mn_name_z(), yzsum_mn_name_x(), border_driving(), set_border_initcond(), getmu(), condensing_species_rate()

subroutine  dustdensity/initialize_dustdensity(f)

Perform any post-parameter-read initialization i.e. calculate derived parameters.

24-nov-02/tony: coded

Parameters:

f (,,*,*) [real,contiguous]

Use :

equationofstate (get_gamma_etc()), sharedvariables (get_shared_variable()), borderprofiles (request_border_driving()), farraymanager (farray_register_global()), general (spline_integral(), itoa()), mpicomm (mpibcast()), special (set_init_parameters())

Call to:

get_gamma_etc(), fatal_error(), init_nd(), warning(), set_init_parameters(), request_border_driving(), sinwave_phase(), blob(), hat3d(), initial_condition_nd(), dot2_mn(), multmv_mn(), del6(), dot_mn(), identify_bcs(), cond_spec_cond(), cond_spec_nucl(), del2fj(), dot2fj(), der6(), special_calc_dustdensity(), calc_diagnostics_dustdensity(), max_mn_name(), xysum_mn_name_z(), yzsum_mn_name_x(), border_driving(), set_border_initcond(), getmu(), condensing_species_rate()

subroutine  dustdensity/init_nd(f)

Initialise dust density; called from start.f90.

7-nov-01/wolf: coded

28-jun-02/axel: added isothermal

Parameters:

f (,,*,*) [real,contiguous]

Use :

equationofstate (cs20()), initcond (hat3d(), sinwave_phase(), posnoise()), initialcondition (initial_condition_nd()), sharedvariables (get_shared_variable()), general (notanumber()), sub (blob())

Call to:

sinwave_phase(), blob(), hat3d(), fatal_error(), initial_condition_nd(), dot2_mn(), multmv_mn(), del6(), warning(), dot_mn(), identify_bcs(), cond_spec_cond(), cond_spec_nucl(), del2fj(), dot2fj(), der6(), special_calc_dustdensity(), calc_diagnostics_dustdensity(), max_mn_name(), xysum_mn_name_z(), yzsum_mn_name_x(), border_driving(), set_border_initcond(), getmu(), condensing_species_rate()

subroutine  dustdensity/pencil_criteria_dustdensity()

All pencils that the Dustdensity module depends on are specified here.

20-11-04/anders: coded

Call to:

fatal_error(), dot2_mn(), multmv_mn(), del6(), warning(), dot_mn(), identify_bcs(), cond_spec_cond(), cond_spec_nucl(), del2fj(), dot2fj(), der6(), special_calc_dustdensity(), calc_diagnostics_dustdensity(), max_mn_name(), xysum_mn_name_z(), yzsum_mn_name_x(), border_driving(), set_border_initcond(), getmu(), condensing_species_rate()

subroutine  dustdensity/pencil_interdep_dustdensity(lpencil_in)

Interdependency among pencils provided by the Dustdensity module is specified here.

20-11-04/anders: coded

Parameters:

lpencil_in (npencils) [logical]

Call to:

fatal_error(), dot2_mn(), multmv_mn(), del6(), warning(), dot_mn(), identify_bcs(), cond_spec_cond(), cond_spec_nucl(), del2fj(), dot2fj(), der6(), special_calc_dustdensity(), calc_diagnostics_dustdensity(), max_mn_name(), xysum_mn_name_z(), yzsum_mn_name_x(), border_driving(), set_border_initcond(), getmu(), condensing_species_rate()

subroutine  dustdensity/calc_pencils_dustdensity(f, p)

Calculate Dustdensity pencils. Most basic pencils should come first, as others may depend on them.

13-nov-04/anders: coded

Parameters:
  • f (,,*,*) [real,inout,contiguous]

  • p [pencil_case,inout]

Use :

sub, general (spline_integral())

Call to:

dot2_mn(), multmv_mn(), del6(), warning(), dot_mn(), identify_bcs(), cond_spec_cond(), fatal_error(), cond_spec_nucl(), del2fj(), dot2fj(), der6(), special_calc_dustdensity(), calc_diagnostics_dustdensity(), max_mn_name(), xysum_mn_name_z(), yzsum_mn_name_x(), border_driving(), set_border_initcond(), getmu(), condensing_species_rate()

subroutine  dustdensity/dndmd_dt(f, df, p)

continuity equation calculate dnd/dt = - u.gradnd - nd*divud

7-jun-02/axel: incoporated from subroutine pde

Parameters:
  • f (,,*,*) [real,in,contiguous]

  • df (,,*,*) [real,inout,contiguous]

  • p [pencil_case,in]

Use :

sub (identify_bcs(), dot_mn(), del2fj(), dot2fj()), special (special_calc_dustdensity()), general (spline_integral()), deriv (der6()), chemistry (cond_spec_cond(), cond_spec_nucl())

Call to:

identify_bcs(), cond_spec_cond(), fatal_error(), cond_spec_nucl(), del2fj(), dot2fj(), der6(), dot_mn(), special_calc_dustdensity(), calc_diagnostics_dustdensity(), max_mn_name(), xysum_mn_name_z(), yzsum_mn_name_x(), border_driving(), set_border_initcond(), getmu(), condensing_species_rate()

subroutine  dustdensity/calc_diagnostics_dustdensity(f, p)

Diagnostic output

Parameters:
  • f (,,*,*) [real,contiguous]

  • p [pencil_case]

Use :

diagnostics

Call to:

max_mn_name(), xysum_mn_name_z(), yzsum_mn_name_x(), border_driving(), set_border_initcond(), getmu(), condensing_species_rate(), fatal_error()

subroutine  dustdensity/dustdensity_before_boundary(f)

Check for dust grain mass interval overflows and redistribute mdbins.

Parameters:

f (,,*,*) [real,contiguous]

Call to:

fatal_error()

subroutine  dustdensity/dustdensity_after_boundary(f)
Parameters:

f (,,*,*) [real,contiguous]

Use :

general (keep_compiler_quiet())

Call to:

fatal_error()

subroutine  dustdensity/read_dustdensity_init_pars(iomsg)
Parameters:

iomsg [character,out]

Use :

file_io (parallel_unit())

Call to:

fatal_error()

subroutine  dustdensity/write_dustdensity_init_pars(unit)
Parameters:

unit [integer,in]

Call to:

fatal_error()

subroutine  dustdensity/read_dustdensity_run_pars(iomsg)
Parameters:

iomsg [character,out]

Use :

file_io (parallel_unit())

Call to:

fatal_error()

subroutine  dustdensity/write_dustdensity_run_pars(unit)
Parameters:

unit [integer,in]

Call to:

fatal_error()

subroutine  dustdensity/rprint_dustdensity(lreset[, lwrite])

Reads and registers print parameters relevant for dust density.

3-may-02/axel: coded

Parameters:
  • lreset [logical]

  • lwrite [logical]

Use :

diagnostics (parse_name()), farraymanager (farray_index_append()), general (itoa(), loptest(), get_species_nr())

Call to:

fatal_error()

subroutine  dustdensity/get_slices_dustdensity(f, slices)

Write slices for animation of Dustdensity variables.

26-jul-06/tony: coded

Parameters:
Use :

slices_methods (assign_slices_scal())

Call to:

fatal_error()

subroutine  dustdensity/impose_dustdensity_floor(f)

Impose a minimum density by setting all lower densities to the minimum value (density_floor). Useful for debugging purposes.

Parameters:

f (,,*,*) [real,inout,contiguous]

Call to:

fatal_error()

subroutine  dustdensity/pushpars2c(p_par)
Parameters:

p_par (200) [integer]

Use :

syscalls (copy_addr()), general (string_to_enum())