magnetic_old

Module

Description

$Id$

This modules deals with all aspects of magnetic fields; if no magnetic fields are invoked, a corresponding replacement dummy routine is used instead which absorbs all the calls to the magnetically relevant subroutines listed in here.

** 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 :: lmagnetic = .true.

MVAR CONTRIBUTION 3 MAUX CONTRIBUTION 0

PENCILS PROVIDED aa(3); a2; aij(3,3); bb(3); bbb(3); ab; ua; uxb(3); exa(3) PENCILS PROVIDED b2; bij(3,3); del2a(3); graddiva(3); jj(3) PENCILS PROVIDED j2; jb; va2; jxb(3); jxbr(3); jxbr2; ub; uxb(3); uxb2 PENCILS PROVIDED uxj(3); beta; uga(3); djuidjbi; jo PENCILS PROVIDED ujxb; oxu(3); oxuxb(3); jxbxb(3); jxbrxb(3) PENCILS PROVIDED glnrhoxb(3); del4a(3); del6a(3); oxj(3); diva PENCILS PROVIDED jij(3,3); sj; ss12 PENCILS PROVIDED etava, etaj, etaj2, etajrho PENCILS PROVIDED cosjb,jparallel;jperp PENCILS PROVIDED cosub

Quick access

Types:: magspectra
Variables:: a_relaxprofile, a_relprof, aakim, aakre, aamxy, aamxz, aamz, ab_r, ab_xy, ab_xy2, ab_xy3, ab_xy4, ab_xz, ab_xz2, ab_yz, abc_a, abc_b, abc_c, abcaa, alev, alp_aniso, alpf_mhd, alpha2_inflation, alpha_cutoff_down, alpha_cutoff_up, alpha_dx1, alpha_effect, alpha_eps, alpha_equator, alpha_equator_gap, alpha_gap_step, alpha_inflation, alpha_input, alpha_quenching, alpha_x, alpha_x1, alpha_xprofile, alpha_y, alpha_yprofile, alphassm, ambipolar_diffusion, amp_relprof, ampl_ax, ampl_ay, ampl_az, ampl_b0, ampl_beltrami, ampl_efield, ampl_fcont_aa, ampl_ff, ampl_uy, amplaa, amplaa2, amplaab, amplaaj, angle_dip, aps_xy, aps_xz, aps_xz2, aps_yz, arms, axisr1, axisr2, axisr3, b0_ext, b0_ext_z, b0_ext_z_h, b0_magfric, b1_ext, b2_diamag, b2_r, b2_xy, b2_xy2, b2_xy3, b2_xy4, b2_xz, b2_xz2, b2_yz, b2min, b_ext11, b_ext2, b_ext21, b_ext_inv, b_ext_tmp, battery_term, bb_r, bb_sph_r, bb_sph_xy, bb_sph_xy2, bb_sph_xy3, bb_sph_xy4, bb_sph_xz, bb_sph_xz2, bb_sph_yz, bb_xy, bb_xy2, bb_xy3, bb_xy4, bb_xz, bb_xz2, bb_yz, beta1_inflation, beta1_r, beta1_xy, beta1_xy2, beta1_xy3, beta1_xy4, beta1_xz, beta1_xz2, beta1_yz, beta2_inflation, beta_xy, beta_xy2, beta_xy3, beta_xy4, beta_xz, beta_yz, betamin_jxb, bij_0d_test, bmz, bmz_beltrami_phase, borderaa, boxchargedensity, boxcurrentux, boxomega, boxshear, brms, brms_target, bthresh, bthresh_per_brms, bthresh_scl, by_left, by_right, bz_left, bz_right, bz_stratified, c_light21, calc_0d_diagnostics_magnetic, calc_1d_diagnostics_magnetic, calc_2d_diagnostics_magnetic, calc_aaxyaver, calc_magnetic_slope_limited, calc_pencils_magnetic, calc_pencils_magnetic_pencpar, calc_pencils_magnetic_std, cmin, coefaa, coefbb, compk_aa, conductivity, coskz, cosx, cosz, cutoff_aa, d_smag, daareset, dalpha_x, dalpha_y, damax, damp, delta_effect, detadz, diagnostic_magnetic, diffus_eta, diffus_eta2, diffus_eta3, dipole_moment, displacement_gun, dispr1, dispr2, dispr3, div_sld_magn, domega_x, domega_y, dtcor, edge_moving_layer, eekim, eekre, ekman_friction_aa, ell_jj, enum_ambipolar_diffusion, enum_borderaa, enum_div_sld_magn, enum_iforcing_continuous_aa, enum_ihall_term, enum_rdep_profile, enum_tdep_eta_type, epol, eps_quench, epsi_dip, epsilonaa, eta1, eta1_aniso, eta1_aniso_d, eta1_aniso_r, eta1_aniso_ratio, eta_ampl, eta_aniso_bb, eta_aniso_hyper3, eta_anom, eta_anom_thresh, eta_bb, eta_cspeed, eta_ext, eta_huge, eta_hyper2, eta_hyper3, eta_hyper3_mesh, eta_j, eta_j2, eta_jrho, eta_jump, eta_jump0, eta_jump1, eta_jump2, eta_jump_shock, eta_max, eta_min, eta_mn, eta_out, eta_power_x, eta_power_z, eta_r, eta_r0, eta_r1, eta_rdep, eta_rwidth, eta_rwidth0, eta_rwidth1, eta_shock, eta_shock2, eta_shock_jump1, eta_sld_thresh, eta_smag, eta_spitzer, eta_tdep, eta_tdep_ascale_power, eta_tdep_exponent, eta_tdep_loverride_ee, eta_tdep_t0, eta_tdep_toffset, eta_total, eta_va, eta_width, eta_width_shock, eta_x, eta_x0, eta_x1, eta_xshock, eta_xtdep, eta_xwidth, eta_xwidth0, eta_xwidth1, eta_xy, eta_xy_max, eta_xy_profile, eta_y, eta_y0, eta_y1, eta_ywidth, eta_z, eta_z0, eta_z1, eta_zdep, eta_zdep_coeff, eta_zdep_prof, eta_zgrid, eta_zshock, eta_zwidth, eta_zwidth2, etab, etaj20, exp_epspb, fac_sld_magn, feta_ztdep, fluxtube_border_width, fmax, forcing_continuous, forcing_continuous_aa_amplfact, forcing_continuous_aa_phasefact, forcing_rhs, fres, fring1, fring2, fring3, fring_profile, gamma_epspb, get_eta_t_and_xtdep, get_eta_zdep, get_resistivity, get_resistivity_implicit, geta, geta_r, geta_x, geta_xy, geta_y, geta_z, h_sld_magn, h_slope_limited, hall_tdep_exponent, hall_tdep_t0, hall_term, hall_zdep_exponent, height_eta, hhall, iacou, iacoux, iacouy, iacouz, iakx, iakxim, iaky, iakyim, iakz, iakzim, ibij, icurlb, icurlbx, icurlby, icurlbz, idiag_a2b2m, idiag_a2m, idiag_a2mz, idiag_aa2m, idiag_ab_int, idiag_abm, idiag_abmh, idiag_abmn, idiag_abms, idiag_abmxy, idiag_abmz, idiag_abph1mz, idiag_abph2mz, idiag_abph3mz, idiag_abrms, idiag_abumx, idiag_abumy, idiag_abumz, idiag_abuxmz, idiag_abuymz, idiag_abuzmz, idiag_acbm, idiag_aem, idiag_ajm, idiag_akxpt, idiag_amax, idiag_ambmz, idiag_ambmzh, idiag_ambmzn, idiag_ambmzs, idiag_aphi2m, idiag_apmphi, idiag_apmr, idiag_armphi, idiag_armr, idiag_arms, idiag_axm, idiag_axmxy, idiag_axmxz, idiag_axmz, idiag_axp2, idiag_axpt, idiag_ay2mz, idiag_aybxmz, idiag_aybym2, idiag_aym, idiag_aymxy, idiag_aymxz, idiag_aymz, idiag_ayp2, idiag_aypt, idiag_az2m, idiag_azm, idiag_azmid_max, idiag_azmid_min, idiag_azmid_ymax, idiag_azmid_ymin, idiag_azmphi, idiag_azmr, idiag_azmxy, idiag_azmxz, idiag_azmz, idiag_azp2, idiag_azpt, idiag_b12m, idiag_b1b23m, idiag_b1b32m, idiag_b1m, idiag_b2b13m, idiag_b2b31m, idiag_b2divum, idiag_b2m, idiag_b2mmx, idiag_b2mphi, idiag_b2mr, idiag_b2mx, idiag_b2mxz, idiag_b2mz, idiag_b2ruzm, idiag_b2sphm, idiag_b2tm, idiag_b2uzm, idiag_b3b12m, idiag_b3b21m, idiag_b4m, idiag_b6m, idiag_b8m, idiag_bbxmax, idiag_bbxmz, idiag_bbymax, idiag_bbymz, idiag_bbzmax, idiag_bbzmz, idiag_bcurlfmz, idiag_bdel2amz, idiag_beta1m, idiag_beta1max, idiag_beta1mxy, idiag_beta1mz, idiag_beta2mx, idiag_beta2mz, idiag_betam, idiag_betamax, idiag_betamin, idiag_betamx, idiag_betamz, idiag_bf2m, idiag_bf2mz, idiag_bf4m, idiag_bfrms, idiag_bij2m, idiag_bij_cov_diffmax, idiag_bjtm, idiag_bm, idiag_bm2, idiag_bmax, idiag_bmin, idiag_bmx, idiag_bmxy_rms, idiag_bmy, idiag_bmz, idiag_bmza2, idiag_bmzph, idiag_bmzphe, idiag_bmzs2, idiag_bp2mphi, idiag_bpbzmphi, idiag_bpcmphi, idiag_bphi2m, idiag_bpmphi, idiag_bpmr, idiag_bpol2m, idiag_bpsmphi, idiag_br2mphi, idiag_brbpmphi, idiag_brbzmphi, idiag_brcmphi, idiag_bresrms, idiag_brmphi, idiag_brmr, idiag_brms, idiag_brmsh, idiag_brmsn, idiag_brmss, idiag_brmsx, idiag_brmsz, idiag_brsmphi, idiag_brsphmphi, idiag_bthmphi, idiag_butm, idiag_bx1mxz, idiag_bx2m, idiag_bx2mx, idiag_bx2mxy, idiag_bx2mxz, idiag_bx2my, idiag_bx2mz, idiag_bx2ph1mz, idiag_bx2ph2mz, idiag_bx2ph3mz, idiag_bx2rmz, idiag_bx2rph1mz, idiag_bx2rph2mz, idiag_bx2rph3mz, idiag_bx3m, idiag_bx4m, idiag_bxbym, idiag_bxbymx, idiag_bxbymxy, idiag_bxbymxz, idiag_bxbymy, idiag_bxbymz, idiag_bxbypt, idiag_bxbzm, idiag_bxbzmx, idiag_bxbzmxy, idiag_bxbzmxz, idiag_bxbzmy, idiag_bxbzmz, idiag_bxm, idiag_bxmax, idiag_bxmin, idiag_bxmx, idiag_bxmxy, idiag_bxmxz, idiag_bxmy, idiag_bxmz, idiag_bxp2, idiag_bxph1mz, idiag_bxph2mz, idiag_bxph3mz, idiag_bxpt, idiag_by1mxz, idiag_by2m, idiag_by2mx, idiag_by2mxy, idiag_by2mxz, idiag_by2my, idiag_by2mz, idiag_by2ph1mz, idiag_by2ph2mz, idiag_by2ph3mz, idiag_by2rmz, idiag_by2rph1mz, idiag_by2rph2mz, idiag_by2rph3mz, idiag_by3m, idiag_by4m, idiag_bybzm, idiag_bybzmx, idiag_bybzmxy, idiag_bybzmxz, idiag_bybzmy, idiag_bybzmz, idiag_bybzpt, idiag_bym, idiag_bymax, idiag_bymin, idiag_bymx, idiag_bymxy, idiag_bymxz, idiag_bymy, idiag_bymz, idiag_byp2, idiag_byph1mz, idiag_byph2mz, idiag_byph3mz, idiag_bypt, idiag_bz1mxz, idiag_bz2m, idiag_bz2mphi, idiag_bz2mx, idiag_bz2mxy, idiag_bz2mxz, idiag_bz2my, idiag_bz2mz, idiag_bz2ph1mz, idiag_bz2ph2mz, idiag_bz2ph3mz, idiag_bz2rmz, idiag_bz2rph1mz, idiag_bz2rph2mz, idiag_bz2rph3mz, idiag_bz3m, idiag_bz4m, idiag_bzaymz, idiag_bzbxpt, idiag_bzcmphi, idiag_bzdivamz, idiag_bzlammz, idiag_bzmax, idiag_bzmin, idiag_bzmphi, idiag_bzmr, idiag_bzmx, idiag_bzmxy, idiag_bzmxz, idiag_bzmy, idiag_bzmz, idiag_bzp2, idiag_bzph1mz, idiag_bzph2mz, idiag_bzph3mz, idiag_bzpt, idiag_bzsmphi, idiag_bzuamz, idiag_coshjbm, idiag_cosjbm, idiag_cosubm, idiag_d6abmz, idiag_d6amz1, idiag_d6amz2, idiag_d6amz3, idiag_dbx2m, idiag_dbxdxmxy, idiag_dbxdymxy, idiag_dbxdzmxy, idiag_dbxm, idiag_dbxmax, idiag_dby2m, idiag_dbydxmxy, idiag_dbydymxy, idiag_dbydzmxy, idiag_dbym, idiag_dbymax, idiag_dbz2m, idiag_dbzdxmxy, idiag_dbzdymxy, idiag_dbzdzmxy, idiag_dbzm, idiag_dbzmax, idiag_dexbmx, idiag_dexbmy, idiag_dexbmz, idiag_divamz, idiag_divarms, idiag_divbmax, idiag_divbrms, idiag_dive2m, idiag_divee2m, idiag_djuidjbim, idiag_dobrms, idiag_dtb, idiag_dtbr, idiag_dteta3, idiag_dtfr, idiag_dthr, idiag_dubrms, idiag_e2m, idiag_e3xamz1, idiag_e3xamz2, idiag_e3xamz3, idiag_eem2, idiag_eem3, idiag_eem4, idiag_ekxpt, idiag_emag, idiag_embmz, idiag_emxamz3, idiag_epsad, idiag_epsm, idiag_epsm2, idiag_epsm3, idiag_epsm4, idiag_epsm_les, idiag_epsmmz, idiag_eta_tdep, idiag_etaaniso, idiag_etaanisobb, idiag_etaj2max, idiag_etajmax, idiag_etajrhomax, idiag_etasmagm, idiag_etasmagmax, idiag_etasmagmin, idiag_etatotalmx, idiag_etatotalmxy, idiag_etatotalmz, idiag_etavamax, idiag_exabot, idiag_examx, idiag_examxy1, idiag_examxy2, idiag_examxy3, idiag_examy, idiag_examz, idiag_examz1, idiag_examz2, idiag_examz3, idiag_exatop, idiag_exatotalmx, idiag_exatotalmy, idiag_exatotalmz, idiag_exatotalmz1, idiag_exatotalmz2, idiag_exatotalmz3, idiag_exaym2, idiag_exjm2, idiag_exjmx, idiag_exjmy, idiag_exjmz, idiag_exmxy, idiag_exmxz, idiag_exp2, idiag_expt, idiag_eymxy, idiag_eymxz, idiag_eyp2, idiag_eypt, idiag_ezmxy, idiag_ezmxz, idiag_ezp2, idiag_ezpt, idiag_fbm, idiag_fxbxm, idiag_gb2m, idiag_glamam, idiag_glambm, idiag_gpxbm, idiag_hjbm, idiag_hjparallelm, idiag_hjperpm, idiag_hjrms, idiag_j2b2m, idiag_j2m, idiag_j2mx, idiag_j2mz, idiag_jb_int, idiag_jbm, idiag_jbmh, idiag_jbmn, idiag_jbmphi, idiag_jbms, idiag_jbmx, idiag_jbmxy, idiag_jbmz, idiag_jbph1mz, idiag_jbph2mz, idiag_jbph3mz, idiag_jbrms, idiag_jbtm, idiag_jdel2am, idiag_jdel2amz, idiag_jem, idiag_jfm, idiag_jh2m1, idiag_jm, idiag_jm2, idiag_jmax, idiag_jmbmz, idiag_jmin, idiag_jmx, idiag_jmy, idiag_jmz, idiag_jparallelm, idiag_jperpm, idiag_jrms, idiag_jutm, idiag_jx2m, idiag_jx2m1, idiag_jx2m2, idiag_jx2m3, idiag_jx4m, idiag_jxbm, idiag_jxbmx, idiag_jxbmy, idiag_jxbmz, idiag_jxbpmphi, idiag_jxbr2m, idiag_jxbrmax, idiag_jxbrmphi, idiag_jxbrms, idiag_jxbrqm, idiag_jxbrxm, idiag_jxbrxmx, idiag_jxbrxmy, idiag_jxbrxmz, idiag_jxbrym, idiag_jxbrymx, idiag_jxbrymy, idiag_jxbrymz, idiag_jxbrzm, idiag_jxbrzmx, idiag_jxbrzmy, idiag_jxbrzmz, idiag_jxbxbm, idiag_jxbxm, idiag_jxbym, idiag_jxbzm, idiag_jxbzmphi, idiag_jxmax, idiag_jxmxy, idiag_jxmxz, idiag_jxmz, idiag_jxp2, idiag_jxph1mz, idiag_jxph2mz, idiag_jxph3mz, idiag_jxpt, idiag_jy2m, idiag_jy2m1, idiag_jy2m2, idiag_jy2m3, idiag_jy4m, idiag_jybxm, idiag_jybym, idiag_jybzm, idiag_jym, idiag_jymax, idiag_jymxy, idiag_jymxz, idiag_jymz, idiag_jyp2, idiag_jyph1mz, idiag_jyph2mz, idiag_jyph3mz, idiag_jypt, idiag_jz2m, idiag_jz3m, idiag_jz4m, idiag_jzbxm, idiag_jzbym, idiag_jzbzm, idiag_jzm, idiag_jzmax, idiag_jzmxy, idiag_jzmxz, idiag_jzmz, idiag_jzp2, idiag_jzph1mz, idiag_jzph2mz, idiag_jzph3mz, idiag_jzpt, idiag_km0em, idiag_km1em, idiag_kmz, idiag_kx_aa, idiag_logbm, idiag_magfricmax, idiag_mflux_x, idiag_mflux_y, idiag_mflux_z, idiag_obm, idiag_obmz, idiag_oxuxbm, idiag_phibmx, idiag_phibmy, idiag_phibmz, idiag_poynxmxy, idiag_poynymxy, idiag_poynzmxy, idiag_poynzmz, idiag_poynzph1mz, idiag_poynzph2mz, idiag_poynzph3mz, idiag_rmmz, idiag_rmrms, idiag_sigma, idiag_sijbibjm, idiag_stokesimxy, idiag_stokesq1mxy, idiag_stokesqmxy, idiag_stokesu1mxy, idiag_stokesumxy, idiag_uabxmz, idiag_uabymz, idiag_uabzmz, idiag_uam, idiag_uamz, idiag_ubbzm, idiag_ubgbpm, idiag_ubm, idiag_ubmxy, idiag_ubmz, idiag_ubtm, idiag_udotxbm, idiag_ugb22m, idiag_ujm, idiag_ujmz, idiag_ujtm, idiag_ujxbm, idiag_ujxbmz, idiag_uxbcmx, idiag_uxbcmy, idiag_uxbdotm, idiag_uxbm, idiag_uxbmx, idiag_uxbmy, idiag_uxbmz, idiag_uxbpmphi, idiag_uxbrmphi, idiag_uxbrms, idiag_uxbsmx, idiag_uxbsmy, idiag_uxbxm, idiag_uxbxmz, idiag_uxbym, idiag_uxbymz, idiag_uxbzm, idiag_uxbzmphi, idiag_uxbzmz, idiag_uxdxuxbm, idiag_uxjm, idiag_uxjxm, idiag_uxjym, idiag_uxjzm, idiag_uybxm, idiag_uybxmxz, idiag_uybxmz, idiag_uybym, idiag_uybymz, idiag_uybzm, idiag_uybzmxz, idiag_uybzmz, idiag_uyjxm, idiag_uyjym, idiag_uyjzm, idiag_uzbxm, idiag_uzbxmz, idiag_uzbym, idiag_uzbymz, idiag_uzbzm, idiag_uzbzmz, idiag_uzjxm, idiag_uzjym, idiag_uzjzm, idiag_va23rms, idiag_va2m, idiag_vam, idiag_vamax, idiag_vamin, idiag_vamxz, idiag_varms, idiag_vmagfricmax, idiag_vmagfricmz, idiag_vmagfricrms, idiag_wl2d, idiag_wl3d, idiag_wl3d2, idiva, iekx, iekxim, ieky, iekyim, iekz, iekzim, iforcing_cont_aa, iforcing_continuous_aa, ihall_term, ilam, imp_alpha0, imp_halpha, inclaa, indzav, inertial_length, initaa, initaa2, initaak, initeek, initpower2_aa, initpower_aa, iresistivity, iring1, iring2, iring3, islope_limiter, iua, ivid_ab, ivid_aps, ivid_b2, ivid_bb, ivid_bb_sph, ivid_beta1, ivid_j2, ivid_jb, ivid_jj, ivid_poynting, izav_start, j2_r, j2_xy, j2_xy2, j2_xy3, j2_xy4, j2_xz, j2_xz2, j2_yz, j_ext, j_ext_quench, jb_r, jb_xy, jb_xy2, jb_xy3, jb_xy4, jb_xz, jb_xz2, jb_yz, jj_xy, jj_xy2, jj_xy3, jj_xy4, jj_xz, jj_xz2, jj_yz, k1_ff, k1_ff_mag, k1x_ff, k1y_ff, k1z_ff, k_relprof, keplerian_gauge, kern_jjsmooth, kgaussian_aa, kinflow, kpeak_aa, ksign, kx_aa, kx_aa2, kx_ax, kx_ay, kx_az, ky_aa, ky_aa2, ky_ax, ky_ay, ky_az, kz_aa, kz_aa2, kz_ax, kz_ay, kz_az, l2d_aa, la_relprof_global, laa_as_aux, laa_cou_as_aux, ladd_bb_init, ladd_disp_current_from_aux, ladd_efield, ladd_global_field, ladvective_gauge, ladvective_gauge2, lalfven_as_aux, lalpha_inflation, lalpha_phi_equation, lambipolar_diffusion, lambipolar_strong_coupling, lax_ext_global, lay_ext_global, laz_ext_global, lb_ext_in_comaux, lb_ext_pot, lback_to_unscaled, lbb_as_aux, lbb_as_comaux, lbb_sph_as_aux, lbbt_as_aux, lbdivu_as_aux, lbeta_as_aux, lbeta_inflation, lbext_curvilinear, lbext_moving_layer, lbgu_as_aux, lbij_as_aux, lbij_test, lboris_correction, lbraginsky, lbx_ext_global, lby_ext_global, lbz_ext_global, lcalc_aamean, lcalc_aameanz, lcheck_positive_va2, lcoulomb, lcoulomb_apply, lcovariant_magnetic, lcurlb_as_aux, ldiamagnetism, learly_set_el_pencil, lee_as_aux, lelectron_inertia, lemf, letasmag_as_aux, lexplicit_resistivity, lfactors_aa, lforce_free_test, lforcing_cont_aa, lforcing_cont_aa_local, lfreeze_aext, lfreeze_aint, lfrozen_bb_bot, lfrozen_bb_top, lgauss, lhalox, lhubble_magnetic, lignore_1rho_in_lorentz, lignore_bext_in_b2, limiter_fact, limp_alpha, limplicit_resistivity, linduction, linertial_2, lj2t_as_aux, lj_ext, ljbt_as_aux, ljj_as_aux, ljj_as_comaux, ljjt_as_aux, lkeplerian_gauge, lkinematic, llambda_aa, llimiter, llocal_friction, llorentz_rhoref, llorentzforce, lmagnetic_slope_limited, lmagneto_friction, lmean_friction, lmeanfield_jxb, lmeanfield_jxb_with_va2, lmeanfield_noalpm, lminus_mode, lncr_anticorrelated, lncr_correlated, lneutralion_heat, lno_eta_tdep, lno_noise_aa, lno_ohmic_heat_bound_z, lno_second_ampl_aa, lnoinduction, lnonzero_eta, lnorm_aa_kk, lohmic_heat, lomega_effect, loverride_ee, loverride_ee2, loverride_ee_decide, lpolarization_basis, lpress_equil, lpress_equil_alt, lpress_equil_via_ss, lpropagate_borderaa, lquench_eta_aniso, lquench_inflation, lrandom_ampl_aa, lreheating, lrelaxprof_glob_scaled, lremove_mean_emf, lremove_meanax, lremove_meanaxy, lremove_meanaxz, lremove_meanay, lremove_meanaz, lremove_volume_average, lreset_aa, lreset_vart_only_at_start, lresi_anomalous, lresi_cspeed, lresi_dep, lresi_dust, lresi_eta_aniso, lresi_eta_const, lresi_eta_proptouz, lresi_eta_shock, lresi_eta_shock2, lresi_eta_shock_perp, lresi_eta_shock_profr, lresi_eta_shock_profz, lresi_eta_tdep, lresi_eta_tdep_t0_norm, lresi_eta_xtdep, lresi_eta_ztdep, lresi_etaj, lresi_etaj2, lresi_etajrho, lresi_etass, lresi_etava, lresi_hyper2, lresi_hyper2_tdep, lresi_hyper3, lresi_hyper3_aniso, lresi_hyper3_csmesh, lresi_hyper3_mesh, lresi_hyper3_polar, lresi_hyper3_strict, lresi_hyper3_tdep, lresi_magfield, lresi_rdep, lresi_shell, lresi_smagorinsky, lresi_smagorinsky_cross, lresi_smagorinsky_nusmag, lresi_spitzer, lresi_sqrtrhoeta_const, lresi_vaspeed, lresi_xdep, lresi_xydep, lresi_ydep, lresi_zdep, lresis_const, lresis_hyper3_mesh, lresis_shock, lresis_zdep, lresistivity, lrhs_max, lrun_initaa, lset_axay_zero, lshearboxcurrent, lskip_projection_aa, lsld_bb, lsmooth_jj, lsquash_aa, lswitch_sign_e2, ltime_integrals_always, lua_as_aux, lugb_as_aux, lupw_aa, luse_bext_in_b2, luse_bgb_as_jxb, lvacuum, lvart_in_shear_frame, lxmask_mag, lymask_mag, magnetic_xaver_range, magnetic_yaver_range, magnetic_zaver_range, magnetosonic_x, magspectra, maxdiffus_eta, maxdiffus_eta3, meanfield_bp, meanfield_bs, meanfield_etab, meanfield_kf, meanfield_molecular_eta, meanfield_qp, meanfield_qs, mesh_hyper_resistivity, mix_factor, mu012, mu_ext_pot, mu_r, n_modes_aa, naareset, nbvec, nbvecmax, ncr_quench, ncutoff_aa, nexp_aa, nfact_aa, nlf_sld_magn, no_ohmic_heat_z0, no_ohmic_heat_zwidth, non_ffree_factor, nquench, nresi_max, nrings, nu_ni, nu_ni1, numag, nzav, omega_bz_ext, omega_dx1, omega_tmp0, omega_x, omega_x1, omega_xprofile, omega_y, omega_yc1, omega_yc2, omega_yprofile, pertaa, pertamplaa, phase_aa, phase_ax, phase_ay, phase_az, phase_beltrami, phasex_aa, phasey_aa, phasez_aa, phix, phiy, phiz, pm_smag1, poynting_r, poynting_xy, poynting_xy2, poynting_xy3, poynting_xy4, poynting_xz, poynting_xz2, poynting_yz, qexp_aa, quench_aniso, r12, r_dip, r_inner, r_outer, radrfp, rdep_profile, reinitialize_aa, relhel_aa, rescale_aa, rescaling_fraction, rfprad, rfpradb, rfpradj, rhomin_jxb, rhoref, rhoref1, rm_ext, rm_int, rmode, rnoise_ext, rnoise_int, robflow_aa, rring1, rring2, rring3, scl_uxb_in_ohm, set_advec_va2, sheet_hyp, sheet_position, sheet_thickness, sigma_ratio, sigma_t1, sigma_t2, sinkz, sinx, sinz, source_zav, speed_moving_layer, ssmax, swirl, t0_bext, t1_sigma, t2_sigma, t_bext, taareset, tau_aa_exterior, tau_jj, tau_relprof, tau_relprof1, tau_remove_meanaxy, tauad, tdep_eta_type, th_spot, two_step_factor, ux_const, uxbb, va2max_boris, va2max_jxb, va2power_jxb, va_min, varms, vcrit_anom, w_sldchar_mag, widthaa, wr1, wr2, wr3, wresistivity, x0aa, xdep_profile, xmask_mag, xyaverages_magnetic, xzaverages_magnetic, y0aa, ydep_profile, ymask_mag, yzaverages_magnetic, z0_gaussian, z0aa, z1_aa, z2_aa, z_allprocs, z_surface, zbot_moving_layer, zdep_profile, zmask_mag, zmode, znoise_ext, znoise_int, ztop_moving_layer
Routines:: alfven_rphi(), alfven_rz(), alfven_x(), alfven_xy(), alfven_xz(), alfven_y(), alfven_z(), alfven_zconst(), alfvenz_bell(), alfvenz_rot(), alfvenz_rot_shear(), bb_unitvec_shock(), bc_frozen_in_bb_z(), beltrami_phase(), calc_ambmz(), calc_ambmzh(), calc_bmx(), calc_bmxy_rms(), calc_bmy(), calc_bmz(), calc_bmz_beltrami_phase(), calc_bmza2(), calc_bmzs2(), calc_bthresh(), calc_diagnostic_auxiliaries_magnetic(), calc_diagnostics_magnetic(), calc_embmz(), calc_emxamz3(), calc_jmbmz(), calc_jmx(), calc_jmy(), calc_jmz(), calc_mfield(), calc_pencils_magnetic(), calc_tau_aa_exterior(), calculate_vars_magnetic(), correct_lorentz_force(), curflux(), curflux_ds(), daa_dt(), df_diagnos_magnetic(), diamagnetism(), dynamical_resistivity(), eta_shell(), eta_xdep(), eta_xy_dep(), eta_ydep(), eta_zdep(), expand_shands_magnetic(), force_free_jet(), geo_benchmark_b(), get_bext(), get_slices_magnetic(), helflux(), init_aa(), initialize_magnetic(), input_persist_magnetic(), input_persist_magnetic_id(), input_persistent_magnetic(), magnetic_after_boundary(), magnetic_after_mn(), magnetic_after_timestep(), magnetic_before_boundary(), magnetic_calc_spectra(), output_persistent_magnetic(), pencil_criteria_magnetic(), pencil_interdep_magnetic(), piecew_dipole_aa(), read_magnetic_init_pars(), read_magnetic_run_pars(), register_magnetic(), remove_mean_emf(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical(), rescaling_magnetic(), rprint_magnetic(), set_ambipolar_diffusion(), set_border_magnetic(), sine_avoid_boundary(), split_update_magnetic(), time_integrals_magnetic(), torus_test(), update_char_vel_magnetic(), write_magnetic_init_pars(), write_magnetic_run_pars()

Needed modules

cparam: $Id$
cdata: $Id$
general (keep_compiler_quiet(), loptest()): $Id$
magnetic_meanfield: $Id$
messages (fatal_error(), inevitably_fatal_error(), warning(), svn_id(), timing()): $Id$

Variables

magnetic/aamz (mz,3) [real]
magnetic/alphassm [real,optional/default=0.0]
magnetic/ampl_b0 [real,optional/default=0.0]
magnetic/ampl_beltrami [real,optional/default=0.0]
magnetic/ampl_fcont_aa [real,optional/default=1.0]
magnetic/ampl_ff [real,optional/default=1.0]
magnetic/amplaa (ninit) [real,optional/default=0.0]
magnetic/b1_ext (3) [real]
magnetic/b2_xy (nx,ny) [real,target]
magnetic/b2_xy2 (nx,ny) [real,target]
magnetic/b2_xy3 (nx,ny) [real,target]
magnetic/b2_xy4 (nx,ny) [real,target]
magnetic/b2_xz (nx,nz) [real,target]
magnetic/b2_yz (ny,nz) [real,target]
magnetic/b_ext (3) [real,optional/default=(/0.0,0.0,0.0/)]
magnetic/b_ext11 [real]
magnetic/b_ext21 [real]
magnetic/b_ext_inv (3) [real]
magnetic/b_ext_tmp (3) [real]
magnetic/bb_xy (nx,ny,3) [real,target]
magnetic/bb_xy2 (nx,ny,3) [real,target]
magnetic/bb_xy3 (nx,ny,3) [real,target]
magnetic/bb_xy4 (nx,ny,3) [real,target]
magnetic/bb_xz (nx,nz,3) [real,target]
magnetic/bb_yz (ny,nz,3) [real,target]
magnetic/bbmz (nz,3) [real]
magnetic/beta_xy (nx,ny) [real,target]
magnetic/beta_xy2 (nx,ny) [real,target]
magnetic/beta_xy3 (nx,ny) [real,target]
magnetic/beta_xy4 (nx,ny) [real,target]
magnetic/beta_xz (nx,nz) [real,target]
magnetic/beta_yz (ny,nz) [real,target]
magnetic/bmz [real,optional/default=0]
magnetic/bmz_beltrami_phase [real,optional/default=0.0]
magnetic/borderaa [character,optional/default='nothing']
magnetic/brms [real,optional/default=0.0]
magnetic/brms_target [real,optional/default=1.0]
magnetic/bthresh [real,optional/default=0.0]
magnetic/bthresh_per_brms [real,optional/default=0.0]
magnetic/bthresh_scl [real,optional/default=1.0]
magnetic/by_left [real,optional/default=0.0]
magnetic/by_right [real,optional/default=0.0]
magnetic/bz_left [real,optional/default=0.0]
magnetic/bz_right [real,optional/default=0.0]
magnetic/center1_x [real,optional/default=0.0]
magnetic/center1_y [real,optional/default=0.0]
magnetic/center1_z [real,optional/default=0.0]
magnetic/coefaa (3) [complex,optional/default=(/0.0,0.0,0.0/)]
magnetic/coefbb (3) [complex,optional/default=(/0.0,0.0,0.0/)]
magnetic/coskz (mz) [real]
magnetic/cutoff_aa [real,optional/default=0.0]
magnetic/d_smag [real,optional/default=0.17]
magnetic/daareset [real,optional/default=0.0]
magnetic/epsilonaa [real,optional/default=0.01]
magnetic/eta [real,optional/default=0.0]
magnetic/eta1 [real,optional/default=0.0]
magnetic/eta_aniso_hyper3 (3) [real]
magnetic/eta_anom [real,optional/default=0.0]
magnetic/eta_ext [real,optional/default=0.0]
magnetic/eta_hyper2 [real,optional/default=0.0]
magnetic/eta_hyper3 [real,optional/default=0.0]
magnetic/eta_int [real,optional/default=0.0]
magnetic/eta_j [real,optional/default=0.0]
magnetic/eta_j2 [real,optional/default=0.0]
magnetic/eta_jrho [real,optional/default=0.0]
magnetic/eta_jump [real,optional/default=0.0]
magnetic/eta_min [real,optional/default=0.0]
magnetic/eta_out [real,optional/default=0.0]
magnetic/eta_shock [real,optional/default=0.0]
magnetic/eta_spitzer [real,optional/default=0.0]
magnetic/eta_va [real,optional/default=0.0]
magnetic/eta_width [real,optional/default=0.0]
magnetic/eta_xy (mx,my) [real]
magnetic/eta_xy_max [real,optional/default=1.0]
magnetic/eta_xy_profile [character,optional/default='schnack89']
magnetic/eta_z (mz) [real]
magnetic/eta_z0 [real,optional/default=1.0]
magnetic/eta_z1 [real,optional/default=1.0]
magnetic/etaj20 [real,optional/default=0.0]
magnetic/fluxtube_border_width [real,optional/default=impossible]
magnetic/forcing_continuous_aa_amplfact [real,optional/default=1.0]
magnetic/forcing_continuous_aa_phasefact [real,optional/default=1.0]
magnetic/geta_xy (mx,my,3) [real]
magnetic/geta_z (mz,3) [real]
magnetic/hall_term [real,optional/default=0.0]
magnetic/height_eta [real,optional/default=0.0]
magnetic/idiag_ab_int [integer,optional/default=0]

DIAG_DOC: $intAvcdotBv;dV$
magnetic/idiag_abm [integer,optional/default=0]

DIAG_DOC: $left<AvcdotBvright>$
magnetic/idiag_abmh [integer,optional/default=0]

DIAG_DOC: $left<AvcdotBvright>$ (temp)
magnetic/idiag_abmn [integer,optional/default=0]

DIAG_DOC: $left<AvcdotBvright>$ (north)
magnetic/idiag_abms [integer,optional/default=0]

DIAG_DOC: $left<AvcdotBvright>$ (south)
magnetic/idiag_abmz [integer,optional/default=0]

DIAG_DOC: $left<AvcdotBvright>|_{xy}$
magnetic/idiag_abrms [integer,optional/default=0]

DIAG_DOC: $left<(AvcdotBv)^2right>^{1/2}$
magnetic/idiag_abumx [integer,optional/default=0]

DIAG_DOC: $left<u_xAvcdotBvright>$
magnetic/idiag_abumy [integer,optional/default=0]

DIAG_DOC: $left<u_yAvcdotBvright>$
magnetic/idiag_abumz [integer,optional/default=0]

DIAG_DOC: $left<u_zAvcdotBvright>$
magnetic/idiag_abuxmz [integer,optional/default=0]

DIAG_DOC: $left<(Av cdot Bv) u_x right>_{xy}$
magnetic/idiag_abuymz [integer,optional/default=0]

DIAG_DOC: $left<(Av cdot Bv) u_y right>_{xy}$
magnetic/idiag_abuzmz [integer,optional/default=0]

DIAG_DOC: $left<(Av cdot Bv) u_z right>_{xy}$
magnetic/idiag_ajm [integer,optional/default=0]

DIAG_DOC: $left<jvcdotAvright>$
magnetic/idiag_amax [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_ambmz [integer,optional/default=0]

DIAG_DOC: $left<left<Avright>_{xy}cdotleft<Bvright>_{xy}right>$
magnetic/idiag_ambmzh [integer,optional/default=0]

DIAG_DOC: $left<left<Avright>_{xy}cdotleft<Bvright>_{xy}right>$
magnetic/idiag_ambmzn [integer,optional/default=0]

DIAG_DOC: $left<left<Avright>_{xy}cdotleft<Bvright>_{xy}right>$
magnetic/idiag_ambmzs [integer,optional/default=0]

DIAG_DOC: $left<left<Avright>_{xy}cdotleft<Bvright>_{xy}right>$
magnetic/idiag_apmphi [integer,optional/default=0]

PHIAVG_DOC:
magnetic/idiag_apmr [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_armphi [integer,optional/default=0]

PHIAVG_DOC:
magnetic/idiag_armr [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_arms [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_axm [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_axmxz [integer,optional/default=0]

DIAG_DOC: $left< A_x right>_{xz}$
magnetic/idiag_axmz [integer,optional/default=0]

DIAG_DOC: $left<{cal A}_xright>_{xy}$
magnetic/idiag_aybym2 [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_aym [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_aymxz [integer,optional/default=0]

DIAG_DOC: $left< A_y right>_{xz}$
magnetic/idiag_aymz [integer,optional/default=0]

DIAG_DOC: $left<{cal A}_yright>_{xy}$
magnetic/idiag_azm [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_azmphi [integer,optional/default=0]

PHIAVG_DOC:
magnetic/idiag_azmr [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_azmxz [integer,optional/default=0]

DIAG_DOC: $left< A_z right>_{xz}$
magnetic/idiag_azmz [integer,optional/default=0]

DIAG_DOC: $left<{cal A}_zright>_{xy}$
magnetic/idiag_b1b32m [integer,optional/default=0]

DIAG_DOC: $left<B_1 B_{3,2} right>$
magnetic/idiag_b1m [integer,optional/default=0]

DIAG_DOC: $left<$\Bv$right>$
magnetic/idiag_b2b13m [integer,optional/default=0]

DIAG_DOC: $left<B_2 B_{1,3} right>$
magnetic/idiag_b2divum [integer,optional/default=0]

DIAG_DOC: $left<Bv^2nablacdotuvright>$
magnetic/idiag_b2m [integer,optional/default=0]

DIAG_DOC: $left<Bv^2right>$
magnetic/idiag_b2mphi [integer,optional/default=0]

PHIAVG_DOC: $left<Bv^2right>_varphi$
magnetic/idiag_b2mr [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_b2mz [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_b2ruzm [integer,optional/default=0]

DIAG_DOC: $left<Bv^2rho u_zright>$
magnetic/idiag_b2tm [integer,optional/default=0]

DIAG_DOC: $left<bv(t)cdotint_0^tbv(t’)
magnetic/idiag_b2uzm [integer,optional/default=0]

DIAG_DOC: $left<Bv^2u_zright>$
magnetic/idiag_b3b21m [integer,optional/default=0]

DIAG_DOC: $left<B_3 B_{2,1} right>$
magnetic/idiag_bbxmz [integer,optional/default=0]

DIAG_DOC: $left<{cal B}’_xright>_{xy}$
magnetic/idiag_bbymz [integer,optional/default=0]

DIAG_DOC: $left<{cal B}’_yright>_{xy}$
magnetic/idiag_bbzmz [integer,optional/default=0]

DIAG_DOC: $left<{cal B}’_zright>_{xy}$
magnetic/idiag_bcosphz [integer,optional/default=0]

DIAG_DOC: cosine of phase of a Beltrami field
magnetic/idiag_beta1m [integer,optional/default=0]

DIAG_DOC: $left<Bv^2/(2mu_0 p)right>$
magnetic/idiag_beta1max [integer,optional/default=0]

DIAG_DOC: $max[Bv^2/(2mu_0 p)]$
magnetic/idiag_bjtm [integer,optional/default=0]

DIAG_DOC: $left<bv(t)cdotint_0^tjv(t’)
magnetic/idiag_bm2 [integer,optional/default=0]

DIAG_DOC: $max(Bv^2)$
magnetic/idiag_bmax [integer,optional/default=0]

DIAG_DOC: $max($\Bv$)$
magnetic/idiag_bmx [integer,optional/default=0]

DIAG_DOC: $left<left<Bvright>_{yz}^2
magnetic/idiag_bmxy_rms [integer,optional/default=0]

DIAG_DOC: $sqrt{[left<b_xright>_z(x,y)]^2 +
magnetic/idiag_bmy [integer,optional/default=0]

DIAG_DOC: $left<left<Bvright>_{xz}^2
magnetic/idiag_bmz [integer,optional/default=0]

DIAG_DOC: $left<left<Bvright>_{xy}^2
magnetic/idiag_bmza2 [integer,optional/default=0]

DIAG_DOC: $left<left<Bv_Aright>_{xy}^2right>$
magnetic/idiag_bmzph [integer,optional/default=0]

DIAG_DOC: Phase of a Beltrami field
magnetic/idiag_bmzphe [integer,optional/default=0]

DIAG_DOC: Error of phase of a Beltrami field
magnetic/idiag_bmzs2 [integer,optional/default=0]

DIAG_DOC: $left<left<Bv_Sright>_{xy}^2right>$
magnetic/idiag_bpmphi [integer,optional/default=0]

PHIAVG_DOC: $left<B_varphiright>_varphi$
magnetic/idiag_bpmr [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_brbpmr [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bresrms [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_brmphi [integer,optional/default=0]

PHIAVG_DOC: $left<B_varpiright>_varphi$
magnetic/idiag_brmr [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_brms [integer,optional/default=0]

DIAG_DOC: $left<Bv^2right>^{1/2}$
magnetic/idiag_brmsh [integer,optional/default=0]
magnetic/idiag_brmsn [integer,optional/default=0]
magnetic/idiag_brmss [integer,optional/default=0]
magnetic/idiag_brsphmphi [integer,optional/default=0]

PHIAVG_DOC: $left<B_rright>_varphi$
magnetic/idiag_bsinphz [integer,optional/default=0]

DIAG_DOC: sine of phase of a Beltrami field
magnetic/idiag_bthmphi [integer,optional/default=0]

PHIAVG_DOC: $left<B_varthetaright>_varphi$
magnetic/idiag_bx2m [integer,optional/default=0]

DIAG_DOC: $left<B_x^2right>$
magnetic/idiag_bx2mx [integer,optional/default=0]

DIAG_DOC: $left<B_x^2right>_{yz}$
magnetic/idiag_bx2mxy [integer,optional/default=0]

DIAG_DOC: $left< B_x^2 right>_{z}$
magnetic/idiag_bx2mxz [integer,optional/default=0]

DIAG_DOC: $left< B_x^2 right>_{xz}$
magnetic/idiag_bx2my [integer,optional/default=0]

DIAG_DOC: $left< B_x^2 right>_{xz}$
magnetic/idiag_bx2mz [integer,optional/default=0]

DIAG_DOC: $left< B_x^2 right>_{xy}$
magnetic/idiag_bx2rmz [integer,optional/default=0]

DIAG_DOC: $left< B_x^2/varrho right>_{xy}$
magnetic/idiag_bxbym [integer,optional/default=0]

DIAG_DOC: $left<B_x B_yright>$
magnetic/idiag_bxbymx [integer,optional/default=0]

DIAG_DOC: $left<B_x B_yright>_{yz}$
magnetic/idiag_bxbymxy [integer,optional/default=0]

DIAG_DOC: $left< B_x B_y right>_{z}$
magnetic/idiag_bxbymxz [integer,optional/default=0]

DIAG_DOC: $left< B_x B_y right>_{y}$
magnetic/idiag_bxbymy [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bxbymz [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bxbzm [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bxbzmxy [integer,optional/default=0]

DIAG_DOC: $left< B_x B_z right>_{z}$
magnetic/idiag_bxbzmxz [integer,optional/default=0]

DIAG_DOC: $left< B_x B_z right>_{y}$
magnetic/idiag_bxbzmy [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bxbzmz [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bxm [integer,optional/default=0]

DIAG_DOC: $left<left<Bright>_{yz}^2right>_x^{1/2}$
magnetic/idiag_bxmax [integer,optional/default=0]

DIAG_DOC: $max(|B_x|)$
magnetic/idiag_bxmin [integer,optional/default=0]

DIAG_DOC: $min(|B_x|)$
magnetic/idiag_bxmx [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bxmxy [integer,optional/default=0]

DIAG_DOC: $left< B_x right>_{xy}$
magnetic/idiag_bxmxz [integer,optional/default=0]

DIAG_DOC: $left< B_x right>_{xz}$
magnetic/idiag_bxmy [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bxmz [integer,optional/default=0]

DIAG_DOC: $left<{cal B}_xright>_{xy}$
magnetic/idiag_bxpt [integer,optional/default=0]

DIAG_DOC: $B_x(x_0,y_0,z_0,t)$
magnetic/idiag_by2m [integer,optional/default=0]

DIAG_DOC: $left<B_y^2right>$
magnetic/idiag_by2mx [integer,optional/default=0]

DIAG_DOC: $left<B_y^2right>_{yz}$
magnetic/idiag_by2mxy [integer,optional/default=0]

DIAG_DOC: $left< B_y^2 right>_{z}$
magnetic/idiag_by2mxz [integer,optional/default=0]

DIAG_DOC: $left< B_y^2 right>_{xz}$
magnetic/idiag_by2my [integer,optional/default=0]

DIAG_DOC: $left< B_y^2 right>_{xz}$
magnetic/idiag_by2mz [integer,optional/default=0]

DIAG_DOC: $left< B_y^2 right>_{xy}$
magnetic/idiag_by2rmz [integer,optional/default=0]

DIAG_DOC: $left< B_y^2/varrho right>_{xy}$
magnetic/idiag_bybzm [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bybzmxy [integer,optional/default=0]

DIAG_DOC: $left< B_y B_z right>_{z}$
magnetic/idiag_bybzmxz [integer,optional/default=0]

DIAG_DOC: $left< B_y B_z right>_{y}$
magnetic/idiag_bybzmy [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bybzmz [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bym [integer,optional/default=0]

DIAG_DOC: $left<left<Bright>_{zx}^2right>_y^{1/2}$
magnetic/idiag_bymax [integer,optional/default=0]

DIAG_DOC: $max(|B_y|)$
magnetic/idiag_bymin [integer,optional/default=0]

DIAG_DOC: $min(|B_y|)$
magnetic/idiag_bymx [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bymxy [integer,optional/default=0]

DIAG_DOC: $left< B_y right>_{xy}$
magnetic/idiag_bymxz [integer,optional/default=0]

DIAG_DOC: $left< B_y right>_{xz}$
magnetic/idiag_bymy [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bymz [integer,optional/default=0]

DIAG_DOC: $left<{cal B}_yright>_{xy}$
magnetic/idiag_bypt [integer,optional/default=0]

DIAG_DOC: $B_y(x_0,y_0,z_0,t)$
magnetic/idiag_bz2m [integer,optional/default=0]

DIAG_DOC: $left<B_z^2right>$
magnetic/idiag_bz2mx [integer,optional/default=0]

DIAG_DOC: $left<B_z^2right>_{yz}$
magnetic/idiag_bz2mxy [integer,optional/default=0]

DIAG_DOC: $left< B_z^2 right>_{z}$
magnetic/idiag_bz2mxz [integer,optional/default=0]

DIAG_DOC: $left< B_z^2 right>_{xz}$
magnetic/idiag_bz2my [integer,optional/default=0]

DIAG_DOC: $left< B_z^2 right>_{xz}$
magnetic/idiag_bz2mz [integer,optional/default=0]

DIAG_DOC: $left< B_z^2 right>_{xy}$
magnetic/idiag_bz2rmz [integer,optional/default=0]

DIAG_DOC: $left< B_z^2/varrho right>_{xy}$
magnetic/idiag_bzm [integer,optional/default=0]

DIAG_DOC: $left<left<Bright>_{xy}^2right>_z^{1/2}$
magnetic/idiag_bzmax [integer,optional/default=0]

DIAG_DOC: $max(|B_z|)$
magnetic/idiag_bzmin [integer,optional/default=0]

DIAG_DOC: $min(|B_z|)$
magnetic/idiag_bzmphi [integer,optional/default=0]

PHIAVG_DOC: $left<B_zright>_varphi$
magnetic/idiag_bzmr [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bzmx [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bzmxy [integer,optional/default=0]

DIAG_DOC: $left< B_z right>_{xy}$
magnetic/idiag_bzmxz [integer,optional/default=0]

DIAG_DOC: $left< B_z right>_{xz}$
magnetic/idiag_bzmy [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_bzmz [integer,optional/default=0]

DIAG_DOC: $left<{cal B}_zright>_{xy}$
magnetic/idiag_bzpt [integer,optional/default=0]

DIAG_DOC: $B_z(x_0,y_0,z_0,t)$
magnetic/idiag_cosjbm [integer,optional/default=0]

DIAG_DOC: $left<JvcdotBv/($\Jv$,|Bv|)right>$
magnetic/idiag_cosubm [integer,optional/default=0]

DIAG_DOC: $left<UvcdotBv/($\Uv$,|Bv|)right>$
magnetic/idiag_dexbmx [integer,optional/default=0]

DIAG_DOC: $left<nablatimesEvtimesBvright>|_x$
magnetic/idiag_dexbmy [integer,optional/default=0]

DIAG_DOC: $left<nablatimesEvtimesBvright>|_y$
magnetic/idiag_dexbmz [integer,optional/default=0]

DIAG_DOC: $left<nablatimesEvtimesBvright>|_z$
magnetic/idiag_djuidjbim [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_dobrms [integer,optional/default=0]

DIAG_DOC: $left<(boldsymbol{omega}-Bv)^2right>^{1/2}$
magnetic/idiag_dtb [integer,optional/default=0]

DIAG_DOC: $delta t / [c_{delta t},delta x
magnetic/idiag_dteta [integer,optional/default=0]

DIAG_DOC: $delta t/[c_{delta t,{rm v}},
magnetic/idiag_dubrms [integer,optional/default=0]

DIAG_DOC: $left<(uv-Bv)^2right>^{1/2}$
magnetic/idiag_emag [integer,optional/default=0]

DIAG_DOC: $int_V{1over2mu_0}Bv^2, dV$
magnetic/idiag_embmz [integer,optional/default=0]

DIAG_DOC: $left<left<Evright>_{xy}cdotleft<Bvright>_{xy}
magnetic/idiag_emxamz3 [integer,optional/default=0]

DIAG_DOC: $left<left<Evright>_{xy}timesleft<Avright>_{xy}
magnetic/idiag_epsad [integer,optional/default=0]

DIAG_DOC: $left<rho^{-1} t_{rm AD} (vec{J}timesvec{B})^2right>$ (heating by ion-neutrals friction)
magnetic/idiag_epsm [integer,optional/default=0]

DIAG_DOC: $left<2etamu_0jv^2right>$
magnetic/idiag_epsm_les [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_etaj2max [integer,optional/default=0]

DIAG_DOC: Max of artificial resistivity
magnetic/idiag_etajmax [integer,optional/default=0]

DIAG_DOC: Max of artificial resistivity
magnetic/idiag_etajrhomax [integer,optional/default=0]

DIAG_DOC: Max of artificial resistivity
magnetic/idiag_etasmagm [integer,optional/default=0]

DIAG_DOC: Mean of Smagorinsky resistivity
magnetic/idiag_etasmagmax [integer,optional/default=0]

DIAG_DOC: Max of Smagorinsky resistivity
magnetic/idiag_etasmagmin [integer,optional/default=0]

DIAG_DOC: Min of Smagorinsky resistivity
magnetic/idiag_etatotalmx [integer,optional/default=0]

DIAG_DOC: $left<etaright>_{yz}$
magnetic/idiag_etatotalmz [integer,optional/default=0]

DIAG_DOC: $left<etaright>_{xy}$
magnetic/idiag_etavamax [integer,optional/default=0]

DIAG_DOC: Max of artificial resistivity
magnetic/idiag_exabot [integer,optional/default=0]

DIAG_DOC: $intEvtimesAv,dS|_{rm bot}$
magnetic/idiag_examx [integer,optional/default=0]

DIAG_DOC: $left<EvtimesAvright>|_x$
magnetic/idiag_examy [integer,optional/default=0]

DIAG_DOC: $left<EvtimesAvright>|_y$
magnetic/idiag_examz [integer,optional/default=0]

DIAG_DOC: $left<EvtimesAvright>|_z$
magnetic/idiag_examz1 [integer,optional/default=0]

DIAG_DOC: $left<EvtimesAvright>_{xy}|_x$
magnetic/idiag_examz2 [integer,optional/default=0]

DIAG_DOC: $left<EvtimesAvright>_{xy}|_y$
magnetic/idiag_examz3 [integer,optional/default=0]

DIAG_DOC: $left<EvtimesAvright>_{xy}|_z$
magnetic/idiag_exatop [integer,optional/default=0]

DIAG_DOC: $intEvtimesAv,dS|_{rm top}$
magnetic/idiag_exaym2 [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_exjm2 [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_exjmx [integer,optional/default=0]

DIAG_DOC: $left<EvtimesJvright>|_x$
magnetic/idiag_exjmy [integer,optional/default=0]

DIAG_DOC: $left<EvtimesJvright>|_y$
magnetic/idiag_exjmz [integer,optional/default=0]

DIAG_DOC: $left<EvtimesJvright>|_z$
magnetic/idiag_exmxz [integer,optional/default=0]

DIAG_DOC: $left<{cal E}_xright>_{y}$
magnetic/idiag_exmz [integer,optional/default=0]

DIAG_DOC: $left<{cal E}_xright>_{xy}$
magnetic/idiag_expt [integer,optional/default=0]

DIAG_DOC: ${cal E}_x(x_0,y_0,z_0,t)$
magnetic/idiag_eymxz [integer,optional/default=0]

DIAG_DOC: $left<{cal E}_yright>_{y}$
magnetic/idiag_eymz [integer,optional/default=0]

DIAG_DOC: $left<{cal E}_yright>_{xy}$
magnetic/idiag_eypt [integer,optional/default=0]

DIAG_DOC: ${cal E}_y(x_0,y_0,z_0,t)$
magnetic/idiag_ezmxz [integer,optional/default=0]

DIAG_DOC: $left<{cal E}_zright>_{y}$
magnetic/idiag_ezmz [integer,optional/default=0]

DIAG_DOC: $left<{cal E}_zright>_{xy}$
magnetic/idiag_ezpt [integer,optional/default=0]

DIAG_DOC: ${cal E}_z(x_0,y_0,z_0,t)$
magnetic/idiag_fbm [integer,optional/default=0]

DIAG_DOC: $left<fvcdotBvright>$
magnetic/idiag_fxbxm [integer,optional/default=0]

DIAG_DOC: $left<f_x B_xright>$
magnetic/idiag_gpxbm [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_j2m [integer,optional/default=0]

DIAG_DOC: $left<jv^2right>$
magnetic/idiag_jb_int [integer,optional/default=0]

DIAG_DOC: $intjvcdotBv;dV$
magnetic/idiag_jbm [integer,optional/default=0]

DIAG_DOC: $left<jvcdotBvright>$
magnetic/idiag_jbmh [integer,optional/default=0]

DIAG_DOC: $left<AvcdotBvright>$ (temp)
magnetic/idiag_jbmn [integer,optional/default=0]

DIAG_DOC: $left<AvcdotBvright>$ (north)
magnetic/idiag_jbmphi [integer,optional/default=0]

PHIAVG_DOC: $left<JvcdotBvright>_varphi$
magnetic/idiag_jbms [integer,optional/default=0]

DIAG_DOC: $left<AvcdotBvright>$ (south)
magnetic/idiag_jbmz [integer,optional/default=0]

DIAG_DOC: $left<JvcdotBvright>|_{xy}$
magnetic/idiag_jbtm [integer,optional/default=0]

DIAG_DOC: $left<jv(t)cdotint_0^tbv(t’)
magnetic/idiag_jfm [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_jm2 [integer,optional/default=0]

DIAG_DOC: $max(jv^2)$
magnetic/idiag_jmax [integer,optional/default=0]

DIAG_DOC: $max($\jv$)$
magnetic/idiag_jmbmz [integer,optional/default=0]

DIAG_DOC: $left<left<Jvright>_{xy}cdotleft<Bvright>_{xy}
magnetic/idiag_jmx [integer,optional/default=0]

DIAG_DOC: $left<left<Jvright>_{yz}^2
magnetic/idiag_jmy [integer,optional/default=0]

DIAG_DOC: $left<left<Jvright>_{xz}^2
magnetic/idiag_jmz [integer,optional/default=0]

DIAG_DOC: $left<left<Jvright>_{xy}^2
magnetic/idiag_jparallelm [integer,optional/default=0]

DIAG_DOC: Mean value of the component
magnetic/idiag_jperpm [integer,optional/default=0]

DIAG_DOC: Mean value of the component
magnetic/idiag_jrms [integer,optional/default=0]

DIAG_DOC: $left<jv^2right>^{1/2}$
magnetic/idiag_jxbm [integer,optional/default=0]

DIAG_DOC: $left<jvtimesBvright>cdotBv_0/B_0^2$
magnetic/idiag_jxbmx [integer,optional/default=0]

DIAG_DOC: $left<(jvtimesBv)_xright>$
magnetic/idiag_jxbmy [integer,optional/default=0]

DIAG_DOC: $left<(jvtimesBv)_yright>$
magnetic/idiag_jxbmz [integer,optional/default=0]

DIAG_DOC: $left<(jvtimesBv)_zright>$
magnetic/idiag_jxbpmphi [integer,optional/default=0]

PHIAVG_DOC:
magnetic/idiag_jxbr2m [integer,optional/default=0]

DIAG_DOC: $left<(JvtimesBv/rho)^2right>$
magnetic/idiag_jxbrmphi [integer,optional/default=0]

PHIAVG_DOC:
magnetic/idiag_jxbrxm [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_jxbrxmx [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_jxbrxmy [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_jxbrxmz [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_jxbrym [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_jxbrymx [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_jxbrymy [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_jxbrymz [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_jxbrzm [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_jxbrzmx [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_jxbrzmy [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_jxbrzmz [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_jxbxbm [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_jxbzmphi [integer,optional/default=0]

PHIAVG_DOC:
magnetic/idiag_jxmxy [integer,optional/default=0]

DIAG_DOC: $left< J_x right>_{xy}$
magnetic/idiag_jxmz [integer,optional/default=0]

DIAG_DOC: $left<{cal J}_xright>_{xy}$
magnetic/idiag_jymxy [integer,optional/default=0]

DIAG_DOC: $left< J_y right>_{xy}$
magnetic/idiag_jymz [integer,optional/default=0]

DIAG_DOC: $left<{cal J}_yright>_{xy}$
magnetic/idiag_jzmxy [integer,optional/default=0]

DIAG_DOC: $left< J_z right>_{xy}$
magnetic/idiag_jzmz [integer,optional/default=0]

DIAG_DOC: $left<{cal J}_zright>_{xy}$
magnetic/idiag_kmz [integer,optional/default=0]

DIAG_DOC: $left<left<Jvright>_{xy}cdotleft<Bvright>_{xy}right>/
magnetic/idiag_kx_aa [integer,optional/default=0]

DIAG_DOC: $k_x$
magnetic/idiag_mflux_x [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_mflux_y [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_mflux_z [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_oxuxbm [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_phibmx [integer,optional/default=0]

DIAG_DOC: $left<phiBvright>|_x$
magnetic/idiag_phibmy [integer,optional/default=0]

DIAG_DOC: $left<phiBvright>|_y$
magnetic/idiag_phibmz [integer,optional/default=0]

DIAG_DOC: $left<phiBvright>|_z$
magnetic/idiag_rmrms [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_uabxmz [integer,optional/default=0]

DIAG_DOC: $left<(uv cdot Av) B_x right>_{xy}$
magnetic/idiag_uabymz [integer,optional/default=0]

DIAG_DOC: $left<(uv cdot Av) B_y right>_{xy}$
magnetic/idiag_uabzmz [integer,optional/default=0]

DIAG_DOC: $left<(uv cdot Av) B_z right>_{xy}$
magnetic/idiag_uam [integer,optional/default=0]

DIAG_DOC: $left<uvcdotAvright>$
magnetic/idiag_uamz [integer,optional/default=0]

DIAG_DOC: $left<uvcdotAvright>|_{xy}$
magnetic/idiag_ubbzm [integer,optional/default=0]

DIAG_DOC: $left<(uvcdotBv)B_zright>$
magnetic/idiag_ubm [integer,optional/default=0]

DIAG_DOC: $left<uvcdotBvright>$
magnetic/idiag_udotxbm [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_ujm [integer,optional/default=0]

DIAG_DOC: $left<uvcdotJvright>$
magnetic/idiag_ujxbm [integer,optional/default=0]

DIAG_DOC: $left<uvcdot(JvtimesBvright>$
magnetic/idiag_uxbcmx [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_uxbcmy [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_uxbdotm [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_uxbm [integer,optional/default=0]

DIAG_DOC: $left<uvtimesBvright>cdotBv_0/B_0^2$
magnetic/idiag_uxbmx [integer,optional/default=0]

DIAG_DOC: $left<(uvtimesBv)_xright>$
magnetic/idiag_uxbmy [integer,optional/default=0]

DIAG_DOC: $left<(uvtimesBv)_yright>$
magnetic/idiag_uxbmz [integer,optional/default=0]

DIAG_DOC: $left<(uvtimesBv)_zright>$
magnetic/idiag_uxbpmphi [integer,optional/default=0]

PHIAVG_DOC:
magnetic/idiag_uxbrmphi [integer,optional/default=0]

PHIAVG_DOC:
magnetic/idiag_uxbrms [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_uxbsmx [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_uxbsmy [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_uxbxm [integer,optional/default=0]

DIAG_DOC: $left<u_xB_xright>$
magnetic/idiag_uxbzmphi [integer,optional/default=0]

PHIAVG_DOC:
magnetic/idiag_uxdxuxbm [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_uxjm [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_uybym [integer,optional/default=0]

DIAG_DOC: $left<u_yB_yright>$
magnetic/idiag_uzbzm [integer,optional/default=0]

DIAG_DOC: $left<u_zB_zright>$
magnetic/idiag_va2m [integer,optional/default=0]

DIAG_DOC:
magnetic/idiag_vamax [integer,optional/default=0]

DIAG_DOC: $max(Bv^2/varrho)^{1/2}$
magnetic/idiag_varms [integer,optional/default=0]

DIAG_DOC: $left<Bv^2/varrhoright>^{1/2}$
magnetic/iforcing_continuous_aa [character,optional/default='fixed_swirl']
magnetic/inclaa [real,optional/default=0.0]
magnetic/inertial_length [real,optional/default=0.0]
magnetic/initaa (ninit) [character,optional/default='nothing']
magnetic/initpower_aa [real,optional/default=0.0]
magnetic/iresistivity (4) [character,optional/default='']
magnetic/jb_xy (nx,ny) [real,target]
magnetic/jb_xy2 (nx,ny) [real,target]
magnetic/jb_xy3 (nx,ny) [real,target]
magnetic/jb_xy4 (nx,ny) [real,target]
magnetic/jb_xz (nx,nz) [real,target]
magnetic/jb_yz (ny,nz) [real,target]
magnetic/jj_xy (nx,ny,3) [real,target]
magnetic/jj_xy2 (nx,ny,3) [real,target]
magnetic/jj_xy3 (nx,ny,3) [real,target]
magnetic/jj_xy4 (nx,ny,3) [real,target]
magnetic/jj_xz (nx,nz,3) [real,target]
magnetic/jj_yz (ny,nz,3) [real,target]
magnetic/jjmz (nz,3) [real]
magnetic/k1_ff [real,optional/default=1.0]
magnetic/k1x_ff [real,optional/default=1.0]
magnetic/k1y_ff [real,optional/default=1.0]
magnetic/k1z_ff [real,optional/default=1.0]
magnetic/kx_aa (ninit) [real,optional/default=1.0]
magnetic/ky_aa (ninit) [real,optional/default=1.0]
magnetic/kz_aa (ninit) [real,optional/default=1.0]
magnetic/ladvective_gauge [logical,optional/default=.false.]
magnetic/lb_ext_pot [logical,optional/default=.false.]
magnetic/lbb_as_aux [logical,optional/default=.false.]
magnetic/lbbt_as_aux [logical,optional/default=.false.]
magnetic/lbext_curvilinear [logical,optional/default=.true.]
magnetic/lcalc_aamean [logical,optional/default=.false.]
magnetic/lcheck_positive_va2 [logical,optional/default=.false.]
magnetic/lelectron_inertia [logical,optional/default=.false.]
magnetic/lforce_free_test [logical,optional/default=.false.]
magnetic/lforcing_cont_aa [logical,optional/default=.false.]
magnetic/lforcing_cont_aa_local [logical,optional/default=.false.]
magnetic/lfreeze_aext [logical,optional/default=.false.]
magnetic/lfreeze_aint [logical,optional/default=.false.]
magnetic/lfrozen_bb_bot (3) [logical,target/optional/default=(/.false.,.false.,.false./)]
magnetic/lfrozen_bb_top (3) [logical,target/optional/default=(/.false.,.false.,.false./)]
magnetic/lgauss [logical,optional/default=.false.]
magnetic/lhalox [logical,optional/default=.false.]
magnetic/linduction [logical,optional/default=.true.]
magnetic/linertial_2 [real]
magnetic/ljj_as_aux [logical,optional/default=.false.]
magnetic/ljjt_as_aux [logical,optional/default=.false.]
magnetic/lkinematic [logical,optional/default=.false.]
magnetic/llambda_aa [real,optional/default=0.0]
magnetic/llorentzforce [logical,optional/default=.true.]
magnetic/lmean_friction [logical,optional/default=.false.]
magnetic/lmeanfield_theory [logical,optional/default=.false.]
magnetic/lneutralion_heat [logical,optional/default=.true.]
magnetic/lohmic_heat [logical,optional/default=.true.]
magnetic/lpress_equil [logical,optional/default=.false.]
magnetic/lpress_equil_alt [logical,optional/default=.false.]
magnetic/lpress_equil_via_ss [logical,optional/default=.false.]
magnetic/lremove_mean_emf [logical,optional/default=.false.]
magnetic/lreset_aa [logical,optional/default=.false.]
magnetic/lresi_anomalous [logical,optional/default=.false.]
magnetic/lresi_dust [logical,optional/default=.false.]
magnetic/lresi_eta_const [logical,optional/default=.false.]
magnetic/lresi_eta_shock [logical,optional/default=.false.]
magnetic/lresi_eta_shock_perp [logical,optional/default=.false.]
magnetic/lresi_etaj [logical,optional/default=.false.]
magnetic/lresi_etaj2 [logical,optional/default=.false.]
magnetic/lresi_etajrho [logical,optional/default=.false.]
magnetic/lresi_etass [logical,optional/default=.false.]
magnetic/lresi_etava [logical,optional/default=.false.]
magnetic/lresi_hyper2 [logical,optional/default=.false.]
magnetic/lresi_hyper3 [logical,optional/default=.false.]
magnetic/lresi_hyper3_aniso [logical,optional/default=.false.]
magnetic/lresi_hyper3_polar [logical,optional/default=.false.]
magnetic/lresi_hyper3_strict [logical,optional/default=.false.]
magnetic/lresi_shell [logical,optional/default=.false.]
magnetic/lresi_smagorinsky [logical,optional/default=.false.]
magnetic/lresi_smagorinsky_cross [logical,optional/default=.false.]
magnetic/lresi_spitzer [logical,optional/default=.false.]
magnetic/lresi_sqrtrhoeta_const [logical,optional/default=.false.]
magnetic/lresi_xydep [logical,optional/default=.false.]
magnetic/lresi_zdep [logical,optional/default=.false.]
magnetic/lrun_initaa [logical,optional/default=.false.]
magnetic/lupw_aa [logical,optional/default=.false.]
magnetic/luse_bext_in_b2 [logical,optional/default=.false.]
magnetic/lweyl_gauge [logical,optional/default=.false.]
magnetic/mu012 [real,optional/default=0.5]

(=1/2mu0)
magnetic/mu_ext_pot [real,optional/default=-0.5]
magnetic/mu_r [real,optional/default=-0.5]

(still needed for backwards compatibility)
magnetic/n_modes_aa [integer,optional/default=1]
magnetic/naareset [integer]
magnetic/nbvec [integer]
magnetic/nbvecmax [integer,optional/default=nx*ny*nz/4]
magnetic/nresi_max [integer,parameter=4]
magnetic/nrings [integer,optional/default=2]
magnetic/nu_ni [real,optional/default=0.0]
magnetic/nu_ni1 [real]
magnetic/omega_bz_ext [real,optional/default=0.0]
magnetic/phase_beltrami [real,optional/default=0.0]
magnetic/phasex_aa (ninit) [real,optional/default=0.0]
magnetic/phasey_aa (ninit) [real,optional/default=0.0]
magnetic/phasez_aa (ninit) [real,optional/default=0.0]
magnetic/radius [real,optional/default=0.1]
magnetic/reinitialize_aa [logical,optional/default=.false.]
magnetic/relhel_aa [real,optional/default=1.0]
magnetic/rescale_aa [real,optional/default=0.0]
magnetic/rescaling_fraction [real,optional/default=1.0]
magnetic/rhomin_jxb [real,optional/default=0.0]
magnetic/rm_ext [real,optional/default=0.0]
magnetic/rm_int [real,optional/default=0.0]
magnetic/rmode [real,optional/default=1.0]
magnetic/rnoise_ext [real,optional/default=impossible]
magnetic/rnoise_int [real,optional/default=impossible]
magnetic/sigma_ratio [real,optional/default=1.0]
magnetic/sinkz (mz) [real]
magnetic/swirl [real,optional/default=1.0]
magnetic/taareset [real,optional/default=0.0]
magnetic/tau_aa_exterior [real,optional/default=0.0]
magnetic/uxbb (nx,3) [real]
magnetic/va2max_jxb [real,optional/default=0.0]
magnetic/va2power_jxb [integer,optional/default=5]
magnetic/vcrit_anom [real,optional/default=1.0]
magnetic/widthaa [real,optional/default=0.5]
magnetic/wresistivity [real,optional/default=0.01]
magnetic/x0aa [real,optional/default=0.0]
magnetic/z0aa [real,optional/default=0.0]
magnetic/z_surface [real,optional/default=0.0]
magnetic/zdep_profile [character,optional/default='fs']
magnetic/zmode [real,target/optional/default=1.0]

(temporary)

Subroutines and functions

subroutine magnetic/register_magnetic()

Initialise variables which should know that we solve for the vector potential: iaa, etc; increase nvar accordingly

1-may-02/wolf: coded

Use :: farraymanager (farray_register_pde())
Called from:: register_modules()
Call to:: svn_id(), eta_xy_dep(), eta_zdep(), fatal_error(), warning(), request_border_driving(), init_aa(), acosy(), modev(), modeb(), sph_constb(), const_lou(), power_randomphase(), random_isotropic_ks(), beltrami(), bessel_x(), bessel_az_x(), wave_uu(), diffrot(), htube(), htube_x(), htube_erf(), hfluxlayer(), hfluxlayer_y(), vfluxlayer(), magsupport(), arcade_x(), halfcos_x(), uniform_x(), uniform_y(), uniform_z(), uniform_phi(), phi_comp_over_r(), vfield(), vfield2(), bipolar(), bipolar_restzero(), vecpatternxy(), bjump(), sinxsinz(), sin2x_sin2y_cosz(), cosx_cosy_cosz(), sinx_siny_cosz(), x_siny_cosz(), x1_siny_cosz(), cosx_siny_cosz(), cosy_sinz(), x3_cosy_cosz(), mdi_init(), cosx_coscosy_cosz(), crazy(), strange(), sinwave(), coswave(), sinwave_phase(), alfven_x(), alfven_y(), alfven_z(), alfven_xy(), alfven_xz(), alfven_rphi(), alfven_zconst(), alfven_rz(), alfvenz_rot(), alfvenz_rot_shear(), sine_avoid_boundary(), piecew_dipole_aa(), ferriere_uniform_x(), ferriere_uniform_y(), robertsflow(), force_free_jet(), geo_benchmark_b(), torus_test(), initial_condition_aa(), boundconds_x(), initiate_isendrcv_bdry(), finalize_isendrcv_bdry(), boundconds_y(), boundconds_z(), curl(), dot2_mn(), stop_it(), gij(), div_mn(), curl_mn(), dot_mn(), cross_mn(), gij_etc(), multsv_mn(), del4v(), del6v(), timing(), identify_bcs(), der6(), del6fjv(), eta_shell(), special_calc_magnetic(), erfunc(), calc_tau_aa_exterior(), set_border_magnetic(), daa_dt_meanfield(), max_mn_name(), integrate_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), save_name(), helflux(), curflux(), dot_mn_sv(), yzsum_mn_name_x(), xzsum_mn_name_y(), xysum_mn_name_z(), phizsum_mn_name_r(), yzintegrate_mn_name_x(), xzintegrate_mn_name_y(), xyintegrate_mn_name_z(), phisum_mn_name_rz(), ysum_mn_name_xz(), calc_bthresh(), vecout(), dot_mn_sm(), der_z(), der2_z(), set_border_initcond(), border_driving(), der_step(), read_snaptime(), update_snaptime(), tridag(), surf_mn_name(), calc_bmx(), calc_bmy(), calc_bmz(), calc_bmzs2(), calc_bmza2(), calc_jmx(), calc_jmy(), calc_jmz(), calc_emxamz3(), calc_embmz(), calc_ambmz(), calc_ambmzh(), calc_jmbmz(), calc_bmxy_rms(), calc_bmz_beltrami_phase(), save_name_halfz(), correct_lorentz_force(), gamma_function(), hypergeometric2f1(), get_radial_distance(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical()

subroutine magnetic/initialize_magnetic(f)

Perform any post-parameter-read initialization

24-nov-02/tony: dummy routine - nothing to do at present 20-may-03/axel: reinitialize_aa added

Parameters:: f (mx,my,mz,mfarray) [real]
Use :: magnetic_meanfield (initialize_magnetic_mf()), borderprofiles (request_border_driving()), farraymanager, sharedvariables (put_shared_variable()), equationofstate (cs0())
Called from:: register_modules(), initialize_modules()
Call to:: eta_xy_dep(), eta_zdep(), fatal_error(), warning(), request_border_driving(), init_aa(), acosy(), modev(), modeb(), sph_constb(), const_lou(), power_randomphase(), random_isotropic_ks(), beltrami(), bessel_x(), bessel_az_x(), wave_uu(), diffrot(), htube(), htube_x(), htube_erf(), hfluxlayer(), hfluxlayer_y(), vfluxlayer(), magsupport(), arcade_x(), halfcos_x(), uniform_x(), uniform_y(), uniform_z(), uniform_phi(), phi_comp_over_r(), vfield(), vfield2(), bipolar(), bipolar_restzero(), vecpatternxy(), bjump(), sinxsinz(), sin2x_sin2y_cosz(), cosx_cosy_cosz(), sinx_siny_cosz(), x_siny_cosz(), x1_siny_cosz(), cosx_siny_cosz(), cosy_sinz(), x3_cosy_cosz(), mdi_init(), cosx_coscosy_cosz(), crazy(), strange(), sinwave(), coswave(), sinwave_phase(), alfven_x(), alfven_y(), alfven_z(), alfven_xy(), alfven_xz(), alfven_rphi(), alfven_zconst(), alfven_rz(), alfvenz_rot(), alfvenz_rot_shear(), sine_avoid_boundary(), piecew_dipole_aa(), ferriere_uniform_x(), ferriere_uniform_y(), robertsflow(), force_free_jet(), geo_benchmark_b(), torus_test(), initial_condition_aa(), boundconds_x(), initiate_isendrcv_bdry(), finalize_isendrcv_bdry(), boundconds_y(), boundconds_z(), curl(), dot2_mn(), stop_it(), gij(), div_mn(), curl_mn(), dot_mn(), cross_mn(), gij_etc(), multsv_mn(), del4v(), del6v(), timing(), identify_bcs(), der6(), del6fjv(), eta_shell(), special_calc_magnetic(), erfunc(), calc_tau_aa_exterior(), set_border_magnetic(), daa_dt_meanfield(), max_mn_name(), integrate_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), save_name(), helflux(), curflux(), dot_mn_sv(), yzsum_mn_name_x(), xzsum_mn_name_y(), xysum_mn_name_z(), phizsum_mn_name_r(), yzintegrate_mn_name_x(), xzintegrate_mn_name_y(), xyintegrate_mn_name_z(), phisum_mn_name_rz(), ysum_mn_name_xz(), calc_bthresh(), vecout(), dot_mn_sm(), der_z(), der2_z(), set_border_initcond(), border_driving(), der_step(), read_snaptime(), update_snaptime(), tridag(), surf_mn_name(), calc_bmx(), calc_bmy(), calc_bmz(), calc_bmzs2(), calc_bmza2(), calc_jmx(), calc_jmy(), calc_jmz(), calc_emxamz3(), calc_embmz(), calc_ambmz(), calc_ambmzh(), calc_jmbmz(), calc_bmxy_rms(), calc_bmz_beltrami_phase(), save_name_halfz(), correct_lorentz_force(), gamma_function(), hypergeometric2f1(), get_radial_distance(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical()

subroutine magnetic/init_aa(f)

initialise magnetic field; called from start.f90 AB: maybe we should here call different routines (such as rings) AB: and others, instead of accummulating all this in a huge routine. We have an init parameter (initaa) to stear magnetic i.c. independently.

7-nov-2001/wolf: coded

Parameters:: f (mx,my,mz,mfarray) [real]
Use :: equationofstate, farraymanager, gravity (gravz(), z1(), z2()), initcond, boundcond, initialcondition (initial_condition_aa()), mpicomm, sharedvariables, sub
Called from:: register_magnetic(), initialize_magnetic(), start
Call to:: acosy(), modev(), modeb(), sph_constb(), const_lou(), power_randomphase(), random_isotropic_ks(), beltrami(), bessel_x(), bessel_az_x(), wave_uu(), diffrot(), htube(), htube_x(), htube_erf(), hfluxlayer(), hfluxlayer_y(), vfluxlayer(), magsupport(), arcade_x(), halfcos_x(), uniform_x(), uniform_y(), uniform_z(), uniform_phi(), phi_comp_over_r(), vfield(), vfield2(), bipolar(), bipolar_restzero(), vecpatternxy(), bjump(), sinxsinz(), sin2x_sin2y_cosz(), cosx_cosy_cosz(), sinx_siny_cosz(), x_siny_cosz(), x1_siny_cosz(), cosx_siny_cosz(), cosy_sinz(), x3_cosy_cosz(), mdi_init(), cosx_coscosy_cosz(), crazy(), strange(), sinwave(), coswave(), sinwave_phase(), alfven_x(), alfven_y(), alfven_z(), alfven_xy(), alfven_xz(), alfven_rphi(), alfven_zconst(), alfven_rz(), alfvenz_rot(), alfvenz_rot_shear(), sine_avoid_boundary(), piecew_dipole_aa(), ferriere_uniform_x(), ferriere_uniform_y(), robertsflow(), force_free_jet(), geo_benchmark_b(), torus_test(), fatal_error(), initial_condition_aa(), boundconds_x(), initiate_isendrcv_bdry(), finalize_isendrcv_bdry(), boundconds_y(), boundconds_z(), curl(), dot2_mn(), stop_it(), gij(), div_mn(), curl_mn(), dot_mn(), cross_mn(), gij_etc(), multsv_mn(), del4v(), del6v(), timing(), identify_bcs(), der6(), del6fjv(), eta_shell(), special_calc_magnetic(), erfunc(), calc_tau_aa_exterior(), set_border_magnetic(), daa_dt_meanfield(), max_mn_name(), integrate_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), save_name(), helflux(), curflux(), dot_mn_sv(), yzsum_mn_name_x(), xzsum_mn_name_y(), xysum_mn_name_z(), phizsum_mn_name_r(), yzintegrate_mn_name_x(), xzintegrate_mn_name_y(), xyintegrate_mn_name_z(), phisum_mn_name_rz(), ysum_mn_name_xz(), calc_bthresh(), vecout(), dot_mn_sm(), der_z(), der2_z(), set_border_initcond(), border_driving(), der_step(), read_snaptime(), update_snaptime(), tridag(), surf_mn_name(), calc_bmx(), calc_bmy(), calc_bmz(), calc_bmzs2(), calc_bmza2(), calc_jmx(), calc_jmy(), calc_jmz(), calc_emxamz3(), calc_embmz(), calc_ambmz(), calc_ambmzh(), calc_jmbmz(), calc_bmxy_rms(), calc_bmz_beltrami_phase(), save_name_halfz(), correct_lorentz_force(), gamma_function(), hypergeometric2f1(), get_radial_distance(), warning(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical()

subroutine magnetic/pencil_criteria_magnetic()

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

19-nov-04/anders: coded

Use :: mpicomm (stop_it())
Called from:: register_modules(), initialize_modules(), finalize_modules(), choose_pencils()
Call to:: stop_it(), dot2_mn(), gij(), div_mn(), curl_mn(), fatal_error(), dot_mn(), cross_mn(), gij_etc(), multsv_mn(), del4v(), del6v(), timing(), identify_bcs(), der6(), del6fjv(), eta_shell(), special_calc_magnetic(), erfunc(), calc_tau_aa_exterior(), set_border_magnetic(), daa_dt_meanfield(), max_mn_name(), integrate_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), save_name(), helflux(), curflux(), dot_mn_sv(), yzsum_mn_name_x(), xzsum_mn_name_y(), xysum_mn_name_z(), phizsum_mn_name_r(), yzintegrate_mn_name_x(), xzintegrate_mn_name_y(), xyintegrate_mn_name_z(), phisum_mn_name_rz(), ysum_mn_name_xz(), calc_bthresh(), vecout(), dot_mn_sm(), der_z(), der2_z(), set_border_initcond(), border_driving(), der_step(), read_snaptime(), update_snaptime(), tridag(), surf_mn_name(), calc_bmx(), calc_bmy(), calc_bmz(), calc_bmzs2(), calc_bmza2(), calc_jmx(), calc_jmy(), calc_jmz(), calc_emxamz3(), calc_embmz(), calc_ambmz(), calc_ambmzh(), calc_jmbmz(), calc_bmxy_rms(), calc_bmz_beltrami_phase(), save_name_halfz(), correct_lorentz_force(), gamma_function(), hypergeometric2f1(), get_radial_distance(), warning(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical(), curl()

subroutine magnetic/pencil_interdep_magnetic(lpencil_in)

Interdependency among pencils from the Magnetic module is specified here.

19-nov-04/anders: coded

Parameters:: lpencil_in (npencils) [logical]
Called from:: register_modules(), initialize_modules(), finalize_modules(), choose_pencils()
Call to:: dot2_mn(), gij(), div_mn(), curl_mn(), fatal_error(), dot_mn(), cross_mn(), gij_etc(), multsv_mn(), del4v(), del6v(), timing(), identify_bcs(), der6(), del6fjv(), eta_shell(), stop_it(), special_calc_magnetic(), erfunc(), calc_tau_aa_exterior(), set_border_magnetic(), daa_dt_meanfield(), max_mn_name(), integrate_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), save_name(), helflux(), curflux(), dot_mn_sv(), yzsum_mn_name_x(), xzsum_mn_name_y(), xysum_mn_name_z(), phizsum_mn_name_r(), yzintegrate_mn_name_x(), xzintegrate_mn_name_y(), xyintegrate_mn_name_z(), phisum_mn_name_rz(), ysum_mn_name_xz(), calc_bthresh(), vecout(), dot_mn_sm(), der_z(), der2_z(), set_border_initcond(), border_driving(), der_step(), read_snaptime(), update_snaptime(), tridag(), surf_mn_name(), calc_bmx(), calc_bmy(), calc_bmz(), calc_bmzs2(), calc_bmza2(), calc_jmx(), calc_jmy(), calc_jmz(), calc_emxamz3(), calc_embmz(), calc_ambmz(), calc_ambmzh(), calc_jmbmz(), calc_bmxy_rms(), calc_bmz_beltrami_phase(), save_name_halfz(), correct_lorentz_force(), gamma_function(), hypergeometric2f1(), get_radial_distance(), warning(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical(), curl()

subroutine magnetic/calc_pencils_magnetic(f, p)

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

19-nov-04/anders: coded

Parameters:

f (mx,my,mz,mfarray) [real,inout]
p [pencil_case,inout] :: consistency check…

Use :

sub, diagnostics (sum_mn_name()), sharedvariables (put_shared_variable())

Called from:

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

Call to:

dot2_mn(), gij(), div_mn(), curl_mn(), fatal_error(), dot_mn(), cross_mn(), gij_etc(), multsv_mn(), del4v(), del6v(), timing(), identify_bcs(), der6(), del6fjv(), eta_shell(), stop_it(), special_calc_magnetic(), erfunc(), calc_tau_aa_exterior(), set_border_magnetic(), daa_dt_meanfield(), max_mn_name(), integrate_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), save_name(), helflux(), curflux(), dot_mn_sv(), yzsum_mn_name_x(), xzsum_mn_name_y(), xysum_mn_name_z(), phizsum_mn_name_r(), yzintegrate_mn_name_x(), xzintegrate_mn_name_y(), xyintegrate_mn_name_z(), phisum_mn_name_rz(), ysum_mn_name_xz(), calc_bthresh(), vecout(), dot_mn_sm(), der_z(), der2_z(), set_border_initcond(), border_driving(), der_step(), read_snaptime(), update_snaptime(), tridag(), surf_mn_name(), calc_bmx(), calc_bmy(), calc_bmz(), calc_bmzs2(), calc_bmza2(), calc_jmx(), calc_jmy(), calc_jmz(), calc_emxamz3(), calc_embmz(), calc_ambmz(), calc_ambmzh(), calc_jmbmz(), calc_bmxy_rms(), calc_bmz_beltrami_phase(), save_name_halfz(), correct_lorentz_force(), gamma_function(), hypergeometric2f1(), get_radial_distance(), warning(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical(), curl()

subroutine magnetic/daa_dt(f, df, p)

Magnetic field evolution.

Calculate dA/dt=uxB+3/2 Omega_0 A_y x_dir -eta mu_0 J. Add jxb/rho to momentum equation. Add eta mu_0 j2/rho to entropy equation.

22-nov-01/nils: coded
1-may-02/wolf: adapted for pencil_modular

17-jun-03/ulf: added bx^2, by^2 and bz^2 as separate diagnostics
8-aug-03/axel: introduced B_ext21=1./B_ext**2, and set =1 to avoid div. by 0

12-aug-03/christer: added alpha effect (alpha in the equation above) 26-may-04/axel: ambipolar diffusion added 18-jun-04/axel: Hall term added

Parameters:

f (mx,my,mz,mfarray) [real,inout]
df (mx,my,mz,mvar) [real,inout]
p [pencil_case,inout]

Use :

deriv (der6()), diagnostics, equationofstate (eoscalc()), io (output_pencil()), mpicomm (stop_it()), special (special_calc_magnetic()), sub

Called from:

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

Call to:

timing(), identify_bcs(), der6(), del6fjv(), eta_shell(), fatal_error(), stop_it(), special_calc_magnetic(), erfunc(), calc_tau_aa_exterior(), set_border_magnetic(), daa_dt_meanfield(), max_mn_name(), integrate_mn_name(), sum_mn_name_halfy(), sum_mn_name_halfz(), save_name(), helflux(), curflux(), cross_mn(), dot2_mn(), dot_mn_sv(), yzsum_mn_name_x(), xzsum_mn_name_y(), xysum_mn_name_z(), phizsum_mn_name_r(), yzintegrate_mn_name_x(), xzintegrate_mn_name_y(), xyintegrate_mn_name_z(), phisum_mn_name_rz(), ysum_mn_name_xz(), calc_bthresh(), vecout(), dot_mn_sm(), dot_mn(), der_z(), der2_z(), set_border_initcond(), border_driving(), der_step(), read_snaptime(), update_snaptime(), tridag(), surf_mn_name(), calc_bmx(), calc_bmy(), calc_bmz(), calc_bmzs2(), calc_bmza2(), calc_jmx(), calc_jmy(), calc_jmz(), calc_emxamz3(), calc_embmz(), calc_ambmz(), calc_ambmzh(), calc_jmbmz(), calc_bmxy_rms(), calc_bmz_beltrami_phase(), save_name_halfz(), correct_lorentz_force(), gamma_function(), hypergeometric2f1(), get_radial_distance(), warning(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical(), curl()

subroutine magnetic/time_integrals_magnetic(f, p)

Calculate time_integrals within each pencil (as long as each pencil case p still contains the current data). This routine is now being called at the end of equ.

28-jun-07/axel+mreinhard: coded 24-jun-08/axel: moved call to this routine to the individual pde routines

1-jul-08/axel: moved this part to magnetic

subroutine magnetic/df_diagnos_magnetic(df, p)

calculate diagnostics that involves df Here we calculate <du/dt x b> and <u x db/dt>. The latter is calculated as <divu dai/dt> - <uji daj/dt> This is used in dynamo theory for checking the minimal tau approximation.

10-oct-06/axel: coded

Parameters:

df (mx,my,mz,mvar) [real,in]
p [pencil_case,in]

Use :

diagnostics (sum_mn_name()), sub

Called from:

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

Call to:

cross_mn(), dot_mn_sv(), dot_mn_sm(), dot_mn(), der_z(), der2_z(), set_border_initcond(), fatal_error(), border_driving(), der_step(), stop_it(), read_snaptime(), update_snaptime(), tridag(), surf_mn_name(), calc_bmx(), calc_bmy(), calc_bmz(), calc_bmzs2(), calc_bmza2(), calc_jmx(), calc_jmy(), calc_jmz(), calc_emxamz3(), calc_embmz(), calc_ambmz(), calc_ambmzh(), calc_jmbmz(), calc_bmxy_rms(), calc_bmz_beltrami_phase(), save_name(), save_name_halfz(), correct_lorentz_force(), gamma_function(), hypergeometric2f1(), get_radial_distance(), warning(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical(), curl()

subroutine magnetic/magnetic_after_boundary(f)

Calculate <A>, which is needed for test-field methods.

2-jan-10/axel: adapted from hydro_after_boundary

subroutine magnetic/rescaling_magnetic(f)

Rescale magnetic field by factor rescale_aa,

22-feb-05/axel: coded 10-feb-09/petri: adapted from testfield

subroutine magnetic/sine_avoid_boundary(ampl, f, kr, r0, rn)

Sine field in cylindrical coordinates, used in Armitage 1998

Bz=B0/r * sin(kr*(r-r0))

And 0 outside of the interval r0-rn Code the field and find Aphi through solving the tridiagonal system for

Bz= d/dr Aphi + Aphi/r

-A_(i-1) + A_(i+1) + 2*A_i*dr/r = 2*dr*Bz

05-apr-08/wlad : coded

Parameters:

ampl [real]
f (mx,my,mz,mfarray) [real]
kr [real]
r0 [real]
rn [real]

Use :

general (tridag()), mpicomm (stop_it())

Called from:

register_magnetic(), initialize_magnetic(), init_aa()

Call to:

stop_it(), tridag(), surf_mn_name(), dot_mn(), calc_bmx(), calc_bmy(), calc_bmz(), calc_bmzs2(), calc_bmza2(), calc_jmx(), calc_jmy(), calc_jmz(), calc_emxamz3(), calc_embmz(), calc_ambmz(), calc_ambmzh(), calc_jmbmz(), calc_bmxy_rms(), calc_bmz_beltrami_phase(), save_name(), save_name_halfz(), correct_lorentz_force(), gamma_function(), hypergeometric2f1(), der_step(), fatal_error(), get_radial_distance(), warning(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical(), curl()

subroutine magnetic/curflux_ds(uxb, jj)

current helicity flux (preliminary)

27-nov-03/axel: adapted from helflux

Parameters:

uxb (nx,3) [real,in]
jj (nx,3) [real,in]

Use :

diagnostics

Call to:

surf_mn_name(), dot_mn(), calc_bmx(), calc_bmy(), calc_bmz(), calc_bmzs2(), calc_bmza2(), calc_jmx(), calc_jmy(), calc_jmz(), calc_emxamz3(), calc_embmz(), calc_ambmz(), calc_ambmzh(), calc_jmbmz(), calc_bmxy_rms(), calc_bmz_beltrami_phase(), save_name(), save_name_halfz(), stop_it(), correct_lorentz_force(), gamma_function(), hypergeometric2f1(), der_step(), fatal_error(), get_radial_distance(), warning(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical(), curl()

subroutine magnetic/read_magnetic_init_pars(iomsg)

subroutine magnetic/write_magnetic_init_pars(unit)

subroutine magnetic/read_magnetic_run_pars(iomsg)

subroutine magnetic/write_magnetic_run_pars(unit)

subroutine magnetic/get_slices_magnetic(f, slices)

Write slices for animation of Magnetic variables.

26-jul-06/tony: coded

Parameters:

f (mx,my,mz,mfarray) [real]
slices [slice_data]

Called from:

wvid_prepare(), wvid()

Call to:

calc_bmx(), calc_bmy(), calc_bmz(), calc_bmzs2(), calc_bmza2(), calc_jmx(), calc_jmy(), calc_jmz(), calc_emxamz3(), calc_embmz(), calc_ambmz(), calc_ambmzh(), calc_jmbmz(), calc_bmxy_rms(), calc_bmz_beltrami_phase(), save_name(), save_name_halfz(), stop_it(), correct_lorentz_force(), gamma_function(), hypergeometric2f1(), der_step(), fatal_error(), get_radial_distance(), warning(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical(), curl()

subroutine magnetic/calc_mfield()

calculate mean magnetic field from xy- or z-averages

19-jun-02/axel: moved from print to here: 9-nov-02/axel: corrected bxmy(m,j); it used bzmy instead!

subroutine magnetic/alfven_rphi(b0, f, mode)

Alfven wave propagating on radial direction with field pointing to the phi direction.

Bphi = B0 cos(k r) ==> Az = -1/k B0 sin(k r)

04-oct-06/wlad: coded

Parameters:

b0 [real]
f (mx,my,mz,mfarray) [real]
mode [real]

Called from:

register_magnetic(), initialize_magnetic(), init_aa()

Call to:

stop_it(), correct_lorentz_force(), gamma_function(), hypergeometric2f1(), der_step(), fatal_error(), get_radial_distance(), warning(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical(), curl()

subroutine magnetic/alfven_zconst(f)

Radially variable field pointing in the z direction 4 Balbus-Hawley wavelengths in the vertical direction

Bz=Lz/(8pi)*Omega ==> Aphi = Lz/(8pi) Omega*r/(2-q)

The smoothed case should be general, since it reduces to the non-smoothed for r0_pot=0.

B=C*(r2+r02)^-q ==> Aphi=C/(r*(2-q))*(r2+r02)^(1-q/2)

04-oct-06/wlad: coded

Parameters:: f (mx,my,mz,mfarray) [real]
Use :: gravity (qgshear(), r0_pot()), sharedvariables, mpicomm (stop_it())
Called from:: register_magnetic(), initialize_magnetic(), init_aa()
Call to:: stop_it(), correct_lorentz_force(), gamma_function(), hypergeometric2f1(), der_step(), fatal_error(), get_radial_distance(), warning(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical(), curl()

subroutine magnetic/alfven_rz(b0, f, mode)

Alfven wave propagating on radial direction with field pointing to the z direction.

Bz = B0 cos(k r) ==> Aphi = B0/k sin(k r) + B0/(k2*r)*cos(k r)

04-oct-06/wlad: coded

Parameters:

b0 [real]
f (mx,my,mz,mfarray) [real]
mode [real]

Use :

mpicomm (stop_it())

Called from:

register_magnetic(), initialize_magnetic(), init_aa()

Call to:

stop_it(), gamma_function(), hypergeometric2f1(), der_step(), fatal_error(), get_radial_distance(), warning(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical(), curl()

subroutine magnetic/correct_lorentz_force(f, lfield, const, pblaw)

Correct for the magnetic term in the centrifugal force. The pressure gradient was already corrected in the density and temperature modules

13-nov-08/wlad : coded

Parameters:

f (mx,my,mz,mfarray) [real]
lfield [logical]
const [real]
pblaw [real]

Use :

sub (get_radial_distance())

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(), curflux_ds(), calc_pencils_magnetic(), sine_avoid_boundary(), alfven_rphi(), alfven_zconst()

Call to:

fatal_error(), get_radial_distance(), warning(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical(), curl()

subroutine magnetic/remove_mean_emf(f, df)

Parameters:

f (mx,my,mz,mfarray) [real,in]
df (mx,my,mz,mvar) [real,inout]

Call to:

fatal_error(), remove_mean_emf_cylindrical(), remove_mean_emf_spherical(), curl()

subroutine magnetic/remove_mean_emf_cylindrical(f, df)

Substract mean emf from the radial component of the induction equation. Activated only when large Bz fields and are present keplerian advection. Due to this u_phi x Bz term, the radial component of the magnetic potential develops a divergence that grows linearly in time. Since it is purely divergent, it is okay analytically. But numerically it leads to problems if this divergent grows bigger than the curl, which it does eventually.

This is a cylindrical version of the rtime_phiavg special file.

13-sep-07/wlad: adapted from remove_mean_momenta

Parameters:

f (mx,my,mz,mfarray) [real,in]
df (mx,my,mz,mvar) [real,inout]

Use :

mpicomm (mpiallreduce_sum()), sub (curl())

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(), curflux_ds(), calc_pencils_magnetic(), sine_avoid_boundary(), alfven_rphi(), alfven_zconst(), alfven_rz(), correct_lorentz_force(), remove_mean_emf()

Call to:

curl(), fatal_error()

subroutine magnetic/remove_mean_emf_spherical(f, df)

Substract mean emf from the radial component of the induction equation. Activated only when large Bz fields and keplerian advection are present. Due to this u_phi x Bz term, the radial component of the magnetic potential develops a divergence that grows linearly in time. Since it is purely divergent, it is okay analytically. But numerically it leads to problems if this divergent grows bigger than the curl, which it does eventually.

This is a cylindrical version of the rtime_phiavg special file.

13-sep-07/wlad: adapted from remove_mean_momenta

Parameters:

f (mx,my,mz,mfarray) [real,in]
df (mx,my,mz,mvar) [real,inout]

Use :

mpicomm (mpiallreduce_sum()), sub (curl())

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(), curflux_ds(), calc_pencils_magnetic(), sine_avoid_boundary(), alfven_rphi(), alfven_zconst(), alfven_rz(), correct_lorentz_force(), remove_mean_emf()

Call to:

curl(), fatal_error()

subroutine magnetic/bb_unitvec_shock(f, bb_hat)

Compute unit vector along the magnetic field. Accurate to 2nd order. Tries to avoid division by zero. Taken from http://nuclear.llnl.gov/CNP/apt/apt/aptvunb.html. If anybody knows a more accurate way of doing this, please modify.

16-aug-06/tobi: coded

subroutine magnetic/output_persistent_magnetic(lun)

Write the stored phase and amplitude for the correction of the Beltrami wave forcing

5-apr-08/axel: adapted from output_persistent_forcing

Parameters:: lun [integer]

subroutine magnetic/rprint_magnetic(lreset[, lwrite])

Reads and registers print parameters relevant for magnetic fields.

3-may-02/axel: coded

27-may-02/axel: added possibility to reset list

Parameters:

lreset [logical]
lwrite [logical,in,]

Use :

diagnostics, farraymanager (farray_index_append())

Called from:

register_modules(), initialize_modules(), finalize_modules(), choose_pencils(), rprint_list()