diagnostics_outlog

Module

Description

$Id$

** AUTOMATIC CPARAM.INC GENERATION ************************

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

Quick access

Variables:: allocate_cnames, allocate_phiaverages_names, allocate_phizaverages_names, allocate_sound_data, allocate_xyaverages_names, allocate_xzaverages_names, allocate_yaverages_names, allocate_yzaverages_names, ch1davg, ch2davg, clean_line, da_xy_rel1, da_xz_rel1, da_yz_rel1, dl_y_rel1, dt_save, dvol_rel1, eps_rkf_save, expand_cname, expand_cname_full, expand_cname_short, fnamexy_cap, init_xaver, it_save, ixav_max, lcalculating_phiavg_norm, max_name, max_name_int, max_name_real, maxranks, mnamer, p_phiavg_norm, parse_name, parse_name_s, parse_name_sa, parse_name_v, phiaverages_clean_up, phiavg_norm, phiavg_profile, phizaverages_clean_up, report_undefined, sound_clean_up, sum_mn_name, sum_mn_name_arr2, sum_mn_name_real, sum_mn_name_std, sum_name, sum_name_int, sum_name_real, write_sound_append, xyaverages_clean_up, xzaverages_clean_up, xzsum_mn_name_y_mpar, yaverages_clean_up, yzaverages_clean_up, zaverages_clean_up, zsum_mn_name_xy, zsum_mn_name_xy_mpar, zsum_mn_name_xy_mpar_vec, zsum_mn_name_xy_scal, zsum_mn_name_xy_vec
Routines:: allocate_diagnostic_arrays(), allocate_diagnostic_names(), allocate_fnames(), allocate_phiaverages(), allocate_phizaverages(), allocate_sound(), allocate_vnames(), allocate_xyaverages(), allocate_xzaverages(), allocate_yaverages(), allocate_yzaverages(), allocate_zaverages(), allocate_zaverages_data(), calc_nnames(), calc_phiavg_profile(), diagnostic(), diagnostics_clean_up(), diagnostics_diag_reductions(), diagnostics_init_reduc_pointers(), fnames_clean_up(), fparse_name(), gen_form_legend(), get_average_density(), get_from_fname(), initialize_diagnostic_arrays(), initialize_diagnostics(), initialize_time_integrals(), integrate_mn(), integrate_mn_name(), max_mn_name(), name_is_present(), phiaverages_rz(), phisum_mn_name_rz(), phizaverages_r(), phizsum_mn_name_r(), prep_finalize_thread_diagnos(), prints(), report_undefined_diagnostics(), restore_diagnostic_controls(), save_diagnostic_controls(), save_name(), save_name_halfz(), save_name_sound(), set_type(), sign_masked_xyaver(), sum_lim_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), sum_weighted_name(), surf_mn_name(), trim_averages(), vnames_clean_up(), write_1daverages(), write_1daverages_prepare(), write_2daverages(), write_2daverages_prepare(), write_sound(), xyaverages_z(), xyintegrate_mn_name_z(), xymax_mn_name_z(), xysum_mn_name_z(), xysum_mn_name_z_npar(), xzaverages_y(), xzintegrate_mn_name_y(), xzsum_mn_name_y(), yaverages_xz(), yintegrate_mn_name_xz(), ysum_mn_name_xz(), ysum_mn_name_xz_npar(), yzaverages_x(), yzintegrate_mn_name_x(), yzsum_mn_name_x(), yzsum_mn_name_x_mpar(), zaverages_xy(), zsum_mn_name_xy_arr(), zsum_mn_name_xy_arr2(), zsum_mn_name_xy_mpar_scal()

Needed modules

cdata: $Id$
messages: $Id$
mpicomm: $Id$

Variables

diagnostics/ch2davg [character,private]
diagnostics/clean_line [private]
diagnostics/dvol_rel1 [real,private]
diagnostics/expand_cname [public]
diagnostics/expand_cname_full [private]
diagnostics/expand_cname_short [private]
diagnostics/fnamexy_cap (*,*,*) [real,private/allocatable]
diagnostics/init_xaver [public]
diagnostics/max_name [public]
diagnostics/max_name_int [private]
diagnostics/max_name_real [private]
diagnostics/mnamer [integer,private]
diagnostics/parse_name [public]
diagnostics/parse_name_s [private]
diagnostics/parse_name_sa [private]
diagnostics/parse_name_v [private]
diagnostics/phiaverages_clean_up [private]
diagnostics/phiavg_profile (nrcyl,nx) [real,private/optional/default=0.0]
diagnostics/phizaverages_clean_up [private]
diagnostics/report_undefined [private]
diagnostics/sound_clean_up [private]
diagnostics/sum_mn_name [public]
diagnostics/sum_mn_name_arr2 [private]
diagnostics/sum_mn_name_real [private]
diagnostics/sum_mn_name_std [private]
diagnostics/sum_name [public]
diagnostics/sum_name_int [private]
diagnostics/sum_name_real [private]
diagnostics/write_sound_append [private]
diagnostics/xyaverages_clean_up [private]
diagnostics/xzaverages_clean_up [private]
diagnostics/xzsum_mn_name_y_mpar [private]
diagnostics/yaverages_clean_up [private]
diagnostics/yzaverages_clean_up [private]
diagnostics/zaverages_clean_up [private]
diagnostics/zsum_mn_name_xy [public]
diagnostics/zsum_mn_name_xy_mpar [public]
diagnostics/zsum_mn_name_xy_mpar_vec [private]
diagnostics/zsum_mn_name_xy_scal [private]
diagnostics/zsum_mn_name_xy_vec [private]

Subroutines and functions

subroutine diagnostics/initialize_diagnostics()

Setup variables needed for output of diagnostic quantities and averages.

14-aug-03/axel: added dxy, dyz, and dxz 26-aug-13/MR: removed switch first; moved calculation of dVol_rel1 from diagnostic 31-mar-15/MR: added switch for proper volume averaging

subroutine diagnostics/prints()

Reads and registers print parameters gathered from the different modules and marked in `print.in’.

3-may-02/axel: coded

20-aug-13/MR: changes for accumulating and complex diagnostics 26-aug-13/MR: corrected insertion of imaginary values 12-jan-17/MR: undefined diagnostics suppressed in output

subroutine diagnostics/report_undefined_diagnostics()

subroutine diagnostics/gen_form_legend(fform[, legend])

19-aug-13/MR: outsourced from prints 10-jan-17/MR: added adjustment of fixed point format to output value 12-jan-17/MR: undefined diagnostics suppressed in format

Parameters:

fform [character,out]
legend [character,out,]

Use :

general (safe_character_append(), itoa()), sub (noform())

Called from:

initialize_diagnostics(), initialize_diagnostic_arrays(), prints()

Call to:

noform(), outlog(), reduce_zsum(), read_snaptime(), update_snaptime(), output_average_2d(), trim_average(), fparse_name(), fatal_error(), name_is_present(), in_overlap_mask(), integrate_mn(), warning(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), ysum_mn_name_xz_npar(), zsum_yy(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up()

subroutine diagnostics/write_sound(tout)

Reads and registers “sound” parameters gathered from the different modules and marked in `sound.in’.

3-dec-10/dhruba+joern: coded

10-jan-11/MR: modified

subroutine diagnostics/get_average_density(mass_per_proc, average_density)

Calculation of average density.

1-dec-09/dhruba+piyali: adapted from diagnostic 3-may-12/axel+MR: to divide by box_volume, not nw

Parameters:

mass_per_proc [real,in]
average_density [real,out]

Called from:

register_density(), initialize_density(), init_lnrho(), density_before_boundary(), pencil_criteria_density(), pencil_interdep_density(), density_after_boundary(), dlnrho_dt(), impose_density_floor(), read_density_init_pars(), write_density_init_pars(), read_density_run_pars(), write_density_run_pars(), rprint_density(), get_slices_density(), get_slices_pressure(), get_init_average_density(), calc_pencils_density(), anelastic_after_mn()

Call to:

reduce_zsum(), read_snaptime(), update_snaptime(), output_average_2d(), trim_average(), fparse_name(), fatal_error(), name_is_present(), in_overlap_mask(), integrate_mn(), warning(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), ysum_mn_name_xz_npar(), zsum_yy(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up()

subroutine diagnostics/diagnostic(vname, nlname[, lcomplex])

Finalize calculation of diagnostic quantities (0-D).

2-sep-01/axel: coded

14-aug-03/axel: began adding surface integrals 26-aug-13/MR: moved calculation of dVol_rel1 to initialize_diagnostics

added optional parameter lcomplex for use with imaginary part

Parameters:

vname (nlname) [real,inout]
nlname [integer,in,]
lcomplex [logical,in,]

Use :

general (loptest())

Called from:

pde(), write_diagnostics(), calc_all_module_diagnostic_auxiliaries(), perform_diagnostics(), finalize_diagnostics()

Call to:

reduce_zsum(), read_snaptime(), update_snaptime(), output_average_2d(), trim_average(), fparse_name(), fatal_error(), name_is_present(), in_overlap_mask(), integrate_mn(), warning(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), ysum_mn_name_xz_npar(), zsum_yy(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up()

subroutine diagnostics/initialize_time_integrals(f)

Initialize time_integrals for full chunks.

28-jun-07/axel+mreinhard: coded 29-oct-20/hongzhe: added iuxst, iuyst, iuzst

subroutine diagnostics/xyaverages_z()

Calculate xy-averages (still depending on z) NOTE: these averages depend on z, so after summation in x and y they are still distributed over nprocz CPUs; hence the dimensions of fsumz (and fnamez). In other words: the whole xy-average is present in one and the same fsumz, but the result is not complete on any of the processors before mpireduce_sum has been called. This is simpler than collecting results first in processors with the same ipz and different ipy, and then assemble result from the subset of ipz processors which have ipy=0 back on the root processor.

6-jun-02/axel: coded

subroutine diagnostics/xzaverages_y()

Calculate xz-averages (still depending on y).

12-oct-05/anders: adapted from xyaverages_z

subroutine diagnostics/yzaverages_x()

Calculate yz-averages (still depending on x).

2-oct-05/anders: adapted from xyaverages_z

subroutine diagnostics/phizaverages_r()

Calculate phiz-averages (still depending on r).

29-jan-07/wlad: adapted from yzaverages_x and phiaverages_rz

subroutine diagnostics/yaverages_xz()

Calculate y-averages (still depending on x and z).

7-jun-05/axel: adapted from zaverages_xy

subroutine diagnostics/zaverages_xy()

Calculate z-averages (still depending on x and y).

19-jun-02/axel: coded 25-mar-16/MR: extensions for Yin-Yang grid

subroutine diagnostics/phiaverages_rz()

Calculate azimuthal averages (as functions of r_cyl,z). NOTE: these averages depend on (r and) z, so after summation they are still distributed over nprocz CPUs; hence the dimensions of fsumrz (and fnamerz).

9-dec-02/wolf: coded

subroutine diagnostics/write_1daverages()

Write 1d averages (z-averages, .., i.e. quantities that are only 1d after averaging). These are written every it1d (default it1) timesteps and appended to their individual files.

7-aug-03/wolf: coded

24-Nov-2018/PABourdin: redesigned

subroutine diagnostics/write_2daverages_prepare(lwrite)

Prepare l2davg for writing 2D averages. This needs to be done in the beginning of each time step, so the various routines know that they need to calculate averages.

23-nov-03/axel: adapted from write_2daverages and wvid_prepare

Parameters:: lwrite [logical,in]
Use :: sub (update_snaptime(), read_snaptime())
Called from:: helper_loop(), reload(), gen_output(), timeloop()
Call to:: read_snaptime(), update_snaptime(), output_average_2d(), trim_average(), fparse_name(), fatal_error(), name_is_present(), in_overlap_mask(), integrate_mn(), warning(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), ysum_mn_name_xz_npar(), zsum_yy(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up()

subroutine diagnostics/write_2daverages()

Write 2d averages (z-averages, phi-averages, .., i.e. quantities that are still 2d after averaging) if it is time. In analogy to 3d output to VARN, the time interval between two writes is determined by a parameter (t2davg) in run.in.

7-aug-03/wolf: adapted from wsnap

24-Nov-2018/PABourdin: redesigned

subroutine diagnostics/trim_averages()

Trim averages for times past the current time.

25-apr-16/ccyang: coded 23-Nov-2018/PABourdin: redesigned

function diagnostics/fparse_name(iname, cname, ctest, itest, cform)

Parse name and format of scalar print variable On output, ITEST is set to INAME if CNAME matches CTEST and CFORM is set to the format given as default. E.g. if CTEST=’bmax’ *i.e. we are testing input line CNAME for ‘bmax’, CNAME=’bmax’ will be parsed to ITEST=INAME, CFORM=’(1pe10.2)’, CNAME=’bmax(G5.1)’ to ITEST=INAME, CFORM=’G5.1’, CNAME=’brms’ to ITEST=<unchanged, normally 0>, CFORM=’(1pe10.2)’

return value is iname if ctest matches

-1

if so and itest/=0 at call (indicates multiple occurrence of the same diagnostic in print.in

0 if ctest does not match

4-may-02/axel: coded 6-apr-04/wolf: more liberate format reading

26-feb-13/MRprepared for ignoring multiple occurrences of
diagnostics in print.in

26-aug-13/MR : removed unneeded 0p setting in format, added handling of D formats 27-aug-13/MR : reinstated 0p 10-jan-17/MR : added correction of floating-point formats if not sufficient to hold sign

Parameters:

iname [integer,in]
cname [character,inout]
ctest [character,in]
itest [integer,inout]
cform [character,out]

Return:

fparse_name [integer]

Use :

general (safe_character_assign(), itoa())

Called from:

initialize_diagnostics(), initialize_diagnostic_arrays(), prints(), report_undefined_diagnostics(), gen_form_legend(), write_sound(), get_average_density(), diagnostic(), initialize_time_integrals(), xyaverages_z(), xzaverages_y(), yzaverages_x(), phizaverages_r(), yaverages_xz(), zaverages_xy(), phiaverages_rz(), write_1daverages(), write_1daverages_prepare(), write_2daverages_prepare(), write_2daverages(), trim_averages(), fparse_name(), register_testflow(), initialize_testflow(), init_uutest(), pencil_criteria_testflow(), pencil_interdep_testflow(), read_testflow_init_pars(), write_testflow_init_pars(), read_testflow_run_pars(), write_testflow_run_pars(), duutest_dt(), get_slices_testflow(), testflow_before_boundary(), calc_ltestflow_nonlin_terms(), rprint_testflow()

Call to:

fparse_name(), fatal_error(), name_is_present(), in_overlap_mask(), integrate_mn(), warning(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), ysum_mn_name_xz_npar(), zsum_yy(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up()

subroutine diagnostics/set_type(iname[, lsqrt[, llog10[, lint[, lsum[, lmax[, lmin]]]]]])

Sets the diagnostic type in itype_name.

21-sep-17/MR: coded

Parameters:

iname [integer]
lsqrt [logical]
llog10 [logical]
lint [logical]
lsum [logical]
lmax [logical]
lmin [logical]

Call to:

in_overlap_mask(), integrate_mn(), warning(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/save_name(a, iname)

Sets the value of a (must be treated as real) in fname array

26-may-02/axel: adapted from max_mn_name 20-sep-17/MR: removed setting of itype_name as it is set to “save” by default.

Parameters:

a [real]
iname [integer]

Called from:

register_magnetic(), initialize_magnetic(), init_aa(), pencil_criteria_magnetic(), pencil_interdep_magnetic(), magnetic_after_boundary(), daa_dt(), read_magnetic_init_pars(), write_magnetic_init_pars(), read_magnetic_run_pars(), write_magnetic_run_pars(), get_slices_magnetic(), time_integrals_magnetic(), df_diagnos_magnetic(), rescaling_magnetic(), calc_mfield(), initialize_diagnostics(), initialize_diagnostic_arrays(), prints(), register_hydro(), initialize_hydro(), init_uu(), duu_dt(), calculate_vars_magnetic(), update_char_vel_magnetic(), set_ambipolar_diffusion(), diamagnetism(), set_border_magnetic(), eta_shell(), calc_bthresh(), calc_tau_aa_exterior(), helflux(), curflux_ds(), curflux(), calc_bmx(), calc_bmy(), calc_bmzs2(), calc_bmza2(), calc_bmz(), calc_jmx(), calc_jmy(), calc_jmz(), calc_embmz(), calc_emxamz3(), calc_ambmz(), calc_ambmzh(), calc_jmbmz(), calc_bmxy_rms(), calc_bmz_beltrami_phase(), calc_means_hydro(), pencil_criteria_hydro(), pencil_interdep_hydro(), hydro_before_boundary(), register_particles_breakup(), initialize_particles_breakup(), read_particles_breakup_init_pars(), write_particles_breakup_init_pars(), read_particles_breakup_run_pars(), write_particles_breakup_run_pars(), rprint_particles_breakup(), dbreakup_dt(), initialize_particles_collisions(), particles_collisions_timestep(), particles_collisions_pencils(), particles_collisions_blocks(), register_pscalar(), initialize_pscalar(), init_lncc(), pencil_criteria_pscalar(), pencil_interdep_pscalar(), calc_pencils_pscalar(), calc_diagnostics_pscalar(), dlncc_dt(), read_pscalar_init_pars(), write_pscalar_init_pars(), read_pscalar_run_pars(), write_pscalar_run_pars(), rprint_pscalar(), get_slices_pscalar(), pscalar_before_boundary(), calc_mpscalar(), register_shear(), initialize_shear(), read_shear_init_pars(), write_shear_init_pars(), read_shear_run_pars(), write_shear_run_pars(), shear_before_boundary(), pencil_criteria_shear(), pencil_interdep_shear(), calc_pencils_shear(), shearing(), calc_diagnostics_shear(), advance_shear(), shear_variables(), register_special(), initialize_special(), finalize_special(), init_special(), pencil_criteria_special(), pencil_interdep_special(), calc_pencils_special(), dspecial_dt(), get_slices_special(), special_calc_hydro(), special_calc_density(), special_calc_dustdensity(), special_calc_energy(), special_calc_magnetic(), special_calc_pscalar(), special_calc_particles(), special_calc_chemistry(), special_before_boundary(), special_boundconds(), special_after_timestep(), calc_ode_diagnostics_special(), register_testscalar(), initialize_testscalar(), init_cctest(), pencil_criteria_testscalar(), pencil_interdep_testscalar(), read_testscalar_init_pars(), write_testscalar_init_pars(), read_testscalar_run_pars(), write_testscalar_run_pars(), dcctest_dt(), initialize_training(), register_training(), read_training_run_pars(), write_training_run_pars(), training_after_boundary(), calc_diagnostics_training(), dtraining_dt(), register_viscosity(), initialize_viscosity(), read_viscosity_run_pars(), write_viscosity_run_pars(), rprint_viscosity(), pencil_criteria_viscosity(), pencil_interdep_viscosity(), calc_pencils_viscosity(), viscosity_after_boundary(), calc_pencils_magnetic(), sine_avoid_boundary(), helper_loop(), reload(), gen_output(), timeloop(), initialize_testfield_general(), calc_uxb(), calc_diffusive_part(), calc_diffusive_part_prof_0d(), calc_diffusive_part_prof_1d(), calc_inverse_matrix(), calc_coefficients()

Call to:

in_overlap_mask(), integrate_mn(), warning(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/save_name_sound(a, iname, iscoord)

Lists the value of a (must be treated as real) in fname array

3-Dec-10/dhruba+joern: adapted from max_mn_name

25-aug-13/MR: removed unneeded setting of itype, added test of iname.

Parameters:

a [real]
iname [integer]
iscoord [integer]

Called from:

update_char_vel_magnetic(), set_ambipolar_diffusion(), diamagnetism()

Call to:

in_overlap_mask(), integrate_mn(), warning(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/save_name_halfz(a, iname)

Lists the value of a (must be treated as real) in fname array

16-may-09/axel: adapted from save_name

Parameters:

a (2) [real]
iname [integer]

Called from:

register_magnetic(), initialize_magnetic(), init_aa(), pencil_criteria_magnetic(), pencil_interdep_magnetic(), magnetic_after_boundary(), daa_dt(), read_magnetic_init_pars(), write_magnetic_init_pars(), read_magnetic_run_pars(), write_magnetic_run_pars(), get_slices_magnetic(), time_integrals_magnetic(), df_diagnos_magnetic(), rescaling_magnetic(), calc_mfield(), update_char_vel_magnetic(), set_ambipolar_diffusion(), diamagnetism(), set_border_magnetic(), eta_shell(), calc_bthresh(), calc_tau_aa_exterior(), helflux(), curflux_ds(), curflux(), calc_bmx(), calc_bmy(), calc_bmzs2(), calc_bmza2(), calc_bmz(), calc_jmx(), calc_jmy(), calc_jmz(), calc_embmz(), calc_emxamz3(), calc_ambmz(), calc_ambmzh(), calc_pencils_magnetic(), sine_avoid_boundary()

Call to:

in_overlap_mask(), integrate_mn(), warning(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/max_mn_name(a, iname[, lsqrt[, l_dt[, lneg[, lreciprocal]]]])

Successively calculate maximum of a, which is supplied at each call. Start from zero if lfirstpoint=.true.

1-apr-01/axel+wolf: coded 4-may-02/axel: adapted for fname array

23-jun-02/axel: allows for taking square root in the end 29-jun-12/MR: incorporated test for iname/=0

Parameters:

a (nx) [real]
iname [integer]
lsqrt [logical]
l_dt [logical]
lneg [logical]
lreciprocal [logical]

Called from:

Call to:

in_overlap_mask(), integrate_mn(), warning(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/sum_mn_name_halfy(a, iname)

To calculate averages over half the size of box, useful for simulations which includes equator (either cartesian or spherical).

??-???-??/dhruba: aped from sum_mn_name

Parameters:

a (nx) [real]
iname [integer]

Called from:

register_magnetic(), initialize_magnetic(), init_aa(), pencil_criteria_magnetic(), pencil_interdep_magnetic(), daa_dt(), update_char_vel_magnetic(), set_ambipolar_diffusion(), diamagnetism(), calc_pencils_magnetic()

Call to:

warning(), integrate_mn(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/sum_mn_name_halfz(a, iname)

To calculate averages over half the size of box (this time divided along the z direction), useful for simulations which includes equator (either cartesian or spherical).

7-may-09/dhruba: aped from sum_mn_name_halfy

Parameters:

a (nx) [real]
iname [integer]

Called from:

register_magnetic(), initialize_magnetic(), init_aa(), pencil_criteria_magnetic(), pencil_interdep_magnetic(), daa_dt(), update_char_vel_magnetic(), set_ambipolar_diffusion(), diamagnetism(), calc_pencils_magnetic()

Call to:

warning(), integrate_mn(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/sum_weighted_name(a, weight, iname[, lsqrt])

Succesively calculate the weighted sum of a. The result is divided by the total weight in the diagnostics subroutine.

17-apr-06/anders: coded

Parameters:

a (*) [real]
weight (*) [real]
iname [integer]
lsqrt [logical]

Called from:

register_particles_number(), initialize_particles_number(), init_particles_number(), set_particle_number(), pencil_criteria_par_number(), dnpswarm_dt_pencil()

Call to:

warning(), integrate_mn(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/sum_lim_mn_name(a, iname, p)

Successively calculate integral of a, which is supplied at each call. Just takes values between r_int < r < r_ext The purpose is to compute quantities just inside a cylinder or sphere

2-nov-05/wlad: adapted from sum_mn_name

Parameters:

a (nx) [real]
iname [integer]
p [pencil_case]

Called from:

register_particles_nbody(), initialize_particles_nbody(), pencil_criteria_par_nbody(), pencil_interdep_par_nbody(), calc_pencils_par_nbody(), init_particles_nbody(), dvvp_dt_nbody_pencil(), dxxp_dt_nbody(), dvvp_dt_nbody(), loop_through_nbodies(), loop_through_nbodies_inertial(), loop_through_nbodies_fixstar(), get_evr(), point_par_name(), read_particles_nbody_init_pars(), write_particles_nbody_init_pars(), read_particles_nbody_run_pars(), write_particles_nbody_run_pars(), reset_center_of_mass(), integrate_selfgravity(), bcast_nbodyarray(), particles_nbody_special(), get_totalmass(), get_gravity_field_nbody(), calc_torque()

Call to:

warning(), integrate_mn(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/surf_mn_name(a, iname)

Successively calculate surface integral. This routine assumes that “a” contains the partial result for each pencil, so here we just need to add up the contributions from all processors. Start from zero if lfirstpoint=.true.

14-aug-03/axel: adapted from sum_mn_name 15-feb-13/MR: test of iname incorporated

Parameters:

a [real,in]
iname [integer,in]

Called from:

register_magnetic(), initialize_magnetic(), init_aa(), pencil_criteria_magnetic(), pencil_interdep_magnetic(), magnetic_after_boundary(), daa_dt(), time_integrals_magnetic(), df_diagnos_magnetic(), rescaling_magnetic(), update_char_vel_magnetic(), set_ambipolar_diffusion(), diamagnetism(), set_border_magnetic(), eta_shell(), calc_bthresh(), calc_tau_aa_exterior(), helflux(), curflux_ds(), register_testflow(), initialize_testflow(), init_uutest(), pencil_criteria_testflow(), pencil_interdep_testflow(), read_testflow_init_pars(), write_testflow_init_pars(), read_testflow_run_pars(), write_testflow_run_pars(), duutest_dt(), get_slices_testflow(), testflow_before_boundary(), calc_ltestflow_nonlin_terms(), calc_pencils_magnetic(), sine_avoid_boundary()

Call to:

integrate_mn(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), warning(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/integrate_mn(a, inta)

Successively calculate sum of a, which is supplied at each call. Start from zero if lfirstpoint=.true. ultimately multiply by dv to get the integral. This differs from sum_mn_name by the setting of ilabel_integrate and hence in the behaviour in the final step.

Note, for regular integration (uniform meshes) it is better to use the usual sum_mn_name routine with lint=.true.

1-dec-09/dhruba+piyali:

Parameters:

a (nx) [real,in]
inta [real,inout]

Called from:

register_density(), initialize_density(), init_lnrho(), density_before_boundary(), pencil_criteria_density(), pencil_interdep_density(), density_after_boundary(), dlnrho_dt(), impose_density_floor(), read_density_init_pars(), write_density_init_pars(), read_density_run_pars(), write_density_run_pars(), rprint_density(), get_init_average_density(), calc_pencils_density(), initialize_diagnostics(), initialize_diagnostic_arrays(), prints(), report_undefined_diagnostics(), gen_form_legend(), write_sound(), get_average_density(), diagnostic(), initialize_time_integrals(), xyaverages_z(), xzaverages_y(), yzaverages_x(), phizaverages_r(), yaverages_xz(), zaverages_xy(), phiaverages_rz(), write_1daverages(), write_1daverages_prepare(), write_2daverages_prepare(), write_2daverages(), trim_averages(), fparse_name(), set_type(), save_name(), save_name_sound(), save_name_halfz(), max_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), sum_weighted_name(), sum_lim_mn_name(), surf_mn_name(), integrate_mn(), integrate_mn_name()

Call to:

integrate_mn(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), warning(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/integrate_mn_name(a, iname)

Successively calculate sum of a, which is supplied at each call. Start from zero if lfirstpoint=.true. ultimately multiply by dv to get the integral. This differs from sum_mn_name by the setting of ilabel_integrate and hence in the behaviour in the final step.

Note, for regular integration (uniform meshes) it is better to use the usual sum_mn_name routine with lint=.true.

30-may-03/tony: adapted from sum_mn_name 13-nov-06/tony: modified to handle stretched mesh 31-mar-15/MR: use integrate_mn to avoid doubled code

Parameters:

a (nx) [real]
iname [integer]

Called from:

register_magnetic(), initialize_magnetic(), init_aa(), pencil_criteria_magnetic(), pencil_interdep_magnetic(), daa_dt(), register_density(), initialize_density(), init_lnrho(), density_before_boundary(), pencil_criteria_density(), pencil_interdep_density(), density_after_boundary(), dlnrho_dt(), calc_pencils_density(), register_energy(), initialize_energy(), init_energy(), pencil_criteria_energy(), pencil_interdep_energy(), calc_pencils_energy(), denergy_dt(), calc_diagnostics_energy(), register_hydro(), initialize_hydro(), init_uu(), update_char_vel_magnetic(), set_ambipolar_diffusion(), diamagnetism(), register_particles_nbody(), initialize_particles_nbody(), pencil_criteria_par_nbody(), pencil_interdep_par_nbody(), calc_pencils_par_nbody(), init_particles_nbody(), dvvp_dt_nbody_pencil(), dxxp_dt_nbody(), dvvp_dt_nbody(), loop_through_nbodies(), loop_through_nbodies_inertial(), loop_through_nbodies_fixstar(), get_evr(), point_par_name(), read_particles_nbody_init_pars(), write_particles_nbody_init_pars(), read_particles_nbody_run_pars(), write_particles_nbody_run_pars(), reset_center_of_mass(), integrate_selfgravity(), bcast_nbodyarray(), particles_nbody_special(), get_totalmass(), get_gravity_field_nbody(), calc_torque(), calc_means_hydro(), pencil_criteria_hydro(), pencil_interdep_hydro(), register_pscalar(), initialize_pscalar(), init_lncc(), pencil_criteria_pscalar(), pencil_interdep_pscalar(), calc_pencils_pscalar(), calc_diagnostics_pscalar(), dlncc_dt(), register_viscosity(), initialize_viscosity(), read_viscosity_run_pars(), write_viscosity_run_pars(), rprint_viscosity(), pencil_criteria_viscosity(), pencil_interdep_viscosity(), calc_pencils_viscosity(), viscosity_after_boundary(), calc_viscous_heat(), calc_viscous_force(), calc_diagnostics_viscosity(), calc_pencils_magnetic()

Call to:

integrate_mn(), xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), warning(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/xysum_mn_name_z(a, iname)

3-sep-13/MR: derived from xysum_mn_name_z

Parameters:

a (nx) [real,in]
iname [integer,in]

Use :

cdata (n())

Called from:

Call to:

xysum_mn_name_z_npar(), yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), warning(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/xysum_mn_name_z_npar(a, n, iname)

Successively calculate sum over x,y of a, which is supplied at each call. The result fnamez is z-dependent. Start from zero if lfirstpoint=.true.

5-jun-02/axel: adapted from sum_mn_name 3-sep-13/MR: derived from xysum_mn_name_z, index n now parameter

Parameters:

a (nx) [real,in]
n [integer,in]
iname [integer,in]

Called from:

initialize_diagnostics(), initialize_diagnostic_arrays(), prints(), report_undefined_diagnostics(), gen_form_legend(), write_sound(), get_average_density(), diagnostic(), initialize_time_integrals(), xyaverages_z(), xzaverages_y(), yzaverages_x(), phizaverages_r(), yaverages_xz(), zaverages_xy(), phiaverages_rz(), write_1daverages(), write_1daverages_prepare(), write_2daverages_prepare(), write_2daverages(), trim_averages(), fparse_name(), set_type(), save_name(), save_name_sound(), save_name_halfz(), max_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), sum_weighted_name(), sum_lim_mn_name(), surf_mn_name(), integrate_mn(), integrate_mn_name(), xysum_mn_name_z()

Call to:

yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), warning(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/xzsum_mn_name_y(a, iname)

3-sep-13/MR: derived from xzsum_mn_name_y

Parameters:

a (nx) [real,in]
iname [integer,in]

Use :

cdata (m())

Called from:

register_magnetic(), initialize_magnetic(), init_aa(), pencil_criteria_magnetic(), pencil_interdep_magnetic(), daa_dt(), register_energy(), initialize_energy(), init_energy(), pencil_criteria_energy(), pencil_interdep_energy(), calc_pencils_energy(), denergy_dt(), calc_diagnostics_energy(), update_char_vel_magnetic(), set_ambipolar_diffusion(), diamagnetism(), register_shock(), initialize_shock(), read_shock_run_pars(), write_shock_run_pars(), rprint_shock(), get_slices_shock(), pencil_criteria_shock(), pencil_interdep_shock(), calc_pencils_shock(), register_pscalar(), initialize_pscalar(), init_lncc(), pencil_criteria_pscalar(), pencil_interdep_pscalar(), calc_pencils_pscalar(), calc_diagnostics_pscalar(), dlncc_dt(), register_selfgravity(), initialize_selfgravity(), pencil_criteria_selfgravity(), pencil_interdep_selfgravity(), calc_pencils_selfgravity(), calc_selfpotential(), addselfgrav(), calc_diagnostics_selfgrav(), calc_diagnostics_shock(), calc_pencils_magnetic()

Call to:

yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), warning(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/yzsum_mn_name_x(a, iname)

3-sep-13/MR: derived from yzsum_mn_name_x

Parameters:

a (nx) [real,in]
iname [integer,in]

Use :

cdata (m())

Called from:

register_magnetic(), initialize_magnetic(), init_aa(), pencil_criteria_magnetic(), pencil_interdep_magnetic(), daa_dt(), register_energy(), initialize_energy(), init_energy(), pencil_criteria_energy(), pencil_interdep_energy(), calc_pencils_energy(), denergy_dt(), calc_diagnostics_energy(), update_char_vel_magnetic(), set_ambipolar_diffusion(), diamagnetism(), register_particles(), initialize_particles(), init_particles(), particles_dragforce_stiff(), pencil_criteria_particles(), pencil_interdep_particles(), calc_pencils_particles(), dxxp_dt(), dvvp_dt(), dxxp_dt_pencil(), register_shock(), initialize_shock(), read_shock_run_pars(), write_shock_run_pars(), rprint_shock(), get_slices_shock(), pencil_criteria_shock(), pencil_interdep_shock(), calc_pencils_shock(), calc_diagnostics_particles(), register_pscalar(), initialize_pscalar(), init_lncc(), pencil_criteria_pscalar(), pencil_interdep_pscalar(), calc_pencils_pscalar(), calc_diagnostics_pscalar(), dlncc_dt(), register_selfgravity(), initialize_selfgravity(), pencil_criteria_selfgravity(), pencil_interdep_selfgravity(), calc_pencils_selfgravity(), calc_selfpotential(), addselfgrav(), calc_diagnostics_selfgrav(), calc_diagnostics_shock(), register_viscosity(), initialize_viscosity(), read_viscosity_run_pars(), write_viscosity_run_pars(), rprint_viscosity(), pencil_criteria_viscosity(), pencil_interdep_viscosity(), calc_pencils_viscosity(), viscosity_after_boundary(), calc_viscous_heat(), calc_viscous_force(), calc_diagnostics_viscosity(), calc_pencils_magnetic()

Call to:

yzsum_mn_name_x_mpar(), fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), warning(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/yzsum_mn_name_x_mpar(a, m, iname)

Successively calculate sum over y,z of a, which is supplied at each call. The result fnamex is x-dependent. Start from zero if lfirstpoint=.true.

2-oct-05/anders: adapted from xysum_mn_name_z 3-sep-13/MR: derived from yzsum_mn_name_x, m now parameter

Parameters:

a (nx) [real,in]
m [integer,in]
iname [integer,in]

Called from:

initialize_diagnostics(), initialize_diagnostic_arrays(), prints(), report_undefined_diagnostics(), gen_form_legend(), write_sound(), get_average_density(), diagnostic(), initialize_time_integrals(), xyaverages_z(), xzaverages_y(), yzaverages_x(), phizaverages_r(), yaverages_xz(), zaverages_xy(), phiaverages_rz(), write_1daverages(), write_1daverages_prepare(), write_2daverages_prepare(), write_2daverages(), trim_averages(), fparse_name(), set_type(), save_name(), save_name_sound(), save_name_halfz(), max_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), sum_weighted_name(), sum_lim_mn_name(), surf_mn_name(), integrate_mn(), integrate_mn_name(), xysum_mn_name_z(), xysum_mn_name_z_npar(), xymax_mn_name_z(), xzsum_mn_name_y(), yzsum_mn_name_x()

Call to:

fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), warning(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/xyintegrate_mn_name_z(a, iname)

Integrate over x and y. Apply trapezoidal rule properly in the case of non-periodic boundaries.

18-jun-07/tobi: adapted from xysum_mn_name_z

subroutine diagnostics/xzintegrate_mn_name_y(a, iname)

Integrate over x and z. Apply trapezoidal rule properly in the case of non-periodic boundaries.

18-jun-07/tobi: adapted from xzsum_mn_name_y

subroutine diagnostics/yzintegrate_mn_name_x(a, iname)

Integrate over y and z. Apply trapezoidal rule properly in the case of non-periodic boundaries.

18-jun-07/tobi: adapted from yzsum_mn_name_x

subroutine diagnostics/phizsum_mn_name_r(a, iname)

Successively calculate sum over phi,z of a, which is supplied at each call. Start from zero if lfirstpoint=.true. The fnamer array uses one of its slots in mnamer where we put ones and sum them up in order to get the normalization correct.

29-jan-07/wlad: adapted from yzsum_mn_name_x and phisum_mn_name

Parameters:

a (nx) [real]
iname [integer]

Called from:

register_magnetic(), initialize_magnetic(), init_aa(), pencil_criteria_magnetic(), pencil_interdep_magnetic(), daa_dt(), update_char_vel_magnetic(), set_ambipolar_diffusion(), diamagnetism(), calc_pencils_magnetic()

Call to:

fatal_error(), ysum_mn_name_xz_npar(), zsum_yy(), warning(), get_from_fname(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/ysum_mn_name_xz(a, iname)

3-sep-13/MR: derived from ysum_mn_name_xz

Parameters:

a (nx) [real,in]
iname [integer,in]

Use :

cdata (n())

Called from:

register_magnetic(), initialize_magnetic(), init_aa(), pencil_criteria_magnetic(), pencil_interdep_magnetic(), daa_dt(), register_energy(), initialize_energy(), init_energy(), pencil_criteria_energy(), pencil_interdep_energy(), calc_pencils_energy(), denergy_dt(), calc_diagnostics_energy(), update_char_vel_magnetic(), set_ambipolar_diffusion(), diamagnetism(), register_pscalar(), initialize_pscalar(), init_lncc(), pencil_criteria_pscalar(), pencil_interdep_pscalar(), calc_pencils_pscalar(), calc_diagnostics_pscalar(), dlncc_dt(), calc_pencils_magnetic()

Call to:

ysum_mn_name_xz_npar(), zsum_yy(), warning(), get_from_fname(), fatal_error(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/ysum_mn_name_xz_npar(a, n, iname)

Successively calculate sum over y of a, which is supplied at each call. The result fnamexz is xz-dependent. Start from zero if lfirstpoint=.true.

7-jun-05/axel: adapted from zsum_mn_name_xy 3-sep-13/MR: derived from ysum_mn_name_xz, n now parameter

Parameters:

a (nx) [real,in]
n [integer,in]
iname [integer,in]

Called from:

initialize_diagnostics(), initialize_diagnostic_arrays(), prints(), report_undefined_diagnostics(), gen_form_legend(), write_sound(), get_average_density(), diagnostic(), initialize_time_integrals(), xyaverages_z(), xzaverages_y(), yzaverages_x(), phizaverages_r(), yaverages_xz(), zaverages_xy(), phiaverages_rz(), write_1daverages(), write_1daverages_prepare(), write_2daverages_prepare(), write_2daverages(), trim_averages(), fparse_name(), set_type(), save_name(), save_name_sound(), save_name_halfz(), max_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), sum_weighted_name(), sum_lim_mn_name(), surf_mn_name(), integrate_mn(), integrate_mn_name(), xysum_mn_name_z(), xysum_mn_name_z_npar(), xymax_mn_name_z(), xzsum_mn_name_y(), yzsum_mn_name_x(), yzsum_mn_name_x_mpar(), xyintegrate_mn_name_z(), xzintegrate_mn_name_y(), yzintegrate_mn_name_x(), phizsum_mn_name_r(), ysum_mn_name_xz()

Call to:

zsum_yy(), warning(), get_from_fname(), fatal_error(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/zsum_mn_name_xy_mpar_scal(a, m, iname)

Accumulates contributions to z-sum within an mn-loop which are later merged with an mpi reduce command. In the Yang part of a Yin-Yang grid, the contributions to the sum along a (extended) phi coordinate line of the Yin grid are calculated by use of predetermined weights. Likewise for those coordinate lines of the Yin-phi which lie completely within the Yang grid (i.e. the polar caps).

3-apr-16/MR: derived from zsum_mn_name_xy_mpar; extensions for Yin-Yang grid 7-jun-16/MR: outsourced initialize_zaver_yy, reduce_zsum, zsum_y and

corresponding Yin-Yang specific data to Yinyang_mpi

Parameters:

a (nx) [real,in]
m [integer,in]
iname [integer,in]

Use :

cdata (n()), yinyang_mpi (zsum_yy())

Call to:

zsum_yy(), warning(), get_from_fname(), fatal_error(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/calc_phiavg_profile(p)

Calculate profile for phi-averaging for given pencil.

2-feb-03/wolf: coded

Parameters:: p [pencil_case,in]
Called from:: pde(), write_diagnostics(), calc_all_module_diagnostic_auxiliaries(), perform_diagnostics(), finalize_diagnostics()
Call to:: warning(), get_from_fname(), fatal_error(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/phisum_mn_name_rz(a, iname)

Successively calculate sum over phi of a, which is supplied at each call. Start from zero if lfirstpoint=.true. The fnamerz array has one extra slice in z where we put ones and sum them up in order to get the normalization correct.

2-feb-03/wolf: adapted from xysum_mn_name_z

Parameters:

a (nx) [real]
iname [integer]

Called from:

register_magnetic(), initialize_magnetic(), init_aa(), pencil_criteria_magnetic(), pencil_interdep_magnetic(), daa_dt(), pde(), write_diagnostics(), calc_all_module_diagnostic_auxiliaries(), perform_diagnostics(), finalize_diagnostics(), update_char_vel_magnetic(), set_ambipolar_diffusion(), diamagnetism(), register_viscosity(), initialize_viscosity(), read_viscosity_run_pars(), write_viscosity_run_pars(), rprint_viscosity(), pencil_criteria_viscosity(), pencil_interdep_viscosity(), calc_pencils_viscosity(), viscosity_after_boundary(), calc_viscous_heat(), calc_viscous_force(), calc_diagnostics_viscosity(), calc_pencils_magnetic()

Call to:

get_from_fname(), fatal_error(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

function diagnostics/get_from_fname(iname)

Gets value from fname.

30-oct-09/MR: coded

Parameters:: iname [integer,in]
Return:: get_from_fname [real]
Called from:: initialize_diagnostics(), initialize_diagnostic_arrays(), prints(), report_undefined_diagnostics(), gen_form_legend(), write_sound(), get_average_density(), diagnostic(), initialize_time_integrals(), xyaverages_z(), xzaverages_y(), yzaverages_x(), phizaverages_r(), yaverages_xz(), zaverages_xy(), phiaverages_rz(), write_1daverages(), write_1daverages_prepare(), write_2daverages_prepare(), write_2daverages(), trim_averages(), fparse_name(), set_type(), save_name(), save_name_sound(), save_name_halfz(), max_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), sum_weighted_name(), sum_lim_mn_name(), surf_mn_name(), integrate_mn(), integrate_mn_name(), xysum_mn_name_z(), xysum_mn_name_z_npar(), xymax_mn_name_z(), xzsum_mn_name_y(), yzsum_mn_name_x(), yzsum_mn_name_x_mpar(), xyintegrate_mn_name_z(), xzintegrate_mn_name_y(), yzintegrate_mn_name_x(), phizsum_mn_name_r(), ysum_mn_name_xz(), ysum_mn_name_xz_npar(), yintegrate_mn_name_xz(), zsum_mn_name_xy_arr2(), zsum_mn_name_xy_arr(), zsum_mn_name_xy_mpar_scal(), calc_phiavg_profile(), diagnostics_init_reduc_pointers(), diagnostics_diag_reductions(), phisum_mn_name_rz(), register_testflow(), initialize_testflow(), init_uutest(), pencil_criteria_testflow(), pencil_interdep_testflow(), read_testflow_init_pars(), write_testflow_init_pars(), read_testflow_run_pars(), write_testflow_run_pars(), duutest_dt(), get_slices_testflow(), testflow_before_boundary(), calc_ltestflow_nonlin_terms()
Call to:: fatal_error(), location_in_proc(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/allocate_sound(nnamel)

Allocate the variables needed for “sound”.

3-Dec-10/dhruba+joern: coded 11-jan-11/MR: parameter nnamel added

subroutine diagnostics/allocate_fnames(nnamel)

Allocate arrays needed for diagnostics.

23-mar-10/Bourdin.KIS: copied from allocate_yaverages 11-jan-11/MR: parameter nnamel added 18-aug-13/MR: accumulation of diagnostics enabled 25-aug-13/MR: added allocation of itype_name

Parameters:: nnamel [integer,in]
Called from:: initialize_diagnostic_arrays(), write_1daverages_prepare(), xymax_mn_name_z(), yintegrate_mn_name_xz(), zsum_mn_name_xy_arr2(), zsum_mn_name_xy_arr(), diagnostics_init_reduc_pointers(), diagnostics_diag_reductions(), prep_finalize_thread_diagnos(), calc_nnames(), allocate_diagnostic_names(), allocate_diagnostic_arrays(), helper_loop()
Call to:: fatal_error(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/allocate_vnames(nnamel)

Allocate arrays needed for video slices.

23-mar-10/Bourdin.KIS: copied from allocate_yaverages 11-jan-11/MR: parameter nnamel added

Parameters:: nnamel [integer,in]
Called from:: initialize_diagnostic_arrays(), write_1daverages_prepare(), xymax_mn_name_z(), yintegrate_mn_name_xz(), zsum_mn_name_xy_arr2(), zsum_mn_name_xy_arr(), diagnostics_init_reduc_pointers(), diagnostics_diag_reductions(), prep_finalize_thread_diagnos(), calc_nnames(), allocate_diagnostic_names()
Call to:: fatal_error(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/allocate_xyaverages(nnamel)

Allocate arrays needed for xy-averages.

24-nov-09/anders: copied from allocate_yaverages 11-jan-11/MR: parameter nnamel added

Parameters:: nnamel [integer,in]
Called from:: initialize_diagnostic_arrays(), write_1daverages_prepare(), xymax_mn_name_z(), yintegrate_mn_name_xz(), zsum_mn_name_xy_arr2(), zsum_mn_name_xy_arr(), diagnostics_init_reduc_pointers(), diagnostics_diag_reductions(), prep_finalize_thread_diagnos(), calc_nnames(), allocate_diagnostic_names(), allocate_diagnostic_arrays()
Call to:: fatal_error(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/allocate_xzaverages(nnamel)

Allocate arrays needed for xz-averages.

24-nov-09/anders: copied from allocate_yaverages 11-jan-11/MR: parameter nnamel added

Parameters:: nnamel [integer,in]
Called from:: initialize_diagnostic_arrays(), write_1daverages_prepare(), xymax_mn_name_z(), yintegrate_mn_name_xz(), zsum_mn_name_xy_arr2(), zsum_mn_name_xy_arr(), diagnostics_init_reduc_pointers(), diagnostics_diag_reductions(), prep_finalize_thread_diagnos(), calc_nnames(), allocate_diagnostic_names(), allocate_diagnostic_arrays()
Call to:: fatal_error(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/allocate_yzaverages(nnamel)

Allocate arrays needed for yz-averages.

24-nov-09/anders: copied from allocate_yaverages 11-jan-11/MR: parameter nnamel added

Parameters:: nnamel [integer,in]
Called from:: initialize_diagnostic_arrays(), write_1daverages_prepare(), xymax_mn_name_z(), yintegrate_mn_name_xz(), zsum_mn_name_xy_arr2(), zsum_mn_name_xy_arr(), diagnostics_init_reduc_pointers(), diagnostics_diag_reductions(), prep_finalize_thread_diagnos(), calc_nnames(), allocate_diagnostic_names(), allocate_diagnostic_arrays()
Call to:: fatal_error(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/allocate_phizaverages(nnamel)

Allocate arrays needed for phiz-averages.

24-nov-09/anders: copied from allocate_yaverages 11-jan-11/MR: parameter nnamel added

Parameters:: nnamel [integer,in]
Called from:: initialize_diagnostic_arrays(), write_1daverages_prepare(), xymax_mn_name_z(), yintegrate_mn_name_xz(), zsum_mn_name_xy_arr2(), zsum_mn_name_xy_arr(), diagnostics_init_reduc_pointers(), diagnostics_diag_reductions(), prep_finalize_thread_diagnos(), calc_nnames(), allocate_diagnostic_names(), allocate_diagnostic_arrays()
Call to:: fatal_error(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/allocate_yaverages(nnamel)

Allocate arrays needed for y-averages.

12-aug-09/dhruba: coded 11-jan-11/MR: parameter nnamel added

Parameters:: nnamel [integer,in]
Called from:: initialize_diagnostic_arrays(), write_1daverages_prepare(), xymax_mn_name_z(), yintegrate_mn_name_xz(), zsum_mn_name_xy_arr2(), zsum_mn_name_xy_arr(), diagnostics_init_reduc_pointers(), diagnostics_diag_reductions(), prep_finalize_thread_diagnos(), calc_nnames(), allocate_diagnostic_names(), allocate_diagnostic_arrays()
Call to:: fatal_error(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/allocate_zaverages(nnamel)

Allocate arrays needed for z-averages.

12-aug-09/dhruba: coded 11-jan-11/MR: parameter nnamel added 11-mar-16/MR: outsourced allocation of fnamexy to allocate_zaverages_data

Parameters:: nnamel [integer,in]
Called from:: initialize_diagnostic_arrays(), write_1daverages_prepare(), xymax_mn_name_z(), yintegrate_mn_name_xz(), zsum_mn_name_xy_arr2(), zsum_mn_name_xy_arr(), diagnostics_init_reduc_pointers(), diagnostics_diag_reductions(), prep_finalize_thread_diagnos(), calc_nnames(), allocate_diagnostic_names()
Call to:: fatal_error(), initialize_zaver_yy(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/allocate_zaverages_data(nnamel)

Allocate data array needed for z-averages. Additional fnamexy_cap for collection of z-sums in polar caps of Yang grid.

11-mar-16/MR: outsourced from allocate_zaverages

subroutine diagnostics/allocate_phiaverages(nnamel)

Allocate arrays needed for phi-averages.

24-nov-09/anders: copied from allocate_zaverages 11-jan-11/MR: parameter nnamel=iadd+nnamerz instead of nnamerz

Parameters:: nnamel [integer,in]
Called from:: initialize_diagnostic_arrays(), write_1daverages_prepare(), xymax_mn_name_z(), yintegrate_mn_name_xz(), zsum_mn_name_xy_arr2(), zsum_mn_name_xy_arr(), diagnostics_init_reduc_pointers(), diagnostics_diag_reductions(), prep_finalize_thread_diagnos(), calc_nnames(), allocate_diagnostic_names(), allocate_diagnostic_arrays()
Call to:: fatal_error(), xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/sign_masked_xyaver(quan, idiag, ncount)

Forms sign-masked averaging over xy-planes (only positive values count).

Parameters:

ncount (nz) [integer,inout] ::
holds the number of points which contribute.

28-sep-16/MR: coded
quan (nx) [real,in]
idiag [integer,in]

Use :

mpicomm (mpiallreduce_sum_int(), ixyplane())

Call to:

xysum_mn_name_z(), vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/diagnostics_clean_up()

16-jan-17/MR: coded

Called from:: helper_loop(), reload(), gen_output(), timeloop(), check_tspec_crit_log_interval(), check_tsnap_crit_log_interval(), check_tvid_crit_log_interval(), setup_signal_files(), print_metrics(), run_start(), start
Call to:: vnames_clean_up(), fnames_clean_up(), name_is_present()

subroutine diagnostics/fnames_clean_up()

Deallocate space needed for reading the print.in file.

20-apr-10/Bourdin.KIS: copied from xyaverages_clean_up 25-aug-13/MR: added deallocation of itype_name.

Called from:: initialize_diagnostics(), initialize_diagnostic_arrays(), prints(), report_undefined_diagnostics(), gen_form_legend(), write_sound(), get_average_density(), diagnostic(), initialize_time_integrals(), xyaverages_z(), xzaverages_y(), yzaverages_x(), phizaverages_r(), yaverages_xz(), zaverages_xy(), phiaverages_rz(), write_1daverages(), write_1daverages_prepare(), write_2daverages_prepare(), write_2daverages(), trim_averages(), fparse_name(), set_type(), save_name(), save_name_sound(), save_name_halfz(), max_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), sum_weighted_name(), sum_lim_mn_name(), surf_mn_name(), integrate_mn(), integrate_mn_name(), xysum_mn_name_z(), xysum_mn_name_z_npar(), xymax_mn_name_z(), xzsum_mn_name_y(), yzsum_mn_name_x(), yzsum_mn_name_x_mpar(), xyintegrate_mn_name_z(), xzintegrate_mn_name_y(), yzintegrate_mn_name_x(), phizsum_mn_name_r(), ysum_mn_name_xz(), ysum_mn_name_xz_npar(), yintegrate_mn_name_xz(), zsum_mn_name_xy_arr2(), zsum_mn_name_xy_arr(), zsum_mn_name_xy_mpar_scal(), calc_phiavg_profile(), diagnostics_init_reduc_pointers(), diagnostics_diag_reductions(), phisum_mn_name_rz(), get_from_fname(), allocate_sound(), allocate_fnames(), allocate_vnames(), allocate_xyaverages(), allocate_xzaverages(), allocate_yzaverages(), allocate_phizaverages(), allocate_yaverages(), allocate_zaverages(), allocate_zaverages_data(), allocate_phiaverages(), sign_masked_xyaver(), diagnostics_clean_up(), pc_configtest, pc_distribute, pc_extract, pc_reduce, pc_tecplot, pc_tecplot_solid
Call to:: name_is_present()

subroutine diagnostics/vnames_clean_up()

Deallocate space needed for reading the video.in file.

20-apr-10/Bourdin.KIS: copied from xyaverages_clean_up

Called from:: initialize_diagnostics(), initialize_diagnostic_arrays(), prints(), report_undefined_diagnostics(), gen_form_legend(), write_sound(), get_average_density(), diagnostic(), initialize_time_integrals(), xyaverages_z(), xzaverages_y(), yzaverages_x(), phizaverages_r(), yaverages_xz(), zaverages_xy(), phiaverages_rz(), write_1daverages(), write_1daverages_prepare(), write_2daverages_prepare(), write_2daverages(), trim_averages(), fparse_name(), set_type(), save_name(), save_name_sound(), save_name_halfz(), max_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), sum_weighted_name(), sum_lim_mn_name(), surf_mn_name(), integrate_mn(), integrate_mn_name(), xysum_mn_name_z(), xysum_mn_name_z_npar(), xymax_mn_name_z(), xzsum_mn_name_y(), yzsum_mn_name_x(), yzsum_mn_name_x_mpar(), xyintegrate_mn_name_z(), xzintegrate_mn_name_y(), yzintegrate_mn_name_x(), phizsum_mn_name_r(), ysum_mn_name_xz(), ysum_mn_name_xz_npar(), yintegrate_mn_name_xz(), zsum_mn_name_xy_arr2(), zsum_mn_name_xy_arr(), zsum_mn_name_xy_mpar_scal(), calc_phiavg_profile(), diagnostics_init_reduc_pointers(), diagnostics_diag_reductions(), phisum_mn_name_rz(), get_from_fname(), allocate_sound(), allocate_fnames(), allocate_vnames(), allocate_xyaverages(), allocate_xzaverages(), allocate_yzaverages(), allocate_phizaverages(), allocate_yaverages(), allocate_zaverages(), allocate_zaverages_data(), allocate_phiaverages(), sign_masked_xyaver(), diagnostics_clean_up(), pc_configtest, pc_distribute, pc_extract, pc_reduce, pc_tecplot, pc_tecplot_solid
Call to:: name_is_present()

function diagnostics/name_is_present(ccname, vlabel[, form])

Verify if the string vlabel is present or not in the ccname array, return index, put format in form if requested and present

16-sep-10/dintrans: coded 16-may-12/MR: changed function type in integer to return index of found name;

new optional parameter form for use with print.in where the name vlabel comes with a format specification

Parameters:

ccname (*) [character,inout]
vlabel [character,in]
form [character,inout,]

Return:

name_is_present [integer]

Called from:

initialize_diagnostics(), initialize_diagnostic_arrays(), prints(), report_undefined_diagnostics(), gen_form_legend(), write_sound(), get_average_density(), diagnostic(), initialize_time_integrals(), xyaverages_z(), xzaverages_y(), yzaverages_x(), phizaverages_r(), yaverages_xz(), zaverages_xy(), phiaverages_rz(), write_1daverages(), write_1daverages_prepare(), write_2daverages_prepare(), write_2daverages(), trim_averages(), fparse_name(), set_type(), save_name(), save_name_sound(), save_name_halfz(), max_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), sum_weighted_name(), sum_lim_mn_name(), surf_mn_name(), integrate_mn(), integrate_mn_name(), xysum_mn_name_z(), xysum_mn_name_z_npar(), xymax_mn_name_z(), xzsum_mn_name_y(), yzsum_mn_name_x(), yzsum_mn_name_x_mpar(), xyintegrate_mn_name_z(), xzintegrate_mn_name_y(), yzintegrate_mn_name_x(), phizsum_mn_name_r(), ysum_mn_name_xz(), ysum_mn_name_xz_npar(), yintegrate_mn_name_xz(), zsum_mn_name_xy_arr2(), zsum_mn_name_xy_arr(), zsum_mn_name_xy_mpar_scal(), calc_phiavg_profile(), diagnostics_init_reduc_pointers(), diagnostics_diag_reductions(), phisum_mn_name_rz(), get_from_fname(), allocate_sound(), allocate_fnames(), allocate_vnames(), allocate_xyaverages(), allocate_xzaverages(), allocate_yzaverages(), allocate_phizaverages(), allocate_yaverages(), allocate_zaverages(), allocate_zaverages_data(), allocate_phiaverages(), sign_masked_xyaver(), diagnostics_clean_up(), fnames_clean_up(), vnames_clean_up(), update_char_vel_magnetic(), set_ambipolar_diffusion(), diamagnetism(), set_border_magnetic(), eta_shell(), calc_bthresh(), calc_tau_aa_exterior(), helflux(), curflux(), calc_bmx(), calc_bmy(), calc_bmzs2(), calc_bmza2(), calc_bmz(), calc_jmx(), calc_jmy(), calc_jmz(), calc_embmz(), calc_emxamz3(), calc_ambmz(), calc_ambmzh(), calc_jmbmz(), calc_bmxy_rms(), calc_bmz_beltrami_phase(), alfven_x(), alfven_y(), alfven_z(), alfven_xy(), alfven_xz(), alfvenz_rot(), alfvenz_bell(), alfvenz_rot_shear(), torus_test(), force_free_jet(), piecew_dipole_aa(), geo_benchmark_b(), eta_xy_dep(), eta_zdep(), eta_ydep(), eta_xdep(), input_persist_magnetic_id(), input_persist_magnetic()