├── .failing-gitlab-ci.yml ├── .gitignore ├── Config.pl ├── Configure.options ├── Configure.pl ├── Copyrights ├── REGISTERED_USER_LICENSE.pdf ├── UM.html ├── UM.long ├── UM.short ├── VAC_UM.long └── VAC_UM.short ├── Doc ├── BNFNamingStand.PDF ├── CFG ├── HTML │ ├── CFG │ ├── csem_home.png │ └── index.html ├── Maple │ ├── CharacteristicSpeed_AnisoPiPe.mpl │ └── CharacteristicSpeed_SemirelAnisoPiPe.mpl ├── RELEASENOTES ├── TECPLOT_PRIMER.pdf ├── Tex │ ├── DESIGN.tex │ ├── HEADER.tex │ ├── Makefile │ ├── NewScaling.eps │ ├── STANDARDS.tex │ ├── USERMANUAL.tex │ ├── adapted_grid.eps │ ├── adapted_grid.pdf │ ├── adaptiveblock.eps │ ├── adaptiveblock.pdf │ ├── bibliography.tex │ ├── datanames.tex │ ├── divb.tex │ ├── equations_physics.tex │ ├── face_numbering.eps │ ├── faces.eps │ ├── fivol.eps │ ├── fivol.pdf │ ├── grid.tex │ ├── hardware_software.tex │ ├── installation.tex │ ├── introduction.tex │ ├── ionosphere.tex │ ├── message_pass_equal.eps │ ├── message_pass_equal.pdf │ ├── message_pass_res_change.eps │ ├── message_pass_res_change.pdf │ ├── numerical_schemes.tex │ ├── output.tex │ ├── param.tex │ ├── proc_dims.eps │ ├── proc_dims.pdf │ ├── quickstart.tex │ ├── runs.tex │ ├── slope.eps │ ├── slope.pdf │ ├── source_code.tex │ ├── spatial.tex │ ├── stable.tex │ ├── sverb.sty │ ├── synopticSW.tex │ ├── temporal.tex │ └── visualization.tex └── index.html ├── LICENSE.txt ├── Makefile ├── Makefile.test ├── PARAM.XML ├── Param ├── 2BODYPLOT │ ├── PARAM.in │ └── TestOutput │ │ ├── REFlog_n000001.log │ │ ├── REFy=0.dat.gz │ │ └── REFy=0.out.gz ├── AMR │ ├── PARAM.in │ ├── PARAM.in.sph │ └── TestOutput │ │ ├── log.log │ │ └── log_sph.log ├── ANISOPRESSURE │ ├── PARAM.in.Alfven │ ├── PARAM.in.anisope.fastwave │ ├── PARAM.in.fastwave │ ├── PARAM.in.soundwave │ └── TestOutput │ │ ├── alfven_ref.out │ │ └── fast_anisope.out ├── B0 │ ├── PARAM.in │ └── PARAM.in.sph ├── COMET │ ├── PARAM.in │ └── TestOutput │ │ └── log_n000001.log ├── COMET3FLUIDSPE │ ├── PARAM.in │ ├── PARAM.in.restartread │ ├── PARAM.in.restartsave │ └── TestOutput │ │ └── log_n000060.log ├── CORONA │ ├── AIA_XRT │ ├── B0local.out │ ├── CR1935_WSO.dat │ ├── CR1962_MDI.dat │ ├── CR2077_GNG.dat │ ├── CR2077_MDI_180.dat │ ├── CR2157_HMI.dat │ ├── EIT_SXR_TEC │ ├── MDI_small.dat │ ├── NEWMAGNETOGRAM.in │ ├── PARAM.in.1Dwedge │ ├── PARAM.in.2Dwedge │ ├── PARAM.in.Awsom │ ├── PARAM.in.Awsom.GPU │ ├── PARAM.in.Awsom.large.GPU │ ├── PARAM.in.AwsomChargeState │ ├── PARAM.in.AwsomR │ ├── PARAM.in.STITCH │ ├── PARAM.in.awsom.bvector │ ├── PARAM.in.bvector │ ├── PARAM.in.heat.cond │ ├── README_los │ ├── RadCoolCorona.dat │ ├── RadCoolCorona_6.0.dat │ ├── RadCoolCorona_8.0.dat │ ├── RadCoolPhoto.dat │ ├── RadCoolPhoto_8.0.dat │ ├── TR.dat │ ├── TR8.dat │ ├── TR_1d.dat │ ├── TestOutput │ │ ├── 1d__mhd_awsomfluids.out │ │ ├── 3d__fux_5_n00000004.dat.gz │ │ ├── 3d__phx_7_n00000004.dat.gz │ │ ├── heatcond_2d_ref.out │ │ ├── log.awsom_large_gpu │ │ ├── log_awsomchargestate.log │ │ ├── log_n000000.bvector │ │ ├── log_n000000.stitch │ │ ├── log_n000001.awsom │ │ ├── log_n000001.awsom.bvector │ │ ├── log_n000001.awsom_gpu │ │ ├── log_n000001.awsomr.td │ │ ├── log_n000001.log │ │ ├── los_dem_1_n0000004.out.gz │ │ ├── los_fux_2_n0000004.out.gz │ │ ├── los_nbi_3_n0000004.out.gz │ │ ├── los_phx_4_n0000004.out.gz │ │ ├── pcl_01__n000001.pcl.gz │ │ └── y=0_nbi_6_n00000004.out.gz │ ├── bvector_cms.dat │ ├── bvector_cr2157.dat │ ├── bvector_cr2157_smoothed.dat │ ├── grid_TR.dat │ ├── grid_awsom.dat │ ├── harmonics-svsm.out.gz │ ├── harmonics_2209_adapt01.dat │ ├── los_Aia.dat │ ├── los_Aia_High.dat │ ├── los_Eit_cor.dat │ ├── los_Eit_photo.dat │ ├── los_EuviA.dat │ ├── los_EuviA_cor.dat │ ├── los_EuviB.dat │ ├── los_EuviB_cor.dat │ ├── los_Sxt.dat │ ├── los_SxtChianti.dat │ ├── los_Xrt.dat │ ├── los_phy.dat │ └── los_tbl.dat ├── CRASH │ ├── 1D.in │ ├── 3D.in │ ├── AMR1.in │ ├── AMR2.in │ ├── AMR3.in │ ├── AMR4.in │ ├── GRAY.in │ ├── GRAYCOND.in │ ├── GRAYCONDHYPRE.in │ ├── GRAYHYPRE.in │ ├── ICS_STANDARD.in │ ├── MULTIGROUP.in │ ├── MULTIGROUPCOND.in │ ├── MULTIGROUPCONDHYPRE.in │ ├── MULTIGROUPHYPRE.in │ ├── PARAM.in.eosgodunov │ ├── PARAM.in.graydiffusion │ ├── PARAM.in.laserpackage │ ├── PARAM.in.mixedcell │ ├── PARAM.in.scaling_hydro │ ├── PARAM.in.scaling_radhydro │ ├── RESTART.in │ ├── RZ.in │ ├── START.in │ ├── TABLE_Au.in │ ├── TABLE_Au_EOS.in │ ├── TABLE_Au_OPACITY.in │ ├── TABLE_Ay.in │ ├── TABLE_Ay_EOS.in │ ├── TABLE_Ay_OPACITY.in │ ├── TABLE_Be.in │ ├── TABLE_Be_EOS.in │ ├── TABLE_Be_OPACITY.in │ ├── TABLE_Pl.in │ ├── TABLE_Pl_EOS.in │ ├── TABLE_Pl_OPACITY.in │ ├── TABLE_Xe.in │ ├── TABLE_Xe_EOS.in │ ├── TABLE_Xe_OPACITY.in │ ├── TestOutput │ │ ├── cut_mixedcell.idl │ │ ├── eosgodunov.out │ │ ├── log_graydiffusion.log │ │ ├── log_laserpackage.log │ │ ├── runlog.eosgodunov │ │ ├── runlog.graydiffusion │ │ ├── runlog.hyades2d │ │ ├── runlog.laserpackage │ │ ├── runlog.levelset │ │ └── runlog_restart.hyades2d │ └── initial_lowrie3.dat.gz ├── CURRENT │ ├── PARAM.in │ ├── TestOutput │ │ └── sat_sat_n000000.sat │ └── sat.dat ├── CYLOUTFLOW │ └── PARAM.in ├── EARTH │ ├── PARAM.in │ ├── PARAM.in.2D │ ├── PARAM.in.L1toBC │ ├── PARAM.in.gpu │ ├── PARAM.in.large.gpu │ ├── PARAM.in.spherical │ ├── TestOutput │ │ ├── L1toBC.log │ │ ├── SPHlog_n000001.log │ │ ├── SPHrunlog │ │ ├── earth_cartesian_mag_grid.out.gz │ │ ├── earth_cartesian_ref.log │ │ ├── earth_cartesian_ref.out │ │ └── earth_large_ref.log │ └── imf19980504.dat ├── EUROPA │ ├── PARAM.in.2fluids │ ├── PARAM.in.3fluids │ ├── TestOutput │ │ ├── trj_satellite_2fpe_e4_n00000010.sat.gz │ │ └── trj_satellite_3fpe_e4_n00000100.sat.gz │ └── satellite_e4.dat ├── FIVEMOMENT │ ├── PARAM.in.alfven │ ├── PARAM.in.langmuir │ ├── PARAM.in.light │ ├── PARAM.in.shock │ └── TestOutput │ │ ├── alfven_ref.out │ │ ├── langmuir_ref.out │ │ └── shock_ref.out ├── FLUXEMERGENCE │ ├── Grid.1D │ ├── Grid.3D │ ├── PARAM.in.FluxEmergence │ ├── PARAM.in.FluxEmergence.1D │ ├── PARAM.in.FluxEmergence.3D │ └── TestOutput │ │ └── log_fluxemergence.log ├── GANYMEDE │ ├── Grid │ ├── PARAM.in │ └── TestOutput │ │ └── REFlog_n000001.log ├── GEMRECONNECTION │ ├── PARAM.in.MhdHypPe │ ├── PARAM.in.double.CS.MHDEPIC │ └── PARAM.in.sixmoment ├── HELIOSPHERE │ ├── Arcade │ ├── CME_fluxrope │ ├── EARTH_TRAJ.in │ ├── EARTH_TRAJ.in.CR │ ├── Heliosphere │ └── SATELLITE.in ├── JUPITER │ ├── Cassini_nov18_IMF │ ├── Ideal_Corotation │ ├── Mag_SSSS_0.25Rj │ ├── Mag_SSSS_0.2Rj │ ├── Mag_SSSS_0.4Rj │ ├── Mag_SSSS_0.5Rj │ ├── Mag_SSSS_0.8Rj │ ├── Mag_SSSS_1.0Rj │ ├── Magnetosphere │ ├── North_IMF │ ├── PARAM.in │ ├── PARAM.in.Mag_SSSS │ ├── Parker_IMF │ ├── Real_Corotation │ ├── South_IMF │ └── TestOutput │ │ └── log_n000001.log ├── MAGNETOMETER │ ├── PARAM.in │ ├── TestOutput │ │ ├── bx0_t00000100.out │ │ ├── eqb_t00000100.dat │ │ ├── eqb_t00000100.out │ │ └── magnetometers_n00000000.mag │ └── magin.dat ├── MARS │ ├── Grid │ ├── PARAM.in │ ├── PARAM.in.restartread │ ├── PARAM.in.restartsave │ ├── PARAM.in.ta │ ├── TestOutput │ │ ├── log_n000001.log │ │ └── log_ta.log │ ├── grid_stretch.dat │ ├── imf.dat │ └── marsmgsp.txt ├── MARSFLUIDS │ ├── PARAM.in │ ├── PARAM.in.restartread │ ├── PARAM.in.restartsave │ └── TestOutput │ │ ├── log_n000001.log │ │ └── runlog ├── MERCURY │ ├── PARAM.in │ └── TestOutput │ │ └── REFlog_n000000.log ├── MHDNONCONS │ ├── PARAM.in │ └── TestOutput │ │ └── cut_mhd_t10m.out ├── MOONIMPACT │ ├── PARAM.in │ ├── PARAM.in.restartread │ ├── PARAM.in.restartsave │ └── TestOutput │ │ ├── log_n000000.log │ │ └── log_n000300.log ├── MULTIFLUID │ ├── PARAM.in │ └── TestOutput │ │ ├── cut_t25.out │ │ └── runlog ├── MULTIION │ ├── PARAM.in │ └── TestOutput │ │ ├── 1d_t2.out │ │ └── runlog ├── OUTERHELIO │ ├── ChargeExchangeRate_LS.dat │ ├── ChargeExchangeRate_MT.dat │ ├── ElectronImpactLookup.py │ ├── ElectronImpactRate.dat │ ├── ExchangeRateLookup.py │ ├── Grid │ ├── INTERPOLATE.in │ ├── MagIntensity.dat │ ├── PARAM.in │ ├── PARAM.in.1d │ ├── PARAM.in.2D │ ├── PARAM.in.awsom │ ├── PARAM.in.awsom.restart │ ├── PARAM.in.pui │ ├── PARAM.in.pui.restart │ ├── PARAM.in.pui1d │ ├── PARAM.in.restart │ ├── SourceIntLookup.pro │ ├── TestOutput │ │ ├── interpolation_ref.dat │ │ ├── log_n000000_1d.log │ │ ├── log_n000000_pui1d.log │ │ ├── log_n000001_2D.log │ │ ├── log_n000010.log │ │ ├── log_n000010_awsom.log │ │ ├── log_n000010_pui.log │ │ ├── log_n000040_awsom_restart.log │ │ ├── runlog │ │ ├── runlog_2D │ │ ├── runlog_pui │ │ ├── runlog_pui1d │ │ └── z=0_2D.out.gz │ ├── oh2d_l1_test.dat │ ├── oh2d_stereoa_test.dat │ ├── solarwind2D.dat │ └── testsat.dat ├── REGION │ ├── PARAM.in.2d │ ├── PARAM.in.cyl │ ├── PARAM.in.cyl_lnr │ └── TestOutput │ │ ├── 2d.out.gz │ │ ├── 2d.pdf │ │ ├── cyl.out.gz │ │ └── cyl_lnr.out.gz ├── ROSETTA │ ├── CG_MOC.bdf.gz │ ├── PARAM.in.fluids.all │ ├── PARAM.in.hd │ ├── PARAM.in.hd.test │ └── TestOutput │ │ ├── CGfluids_log_n000000.log │ │ └── hd_log_n000000.log ├── ROTATINGFRAME │ └── PARAM.in ├── RUNGEKUTTA │ └── PARAM.in ├── SATURN │ ├── GRID │ ├── Ideal_Corotation │ ├── Mag_SSSS_0.2Rs │ ├── Mag_SSSS_0.4Rs │ ├── Mag_SSSS_0.8Rs │ ├── Magnetosphere │ ├── Magnetosphere.nominal │ ├── North_IMF │ ├── PARAM.in │ ├── Parker_IMF │ ├── Real_Corotation │ ├── South_IMF │ └── TestOutput │ │ └── log_n000001.log ├── SHOCKTUBE │ ├── PARAM.in │ ├── PARAM.in.FiveMoment │ ├── PARAM.in.HallTest │ ├── PARAM.in.HybridTest │ ├── PARAM.in.Orszag-Tang │ ├── PARAM.in.TwofluidTest │ ├── PARAM.in.alfven_sa │ ├── PARAM.in.bx0 │ ├── PARAM.in.fast_wave │ ├── PARAM.in.fast_wave_2d │ ├── PARAM.in.fast_wave_multiGPU │ ├── PARAM.in.hillvortex │ ├── PARAM.in.hyp │ ├── PARAM.in.hyp.1d │ ├── PARAM.in.mhd_blastwave │ ├── PARAM.in.partsteady │ ├── PARAM.in.rayleigh_taylor │ ├── PARAM.in.rotation │ ├── PARAM.in.roundcube │ ├── PARAM.in.shockramp │ ├── PARAM.in.sphalfven │ ├── PARAM.in.sphericalwedge │ ├── PARAM.in.timewarp1d │ ├── PARAM.in.timewarp2d │ ├── TestOutput │ │ ├── bx0_ref.out │ │ ├── fast_wave_ref.out │ │ ├── fastwave_2d_ref.out │ │ ├── hall_ref.out │ │ ├── mhdhyp_1d_ref.out │ │ ├── mhdhyp_ref.out │ │ ├── mhdhyp_runlog │ │ ├── partsteady_ref.out │ │ ├── partsteady_runlog │ │ ├── runlog │ │ ├── shockramp.out.gz │ │ ├── shockround.log │ │ ├── timewarp1d.out │ │ ├── timewarp2d.out.gz │ │ └── twofluid_ref.out │ ├── briowu │ ├── grid256 │ ├── grid64 │ ├── rotatingM.pro │ └── ryu ├── SIXMOMENT │ ├── PARAM.in.alfven │ ├── PARAM.in.fast │ ├── PARAM.in.light │ └── PARAM.in.shock ├── TITAN │ ├── Grid │ ├── PARAM.in │ ├── PARAM.in.restartread │ ├── PARAM.in.restartsave │ ├── TestOutput │ │ ├── log_n000001.log │ │ └── runlog │ └── TitanInput.tgz ├── VENUS │ ├── Grid │ ├── PARAM.in │ ├── PARAM.in.restartread │ ├── PARAM.in.restartsave │ └── TestOutput │ │ ├── log_n000001.log │ │ └── runlog └── VISCOSITY │ ├── PARAM.in │ └── TestOutput │ └── z=0_ref.outs ├── README.md ├── Scripts ├── FixNonconservative.pl ├── FixUnitVar.pl ├── IDL │ ├── FixEndian.pl │ ├── call_euv.pro │ ├── euv.csh │ └── euv.pro ├── Scaling.job.h2ologin ├── Scaling.job.hera ├── Scaling.job.jaguarpf-ext ├── Scaling.job.pfe ├── Scaling.job.ubgl ├── Scaling.job.yslogin ├── Scaling.pl ├── TEC │ ├── pTEC │ ├── ppA │ ├── ppAr │ └── ppgz ├── cat_sat_log.pl └── shift_satellite.pl ├── TestParam.pl ├── src ├── Makefile ├── Makefile.RULES.all ├── Makefile.def ├── ModAMR.f90 ├── ModAdvance.f90 ├── ModAdvanceExplicit.f90 ├── ModAdvectPoints.f90 ├── ModB0.f90 ├── ModBatlInterface.f90 ├── ModBatsrusMethods.f90 ├── ModBatsrusUtility.f90 ├── ModBlockData.f90 ├── ModBorisCorrection.f90 ├── ModBoundaryGeometry.f90 ├── ModBuffer.f90 ├── ModCalcSource.f90 ├── ModCellBoundary.f90 ├── ModCellGradient.f90 ├── ModCharacteristic.f90 ├── ModChromosphere.f90 ├── ModCleanDivB.f90 ├── ModCoarseAxis.f90 ├── ModConservative.f90 ├── ModConserveFlux.f90 ├── ModConstrainDivB.f90 ├── ModCoronalHeating.f90 ├── ModCurrent.f90 ├── ModElectricField.f90 ├── ModEnergy.f90 ├── ModExtraVariables.f90 ├── ModFaceBoundary.f90 ├── ModFaceFlux.f90 ├── ModFaceGradient.f90 ├── ModFaceValue.f90 ├── ModFieldLineThread.f90 ├── ModFieldLineThreadNew.f90 ├── ModFieldTrace.f90 ├── ModFieldTraceFast.f90 ├── ModFixAxisCells.f90 ├── ModGeometry.f90 ├── ModGroundMagPerturb.f90 ├── ModHallResist.f90 ├── ModHdf5_empty.f90 ├── ModHdf5_orig.f90 ├── ModHeatConduction.f90 ├── ModHeatFluxCollisionless.f90 ├── ModIO.f90 ├── ModIeCoupling.f90 ├── ModImCoupling.f90 ├── ModImplHypre_empty.f90 ├── ModImplHypre_orig.f90 ├── ModImplicit.f90 ├── ModIonElectron.f90 ├── ModLaserHeating.f90 ├── ModLoadBalance.f90 ├── ModLocalTimeStep.f90 ├── ModMain.f90 ├── ModMessagePass.f90 ├── ModMultiFluid.f90 ├── ModMultiIon.f90 ├── ModNodes.f90 ├── ModOptimizeParam_orig.f90 ├── ModPIC.f90 ├── ModPUI.f90 ├── ModParallel.f90 ├── ModPartImplicit.f90 ├── ModPartSteady.f90 ├── ModParticleFieldLine.f90 ├── ModParticleMover.f90 ├── ModParticles.f90 ├── ModPhysicalFlux.f90 ├── ModPhysics.f90 ├── ModPlotBox.f90 ├── ModPlotShell.f90 ├── ModPlotShock.f90 ├── ModPointImplicit.f90 ├── ModProjectDivB.f90 ├── ModRadDiffusion.f90 ├── ModRadiativeCooling.f90 ├── ModRadioWaveImage.f90 ├── ModResistivity.f90 ├── ModRestartFile.f90 ├── ModSatelliteFile.f90 ├── ModSemiImplicit.f90 ├── ModSetInitialCondition.f90 ├── ModSetParameters.f90 ├── ModSingleFluid.f90 ├── ModSize.f90 ├── ModSolarwind.f90 ├── ModSpectrum.f90 ├── ModStreamAligned.f90 ├── ModThreadedLC.f90 ├── ModThreadedLCNew.f90 ├── ModTimeStepControl.f90 ├── ModTimewarp.f90 ├── ModTurbulence.f90 ├── ModUaCoupling.f90 ├── ModUpdateState.f90 ├── ModUpdateStateFast.f90 ├── ModUserDefault.f90 ├── ModUserEmpty.f90 ├── ModUserInterface.f90 ├── ModViscosity.f90 ├── ModWaves.f90 ├── ModWriteLogSatFile.f90 ├── ModWritePlot.f90 ├── ModWritePlotIdl.f90 ├── ModWritePlotLos.f90 ├── ModWritePlotRadioWave.f90 ├── ModWriteProgress.f90 ├── ModWriteTecplot.f90 ├── main.f90 ├── show_git_info.h.orig ├── user_interface.f90 ├── user_module.h └── vector_functions.h ├── srcEquation ├── ModEquationAwsom.f90 ├── ModEquationAwsomAnisoPi.f90 ├── ModEquationAwsomAnisoPiSA.f90 ├── ModEquationAwsomChargeState.f90 ├── ModEquationAwsomFluids.f90 ├── ModEquationAwsomFluidsAnisoPi.f90 ├── ModEquationAwsomSA.f90 ├── ModEquationAwsomSAWdiff.f90 ├── ModEquationAwsomWdiff.f90 ├── ModEquationComet3FluidsPe.f90 ├── ModEquationCometCG3FluidsPe.f90 ├── ModEquationCometCG3FluidsPeHyp.f90 ├── ModEquationCrash.f90 ├── ModEquationCrashTe.f90 ├── ModEquationFiveMoment.f90 ├── ModEquationFiveMomentHyp.f90 ├── ModEquationFourIonFluidPe.f90 ├── ModEquationHd.f90 ├── ModEquationHdCrash.f90 ├── ModEquationHdEos.f90 ├── ModEquationHdEosRad.f90 ├── ModEquationHdRadCrash.f90 ├── ModEquationIonsHNO.f90 ├── ModEquationMarsFluids.f90 ├── ModEquationMarsFluidsPe.f90 ├── ModEquationMarsFluidsSw.f90 ├── ModEquationMhd.f90 ├── ModEquationMhdAnisoP.f90 ├── ModEquationMhdComet.f90 ├── ModEquationMhdCorona.f90 ├── ModEquationMhdCrash.f90 ├── ModEquationMhdEos.f90 ├── ModEquationMhdEosRad.f90 ├── ModEquationMhdEuropa2FluidsPe.f90 ├── ModEquationMhdEuropa3FluidsPe.f90 ├── ModEquationMhdFluidsPeAniso.f90 ├── ModEquationMhdHd.f90 ├── ModEquationMhdHpOp.f90 ├── ModEquationMhdHyp.f90 ├── ModEquationMhdHypPe.f90 ├── ModEquationMhdMars.f90 ├── ModEquationMhdMarsPe.f90 ├── ModEquationMhdNonCons.f90 ├── ModEquationMhdPe.f90 ├── ModEquationMhdPeAniso.f90 ├── ModEquationMhdPeAnisoPi.f90 ├── ModEquationMhdPeSA.f90 ├── ModEquationMhdPw.f90 ├── ModEquationMhdSA.f90 ├── ModEquationMhdSaturn3sp.f90 ├── ModEquationMhdSwIono.f90 ├── ModEquationMhdTitan.f90 ├── ModEquationMhdTitanPe.f90 ├── ModEquationMhdWaves.f90 ├── ModEquationMhdWavesHyp.f90 ├── ModEquationMhdWavesHypPe.f90 ├── ModEquationMultiIon.f90 ├── ModEquationMultiIonPe.f90 ├── ModEquationOuterHelio.f90 ├── ModEquationOuterHelio2PUIPe.f90 ├── ModEquationOuterHelio2d.f90 ├── ModEquationOuterHelioAwsom.f90 ├── ModEquationOuterHelioAwsomPuiBin.f90 ├── ModEquationOuterHelioPUI.f90 ├── ModEquationOuterHelioPUIPe.f90 ├── ModEquationScalar.f90 ├── ModEquationSixMoment.f90 ├── ModEquationSixMomentHyp.f90 ├── ModEquationSwh.f90 ├── ModEquationSwh2PuiPe.f90 ├── ModEquationSwhPui.f90 ├── ModEquationSwhPuiPe.f90 └── ModEquationThreeIonFluidPe.f90 ├── srcInterface ├── GM_couple_ie.f90 ├── GM_couple_ih.f90 ├── GM_couple_im.f90 ├── GM_couple_pc.f90 ├── GM_couple_ps.f90 ├── GM_couple_pt.f90 ├── GM_couple_pw.f90 ├── GM_couple_rb.f90 ├── GM_couple_ua.f90 ├── GM_wrapper.f90 ├── Makefile ├── ModGridDescriptor.f90 ├── read_ih_buffer.f90 └── read_pw_buffer.f90 ├── srcPostProc ├── CON_geopack_internal.f90 ├── ConvertRestart.f90 ├── Makefile ├── Makefile.RULES.all ├── ModUT.f90 ├── PHOTOEXC_input.pro ├── SPECTRUM_input.pro ├── chianti_dblpre.ioneq ├── earth_traj.f90 ├── interpolate_output.f90 ├── select_snapshot.f90 ├── spectrum.f90 ├── test-SPECTRUM.in ├── test-spectrum.ref.gz ├── test-spectrum2-input.info ├── test-spectrum2-input.out ├── test-spectrum2-input.tree ├── time_convert.f90 └── time_to_cr_no.f90 └── srcUser ├── ModUserAwsom.f90 ├── ModUserCCMC.f90 ├── ModUserComet1Sp.f90 ├── ModUserComet3FluidsPe.f90 ├── ModUserComet6Sp.f90 ├── ModUserCometCG.f90 ├── ModUserCometCGfluids.f90 ├── ModUserCometNeutralFluids.f90 ├── ModUserCylOutflow.f90 ├── ModUserEarthXray.f90 ├── ModUserEe.f90 ├── ModUserFluxRope.f90 ├── ModUserGemReconnect.f90 ├── ModUserJupiter.f90 ├── ModUserKelvinHelmholtz.f90 ├── ModUserLogAdvection.f90 ├── ModUserMars.f90 ├── ModUserMarsFluids.f90 ├── ModUserMercury.f90 ├── ModUserMoonImpact.f90 ├── ModUserNonGyro.f90 ├── ModUserOuterHelio.f90 ├── ModUserOuterHelio2d.f90 ├── ModUserPointImplicit.f90 ├── ModUserSaturn.f90 ├── ModUserStretchedDipole.f90 ├── ModUserSwIono.f90 ├── ModUserSwarm.f90 ├── ModUserSwitchback.f90 ├── ModUserTitan.f90 ├── ModUserVenus.f90 ├── ModUserWaveReflection.f90 └── ModUserWaves.f90 /.failing-gitlab-ci.yml: -------------------------------------------------------------------------------- 1 | #image: ubuntu:latest 2 | image: ubuntu:22.04 3 | 4 | variables: 5 | GCC_V: 9 6 | OMPI_ALLOW_RUN_AS_ROOT: 1 7 | OMPI_ALLOW_RUN_AS_ROOT_CONFIRM: 1 8 | 9 | stages: 10 | - build 11 | - test 12 | 13 | install: 14 | stage: build 15 | before_script: 16 | - 'which ssh-agent || ( apt-get update -y && apt-get install openssh-client git -y )' 17 | - apt-get install -y make 18 | - eval $(ssh-agent -s) 19 | - echo "$SSH_PRIVATE_KEY_SHARE" | tr -d '\r' | ssh-add - 20 | - echo "$SSH_PRIVATE_KEY_UTIL" | tr -d '\r' | ssh-add - 21 | - echo "$SSH_PRIVATE_KEY_BATL" | tr -d '\r' | ssh-add - 22 | - mkdir -p ~/.ssh 23 | - chmod 700 ~/.ssh 24 | - gitlab_hostname=$(echo "${CI_REPOSITORY_URL}" | sed -e 's|https\?://gitlab-ci-token:.*@||g' | sed -e 's|/.*||g') 25 | - ssh-keyscan "${gitlab_hostname}" >> ~/.ssh/known_hosts 26 | - chmod 644 ~/.ssh/known_hosts 27 | script: 28 | - perl -e '$a=`find * -type f -size +1M`; die "Files exceeding 1M size limit $a\n" if $a; exit 0' 29 | - ./Config.pl -compiler=gfortran -install 30 | artifacts: 31 | paths: 32 | - share 33 | - util 34 | - srcBATL 35 | expire_in: 1 hour 36 | 37 | validity: 38 | stage: test 39 | before_script: 40 | - apt-get update -y && apt-get install --yes gcc gfortran g++ make 41 | - DEBIAN_FRONTEND=noninteractive apt-get install --yes openmpi-bin libopenmpi-dev 42 | - chmod +x share/Scripts/CheckEmptyFile.sh 43 | script: 44 | - ./Config.pl -O0 -compiler=gfortran -install 45 | - make -j test_hallmhd NP=1 46 | - ./share/Scripts/CheckEmptyFile.sh test_hallmhd.diff 47 | 48 | # vim: softtabstop=2 shiftwidth=2 49 | -------------------------------------------------------------------------------- /.gitignore: -------------------------------------------------------------------------------- 1 | *.mod 2 | *.pyc 3 | *.o 4 | *.a 5 | *.diff 6 | *.exe 7 | *~ 8 | *_orig_ 9 | *.speed 10 | *.protex 11 | *.aux 12 | *.dvi 13 | *.toc 14 | *.xmltex 15 | *# 16 | .DS_Store 17 | Doc/Tex/*.log 18 | Doc/USERMANUAL.pdf 19 | Doc/DESIGN.pdf 20 | Doc/STANDARDS.pdf 21 | 22 | gitinfo.txt 23 | src/show_git_info.h 24 | src/ModOptimizeParam.f90 25 | 26 | run*/ 27 | Makefile.conf 28 | Makefile.def 29 | Makefile.RULES 30 | Makefile.DEPEND 31 | stand_alone.f90 32 | BATL_size.f90 33 | ModUser.f90 34 | ModUser.f90.safe 35 | ModEquation.f90 36 | ModEquation.f90.safe 37 | ModImplHypre.f90 38 | ModHdf5.f90 39 | Tmp_* 40 | util/ 41 | share/ 42 | srcBATL/ 43 | srcUserExtra/ 44 | data 45 | dataCRASH 46 | srcPostProc/*.log 47 | srcPostProc/*.out 48 | srcPostProc/*.dat 49 | srcPostProc/SPECTRUM.in 50 | -------------------------------------------------------------------------------- /Configure.options: -------------------------------------------------------------------------------- 1 | # Copyright (C) 2002 Regents of the University of Michigan, portions used with permission 2 | # For more information, see http://csem.engin.umich.edu/tools/swmf 3 | #^CFG FILE CONFIGURE 4 | # 5 | # Set configure to true to keep the build configurable 6 | # 7 | CONFIGURE OFF 8 | 9 | 10 | # ^CFG IF NOT COPYRIGHT BEGIN 11 | # Inserting copyright 12 | # 13 | COPYRIGHT OFF ^CFG END COPYRIGHT 14 | 15 | # ^CFG IF DOC BEGIN 16 | # Including documentation 17 | # 18 | DOC OFF 19 | DOCHTML OFF NEEDS DOC ^CFG IF DOCHTML 20 | 21 | # ^CFG IF NOT REMOVEDOCTEX BEGIN 22 | # 23 | # Making documentation for the configured code 24 | # 25 | MAKEPDF OFF NEEDS DOC 26 | MAKEHTML OFF NEEDS DOCHTML ^CFG IF DOCHTML 27 | # 28 | # Removing the source Doc/Tex files 29 | # 30 | REMOVEDOCTEX OFF NEEDS DOC 31 | # ^CFG END REMOVEDOCTEX 32 | # ^CFG END DOC 33 | 34 | # ^CFG IF TESTING BEGIN 35 | # Including test suite 36 | # 37 | TESTING OFF ^CFG END TESTING 38 | 39 | # ^CFG IF DEBUGGING BEGIN 40 | # 41 | # Temporary or Debugging parts of source code 42 | # 43 | DEBUGGING OFF ^CFG END DEBUGGING 44 | 45 | # 46 | # Physics modules 47 | # 48 | IONOSPHERE ON #^CFG IF IONOSPHERE 49 | RCM OFF NEEDS RAYTRACE #^CFG IF RCM 50 | 51 | # 52 | # Extra physics 53 | # 54 | SECONDBODY OFF #^CFG IF SECONDBODY 55 | 56 | # 57 | # Numerical schemes 58 | # 59 | IMPLICIT OFF #^CFG IF IMPLICIT 60 | 61 | DIVBDIFFUSE OFF #^CFG IF DIVBDIFFUSE 62 | PROJECTION OFF #^CFG IF PROJECTION 63 | CONSTRAINB OFF NEEDS PROJECTION #^CFG IF CONSTRAINB 64 | 65 | RAYTRACE OFF #^CFG IF RAYTRACE 66 | 67 | SIMPLEBORIS OFF #^CFG IF SIMPLEBORIS 68 | BORISCORR ON #^CFG IF BORISCORR 69 | 70 | RUSANOVFLUX ON #^CFG IF RUSANOVFLUX 71 | LINDEFLUX OFF #^CFG IF LINDEFLUX 72 | AWFLUX OFF #^CFG IF AWFLUX 73 | HLLDFLUX OFF #^CFG IF HLLDFLUX 74 | ROEFLUX OFF #^CFG IF ROEFLUX 75 | DISSFLUX OFF #^CFG IF DISSFLUX 76 | 77 | 78 | 79 | # #^CFG IF DYNAMIC BEGIN 80 | # Dynamic allocation of big arrays 81 | # 82 | DYNAMIC OFF #^CFG END DYNAMIC 83 | -------------------------------------------------------------------------------- /Configure.pl: -------------------------------------------------------------------------------- 1 | #!/usr/bin/perl 2 | # Copyright (C) 2002 Regents of the University of Michigan, portions used with permission 3 | # For more information, see http://csem.engin.umich.edu/tools/swmf 4 | #^CFG FILE CONFIGURE 5 | # Find and run share/Scripts/Configure.pl -c=CFG 6 | use strict; 7 | 8 | my $ConfigureScript = "share/Scripts/Configure.pl"; 9 | my @ConfigureScript = ($ConfigureScript, 10 | "../../$ConfigureScript", 11 | "../../../$ConfigureScript"); 12 | 13 | my $Script; 14 | foreach $Script (@ConfigureScript){ 15 | if(-f $Script){ 16 | my $Command = "$Script -c=CFG ".join(' ',@ARGV); 17 | exec($Command) or die "Could not execute $Command\n"; 18 | } 19 | } 20 | 21 | die "Could not find $ConfigureScript\n"; 22 | -------------------------------------------------------------------------------- /Copyrights/REGISTERED_USER_LICENSE.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Copyrights/REGISTERED_USER_LICENSE.pdf -------------------------------------------------------------------------------- /Copyrights/UM.html: -------------------------------------------------------------------------------- 1 | COPYRIGHT (c) 2002
2 | THE REGENTS OF THE UNIVERSITY OF MICHIGAN
3 | ALL RIGHTS RESERVED 4 | 5 |

6 | PERMISSION IS GRANTED TO A SINGLE USER TO USE THIS SOFTWARE FOR 7 | NONCOMMERCIAL EDUCATION AND RESEARCH PURPOSES. THE USER IS GIVEN 8 | PERMISSION TO INSTALL AND RUN THE SOFTWARE ON MULTIPLE PLATFORMS. NO 9 | PERMISSION IS GIVEN TO COPY OR REDISTRIBUTE THIS SOFTWARE IN ANY FORM 10 | INCLUDING, BUT NOT LIMITED TO, COPYING SECTIONS OF THE SOURCE CODE OR 11 | DOCUMENTATION AND DISTRIBUTING THE CODE TO OTHER USERS. ALTHOUGH 12 | SOURCE CODE IS DISTRIBUTED, USERS ARE RESTRICTED FROM MODIFYING THE 13 | SOURCE CODE IN ANYWAY WITHOUT PRIOR CONSENT OF THE CENTER FOR SPACE 14 | ENVIRONMENT MODELING AT THE UNIVERSITY OF MIGHIGAN. 15 | 16 |

17 | THIS SOFTWARE IS PROVIDED AS IS, WITHOUT REPRESENTATION AS TO ITS 18 | FITNESS FOR ANY PURPOSE, AND WITHOUT WARRANTY OF ANY KIND, EITHER 19 | EXPRESS OR IMPLIED, INCLUDING WITHOUT LIMITATION THE IMPLIED WARRANTIES 20 | OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE. THE REGENTS OF 21 | THE UNIVERSITY OF MICHIGAN SHALL NOT BE LIABLE FOR ANY DAMAGES, 22 | INCLUDING SPECIAL, INDIRECT, INCIDENTAL, OR CONSEQUENTIAL DAMAGES, WITH 23 | RESPECT TO ANY CLAIM ARISING OUT OF OR IN CONNECTION WITH THE USE OF 24 | THE SOFTWARE, EVEN IF IT HAS BEEN OR IS HEREAFTER ADVISED OF THE 25 | POSSIBILITY OF SUCH DAMAGES. 26 | 27 | -------------------------------------------------------------------------------- /Copyrights/UM.long: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Copyrights/UM.long -------------------------------------------------------------------------------- /Copyrights/UM.short: -------------------------------------------------------------------------------- 1 | This code is a copyright protected software. 2 | (c) 2002- University of Michigan 3 | 4 | -------------------------------------------------------------------------------- /Copyrights/VAC_UM.long: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Copyrights/VAC_UM.long -------------------------------------------------------------------------------- /Copyrights/VAC_UM.short: -------------------------------------------------------------------------------- 1 | This code is a copyright protected software. 2 | (c) 2002- University of Michigan 3 | 4 | -------------------------------------------------------------------------------- /Doc/BNFNamingStand.PDF: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Doc/BNFNamingStand.PDF -------------------------------------------------------------------------------- /Doc/CFG: -------------------------------------------------------------------------------- 1 | DOC 2 | -------------------------------------------------------------------------------- /Doc/HTML/CFG: -------------------------------------------------------------------------------- 1 | DOCHTML 2 | -------------------------------------------------------------------------------- /Doc/HTML/csem_home.png: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Doc/HTML/csem_home.png -------------------------------------------------------------------------------- /Doc/HTML/index.html: -------------------------------------------------------------------------------- 1 | 2 | 4 | 5 | BATS-R-US MANUALS 6 | 7 | 8 |

9 |

BATS-R-US MANUALS

10 |
11 | 12 |

DESIGN DOCUMENT

13 | 14 |

USER MANUAL

15 | 16 |

CODING STANDARDS

17 | 18 | 19 | -------------------------------------------------------------------------------- /Doc/TECPLOT_PRIMER.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Doc/TECPLOT_PRIMER.pdf -------------------------------------------------------------------------------- /Doc/Tex/DESIGN.tex: -------------------------------------------------------------------------------- 1 | % Copyright (C) 2002 Regents of the University of Michigan, portions used with permission 2 | % For more information, see http://csem.engin.umich.edu/tools/swmf 3 | \documentclass[twoside,10pt]{book} 4 | 5 | \input HEADER 6 | 7 | \title{\BATSRUS\ Design \\ \large Code Version 7.4.4} 8 | \author{K.C. Hansen, G\'abor T\'oth, Aaron Ridley, Darren DeZeeuw} 9 | 10 | \begin{document} 11 | 12 | \pagestyle{headings} 13 | 14 | \maketitle 15 | 16 | \tableofcontents 17 | 18 | % Files for chapter on equations and physics 19 | %\input{equations_physics.tex} 20 | 21 | % Files for chapter on numerical discretization 22 | \chapter{Numerical Discretization \label{chapter:numerics}} 23 | \input{grid.tex} 24 | \input{temporal.tex} 25 | \input{spatial.tex} 26 | \input{divb.tex} 27 | 28 | % Files for chapter on modifying the source code 29 | %\input{source_code.tex} 30 | 31 | \end{document} 32 | -------------------------------------------------------------------------------- /Doc/Tex/HEADER.tex: -------------------------------------------------------------------------------- 1 | % Copyright (C) 2002 Regents of the University of Michigan, portions used with permission 2 | % For more information, see http://csem.engin.umich.edu/tools/swmf 3 | \usepackage{times} 4 | \usepackage{alltt} 5 | \usepackage{amsmath} 6 | \usepackage{epsfig} 7 | 8 | % use these lengths for a format for an offset margin 9 | % this is useful for binding or for 3 hole punching 10 | %\setlength{\oddsidemargin}{.25 in} 11 | %\setlength{\evensidemargin}{-.55 in} 12 | %\setlength{\textwidth}{6.5 in} 13 | %\setlength{\topmargin}{0 in} 14 | %\setlength{\textheight}{8.5 in} 15 | 16 | % use these lengths for a more uniform margin 17 | % this format is more pleasing for stapling 18 | \setlength{\oddsidemargin}{-.1 in} 19 | \setlength{\evensidemargin}{0.0 in} 20 | \setlength{\textwidth}{6.5 in} 21 | \setlength{\topmargin}{0 in} 22 | \setlength{\textheight}{8.5 in} 23 | 24 | \renewcommand{\deg}{^{\circ}} 25 | \newcommand{\bA}{\mathbf{A}} 26 | \newcommand{\bb}{\mathbf{b}} 27 | \newcommand{\bB}{\mathbf{B}} 28 | \newcommand{\divb}{\nabla\cdot\bb} 29 | \newcommand{\divB}{\nabla\cdot\bB} 30 | \newcommand{\bE}{\mathbf{E}} 31 | \newcommand{\bv}{\mathbf{v}} 32 | \newcommand{\bu}{\mathbf{u}} 33 | \newcommand{\bx}{\mathbf x} 34 | \newcommand{\bn}{\mathbf n} 35 | \newcommand{\bF}{\mathbf F} 36 | \newcommand{\DU}{\Delta U} 37 | \newcommand{\minmod}{\mathrm{minmod}} 38 | \newcommand{\BATSRUS}{BATS-R-US} 39 | -------------------------------------------------------------------------------- /Doc/Tex/Makefile: -------------------------------------------------------------------------------- 1 | # Copyright (C) 2002 Regents of the University of Michigan, portions used with permission 2 | # For more information, see http://csem.engin.umich.edu/tools/swmf 3 | include ../../Makefile.def 4 | include ${DIR}/share/build/Makefile.doc 5 | 6 | help: 7 | @echo " Make manuals for BATSRUS" 8 | @echo "" 9 | @echo "Targets:" 10 | @echo "" 11 | @echo "PDF - produces USERMANUAL.pdf DESIGN.pdf STANDARDS.pdf" 12 | @echo "" 13 | @echo "clean - rm -rf *.aux *.log *.toc *.dvi *~" 14 | @echo "" 15 | @echo "cleanpdf - clean + remove PDF manuals" 16 | 17 | USERMANUAL.pdf: \ 18 | USERMANUAL.tex \ 19 | HEADER.tex \ 20 | introduction.tex \ 21 | quickstart.tex \ 22 | hardware_software.tex \ 23 | installation.tex \ 24 | runs.tex \ 25 | param.tex \ 26 | ../../PARAM.xmltex \ 27 | output.tex \ 28 | visualization.tex \ 29 | synopticSW.tex 30 | 31 | 32 | PDF: USERMANUAL.pdf DESIGN.pdf STANDARDS.pdf 33 | 34 | clean: 35 | rm -rf *.aux *.log *.toc *.dvi *.ind *.idx *~ ../../PARAM.xmltex 36 | 37 | cleanpdf: clean 38 | cd ..; rm -rf USERMANUAL.pdf DESIGN.pdf STANDARDS.pdf 39 | 40 | 41 | -------------------------------------------------------------------------------- /Doc/Tex/STANDARDS.tex: -------------------------------------------------------------------------------- 1 | % Copyright (C) 2002 Regents of the University of Michigan, portions used with permission 2 | % For more information, see http://csem.engin.umich.edu/tools/swmf 3 | \documentclass[twoside,10pt]{book} 4 | 5 | \input HEADER 6 | 7 | \title{STANDARDS} 8 | \author{David Chesney and G\'abor T\'oth} 9 | 10 | \begin{document} 11 | 12 | \pagestyle{headings} 13 | 14 | \maketitle 15 | 16 | \tableofcontents 17 | 18 | % Data naming standard 19 | \input{datanames.tex} 20 | 21 | \end{document} 22 | -------------------------------------------------------------------------------- /Doc/Tex/USERMANUAL.tex: -------------------------------------------------------------------------------- 1 | % Copyright (C) 2002 Regents of the University of Michigan, portions used with permission 2 | % For more information, see http://csem.engin.umich.edu/tools/swmf 3 | \documentclass[twoside,10pt]{book} 4 | 5 | \usepackage{sverb} 6 | 7 | \input HEADER 8 | 9 | \title{\BATSRUS\ and CRASH User Manual} 10 | \author{CSEM and CRASH team} 11 | 12 | \begin{document} 13 | 14 | \pagestyle{headings} 15 | 16 | \maketitle 17 | 18 | \tableofcontents 19 | 20 | % Files for introduction 21 | \input{introduction} 22 | 23 | % Files for chapter in installation, compilation and running the code. 24 | \chapter{Using \BATSRUS\ and CRASH \label{chapter:usage}} 25 | \input{quickstart} 26 | \input{hardware_software} 27 | \input{installation} 28 | \input{runs} 29 | 30 | % Files for chapter on input parameteres 31 | \chapter{Input Parameters \label{chapter:param}} 32 | \input{param.tex} 33 | 34 | \input ../../PARAM.xmltex 35 | 36 | % Files for chapter on Code output 37 | \input{output.tex} 38 | 39 | % Files for chapter on visualization 40 | \chapter{Visualization \label{chapter:visualization}} 41 | \input{visualization.tex} 42 | 43 | % Files for chapter on the synoptic solar wind model 44 | \input{synopticSW.tex} 45 | 46 | %\input{bibliography.tex} 47 | 48 | \end{document} 49 | -------------------------------------------------------------------------------- /Doc/Tex/adapted_grid.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Doc/Tex/adapted_grid.pdf -------------------------------------------------------------------------------- /Doc/Tex/adaptiveblock.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Doc/Tex/adaptiveblock.pdf -------------------------------------------------------------------------------- /Doc/Tex/fivol.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Doc/Tex/fivol.pdf -------------------------------------------------------------------------------- /Doc/Tex/ionosphere.tex: -------------------------------------------------------------------------------- 1 | % Copyright (C) 2002 Regents of the University of Michigan, portions used with permission 2 | % For more information, see http://csem.engin.umich.edu/tools/swmf 3 | \section{The Ionosphere} 4 | 5 | The \BATSRUS\ code is capable of using an inner boundary condition 6 | which is similar to the ionosphere. The field aligned currents near 7 | the inner boundary are used to drive an ionospheric flow, which is 8 | used as a lower boundary in the magnetosphere. 9 | 10 | \begin{itemize} 11 | 12 | \item 13 | Currents are calculated from $\mathbf{\nabla \times B}$ at some 14 | distance away from the inner boundary ({\tt rCurrents}). 15 | Specifically, the code runs through all of the cells in the simulation 16 | to determine whether the particular cell is the first cell (as 17 | compared to it's neighbors) which is outside of {\tt rCurrents}. 18 | This is in {\tt magnetosphere.f90}. 19 | 20 | \item 21 | The code determines whether the cell maps to the northern or southern 22 | hemisphere by testing whether it maps to the northern hemisphere. It 23 | stores the value of the full 3-D current as well as the location in 24 | arrays. This is in {\tt magnetosphere.f90}. 25 | 26 | \item 27 | If {\tt UseFullCurrents} is {\tt .false.}, the code determines the 28 | magnetic field direction in the magnetosphere and takes a dot product 29 | of the 3-D current and the field direction. If not, the 30 | 31 | \end{itemize} 32 | -------------------------------------------------------------------------------- /Doc/Tex/message_pass_equal.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Doc/Tex/message_pass_equal.pdf -------------------------------------------------------------------------------- /Doc/Tex/message_pass_res_change.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Doc/Tex/message_pass_res_change.pdf -------------------------------------------------------------------------------- /Doc/Tex/proc_dims.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Doc/Tex/proc_dims.pdf -------------------------------------------------------------------------------- /Doc/Tex/slope.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Doc/Tex/slope.pdf -------------------------------------------------------------------------------- /Doc/index.html: -------------------------------------------------------------------------------- 1 | 2 | 3 | 5 | 6 | BATS-R-US MANUALS 7 | 8 | 9 |
10 |

BATS-R-US MANUALS

11 |
12 | 13 | 14 | 27 | 28 | 29 | 30 | -------------------------------------------------------------------------------- /LICENSE.txt: -------------------------------------------------------------------------------- 1 | Copyright (C) 2002 Regents of the University of Michigan, 2 | portions used with permission. 3 | 4 | Licensed under the Apache License, Version 2.0 (the "License"); 5 | you may not use this file except in compliance with the License. 6 | You may obtain a copy of the License at 7 | 8 | http://www.apache.org/licenses/LICENSE-2.0 9 | 10 | Unless required by applicable law or agreed to in writing, software 11 | distributed under the License is distributed on an "AS IS" BASIS, 12 | WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 13 | See the License for the specific language governing permissions and 14 | limitations under the License. 15 | -------------------------------------------------------------------------------- /Param/2BODYPLOT/TestOutput/REFy=0.dat.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/2BODYPLOT/TestOutput/REFy=0.dat.gz -------------------------------------------------------------------------------- /Param/2BODYPLOT/TestOutput/REFy=0.out.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/2BODYPLOT/TestOutput/REFy=0.out.gz -------------------------------------------------------------------------------- /Param/AMR/TestOutput/log_sph.log: -------------------------------------------------------------------------------- 1 | Volume averages, fluxes, etc 2 | it t dt rho mx my mz p bx by bz pmin pmax status 3 | 0 0.00000E+00 0.00000E+00 5.00000E+00 -3.34520E-15 0.00000E+00 0.00000E+00 1.72587E-03 0.00000E+00 0.00000E+00 -5.00000E+00 1.72587E-03 1.72587E-03 0.00000E+00 4 | 1 0.00000E+00 0.00000E+00 5.00001E+00 -3.34495E-15 -2.11430E-34 -2.80566E-34 1.72610E-03 8.14538E-18 -1.28811E-17 -5.00206E+00 1.63830E-03 1.81711E-03 0.00000E+00 5 | 2 0.00000E+00 0.00000E+00 5.00003E+00 -3.34468E-15 5.67773E-34 -7.26832E-34 1.72626E-03 1.25022E-17 -5.30397E-18 -5.00410E+00 1.56514E-03 1.90176E-03 0.00000E+00 6 | 3 0.00000E+00 0.00000E+00 5.00006E+00 -3.34441E-15 6.77699E-34 4.30635E-34 1.72640E-03 -6.51631E-17 -1.66696E-17 -5.00612E+00 1.50408E-03 1.97815E-03 0.00000E+00 7 | 4 0.00000E+00 0.00000E+00 5.00008E+00 -3.34412E-15 1.61499E-33 5.34730E-34 1.72657E-03 -3.78855E-17 -2.27313E-17 -5.00812E+00 1.45339E-03 2.04537E-03 0.00000E+00 8 | 5 0.00000E+00 0.00000E+00 5.00012E+00 -3.34382E-15 1.21761E-33 -8.66453E-34 1.72677E-03 1.04564E-16 -1.33357E-16 -5.01012E+00 1.41189E-03 2.10223E-03 0.00000E+00 9 | 6 0.00000E+00 0.00000E+00 5.00017E+00 -3.34351E-15 1.05046E-33 -1.32084E-33 1.72701E-03 -7.57710E-18 -1.51542E-17 -5.01211E+00 1.37768E-03 2.14659E-03 0.00000E+00 10 | 7 0.00000E+00 0.00000E+00 5.00023E+00 -3.34319E-15 1.15110E-33 -2.46719E-33 1.72729E-03 -6.97093E-17 -1.75789E-16 -5.01410E+00 1.34976E-03 2.18265E-03 0.00000E+00 11 | 8 0.00000E+00 0.00000E+00 5.00030E+00 -3.34287E-15 1.91173E-33 -2.41190E-33 1.72762E-03 -1.12141E-16 -7.27402E-17 -5.01609E+00 1.32731E-03 2.21048E-03 0.00000E+00 12 | 9 0.00000E+00 0.00000E+00 5.00038E+00 -3.34254E-15 3.78112E-33 2.24604E-35 1.72799E-03 1.69727E-16 1.09110E-16 -5.01811E+00 1.30968E-03 2.23954E-03 0.00000E+00 13 | 10 0.00000E+00 0.00000E+00 5.00047E+00 -3.34220E-15 3.77680E-33 7.60198E-34 1.72839E-03 -9.09252E-18 -9.09252E-18 -5.02014E+00 1.29639E-03 2.26473E-03 0.00000E+00 14 | -------------------------------------------------------------------------------- /Param/ANISOPRESSURE/PARAM.in.Alfven: -------------------------------------------------------------------------------- 1 | #PLANET 2 | NONE NamePlanet 3 | 4 | #IOUNITS 5 | NONE NameIoUnits 6 | 7 | #NORMALIZATION 8 | NONE NameNormalization 9 | 10 | #NONCONSERVATIVE 11 | T UseNonConservative 12 | 13 | #CHECKGRIDSIZE 14 | 100 nI 15 | 2 nJ 16 | 2 nK 17 | 100 MinBlockAll 18 | 19 | #GRID 20 | 4 nRootBlockX 21 | 1 nRootBlockY 22 | 1 nRootBlockZ 23 | -6.0 xMin 24 | 6.0 xMax 25 | -0.05 yMin 26 | 0.05 yMax 27 | -0.05 zMin 28 | 0.05 zMax 29 | 30 | #UNIFORMSTATE 31 | 1.0 StateVar Rho 32 | 0.0 StateVar Ux 33 | 0.0 StateVar Uy 34 | 0.0 StateVar Uz 35 | 10.0 StateVar Bx 36 | 0.0 StateVar By 37 | 0.0 StateVar Bz 38 | 6.0 StateVar Ppar 39 | 35.3333 StateVar p 40 | 41 | #WAVE 42 | Uy NameVar 43 | 1.0 Width 44 | 0.12 Amplitude 45 | 6.0 LambdaX 46 | -1.0 LambdaY 47 | -1.0 LambdaZ 48 | 0.0 Phase 49 | 50 | #WAVE 51 | By NameVar 52 | 1.0 Width 53 | 0.1 Amplitude 54 | 6.0 LambdaX 55 | -1.0 LambdaY 56 | -1.0 LambdaZ 57 | 0.0 Phase 58 | 59 | #OUTERBOUNDARY 60 | periodic xMinBc 61 | periodic xMaxBc 62 | float yMinBc 63 | float yMaxBc 64 | float zMinBc 65 | float zMaxBc 66 | 67 | #TIMESTEPPING 68 | 2 nStage 69 | 0.8 CflExpl 70 | 71 | #SCHEME 72 | 2 nOrder 73 | Rusanov TypeFlux 74 | mc3 TypeLimiter 75 | 1.5 LimiterBeta 76 | 77 | #SAVELOGFILE 78 | T DoSaveLogfile 79 | raw StringLogfile 80 | 1 DnSaveLogfile 81 | -1. DtSaveLogfile 82 | 83 | #SAVERESTART 84 | F DoSaverestart 85 | 86 | #PLOTFILENAME 87 | second NameMaxTimeUnit 88 | 89 | #SAVEINITIAL 90 | T DoSaveInitial 91 | 92 | #SAVEPLOT 93 | 1 nPlotfile 94 | 1d mhd idl_ascii PlotString 95 | -1 DnSavePlot 96 | 0.1 DtSavePlot 97 | 0.0 DxSavePlot 98 | 99 | #STOP 100 | -1 MaxIteration 101 | 3.0 tSimulationMax 102 | -------------------------------------------------------------------------------- /Param/ANISOPRESSURE/PARAM.in.fastwave: -------------------------------------------------------------------------------- 1 | #DESCRIPTION 2 | Param/ANISOPRESSURE/PARAM.in.fastwave: P anisotropic StringDescription 3 | 4 | #PLANET 5 | NONE NamePlanet 6 | 7 | #IOUNITS 8 | NONE TypeIoUnit 9 | 10 | #NORMALIZATION 11 | NONE TypeNormalization 12 | 13 | #ANISOTROPICPRESSURE 14 | T UseConstantTau 15 | -1.0 TauInstabilitySi 16 | -1.0 TauGlobalSi 17 | 18 | #NONCONSERVATIVE 19 | T UseNonConservative 20 | 21 | #ENTROPY 22 | T UseEntropy 23 | 24 | #GRIDBLOCKALL 25 | 20 MaxBlock 26 | 27 | #GRID 28 | 10 nRootBlock1 29 | 1 nRootBlock2 30 | 1 nRootBlock3 31 | -30.0 xMin 32 | 30.0 xMax 33 | -0.05 yMin 34 | 0.05 yMax 35 | -0.05 zMin 36 | 0.05 zMax 37 | 38 | #UNIFORMSTATE 39 | 1.0 StateVar Rho (left state) 40 | 12.0 StateVar Ux 41 | 0.0 StateVar Uy 42 | 0.0 StateVar Uz 43 | 0.0 StateVar Bx 44 | 30.0 StateVar By 45 | 0.0 StateVar Bz 46 | 14.0 StateVar Ppar 47 | 16.0 StateVar P 48 | 49 | #WAVE 50 | Rho NameVar 51 | 3.0 Width 52 | 0.01 Amplitude 53 | 12.0 LambdaX 54 | -1.0 LambdaY 55 | -1.0 LambdaZ 56 | 0.0 Phase 57 | 58 | #WAVE 59 | P NameVar 60 | 3.0 Width 61 | 0.16 Amplitude 62 | 12.0 LambdaX 63 | -1.0 LambdaY 64 | -1.0 LambdaZ 65 | 0.0 Phase 66 | 67 | #WAVE 68 | Ppar NameVar 69 | 3.0 Width 70 | 0.14 Amplitude 71 | 12.0 LambdaX 72 | -1.0 LambdaY 73 | -1.0 LambdaZ 74 | 0.0 Phase 75 | 76 | #WAVE 77 | By NameVar 78 | 3.0 Width 79 | 0.3 Amplitude 80 | 12.0 LambdaX 81 | -1.0 LambdaY 82 | -1.0 LambdaZ 83 | 0.0 Phase 84 | 85 | #OUTERBOUNDARY 86 | periodic TypeBc1 87 | periodic TypeBc2 88 | float TypeBc3 89 | float TypeBc4 90 | float TypeBc5 91 | float TypeBc6 92 | 93 | #TIMESTEPPING 94 | 2 nStage 95 | 0.8 CflExpl 96 | 97 | #SCHEME 98 | 2 nOrder 99 | Rusanov TypeFlux 100 | mc3 TypeLimiter 101 | 1.5 LimiterBeta 102 | 103 | #SAVELOGFILE 104 | T DoSaveLogfile 105 | raw StringLog 106 | 1 DnSaveLogfile 107 | -1. DtSaveLogfile 108 | 109 | #SAVERESTART 110 | F DoSaveRestart 111 | 112 | #PLOTFILENAME 113 | second NameMaxTimeUnit 114 | 115 | #SAVEINITIAL 116 | T DoSaveInitial 117 | 118 | #SAVEPLOT 119 | 1 nPlotFile 120 | 1d var idl_ascii StringPlot 121 | -1 DnSavePlot 122 | 0.1 DtSavePlot 123 | 0.0 DxSavePlot 124 | {mhd} NameVars 125 | {default} NamePars 126 | 127 | #STOP 128 | -1 MaxIteration 129 | 1.0 tSimulationMax 130 | -------------------------------------------------------------------------------- /Param/ANISOPRESSURE/PARAM.in.soundwave: -------------------------------------------------------------------------------- 1 | #PLANET 2 | NONE NamePlanet 3 | 4 | #IOUNITS 5 | NONE NameIoUnits 6 | 7 | #NORMALIZATION 8 | NONE NameNormalization 9 | 10 | #NONCONSERVATIVE 11 | T UseNonConservative 12 | 13 | #CHECKGRIDSIZE 14 | 100 nI 15 | 2 nJ 16 | 2 nK 17 | 100 MinBlockAll 18 | 19 | #GRID 20 | 12 nRootBlockX 21 | 1 nRootBlockY 22 | 1 nRootBlockZ 23 | -6.0 xMin 24 | 6.0 xMax 25 | -0.05 yMin 26 | 0.05 yMax 27 | -0.05 zMin 28 | 0.05 zMax 29 | 30 | #SHOCKTUBE 31 | 1.0 LeftState rho 32 | 0.0 LeftState Ux 33 | 0.0 LeftState Uy 34 | 0.0 LeftState Uz 35 | 10.0 LeftState Bx 36 | 0.0 LeftState By 37 | 0.0 LeftState Bz 38 | 3.0 LeftState Ppar 39 | 3.0 LeftState p 40 | 1.0 RightState rho 41 | 0.0 RightState Ux 42 | 0.0 RightState Uy 43 | 0.0 RightState Uz 44 | 10.0 RightState Bx 45 | 0.0 RightState By 46 | 0.0 RightState Bz 47 | 3.0 RightState Ppar 48 | 3.0 RightState p 49 | 50 | #WAVE 51 | Rho NameVar 52 | 12.0 Width 53 | 0.02 Amplitude 54 | 6.0 LambdaX 55 | -1.0 LambdaY 56 | -1.0 LambdaZ 57 | 0.0 Phase 58 | 59 | #WAVE 60 | Ppar NameVar 61 | 12.0 Width 62 | 0.18 Amplitude 63 | 6.0 LambdaX 64 | -1.0 LambdaY 65 | -1.0 LambdaZ 66 | 0.0 Phase 67 | 68 | #WAVE 69 | p NameVar 70 | 12.0 Width 71 | 0.1 Amplitude 72 | 6.0 LambdaX 73 | -1.0 LambdaY 74 | -1.0 LambdaZ 75 | 0.0 Phase 76 | 77 | #OUTERBOUNDARY 78 | periodic xMinBc 79 | periodic xMaxBc 80 | float yMinBc 81 | float yMaxBc 82 | float zMinBc 83 | float zMaxBc 84 | 85 | #TIMESTEPPING 86 | 2 nStage 87 | 0.8 CflExpl 88 | 89 | #SCHEME 90 | 2 nOrder 91 | Rusanov TypeFlux 92 | mc3 TypeLimiter 93 | 1.5 LimiterBeta 94 | 95 | #SAVELOGFILE 96 | T DoSaveLogfile 97 | raw StringLogfile 98 | 1 DnSaveLogfile 99 | -1. DtSaveLogfile 100 | 101 | #SAVERESTART 102 | F DoSaverestart 103 | 104 | #PLOTFILENAME 105 | second NameMaxTimeUnit 106 | 107 | #SAVEINITIAL 108 | T DoSaveInitial 109 | 110 | #SAVEPLOT 111 | 1 nPlotfile 112 | cut mhd idl PlotString 113 | -1 DnSavePlot 114 | 0.1 DtSavePlot 115 | -6.0 Coord1MinCut 116 | +6.0 Coord1MaxCut 117 | -0.001 Coord2MinCut 118 | +0.001 Coord2MaxCut 119 | -0.001 Coord3MinCut 120 | +0.001 Coord3MaxCut 121 | 0.0 DxSavePlot 122 | 123 | #STOP 124 | -1 MaxIteration 125 | 6.0 tSimulationMax 126 | -------------------------------------------------------------------------------- /Param/B0/PARAM.in: -------------------------------------------------------------------------------- 1 | #DESCRIPTION 2 | Param/B0/PARAM.in: check equilibrium 3 | 4 | ! Switch off update B0 5 | #IDEALAXES 6 | 7 | ! Should work with conservative energy and non-conservative pressure equation (trivial) 8 | #NONCONSERVATIVE 9 | F UseNonConservative 10 | 11 | ! Large errors if B0 source terms are switched off 12 | B0SOURCE 13 | F UseB0Source 14 | 15 | #SAVERESTART 16 | F DoSaveRestart 17 | 18 | ! Avoid singularity at the origin by starting at zMin=3 19 | #GRID 20 | 2 nRootBlockX 21 | 2 nRootBlockY 22 | 1 nRootBlockZ 23 | -16. xMin 24 | +16. xMax 25 | -16.0 yMin 26 | +16.0 yMax 27 | +3.0 zMin 28 | +19.0 zMax 29 | 30 | #GRIDLEVEL 31 | 2 nLevel 32 | initial NameArea 33 | 34 | #GRIDLEVEL 35 | 3 Resolution 36 | sphere0 NameArea 37 | 6.0 Radius 38 | 39 | #DOAMR 40 | T DoAmr 41 | 6 DnAmr 42 | -1.0 DtAmr 43 | T UseStrictAmr 44 | 45 | #BODY 46 | F UseBody 47 | 48 | #SOLARWIND 49 | 5.0 SwRhoDim [amu/cc] 50 | 25000.0 SwTDim [K] 51 | 0.0 SwUxDim [km/s] 52 | 0.0 SwUyDim [km/s] 53 | 0.0 SwUzDim [km/s] 54 | 30.0 SwBxDim [nT] 55 | 40.0 SwByDim [nT] 56 | 60.0 SwBzDim [nT] 57 | 58 | #OUTERBOUNDARY 59 | fixedb1 TypeBc1 60 | fixedb1 TypeBc2 61 | fixedb1 TypeBc3 62 | fixedb1 TypeBc4 63 | fixedb1 TypeBc5 64 | fixedb1 TypeBc6 65 | 66 | #TIMESTEPPING 67 | 2 nStage 68 | 0.8 CflExlp 69 | 70 | #SCHEME 71 | 2 nOrder 72 | Rusanov TypeFlux 73 | mc3 TypeLimiter 74 | 1.5 BetaLimiter 75 | 76 | #SAVELOGFILE 77 | T DoSaveLogfile 78 | RAW StringLogfile 79 | 1 DnSaveLogfile 80 | -1. DtSaveLogfile 81 | 82 | #SAVEPLOT 83 | 2 nPlotFiles 84 | x=0 FUL idl StringPlot 85 | 1 DnSavePlot 86 | -1. DtSavePlot 87 | -1.0 Dx 88 | y=0 FUL idl StringPlot 89 | 1 DnSavePlot 90 | -1. DtSavePlot 91 | -1.0 Dx 92 | 93 | #STOP 94 | 10 MaxIteration 95 | -1.0 tSimulationMax 96 | -------------------------------------------------------------------------------- /Param/B0/PARAM.in.sph: -------------------------------------------------------------------------------- 1 | #DESCRIPTION 2 | Param/B0/PARAM.in.sph: check equilibrium 3 | 4 | ! Switch off update B0 5 | #IDEALAXES 6 | 7 | TESTXYZ 8 | 1.0 xTest 9 | 0.01 yTest 10 | 3.5 zTest 11 | 12 | TESTVAR 13 | 2 Ux 14 | 15 | TEST 16 | update_states update_states_mhd calc_sources set_b0_reschange 17 | 18 | ! Should work with conservative energy and non-conservative pressure equation (trivial) 19 | NONCONSERVATIVE 20 | F UseNonConservative 21 | 22 | ! Large errors if B0 source terms are switched off 23 | B0SOURCE 24 | F UseB0Source 25 | 26 | #SAVERESTART 27 | F DoSaveRestart 28 | 29 | #GRID 30 | 2 nRootBlockX 31 | 2 nRootBlockY 32 | 1 nRootBlockZ 33 | -16. xMin 34 | +16. xMax 35 | -16.0 yMin 36 | +16.0 yMax 37 | -16.0 zMin 38 | +16.0 zMax 39 | 40 | #GRIDGEOMETRY 41 | spherical TypeGeometry 42 | 43 | #LIMITRADIUS 44 | 3.0 RadiusMin 45 | 16.0 RadiusMax 46 | 47 | #GRIDLEVEL 48 | 1 nLevel 49 | initial NameArea 50 | 51 | #GRIDLEVEL 52 | 3 Resolution 53 | sphere0 NameArea 54 | 6.0 Radius 55 | 56 | #DOAMR 57 | T DoAmr 58 | 6 DnAmr 59 | -1.0 DtAmr 60 | T UseStrictAmr 61 | 62 | #BODY 63 | F UseBody 64 | 65 | #SOLARWIND 66 | 5.0 SwRhoDim [amu/cc] 67 | 25000.0 SwTDim [K] 68 | 0.0 SwUxDim [km/s] 69 | 0.0 SwUyDim [km/s] 70 | 0.0 SwUzDim [km/s] 71 | 30.0 SwBxDim [nT] 72 | 40.0 SwByDim [nT] 73 | 60.0 SwBzDim [nT] 74 | 75 | #OUTERBOUNDARY 76 | float TypeBc1 77 | float TypeBc2 78 | float TypeBc3 79 | float TypeBc4 80 | float TypeBc5 81 | float TypeBc6 82 | 83 | #TIMESTEPPING 84 | 2 nStage 85 | 0.8 CflExlp 86 | 87 | #SCHEME 88 | 2 nOrder 89 | Rusanov TypeFlux 90 | mc3 TypeLimiter 91 | 1.5 BetaLimiter 92 | 93 | #SAVELOGFILE 94 | T DoSaveLogfile 95 | RAW StringLogfile 96 | 1 DnSaveLogfile 97 | -1. DtSaveLogfile 98 | 99 | #SAVEPLOT 100 | 3 nPlotFiles 101 | x=0 FUL idl StringPlot 102 | 1 DnSavePlot 103 | -1. DtSavePlot 104 | -1. Dx 105 | y=0 FUL idl StringPlot 106 | 1 DnSavePlot 107 | -1. DtSavePlot 108 | -1. Dx 109 | z=0 FUL idl StringPlot 110 | 1 DnSavePlot 111 | -1. DtSavePlot 112 | -1. Dx 113 | 114 | #STOP 115 | 10 MaxIteration 116 | -1.0 tSimulationMax 117 | -------------------------------------------------------------------------------- /Param/COMET3FLUIDSPE/TestOutput/log_n000060.log: -------------------------------------------------------------------------------- 1 | Volume averages, fluxes, etc 2 | it year mo dy hr mn sc msc t Rho Pmin Pmax Pe Bx By Bz SwRho SwUx SwUy SwUz SwP H2OpRho H2OpUx H2OpUy H2OpUz H2OpP 3 | 60 2000 01 01 00 00 00 000 0.00000E+00 1.18816E+003 5.17634E+000 1.35993E+002 2.31630E-002 -1.60533E+000 1.90052E+000 -2.23346E-003 6.19891E+001 9.04586E+000 -4.87113E+000 1.23335E+001 1.39537E-002 1.12364E+003 -2.85320E+000 1.10644E+000 4.44492E+000 3.20676E+001 4 | -------------------------------------------------------------------------------- /Param/CORONA/AIA_XRT: -------------------------------------------------------------------------------- 1 | ! Full euv and X-ray images 2 | ! SDO(AIA) - 6 bands of EUV 3 | ! Hinode(XRT) - XRT (three channels) 4 | !------------------------- 5 | ! Position should be in HGI 6 | ! Coordinates, Time as set in PARAM.in 7 | ! file should coinside with the 8 | ! observation time 9 | 10 | #LOOKUPTABLE 11 | AiaXrt NameTable 12 | load TypeAction 13 | Param/los_tbl.dat NameFile 14 | ascii TypeFIle 15 | 16 | 17 | 18 | 19 | #SAVEPLOT 20 | 1 nPlotFiles 21 | los tbl tec StringPlot 22 | 40000 DnSavePlot 23 | -1 DtSavePlot 24 | -4.21 xPos(Earth) 25 | 211.56 yPos(Earth) 26 | -26.91 zPos(Earth) 27 | 0. OffsetAngle 28 | 1.98 rSizeImage 29 | 0. xOffset 30 | 0. yOffset 31 | 0. rOccult 32 | 0. MuLimbDarkening 33 | 512 nPix 34 | AiaXrt NameLosTable 35 | 36 | #END -------------------------------------------------------------------------------- /Param/CORONA/B0local.out: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/CORONA/B0local.out -------------------------------------------------------------------------------- /Param/CORONA/EIT_SXR_TEC: -------------------------------------------------------------------------------- 1 | #SAVEPLOT 2 | 2 nPlotFiles 3 | los euv tec StringPlot 4 | 40000 DnSavePlot 5 | -1 DtSavePlot 6 | #INCLUDE 7 | POS_EARTH 8 | 0. OffsetAngle 9 | 1.98 rSizeImage 10 | 0. xOffset 11 | 0. yOffset 12 | 0. rOccult 13 | 0. MuLimbDarkening 14 | 512 nPix 15 | los sxr tec StringPlot 16 | 40000 DnSavePlot 17 | -1 DtSavePlot 18 | #INCLUDE 19 | POS_EARTH 20 | 0. OffsetAngle 21 | 2.11 rSizeImage 22 | 0. xOffset 23 | 0. yOffset 24 | 0. rOccult 25 | 0. MuLimbDarkening 26 | 512 nPix 27 | -------------------------------------------------------------------------------- /Param/CORONA/MDI_small.dat: -------------------------------------------------------------------------------- 1 | MDI_small.dat 2 | 0 0.0 -2 3 2 3 | 3 1 4 | 1 2044 0.0 5 | n m g h nOrder CR dLon 6 | 0 0 0.00000 0.00000 7 | 1 0 20.00000 0.00000 8 | 1 1 10.00000 0.00000 9 | -------------------------------------------------------------------------------- /Param/CORONA/NEWMAGNETOGRAM.in: -------------------------------------------------------------------------------- 1 | 2 | #NEWMAGNETOGRAM 3 | T UseMagnetogram 4 | 1.0 rMagnetogram 5 | 2.5 rSourceSurface 6 | 0.0 HeightInnerBc 7 | SC/newmf.dat NameMagnetogramFile 8 | 12 nHeaderLine 9 | -1.0 PhiShift 10 | 1.0 UnitB 11 | 12 | #HELIOUPDATEB0 13 | 3600.0 14 | 15 | #INCLUDE 16 | FINALTIME.in 17 | 18 | #END -------------------------------------------------------------------------------- /Param/CORONA/RadCoolCorona_6.0.dat: -------------------------------------------------------------------------------- 1 | Radiative losses of optically thin Coronal plasma, sun_coronal_ext.abund (CHIANTI 5.0): [log10 K] [ergs*cm^3*1e22] [#] 2 | 0 0.0 1 1 1 3 | 41 4 | 5.0 5 | logTe lambdaT ChiantiVersion 6 | 0.4000000E+01 0.9775910E-01 7 | 0.4100002E+01 0.6204106E+00 8 | 0.4199999E+01 0.1821040E-01 9 | 0.4300000E+01 0.1402317E-01 10 | 0.4400001E+01 0.8884599E-01 11 | 0.4500000E+01 0.8374518E-01 12 | 0.4600000E+01 0.1026979E-01 13 | 0.4700000E+01 0.1498370E-01 14 | 0.4800001E+01 0.2397494E-01 15 | 0.4900000E+01 0.3263567E-01 16 | 0.5000000E+01 0.3730172E-01 17 | 0.5100002E+01 0.3606443E-01 18 | 0.5199999E+01 0.3942354E-01 19 | 0.5300000E+01 0.4287060E-01 20 | 0.5400001E+01 0.4254672E-01 21 | 0.5500000E+01 0.2678095E-01 22 | 0.5600000E+01 0.1801587E-01 23 | 0.5700000E+01 0.1726192E-01 24 | 0.5800001E+01 0.1659424E-01 25 | 0.5900000E+01 0.1713749E-01 26 | 0.6000000E+01 0.1876034E-01 27 | 0.6100002E+01 0.1862474E-01 28 | 0.6199999E+01 0.1710604E-01 29 | 0.6300000E+01 0.1240186E-01 30 | 0.6400001E+01 0.7562685E+00 31 | 0.6500000E+01 0.5045575E+00 32 | 0.6600000E+01 0.4130752E+00 33 | 0.6700000E+01 0.4010693E+00 34 | 0.6800001E+01 0.4146723E+00 35 | 0.6900000E+01 0.4200034E+00 36 | 0.7000000E+01 0.4041402E+00 37 | 0.7100002E+01 0.3526888E+00 38 | 0.7199999E+01 0.2736992E+00 39 | 0.7300000E+01 0.2232313E+00 40 | 0.7400001E+01 0.2027298E+00 41 | 0.7500000E+01 0.1990994E+00 42 | 0.7600000E+01 0.2050792E+00 43 | 0.7700000E+01 0.2172907E+00 44 | 0.7800001E+01 0.2331289E+00 45 | 0.7900000E+01 0.2522123E+00 46 | 0.8000000E+01 0.2731431E+00 47 | -------------------------------------------------------------------------------- /Param/CORONA/TR8.dat: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/CORONA/TR8.dat -------------------------------------------------------------------------------- /Param/CORONA/TestOutput/3d__fux_5_n00000004.dat.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/CORONA/TestOutput/3d__fux_5_n00000004.dat.gz -------------------------------------------------------------------------------- /Param/CORONA/TestOutput/3d__phx_7_n00000004.dat.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/CORONA/TestOutput/3d__phx_7_n00000004.dat.gz -------------------------------------------------------------------------------- /Param/CORONA/TestOutput/log.awsom_large_gpu: -------------------------------------------------------------------------------- 1 | Volume averages, fluxes, etc 2 | it t dt Rho Mx My Mz Bx By Bz Ehot Ew Pe E Pmin Pmax 3 | 167 2.50000E+01 1.69765E-001 2.21392E-019 2.71825E-016 -6.10242E-016 -1.43698E-016 -2.31993E-006 3.90559E-007 6.04929E-007 1.46860E-006 3.53025E-005 4.03877E-006 2.00265E-005 1.24899E-006 1.38700E+000 4 | -------------------------------------------------------------------------------- /Param/CORONA/TestOutput/log_awsomchargestate.log: -------------------------------------------------------------------------------- 1 | Volume averages, fluxes, etc 2 | it t dt Rho Mx My Mz Bx By Bz Ehot o1 o2 o3 o4 o5 o6 o7 o8 o9 Ew Pe Ppar E Pmin Pmax 3 | 1 0.00000E+00 0.00000E+000 1.91842E-019 -1.76875E-028 0.00000E+000 -4.50338E-028 4.57159E-022 2.17505E-021 -2.44792E-006 1.48343E-006 5.03826E-024 1.33710E-021 5.89923E-022 8.13395E-025 6.79575E-027 1.12613E-024 3.54996E-022 7.42011E-023 1.67104E-024 4.09990E-005 4.04507E-006 5.40523E-006 1.70390E-005 8.26826E-007 9.67403E-001 4 | 2 0.00000E+00 0.00000E+000 2.00572E-019 -2.65122E-028 -1.70960E-028 5.46153E-026 -8.16388E-022 -2.34067E-021 -4.96672E-006 1.52111E-006 5.22761E-024 1.38784E-021 6.13888E-022 1.79115E-024 3.51226E-025 1.15102E-024 3.54904E-022 7.41773E-023 1.67044E-024 3.96459E-005 4.14084E-006 5.85940E-006 1.78822E-005 9.39702E-007 1.05612E+000 5 | 3 0.00000E+00 0.00000E+000 2.13337E-019 -4.12136E-028 -2.68652E-028 5.79372E-025 -2.79245E-021 -3.62945E-021 -7.54613E-006 1.55820E-006 5.55059E-024 1.47387E-021 6.52498E-022 2.77168E-024 8.25405E-025 1.23104E-024 3.54819E-022 7.41554E-023 1.66990E-024 3.83524E-005 4.25266E-006 6.28869E-006 1.87112E-005 1.03616E-006 1.15049E+000 6 | 4 0.00000E+00 0.00000E+000 2.28931E-019 -5.36350E-028 -2.19806E-028 9.15232E-025 3.21319E-021 -3.65529E-021 -1.01902E-005 1.59128E-006 5.97698E-024 1.58714E-021 7.03078E-022 3.61459E-024 1.39144E-024 1.35314E-024 3.54741E-022 7.41351E-023 1.66941E-024 3.71467E-005 4.37115E-006 6.67824E-006 1.94965E-005 1.11853E-006 1.24353E+000 7 | -------------------------------------------------------------------------------- /Param/CORONA/TestOutput/log_n000000.bvector: -------------------------------------------------------------------------------- 1 | Volume averages, fluxes, etc 2 | it t dt absjxb emag1 ekin 3 | 0 0.00000E+00 0.00000E+000 6.35945E+002 0.00000E+000 0.00000E+000 4 | 1 0.00000E+00 0.00000E+000 8.14357E+002 2.26328E-001 8.16645E-003 5 | 2 0.00000E+00 0.00000E+000 9.97740E+002 7.37194E-001 7.42416E-003 6 | 3 0.00000E+00 0.00000E+000 1.14625E+003 1.34299E+000 7.22414E-003 7 | 4 0.00000E+00 0.00000E+000 1.26566E+003 1.95367E+000 7.27579E-003 8 | 5 0.00000E+00 0.00000E+000 1.36094E+003 2.52574E+000 7.38155E-003 9 | 6 0.00000E+00 0.00000E+000 1.43693E+003 3.04629E+000 7.47062E-003 10 | 7 0.00000E+00 0.00000E+000 1.49806E+003 3.52028E+000 7.51549E-003 11 | 8 0.00000E+00 0.00000E+000 1.54771E+003 3.95132E+000 7.52128E-003 12 | 9 0.00000E+00 0.00000E+000 1.58970E+003 4.35261E+000 7.49166E-003 13 | 10 0.00000E+00 0.00000E+000 1.62579E+003 4.72922E+000 7.43765E-003 14 | 11 0.00000E+00 0.00000E+000 1.65738E+003 5.08568E+000 7.37021E-003 15 | 12 0.00000E+00 0.00000E+000 1.68538E+003 5.42491E+000 7.29557E-003 16 | 13 0.00000E+00 0.00000E+000 1.71033E+003 5.74871E+000 7.21838E-003 17 | 14 0.00000E+00 0.00000E+000 1.73286E+003 6.05735E+000 7.13919E-003 18 | 15 0.00000E+00 0.00000E+000 1.75349E+003 6.35365E+000 7.06273E-003 19 | 16 0.00000E+00 0.00000E+000 1.77252E+003 6.63933E+000 6.98847E-003 20 | 17 0.00000E+00 0.00000E+000 1.79001E+003 6.91499E+000 6.91499E-003 21 | 18 0.00000E+00 0.00000E+000 1.80613E+003 7.18151E+000 6.84280E-003 22 | 19 0.00000E+00 0.00000E+000 1.82109E+003 7.43972E+000 6.77285E-003 23 | 20 0.00000E+00 0.00000E+000 1.83505E+003 7.69015E+000 6.70491E-003 24 | -------------------------------------------------------------------------------- /Param/CORONA/TestOutput/log_n000001.awsom.bvector: -------------------------------------------------------------------------------- 1 | Volume averages, fluxes, etc 2 | it t dt absjxb emag1 ekin 3 | 1 0.00000E+00 0.00000E+000 2.65030E+003 2.55741E-003 9.04846E-006 4 | 2 0.00000E+00 0.00000E+000 2.56204E+003 2.72439E-002 9.54488E-006 5 | 3 0.00000E+00 0.00000E+000 2.45385E+003 1.02502E-001 1.02085E-005 6 | 4 0.00000E+00 0.00000E+000 2.33602E+003 2.52529E-001 1.09695E-005 7 | 5 0.00000E+00 0.00000E+000 2.21960E+003 4.93787E-001 1.17824E-005 8 | 6 0.00000E+00 0.00000E+000 2.11154E+003 8.35328E-001 1.26185E-005 9 | 7 0.00000E+00 0.00000E+000 2.01418E+003 1.28001E+000 1.34605E-005 10 | 8 0.00000E+00 0.00000E+000 1.93070E+003 1.82627E+000 1.42993E-005 11 | 9 0.00000E+00 0.00000E+000 1.86926E+003 2.46997E+000 1.51307E-005 12 | 10 0.00000E+00 0.00000E+000 1.81637E+003 3.20503E+000 1.59533E-005 13 | 11 0.00000E+00 0.00000E+000 1.77223E+003 4.02464E+000 1.67674E-005 14 | 12 0.00000E+00 0.00000E+000 1.73943E+003 4.92081E+000 1.75741E-005 15 | 13 0.00000E+00 0.00000E+000 1.71697E+003 5.88450E+000 1.83747E-005 16 | 14 0.00000E+00 0.00000E+000 1.70306E+003 6.90620E+000 1.91706E-005 17 | 15 0.00000E+00 0.00000E+000 1.69876E+003 7.97608E+000 1.99631E-005 18 | 16 0.00000E+00 0.00000E+000 1.70428E+003 9.08381E+000 2.07539E-005 19 | 17 0.00000E+00 0.00000E+000 1.71441E+003 1.02187E+001 2.15442E-005 20 | 18 0.00000E+00 0.00000E+000 1.73500E+003 1.13709E+001 2.23351E-005 21 | 19 0.00000E+00 0.00000E+000 1.76482E+003 1.25310E+001 2.31288E-005 22 | 20 0.00000E+00 0.00000E+000 1.80503E+003 1.36898E+001 2.39276E-005 23 | -------------------------------------------------------------------------------- /Param/CORONA/TestOutput/los_dem_1_n0000004.out.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/CORONA/TestOutput/los_dem_1_n0000004.out.gz -------------------------------------------------------------------------------- /Param/CORONA/TestOutput/los_fux_2_n0000004.out.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/CORONA/TestOutput/los_fux_2_n0000004.out.gz -------------------------------------------------------------------------------- /Param/CORONA/TestOutput/los_nbi_3_n0000004.out.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/CORONA/TestOutput/los_nbi_3_n0000004.out.gz -------------------------------------------------------------------------------- /Param/CORONA/TestOutput/los_phx_4_n0000004.out.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/CORONA/TestOutput/los_phx_4_n0000004.out.gz -------------------------------------------------------------------------------- /Param/CORONA/TestOutput/pcl_01__n000001.pcl.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/CORONA/TestOutput/pcl_01__n000001.pcl.gz -------------------------------------------------------------------------------- /Param/CORONA/TestOutput/y=0_nbi_6_n00000004.out.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/CORONA/TestOutput/y=0_nbi_6_n00000004.out.gz -------------------------------------------------------------------------------- /Param/CORONA/harmonics-svsm.out.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/CORONA/harmonics-svsm.out.gz -------------------------------------------------------------------------------- /Param/CRASH/1D.in: -------------------------------------------------------------------------------- 1 | #OUTERBOUNDARY 2 | float xMinBc 3 | float xMaxBc 4 | reflect yMinBc 5 | reflect yMaxBc 6 | reflect zMinBc 7 | reflect zMaxBc 8 | 9 | -------------------------------------------------------------------------------- /Param/CRASH/3D.in: -------------------------------------------------------------------------------- 1 | #OUTERBOUNDARY 2 | float xMinBc 3 | float xMaxBc 4 | reflect yMinBc 5 | float yMaxBc 6 | reflect zMinBc 7 | float zMaxBc 8 | 9 | #GRIDSYMMETRY 10 | F IsMirrorX 11 | T IsMirrorY 12 | T IsMirrorZ 13 | 14 | -------------------------------------------------------------------------------- /Param/CRASH/AMR1.in: -------------------------------------------------------------------------------- 1 | #USERINPUTBEGIN ----------------------------------- 2 | 3 | #USERAMR 4 | 20.0 RhoMinAmr 5 | 2500.0 xMaxAmr 6 | 1.0 BetaProlong 7 | 8 | #USERINPUTEND ------------------------------------- 9 | 10 | #AMRINITPHYSICS 11 | 2 nRefineLevelIC 12 | 13 | #AMRLEVELS 14 | 0 MinBlockLevel 15 | 1 MaxBlockLevel 16 | 17 | #DOAMR 18 | T DoAmr 19 | 1 DnAmr 20 | -1 DtAmr 21 | T IsStrictAmr 22 | 23 | #AMRCRITERIA 24 | 1 nRefineCrit 25 | user TypeRefine 26 | 0.5 CoarsenLimit 27 | 0.5 RefineLimit 28 | -------------------------------------------------------------------------------- /Param/CRASH/AMR2.in: -------------------------------------------------------------------------------- 1 | #USERINPUTBEGIN ---------------------- 2 | 3 | #USERAMR 4 | 20.0 RhoMinAmr 5 | 2500.0 xMaxAmr 6 | 1.0 BetaProlong 7 | 8 | #USERINPUTEND ---------------------- 9 | 10 | #AMRINITPHYSICS 11 | 3 nRefineLevelIC 12 | 13 | #AMRLEVELS 14 | 0 MinBlockLevel 15 | 2 MaxBlockLevel 16 | 17 | #DOAMR 18 | T DoAmr 19 | 1 DnAmr 20 | -1 DtAmr 21 | T IsStrictAmr 22 | 23 | #AMRCRITERIA 24 | 1 nRefineCrit 25 | user TypeRefine 26 | 0.5 LimitCoarsen 27 | 0.5 LimitRefine 28 | -------------------------------------------------------------------------------- /Param/CRASH/AMR3.in: -------------------------------------------------------------------------------- 1 | #USERINPUTBEGIN ---------------------- 2 | 3 | #USERAMR 4 | 20.0 RhoMinAmr 5 | 2500.0 xMaxAmr 6 | 1.0 BetaProlong 7 | 8 | #USERINPUTEND ---------------------- 9 | 10 | #AMRINITPHYSICS 11 | 4 nRefineLevelIC 12 | 13 | #AMRLEVELS 14 | 0 MinBlockLevel 15 | 3 MaxBlockLevel 16 | 17 | #DOAMR 18 | T DoAmr 19 | 1 DnAmr 20 | -1 DtAmr 21 | T IsStrictAmr 22 | 23 | #AMRCRITERIA 24 | 1 nRefineCrit 25 | user TypeRefine 26 | 0.5 LimitCoarsen 27 | 0.5 LimitRefine 28 | -------------------------------------------------------------------------------- /Param/CRASH/AMR4.in: -------------------------------------------------------------------------------- 1 | #USERINPUTBEGIN ---------------------- 2 | 3 | #USERAMR 4 | 20.0 RhoMinAmr 5 | 2500.0 xMaxAmr 6 | 1.0 BetaProlong 7 | 8 | #USERINPUTEND ---------------------- 9 | 10 | #AMRINITPHYSICS 11 | 5 nRefineLevelIC 12 | 13 | #AMRLEVELS 14 | 0 MinBlockLevel 15 | 4 MaxBlockLevel 16 | 17 | #DOAMR 18 | T DoAmr 19 | 1 DnAmr 20 | -1 DtAmr 21 | T IsStrictAmr 22 | 23 | #AMRCRITERIA 24 | 1 nRefineCrit 25 | user TypeRefine 26 | 0.5 LimitCoarsen 27 | 0.5 LimitRefine 28 | -------------------------------------------------------------------------------- /Param/CRASH/GRAY.in: -------------------------------------------------------------------------------- 1 | #LOOKUPTABLE 2 | Opacity(rho,T) NameTable 3 | load NameCommand 4 | table_opacities.out NameFile 5 | ascii TypeFile 6 | 7 | #SEMIIMPLICIT 8 | T UseSemiImplicit 9 | radiation TypeSemiImplicit 10 | 11 | #SEMIKRYLOV 12 | CG TypeKrylov 13 | 1.0e-8 ErrorMaxKrylov 14 | 900 MaxMatvecKrylov 15 | 16 | #RADIATION 17 | T UseRadDiffusion 18 | T UseRadFluxLimiter 19 | larsen TypeRadFluxLimiter 20 | 300.0 TradMinSi 21 | 22 | 23 | -------------------------------------------------------------------------------- /Param/CRASH/GRAYCOND.in: -------------------------------------------------------------------------------- 1 | #LOOKUPTABLE 2 | Opacity(rho,T) NameTable 3 | load NameCommand 4 | table_opacities.out NameFile 5 | ascii TypeFile 6 | 7 | #SEMIIMPLICIT 8 | T UseSemiImplicit 9 | radcond TypeSemiImplicit 10 | 11 | #SEMIKRYLOV 12 | CG TypeKrylov 13 | 1.0e-8 ErrorMaxKrylov 14 | 900 MaxMatvecKrylov 15 | 16 | #SEMIPRECONDITIONER 17 | T DoPrecond 18 | DILU TypePrecond 19 | 20 | #RADIATION 21 | T UseRadDiffusion 22 | T UseRadFluxLimiter 23 | larsen TypeRadFluxLimiter 24 | 300.0 TradMinSi 25 | 26 | #ELECTRONPRESSURE 27 | 1.1e5 PeMinSi 28 | 29 | 30 | -------------------------------------------------------------------------------- /Param/CRASH/GRAYCONDHYPRE.in: -------------------------------------------------------------------------------- 1 | #LOOKUPTABLE 2 | Opacity(rho,T) NameTable 3 | load NameCommand 4 | table_opacities.out NameFile 5 | ascii TypeFile 6 | 7 | #SEMIIMPLICIT 8 | T UseSemiImplicit 9 | radcond TypeSemiImplicit 10 | 11 | #SEMIKRYLOV 12 | GMRES TypeKrylov 13 | 1.0e-8 ErrorMaxKrylov 14 | 100 MaxMatvecKrylov 15 | 16 | #SEMIPRECONDITIONER 17 | T DoPrecond 18 | HYPRE TypePrecond 19 | 20 | #RADIATION 21 | T UseRadDiffusion 22 | T UseRadFluxLimiter 23 | larsen TypeRadFluxLimiter 24 | 300.0 TradMinSi 25 | 26 | #ELECTRONPRESSURE 27 | 1.1e5 PeMinSi 28 | 29 | #TEST 30 | krylov StringTest 31 | 32 | -------------------------------------------------------------------------------- /Param/CRASH/GRAYHYPRE.in: -------------------------------------------------------------------------------- 1 | #LOOKUPTABLE 2 | Opacity(rho,T) NameTable 3 | load NameCommand 4 | table_opacities.out NameFile 5 | ascii TypeFile 6 | 7 | #SEMIIMPLICIT 8 | T UseSemiImplicit 9 | radiation TypeSemiImplicit 10 | 11 | #SEMIKRYLOV 12 | GMRES TypeKrylov 13 | 1.0e-8 ErrorMaxKrylov 14 | 100 MaxMatvecKrylov 15 | 16 | #SEMIPRECONDITIONER 17 | T DoPrecond 18 | HYPRE TypePrecond 19 | 20 | #RADIATION 21 | T UseRadDiffusion 22 | T UseRadFluxLimiter 23 | larsen TypeRadFluxLimiter 24 | 300.0 TradMinSi 25 | 26 | #TEST 27 | krylov StringTest 28 | 29 | -------------------------------------------------------------------------------- /Param/CRASH/ICS_STANDARD.in: -------------------------------------------------------------------------------- 1 | #USERINPUTBEGIN ----------------------------------- 2 | 3 | #STATEDEFINITION 4 | rho p StringStateVar 5 | 7 nMaterialState 6 | M1 0.0065 0.00012 Name, State 7 | M21 6.5e-6 1.7491328e-06 Name, State 8 | M22 1.85 0.49783012 Name, State 9 | M31 1.43 0.35375 Name, State 10 | M32 0.0065 0.01 Name, State 11 | M4 20.0 0.24625032 Name, State 12 | M5 1.13 0.41057944 Name, State 13 | 14 | #STATEINTERFACE 15 | 11 nSegment 16 | M21 M22 0.0, 0.0, 0.0, 1000.0 Material1 Material2 x1 y1 x2 y2 17 | M22 M1 20.0, 0.0, 20.0, 287.5 Material1 Material2 x1 y1 x2 y2 18 | M22 M4 20.0, 287.5, 20.0, 1000.0 Material1 Material2 x1 y1 x2 y2 19 | M4 M1 20.0, 287.5, 100.0, 287.5 Material1 Material2 x1 y1 x2 y2 20 | M4 M5 100.0, 287.5, 100.0, 1000.0 Material1 Material2 x1 y1 x2 y2 21 | M5 M1 100.0, 287.5, 300.0, 287.5 Material1 Material2 x1 y1 x2 y2 22 | M5 M31 300.0, 287.5, 300.0, 312.5 Material1 Material2 x1 y1 x2 y2 23 | M5 M31 300.0, 312.5, 900.0, 312.5 Material1 Material2 x1 y1 x2 y2 24 | M5 M32 900.0, 312.5, 900.0, 1000.0 Material1 Material2 x1 y1 x2 y2 25 | M32 M31 900.0, 312.5, 5000.0, 312.5 Material1 Material2 x1 y1 x2 y2 26 | M31 M1 300.0, 287.5, 5000.0, 287.5 Material1 Material2 x1 y1 x2 y2 27 | 28 | 29 | STATEDEFINITION 30 | rho p StringStateVar 31 | 6 nMaterialState 32 | Xe 0.0065 0.00012 Name, State 33 | Be1 6.5e-6 1.7491328e-06 Name, State 34 | Be2 1.85 0.49783012 Name, State 35 | Pl 1.43 0.35375 Name, State 36 | Au 20.0 0.24625032 Name, State 37 | Ay 1.13 0.41057944 Name, State 38 | 39 | STATEINTERFACE 40 | 10 nSegment 41 | Be1 Be2 0.0, 0.0, 0.0, 1000.0 Material1 Material2 x1 y1 x2 y2 42 | Be2 Xe 20.0, 0.0, 20.0, 287.5 Material1 Material2 x1 y1 x2 y2 43 | Be2 Au 20.0, 287.5, 20.0, 1000.0 Material1 Material2 x1 y1 x2 y2 44 | Au Xe 20.0, 287.5, 100.0, 287.5 Material1 Material2 x1 y1 x2 y2 45 | Au Ay 100.0, 287.5, 100.0, 1000.0 Material1 Material2 x1 y1 x2 y2 46 | Ay Xe 100.0, 287.5, 300.0, 287.5 Material1 Material2 x1 y1 x2 y2 47 | Ay Pl 300.0, 287.5, 300.0, 312.5 Material1 Material2 x1 y1 x2 y2 48 | Ay Pl 300.0, 312.5, 900.0, 312.5 Material1 Material2 x1 y1 x2 y2 49 | Ay Pl 900.0, 312.5, 900.0, 1000.0 Material1 Material2 x1 y1 x2 y2 50 | Pl Xe 300.0, 287.5, 5000.0, 287.5 Material1 Material2 x1 y1 x2 y2 51 | 52 | #USERINPUTEND ------------------------------------- 53 | -------------------------------------------------------------------------------- /Param/CRASH/MULTIGROUP.in: -------------------------------------------------------------------------------- 1 | #USERINPUTBEGIN ---------------------- 2 | 3 | #OPACITY 4 | T UseExcitation 5 | F UseCoulombCorrection 6 | T DoStateElimination 7 | 8 | #GROUPRANGE 9 | 0.1 EnergyPhotonMin 10 | 20000.0 EnergyPhotonMax 11 | 12 | #USERINPUTEND ---------------------- 13 | 14 | #SEMIIMPLICIT 15 | T UseSemiImplicit 16 | radiation TypeSemiImplicit 17 | 18 | #SEMIKRYLOV 19 | CG TypeKrylov 20 | 1.0e-5 ErrorMaxKrylov 21 | 200 MaxMatvecKrylov 22 | 23 | #RADIATION 24 | T UseRadDiffusion 25 | T UseRadFluxLimiter 26 | larsen TypeRadFluxLimiter 27 | 300.0 TradMinSi 28 | 29 | -------------------------------------------------------------------------------- /Param/CRASH/MULTIGROUPCOND.in: -------------------------------------------------------------------------------- 1 | #USERINPUTBEGIN ----------------------------------- 2 | 3 | #OPACITY 4 | T UseExcitation 5 | F UseCoulombCorrection 6 | T DoStateElimination 7 | 8 | #GROUPRANGE 9 | 0.1 EnergyPhotonMin 10 | 20000.0 EnergyPhotonMax 11 | 12 | #USERINPUTEND ------------------------------------- 13 | 14 | #SEMIIMPLICIT 15 | T UseSemiImplicit 16 | radcond TypeSemiImplicit 17 | 18 | #SEMIKRYLOV 19 | CG TypeKrylov 20 | 1.0e-5 ErrorMaxKrylov 21 | 200 MaxMatvecKrylov 22 | 23 | #RADIATION 24 | T UseRadDiffusion 25 | T UseRadFluxLimiter 26 | larsen TypeRadFluxLimiter 27 | 300.0 TradMinSi 28 | 29 | #ELECTRONPRESSURE 30 | 1.0 PeMinSi 31 | -------------------------------------------------------------------------------- /Param/CRASH/MULTIGROUPCONDHYPRE.in: -------------------------------------------------------------------------------- 1 | #USERINPUTBEGIN ---------------------- 2 | 3 | #OPACITY 4 | T UseExcitation 5 | F UseCoulombCorrection 6 | T DoStateElimination 7 | 8 | #GROUPRANGE 9 | 0.1 EnergyPhotonMin 10 | 20000.0 EnergyPhotonMax 11 | 12 | #USERINPUTEND ---------------------- 13 | 14 | #SEMIIMPLICIT 15 | T UseSemiImplicit 16 | radcond TypeSemiImplicit 17 | 18 | #SEMIKRYLOV 19 | GMRES TypeKrylov 20 | 1.0e-8 ErrorMaxKrylov 21 | 100 MaxMatvecKrylov 22 | 23 | #SEMIPRECONDITIONER 24 | T DoPrecond 25 | HYPRE TypePrecond 26 | 27 | #RADIATION 28 | T UseRadDiffusion 29 | T UseRadFluxLimiter 30 | larsen TypeRadFluxLimiter 31 | 300.0 TradMinSi 32 | 33 | #ELECTRONPRESSURE 34 | 1.0 PeMinSi 35 | -------------------------------------------------------------------------------- /Param/CRASH/MULTIGROUPHYPRE.in: -------------------------------------------------------------------------------- 1 | #USERINPUTBEGIN ---------------------- 2 | 3 | #OPACITY 4 | T UseExcitation 5 | F UseCoulombCorrection 6 | T DoStateElimination 7 | 8 | #GROUPRANGE 9 | 0.1 EnergyPhotonMin 10 | 20000.0 EnergyPhotonMax 11 | 12 | #USERINPUTEND ---------------------- 13 | 14 | #SEMIIMPLICIT 15 | T UseSemiImplicit 16 | radiation TypeSemiImplicit 17 | 18 | #SEMIKRYLOV 19 | gmres TypeKrylov 20 | 1.0e-8 ErrorMaxKrylov 21 | 100 MaxMatvecKrylov 22 | 23 | #SEMIPRECONDITIONER 24 | T DoPrecond 25 | HYPRE TypePrecond 26 | 27 | #RADIATION 28 | T UseRadDiffusion 29 | T UseRadFluxLimiter 30 | larsen TypeRadFluxLimiter 31 | 300.0 TradMinSi 32 | 33 | -------------------------------------------------------------------------------- /Param/CRASH/PARAM.in.eosgodunov: -------------------------------------------------------------------------------- 1 | #DESCRIPTION 2 | Param/CRASH/PARAM.in.eosgodunov StringDescription 3 | 4 | #INCLUDE 5 | Param/CRASH/START.in 6 | 7 | #SAVEINITIAL 8 | F DoSaveInitial 9 | 10 | #USERSWITCH 11 | +update StringSwitch 12 | 13 | #CHECKGRIDSIZE 14 | 10 nI 15 | 1 nJ 16 | 1 nK 17 | 8 MinBlockAll 18 | 19 | #GRID 20 | 8 nRootBlock1 21 | 1 nRootBlock2 22 | 1 nRootBlock3 23 | 0.0 xMin 24 | 5000.0 xMax 25 | -1.0 yMin 26 | 1.0 yMax 27 | -1.0 zMin 28 | 1.0 zMax 29 | 30 | #SHOCKTUBE 31 | 100. LeftState Rho !~30 times the normal density 32 | -3.0e4 LeftState Ux !Piston speed 33 | 0.0 LeftState Uy 34 | 0.0 LeftState Uz 35 | 0.0 LeftState Eint 36 | 2.04e+7 LeftState P !3000 K of temperature 37 | 100. RightState rho 38 | -3.0e4 RightState Ux 39 | 0.0 RightState Uy 40 | 0.0 RightState Uz 41 | 0.0 RightState Eint 42 | 2.04e+7 RightState P 43 | 44 | Reflective boundary at the left 45 | #OUTERBOUNDARY 46 | reflect TypeBc1 47 | float TypeBc2 48 | 49 | #SCHEME 50 | 2 nOrder 51 | Godunov TypeFlux 52 | beta TypeLimiter 53 | 1.2 LimiterBeta 54 | 55 | #SAVERESTART 56 | F DoSaveRestart 57 | 58 | #SAVEPLOT 59 | 1 nPlotFile 60 | cut HD idl_ascii StringPlot 61 | -1 DnSavePlot 62 | 1.0e-8 DtSavePlot 63 | 0.0 Coord1MinCut 64 | 5.0e3 Coord1MaxCut 65 | -1.0e-7 Coord2MinCut 66 | +1.0e-7 Coord2MaxCut 67 | -1.0e-7 Coord3MinCut 68 | +1.0e-7 Coord3MaxCut 69 | 0.0 DxSavePlot 70 | 71 | #STOP 72 | -1 MaxIteration 73 | 5.0e-7 tSimulationMax 74 | 75 | -------------------------------------------------------------------------------- /Param/CRASH/RESTART.in: -------------------------------------------------------------------------------- 1 | #INCLUDE 2 | restartIN/restart.H 3 | 4 | #LIMITER 5 | F UseLogRhoLimiter 6 | F UseLogPLimiter 7 | T UseRhoRatioLimiter 8 | EInt NameVarLimitRatio 9 | 10 | #RESCHANGE 11 | T UseAccurateResChange 12 | 13 | #PLOTFILENAME 14 | nanosecond NameMaxTimeUnit 15 | -------------------------------------------------------------------------------- /Param/CRASH/RZ.in: -------------------------------------------------------------------------------- 1 | #GRIDGEOMETRY 2 | rz TypeGeometry 3 | 4 | #OUTERBOUNDARY 5 | float TypeBc1 6 | float TypeBc2 7 | reflect TypeBc3 8 | float TypeBc4 9 | -------------------------------------------------------------------------------- /Param/CRASH/START.in: -------------------------------------------------------------------------------- 1 | #PLANET 2 | NONE NamePlanet 3 | 4 | #IOUNITS 5 | SI TypeIoUnit 6 | 7 | #NORMALIZATION 8 | READ TypeNormalization 9 | 1.0e-6 No2SiUnitX Micron 10 | 1.0e+3 No2SiUnitU km/s 11 | 1000.0 No2SiUnitRho g/cm^3 12 | 13 | #LIMITER 14 | F UseLogRhoLimiter 15 | F UseLogPLimiter 16 | T UseRhoRatioLimiter 17 | EInt NameVarLimitRatio 18 | 19 | #RESCHANGE 20 | T UseAccurateReschange 21 | 22 | #SAVEINITIAL 23 | T DoSaveInitial 24 | 25 | #PLOTFILENAME 26 | nanosecond NameMaxTimeUnit 27 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Au.in: -------------------------------------------------------------------------------- 1 | #INCLUDE 2 | Param/CRASH/TABLE_Au_EOS.in 3 | 4 | #INCLUDE 5 | Param/CRASH/TABLE_Au_OPACITY.in 6 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Au_EOS.in: -------------------------------------------------------------------------------- 1 | #LOOKUPTABLE 2 | Au_eos NameTable 3 | use param NameCommand 4 | Tables/Au_eos_CRASH.dat NameFile 5 | real8 TypeFile 6 | Mass NameTableParam 7 | 196.96656 TableParam [amu] 8 | CRASH EOS for Au StringDescription 9 | logTe logNa P E Pe Ee Cv Cve Gamma GammaE TeTi Cond Z Z2 DPOverDRho DPOverDT DPEOverDRho DPEOverDT NameVar 10 | 2 nIndex 11 | 201 nIndex1 12 | 1.0E-3 Index1Min [eV] 13 | 1.0E+2 Index1Max [eV] 14 | 201 nIndex2 15 | 1.0E+24 Index2Min [m-3] 16 | 1.2E+29 Index2Max [m-3] 17 | 18 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Au_OPACITY.in: -------------------------------------------------------------------------------- 1 | #LOOKUPTABLE 2 | Au_opac NameTable 3 | use param NameCommand 4 | Tables/Au_opac_CRASH.dat NameFile 5 | real8 TypeFile 6 | EvMin EvMax NameTableParam 7 | 0.1 TableParam [eV] 8 | 20000.0 TableParam [eV] 9 | Opacity(rho,Te) for Au StringDescription 10 | logrho logT Planck(30) Ross(30) NameVar 11 | 2 nIndex 12 | 201 nIndex1 13 | 0.3270688818768 Index1Min [kg/m3] 14 | 3.9248265825216E+04 Index1Max [kg/m3] 15 | 201 nIndex2 16 | 2.0E-1 Index2Min [eV] 17 | 1.0E+2 Index2Max [eV] 18 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Ay.in: -------------------------------------------------------------------------------- 1 | #INCLUDE 2 | Param/CRASH/TABLE_Ay_EOS.in 3 | 4 | #INCLUDE 5 | Param/CRASH/TABLE_Ay_OPACITY.in 6 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Ay_EOS.in: -------------------------------------------------------------------------------- 1 | #LOOKUPTABLE 2 | Ay_eos NameTable 3 | use param NameCommand 4 | Tables/Ay_eos_CRASH.dat NameFile 5 | real8 TypeFile 6 | Mass NameTableParam 7 | 6.6745291 TableParam [amu] 8 | CRASH EOS for Ay StringDescription 9 | logTe logNa P E Pe Ee Cv Cve Gamma GammaE TeTi Cond Z Z2 DPOverDRho DPOverDT DPEOverDRho DPEOverDT NameVar 10 | 2 nIndex 11 | 201 nIndex1 12 | 1.0E-3 Index1Min [eV] 13 | 1.0E+2 Index1Max [eV] 14 | 201 nIndex2 15 | 1.0E+24 Index2Min [m-3] 16 | 1.5E+29 Index2Max [m-3] 17 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Ay_OPACITY.in: -------------------------------------------------------------------------------- 1 | #LOOKUPTABLE 2 | Ay_opac NameTable 3 | use param NameCommand 4 | Tables/Ay_opac_CRASH.dat NameFile 5 | real8 TypeFile 6 | EvMin EvMax NameTableParam 7 | 0.1 TableParam [eV] 8 | 20000.0 TableParam [eV] 9 | Opacity(rho,Te) for Ay StringDescription 10 | logrho logT Planck(30) Ross(30) NameVar 11 | 2 nIndex 12 | 201 nIndex1 13 | 1.1083255751072E-02 Index1Min [kg/m3] 14 | 1.6624883626608E+03 Index1Max [kg/m3] 15 | 201 nIndex2 16 | 5.0E-2 Index2Min [eV] 17 | 1.0E+2 Index2Max [eV] 18 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Be.in: -------------------------------------------------------------------------------- 1 | #INCLUDE 2 | Param/CRASH/TABLE_Be_EOS.in 3 | 4 | #INCLUDE 5 | Param/CRASH/TABLE_Be_OPACITY.in 6 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Be_EOS.in: -------------------------------------------------------------------------------- 1 | #LOOKUPTABLE 2 | Be_eos NameTable 3 | use param NameCommand 4 | Tables/Be_eos_CRASH.dat NameFile 5 | real8 TypeFile 6 | Mass NameTableParam 7 | 9.0121823 TableParam [amu] 8 | CRASH EOS for Be StringDescription 9 | logTe logNa P E Pe Ee Cv Cve Gamma GammaE TeTi Cond Z Z2 DPOverDRho DPOverDT DPEOverDRho DPEOverDT NameVar 10 | 2 nIndex 11 | 201 nIndex1 12 | 1.0E-3 Index1Min [eV] 13 | 4.0E+3 Index1Max [eV] 14 | 201 nIndex2 15 | 1.0E+23 Index2Min [m-3] 16 | 2.0E+29 Index2Max [m-3] 17 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Be_OPACITY.in: -------------------------------------------------------------------------------- 1 | #LOOKUPTABLE 2 | Be_opac NameTable 3 | use param NameCommand 4 | Tables/Be_opac_CRASH.dat NameFile 5 | real8 TypeFile 6 | EvMin EvMax NameTableParam 7 | 0.1 TableParam [eV] 8 | 20000.0 TableParam [eV] 9 | Opacity(rho,Te) for Be StringDescription 10 | logrho logT Planck(30) Ross(30) NameVar 11 | 2 nIndex 12 | 201 nIndex1 13 | 0.0014965 Index1Min [kg/m3] 14 | 2993.0 Index1Max [kg/m3] 15 | 201 nIndex2 16 | 3.0E-2 Index2Min [eV] 17 | 4.0E+3 Index2Max [eV] 18 | 19 | 20 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Pl.in: -------------------------------------------------------------------------------- 1 | #INCLUDE 2 | Param/CRASH/TABLE_Pl_EOS.in 3 | 4 | #INCLUDE 5 | Param/CRASH/TABLE_Pl_OPACITY.in 6 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Pl_EOS.in: -------------------------------------------------------------------------------- 1 | #LOOKUPTABLE 2 | Pl_eos NameTable 3 | use param NameCommand 4 | Tables/Pl_eos_CRASH.dat NameFile 5 | real8 TypeFile 6 | Mass NameTableParam 7 | 9.8034439 TableParam [amu] 8 | CRASH EOS for Pl StringDescription 9 | logTe logNa P E Pe Ee Cv Cve Gamma GammaE TeTi Cond Z Z2 DPOverDRho DPOverDT DPEOverDRho DPEOverDT NameVar 10 | 2 nIndex 11 | 201 nIndex1 12 | 1.0E-3 Index1Min [eV] 13 | 1.0E+2 Index1Max [eV] 14 | 201 nIndex2 15 | 1.0E+24 Index2Min [m-3] 16 | 1.5E+29 Index2Max [m-3] 17 | 18 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Pl_OPACITY.in: -------------------------------------------------------------------------------- 1 | #LOOKUPTABLE 2 | Pl_opac NameTable 3 | use param NameCommand 4 | Tables/Pl_opac_CRASH.dat NameFile 5 | real8 TypeFile 6 | EvMin EvMax NameTableParam 7 | 0.1 TableParam [eV] 8 | 20000.0 TableParam [eV] 9 | Opacity(rho,Te) for Pl StringDescription 10 | logrho logT Planck(30) Ross(30) NameVar 11 | 2 nIndex 12 | 201 nIndex1 13 | 1.62789127801646E-02 Index1Min [kg/m3] 14 | 2.44183691702469E+03 Index1Max [kg/m3] 15 | 201 nIndex2 16 | 5.0E-2 Index2Min [eV] 17 | 1.0E+2 Index2Max [eV] 18 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Xe.in: -------------------------------------------------------------------------------- 1 | #INCLUDE 2 | Param/CRASH/TABLE_Xe_EOS.in 3 | 4 | #INCLUDE 5 | Param/CRASH/TABLE_Xe_OPACITY.in 6 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Xe_EOS.in: -------------------------------------------------------------------------------- 1 | #LOOKUPTABLE 2 | Xe_eos NameTable 3 | use param NameCommand 4 | Tables/Xe_eos_CRASH.dat NameFile 5 | real8 TypeFile 6 | Mass NameTableParam 7 | 131.292 TableParam [amu] 8 | CRASH EOS for Xe StringDescription 9 | logTe logNa P E Pe Ee Cv Cve Gamma GammaE TeTi Cond Z Z2 DPOverDRho DPOverDT DPEOverDRho DPEOverDT NameVar 10 | 2 nIndex 11 | 201 nIndex1 12 | 1.0E-2 Index1Min [eV] 13 | 1.0E+3 Index1Max [eV] 14 | 201 nIndex2 15 | 1.0E+24 Index2Min [m-3] 16 | 1.0E+29 Index2Max [m-3] 17 | -------------------------------------------------------------------------------- /Param/CRASH/TABLE_Xe_OPACITY.in: -------------------------------------------------------------------------------- 1 | #LOOKUPTABLE 2 | Xe_opac NameTable 3 | use param NameCommand 4 | Tables/Xe_opac_CRASH.dat NameFile 5 | real8 TypeFile 6 | EvMin EvMax NameTableParam 7 | 0.1 TableParam [eV] 8 | 20000.0 TableParam [eV] 9 | Opacity(rho,Te) for Xe StringDescription 10 | logrho logT Planck(30) Ross(30) NameVar 11 | 2 nIndex 12 | 201 nIndex1 13 | 0.21801430476 Index1Min [kg/m3] 14 | 2.1801430476E+04 Index1Max [kg/m3] 15 | 201 nIndex2 16 | 3.0E-2 Index2Min [eV] 17 | 1.0E+3 Index2Max [eV] 18 | -------------------------------------------------------------------------------- /Param/CRASH/initial_lowrie3.dat.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/CRASH/initial_lowrie3.dat.gz -------------------------------------------------------------------------------- /Param/CURRENT/sat.dat: -------------------------------------------------------------------------------- 1 | Satellite in a box 2 | 3 | #START 4 | 00 3 21 10 45 0 1 -5.0 0.0 0.0 5 | 00 3 21 10 45 10 0 +5.0 0.0 0.0 6 | 00 3 22 23 59 59 0 +8.0 0.0 0.0 7 | -------------------------------------------------------------------------------- /Param/CYLOUTFLOW/PARAM.in: -------------------------------------------------------------------------------- 1 | Config.pl -e=Hd -u=CylOutflow -g=124,512,1 2 | 3 | #DESCRIPTION 4 | Param/CYLOUTFLOW/PARAM.in: Radial Outflow StringDescription 5 | 6 | #PLANET 7 | NONE NamePlanet 8 | 9 | #IOUNITS 10 | NONE TypeIoUnit 11 | 12 | #NORMALIZATION 13 | NONE TypeNormalization 14 | 15 | #GRIDBLOCKALL 16 | 2 MaxBlock 17 | 18 | #GRID 19 | 1 nRootBlock1 20 | 1 nRootBlock2 21 | 1 nRootBlock3 22 | -25.0 xMin 23 | 25.0 xMax 24 | -25.0 yMin 25 | 25.0 yMax 26 | -0.5 zMin 27 | 0.5 zMax 28 | 29 | #GRIDGEOMETRY 30 | cylindrical TypeGeometry 31 | 32 | #LIMITRADIUS 33 | 4.13804 rMin 34 | 23.5433 rMax 35 | 36 | BUMP 37 | Rho NameVar 38 | 10.0 Amplitude 39 | 3.0 WidthX 40 | 3.0 WidthY 41 | -1.0 WidthZ 42 | 5.0 CenterX 43 | 6.0 CenterY 44 | 0.0 CenterZ 45 | 2 nPower 46 | 47 | #USERINPUTBEGIN ----------------------------------- 48 | 49 | #OUTFLOW 50 | 10.0 u0 51 | 0.0 u1 52 | -1.0 Exponent 53 | 54 | #USERINPUTEND ------------------------------------- 55 | 56 | #UPDATEVAR 57 | Rho StringVarUpdate 58 | 59 | #OUTERBOUNDARY 60 | user TypeBc1 rMin 61 | float TypeBc2 rMax 62 | periodic TypeBc3 PhiMin 63 | periodic TypeBc4 PhiMax 64 | 65 | #SCHEME 66 | 2 nOrder 67 | Rusanov TypeFlux 68 | mc3 TypeLimiter 69 | 2.0 LimiterBeta 70 | 71 | #NONCONSERVATIVE 72 | T UseNonConservative 73 | 74 | #SAVELOGFILE 75 | T DoSaveLogfile 76 | raw StringLog 77 | 1 DnSaveLogfile 78 | -1. DtSaveLogfile 79 | 80 | #SAVERESTART 81 | F DoSaveRestart 82 | 83 | #PLOTFILENAME 84 | second NameMaxTimeUnit 85 | 86 | #SAVEINITIAL 87 | T DoSaveInitial 88 | 89 | #SAVEPLOT 90 | 1 nPlotFile 91 | z=0 hd idl StringPlot 92 | -1 DnSavePlot 93 | 0.1 DtSavePlot 94 | 0.0 DxSavePlot 95 | 0.0 DxSavePlot 96 | 97 | #STOP 98 | -1 MaxIteration 99 | 1.0 tSimulationMax 100 | -------------------------------------------------------------------------------- /Param/EARTH/PARAM.in.L1toBC: -------------------------------------------------------------------------------- 1 | #DESCRIPTION 2 | Param/EARTH/PARAM.in.L1toBC: Solar wind propagation 3 | 4 | #COORDSYSTEM 5 | GSE NameCoordSystem 6 | 7 | #B0 8 | F UseB0 9 | 10 | #STARTTIME 11 | 1998 iYear 12 | 5 iMonth 13 | 4 iDay 14 | 0 iHour 15 | 45 iMinute 16 | 0 iSecond 17 | 18 | #SAVERESTART 19 | F DoSaveRestart 20 | 21 | #GRID 22 | 40 nRootBlockX 23 | 1 nRootBlockY 24 | 1 nRootBlockZ 25 | +31. xMin 26 | +235. xMax 27 | -1. yMin 28 | +1. yMax 29 | -1. zMin 30 | +1. zMax 31 | 32 | #GRIDBLOCKALL 33 | 700 MaxBlock 34 | 35 | #SOLARWINDFILE 36 | T UseSolarWindFile 37 | L1.dat NameSolarWindFile 38 | 39 | #OUTERBOUNDARY 40 | outflow TypeBc1 41 | inflow TypeBc2 42 | 43 | #SCHEME 44 | 5 nOrder 45 | Linde TypeFlux 46 | mc3 TypeLimiter 47 | 1.5 BetaLimiter 48 | 49 | ! Propagate Bx with the flow 50 | #DIVB 51 | T UseDivbSource 52 | F UseDivbDiffusion 53 | F UseProjection 54 | F UseConstrainB 55 | 56 | #TESTXYZ 57 | 32.0 xTest 58 | 0.0 yTest 59 | 0.0 zTest 60 | 61 | #SAVELOGFILE 62 | T DoSaveLogfile 63 | VAR date GSM StringLogfile 64 | -1 DnSaveLogfile 65 | 60.0 DtSaveLogfile 66 | BXPNT BYPNT BZPNT UXPNT UYPNT UZPNT RHOPNT TPNT 67 | 68 | #SAVEINITIAL 69 | T DoSaveInitial 70 | 71 | #SAVEPLOT 72 | 1 nPlotFiles 73 | 1d MHD idl_ascii StringPlot 74 | -1 DnSavePlot 75 | 60. DtSavePlot 76 | 0.0 DxSavePlot 77 | 78 | #ENDTIME 79 | 1998 iYear 80 | 5 iMonth 81 | 4 iDay 82 | 1 iHour 83 | 45 iMinute 84 | 0 iSecond 85 | 86 | STOP 87 | -1 MaxIteration 88 | 3600.0 tSimulationMax 89 | 90 | -------------------------------------------------------------------------------- /Param/EARTH/TestOutput/earth_cartesian_mag_grid.out.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/EARTH/TestOutput/earth_cartesian_mag_grid.out.gz -------------------------------------------------------------------------------- /Param/EARTH/TestOutput/earth_cartesian_ref.log: -------------------------------------------------------------------------------- 1 | Volume averages, fluxes, etc 2 | it t dt rho mx my mz p bx by bz pmin pmax dst lon1 lat1 lon2 lat2 status 3 | 0 0.00000E+00 0.00000E+000 3.06000E+000 -3.14802E-015 -1.25907E-017 3.53920E-016 3.96994E-003 7.43992E+000 1.24263E+001 3.60607E+000 3.96994E-003 3.96994E-003 -7.62584E-017 0.00000E+000 0.00000E+000 0.00000E+000 0.00000E+000 0.00000E+000 4 | 25 1.00000E+00 1.45817E-002 3.06000E+000 -3.14801E-015 -1.25906E-017 3.53919E-016 3.96995E-003 7.43991E+000 1.24263E+001 3.60606E+000 2.56663E-003 5.71470E-003 -9.40322E-002 -1.96612E+002 -9.84516E+001 -1.98632E+002 -9.99213E+001 2.10059E-002 5 | -------------------------------------------------------------------------------- /Param/EARTH/TestOutput/earth_large_ref.log: -------------------------------------------------------------------------------- 1 | Volume averages, fluxes, etc 2 | it t dt rho mx my mz p bx by bz pmin pmax dst lon1 lat1 lon2 lat2 status 3 | 0 0.00000E+00 0.00000E+000 3.06000E+000 -3.14802E-015 -1.25907E-017 3.53920E-016 3.96994E-003 7.43992E+000 1.24263E+001 3.60607E+000 3.96994E-003 3.96994E-003 0.00000E+000 0.00000E+000 0.00000E+000 0.00000E+000 0.00000E+000 0.00000E+000 4 | 28 1.00000E+00 1.26887E-002 3.06000E+000 -3.14801E-015 -1.25906E-017 3.53919E-016 3.96994E-003 7.43991E+000 1.24263E+001 3.60606E+000 2.53545E-003 6.14392E-003 -4.29772E-001 -1.96695E+002 -9.84900E+001 -1.98619E+002 -9.99204E+001 2.06010E-002 5 | -------------------------------------------------------------------------------- /Param/EUROPA/TestOutput/trj_satellite_2fpe_e4_n00000010.sat.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/EUROPA/TestOutput/trj_satellite_2fpe_e4_n00000010.sat.gz -------------------------------------------------------------------------------- /Param/EUROPA/TestOutput/trj_satellite_3fpe_e4_n00000100.sat.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/Param/EUROPA/TestOutput/trj_satellite_3fpe_e4_n00000100.sat.gz -------------------------------------------------------------------------------- /Param/FIVEMOMENT/PARAM.in.langmuir: -------------------------------------------------------------------------------- 1 | #DESCRIPTION 2 | Param/FIVEMOMENT/PARAM.in.langmuir: 5-moment test with Langmuir waves StringDescription 3 | 4 | #CHECKGRIDSIZE 5 | 50 nI 6 | 50 nJ 7 | 1 nK 8 | 10 MinBlockAll 9 | 10 | #PLANET 11 | NONE NamePlanet 12 | 13 | #IOUNITS 14 | NONE TypeIoUnit 15 | 16 | #NORMALIZATION 17 | NONE TypeNormalization 18 | 19 | #GRID 20 | 1 nRootBlock1 21 | 1 nRootBlock2 22 | 1 nRootBlock3 23 | -0.5 xMin 24 | +0.5 xMax 25 | -0.5 yMin 26 | 0.5 yMax 27 | -1. zMin 28 | 1. zMax 29 | 30 | 31 | #POINTIMPLICIT 32 | T UsePointImplicit 33 | 0.5 BetaPointImplicit 34 | F IsAsymmetric 35 | F DoNormalizeCell 36 | 37 | #LIGHTSPEED 38 | 1.0 cLightDim 39 | 40 | #PLASMA 41 | 1.0 FluidMass 42 | 1/100.0 FluidMass 43 | 1.0 IonCharge 44 | -1.0 IonCharge 45 | 1.0 ElectronTemperatureRatio 46 | 47 | #UNIFORMSTATE 48 | 1.0 StateVar Rho Left 49 | 0.0 StateVar Ux 50 | 0.0 StateVar Uy 51 | 0.0 StateVar Uz 52 | 0.0 StateVar Bx 53 | 0.0 StateVar By 54 | 0.0 StateVar Bz 55 | 0.0 StateVar Ex 56 | 0.0 StateVar Ey 57 | 0.0 StateVar Ez 58 | 0.0 StateVar HypE 59 | 5e-5 StateVar P 60 | 1.0/100.0 StateVar ElRho 61 | 0.0 StateVar ElUx 62 | 0.0 StateVar ElUy 63 | 0.0 StateVar ElUz 64 | 5.e-5 StateVar ElP 65 | 66 | #SHOCKPOSITION 67 | 0.0 ShockPosition 68 | 1.0 ShockSlope 69 | 70 | #WAVE 71 | ElUx NameVar 72 | 1000.0 Width 73 | 1e-2 Amplitude 74 | 0.707106781186 LambdaX 75 | -1.0 LambdaY 76 | -1.0 LambdaZ 77 | 0.0 Phase 78 | 79 | #OUTERBOUNDARY 80 | periodic TypeBc1 81 | periodic TypeBc2 82 | periodic TypeBc3 83 | periodic TypeBc4 84 | 85 | #NONCONSERVATIVE 86 | T UseNonConservative 87 | 88 | #SCHEME 89 | 2 nOrder 90 | Rusanov TypeFlux 91 | mc TypeLimiter 92 | 1.5 LimiterBeta 93 | 94 | #HYPERBOLICDIVE 95 | 0.5 HypEDecay 96 | 97 | FIXEDTIMESTEP 98 | T UseFixedTimeStep 99 | 1.2566/30.0 100 | 101 | TIMESTEPPING 102 | 2 nStage 103 | 0.4 CflExpl 104 | 105 | #SAVERESTART 106 | F DoSaveRestart 107 | 108 | #SAVELOGFILE 109 | T DoSaveLogfile 110 | raw StringLog 111 | 10 DnSaveLogfile 112 | -1. DtSaveLogfile 113 | 114 | #SAVEINITIAL 115 | T DoSaveInitial 116 | 117 | #SAVEPLOT 118 | 2 nPlotFile 119 | 1d mhd idl_ascii StringPlot 120 | -1 DnSavePlot 121 | 0.1 DtSavePlot 122 | 0. DxSavePlot 123 | z=0 MHD idl StringPlot 124 | -1 DnSavePlot 125 | 0.1 DtSavePlot 126 | 0. DxSavePlot 127 | 128 | #STOP 129 | -1 MaxIteration 130 | 1.0 tSimulationMax 131 | 132 | 1.2566 tSimulationMax 133 | -------------------------------------------------------------------------------- /Param/FIVEMOMENT/PARAM.in.shock: -------------------------------------------------------------------------------- 1 | 2 | #DESCRIPTION 3 | Param/FIVEMOMENT/PARAM.in.shock: Five-moment test StringDescription 4 | 5 | #CHECKGRIDSIZE 6 | 100 nI 7 | 1 nJ 8 | 1 nK 9 | 10 MinBlockAll 10 | 11 | #PLANET 12 | NONE NamePlanet 13 | 14 | #IOUNITS 15 | NONE TypeIoUnit 16 | 17 | #NORMALIZATION 18 | NONE TypeNormalization 19 | 20 | #LIGHTSPEED 21 | 1.0 cLightDim 22 | 23 | #PLASMA 24 | 0.1 FluidMass 25 | 0.001 FluidMass 26 | 1.0 IonCharge 27 | -1.0 IonCharge 28 | 1.0 ElectronTemperatureRatio 29 | 30 | #SHOCKTUBE 31 | 1.0 LeftState Rho 32 | 0.0 LeftState Ux 33 | 0.0 LeftState Uy 34 | 0.0 LeftState Uz 35 | 0.75e-2 LeftState Bx 36 | 1.0e-2 LeftState By 37 | 0.0 LeftState Bz 38 | 0.0 LeftState Ex 39 | 0.0 LeftState Ey 40 | 0.0 LeftState Ez 41 | 0.0 LeftState HypE 42 | 0.5e-4 LeftState P 43 | 1.0/100.0 LeftState eRho 44 | 0.0 LeftState eUx 45 | 0.0 LeftState eUy 46 | 0.0 LeftState eUz 47 | 0.5e-4 LeftState eP 48 | 0.125 RightState Rho 49 | 0.0 RightState Ux 50 | 0.0 RightState Uy 51 | 0.0 RightState Uz 52 | 0.75e-2 RightState Bx 53 | -1.0e-2 RightState By 54 | 0.0 RightState Bz 55 | 0.0 RightState Ex 56 | 0.0 RightState Ey 57 | 0.0 RightState Ez 58 | 0.0 RightState HypE 59 | 0.05e-4 RightState P 60 | 0.125/100.0 RightState eRho 61 | 0.0 RightState eUx 62 | 0.0 RightState eUy 63 | 0.0 RightState eUz 64 | 0.05e-4 RightState eP 65 | 66 | #GRID 67 | 8 nRootBlock1 68 | 1 nRootBlock2 69 | 1 nRootBlock3 70 | -0.5 xMin 71 | +0.5 xMax 72 | -1. yMin 73 | 1. yMax 74 | -1. zMin 75 | 1. zMax 76 | 77 | #OUTERBOUNDARY 78 | float TypeBc1 79 | float TypeBc2 80 | 81 | #SCHEME 82 | 2 nOrder 83 | Rusanov TypeFlux 84 | mc TypeLimiter 85 | 1.5 LimiterBeta 86 | 87 | TIMESTEPPING 88 | 3 nStage 89 | 0.8 CflExpl 90 | 91 | #SAVERESTART 92 | F DoSaveRestart 93 | 94 | #SAVELOGFILE 95 | T DoSaveLogfile 96 | raw StringLog 97 | 10 DnSaveLogfile 98 | -1. DtSaveLogfile 99 | 100 | #SAVEINITIAL 101 | T DoSaveInitial 102 | 103 | #SAVEPLOT 104 | 1 nPlotFile 105 | 1d mhd idl_ascii StringPlot 106 | -1 DnSavePlot 107 | 0.1 DtSavePlot 108 | 0. DxSavePlot 109 | 110 | #STOP 111 | -1 MaxIteration 112 | 2.0 tSimulationMax 113 | -------------------------------------------------------------------------------- /Param/FLUXEMERGENCE/Grid.1D: -------------------------------------------------------------------------------- 1 | #GRID 2 | 2 nRootBlock1 3 | 2 nRootBlock2 4 | 56 nRootBlock3 5 | -5. xMin 6 | 5. xMax 7 | -5. yMin 8 | 5. yMax 9 | -140. zMin 10 | 140. zMax 11 | 12 | #GRIDBLOCKALL 13 | 500 MaxBlock 14 | 15 | #GRIDRESOLUTION 16 | 5./64 Resolution 17 | initial StringShape 18 | 19 | #END 20 | -------------------------------------------------------------------------------- /Param/FLUXEMERGENCE/Grid.3D: -------------------------------------------------------------------------------- 1 | #GRID 2 | 1 nRootBlock1 3 | 1 nRootBlock2 4 | 7 nRootBlock3 5 | -20 xMin 6 | 20 xMax 7 | -20 yMin 8 | 20 yMax 9 | -140 zMin 10 | 140 zMax 11 | 12 | #GRIDBLOCKALL 13 | 500 MaxBlock 14 | 15 | #GRIDRESOLUTION 16 | 35./256 Resolution 17 | initial StringShape 18 | 19 | #GRIDLEVEL 20 | 2 nLevel 21 | all StringShape 22 | 23 | GRIDLEVEL 24 | 3 nLevel 25 | box NameArea 26 | -20 xMinBox 27 | -20 yMinBox 28 | -60 zMinBox 29 | 20 xMaxBox 30 | 20 yMaxBox 31 | 20 zMaxBox 32 | 33 | GRIDLEVEL 34 | 4 nLevel 35 | box NameArea 36 | -20 xMinBox 37 | -20 yMinBox 38 | -10 zMinBox 39 | 20 xMaxBox 40 | 20 yMaxBox 41 | 0 zMaxBox 42 | 43 | #END 44 | -------------------------------------------------------------------------------- /Param/FLUXEMERGENCE/TestOutput/log_fluxemergence.log: -------------------------------------------------------------------------------- 1 | Volume averages, fluxes, etc 2 | it t rho Ux Uy Uz E Pmin Pmax 3 | 11 4.47557E+00 4.61419E-004 -8.03507E-010 -2.26747E-010 -5.59218E-001 9.82691E+009 1.87158E-004 5.40250E+010 4 | 12 4.92173E+00 4.61422E-004 -1.59988E-009 -4.53107E-010 -6.14793E-001 9.82700E+009 1.87158E-004 5.40259E+010 5 | 13 5.36765E+00 4.61424E-004 -2.38998E-009 -6.79019E-010 -6.70316E-001 9.82707E+009 1.87158E-004 5.40266E+010 6 | 14 5.81331E+00 4.61426E-004 -3.17439E-009 -9.04459E-010 -7.25782E-001 9.82715E+009 1.87158E-004 5.40269E+010 7 | 15 6.25876E+00 4.61428E-004 -3.95371E-009 -1.12943E-009 -7.81196E-001 9.82723E+009 1.87158E-004 5.40269E+010 8 | 16 6.70399E+00 4.61430E-004 -4.72837E-009 -1.35384E-009 -8.36555E-001 9.82730E+009 1.87158E-004 5.40264E+010 9 | 17 7.14900E+00 4.61432E-004 -5.49869E-009 -1.57764E-009 -8.91856E-001 9.82737E+009 1.87158E-004 5.40256E+010 10 | 18 7.59380E+00 4.61435E-004 -6.26494E-009 -1.80070E-009 -9.47100E-001 9.82744E+009 1.87158E-004 5.40243E+010 11 | 19 8.03840E+00 4.61437E-004 -7.02731E-009 -2.02300E-009 -1.00229E+000 9.82751E+009 1.87158E-004 5.40226E+010 12 | 20 8.48278E+00 4.61439E-004 -7.78594E-009 -2.24450E-009 -1.05742E+000 9.82757E+009 1.87158E-004 5.40205E+010 13 | -------------------------------------------------------------------------------- /Param/GANYMEDE/Grid: -------------------------------------------------------------------------------- 1 | ---------------grid structure----------------- 2 | #GRIDGEOMETRY 3 | spherical_lnr TypeGeometry 4 | 5 | #GRID 6 | 10 nRootBlock1 7 | 8 nRootBlock2 8 | 8 nRootBlock3 9 | -100 xMin 10 | 100 xMax 11 | -100 yMin 12 | 100 yMax 13 | -100 zMin 14 | 100 zMax 15 | 16 | #GRIDBLOCKALL 17 | 700 MaxBlock 18 | 19 | #LIMITRADIUS 20 | 0.5 rMin 21 | 174.0 rMax > (sqrt(100^2+100^2+100^2)) 22 | 23 | GRIDLEVEL 24 | 2 nLevelArea 25 | initial NameArea 26 | 27 | GRIDLEVEL 28 | 2 nLevelArea 29 | box_gen TypeRegion 30 | 1.1 rmin 31 | 0.0 LonMin 32 | -65.0 LatMin 33 | 4.0 rmax 34 | 360.0 LonMax 35 | 65.0 LatMax 36 | 37 | GRIDLEVEL 38 | 1 nLevelArea 39 | box_gen TypeRegion 40 | 5.0 rmin 41 | -70.0 LonMin 42 | 35.0 LatMin 43 | 15.0 rmax 44 | 70.0 LonMax 45 | 65.0 Latmax 46 | 47 | GRIDLEVEL 48 | 1 nLevelArea 49 | box_gen TypeRegion 50 | 5.0 rmin 51 | -70.0 LonMin 52 | -35.0 LatMin 53 | 15.0 rmax 54 | 70.0 LonMax 55 | -65.0 LatMax 56 | 57 | -------------------------------------------------------------------------------- /Param/GANYMEDE/TestOutput/REFlog_n000001.log: -------------------------------------------------------------------------------- 1 | Volume averages, fluxes, etc 2 | it t dt Rho Mx My Mz Bx By Bz Hyp Pe E Pmin Pmax 3 | 1 0.00000E+00 0.00000E+000 5.60002E+001 1.31132E-014 1.17543E-024 -1.51660E-023 -1.00000E+001 -6.00000E+000 -8.60000E+001 2.45992E-010 2.00000E-001 9.31481E-009 1.15000E-001 3.61415E+000 4 | 2 0.00000E+00 0.00000E+000 5.60002E+001 1.31132E-014 2.24292E-024 -3.10428E-023 -1.00000E+001 -6.00000E+000 -8.60000E+001 9.18896E-010 2.00000E-001 9.31481E-009 1.15000E-001 3.62507E+000 5 | 3 0.00000E+00 0.00000E+000 5.60002E+001 1.31132E-014 3.25511E-024 -4.78043E-023 -1.00000E+001 -6.00000E+000 -8.60000E+001 1.88562E-009 2.00000E-001 9.31481E-009 1.15000E-001 3.62665E+000 6 | 4 0.00000E+00 0.00000E+000 5.60002E+001 1.31132E-014 4.24307E-024 -6.55238E-023 -1.00000E+001 -6.00000E+000 -8.60000E+001 3.02697E-009 2.00000E-001 9.31481E-009 1.15000E-001 3.61625E+000 7 | 5 0.00000E+00 0.00000E+000 5.60002E+001 1.31132E-014 5.22451E-024 -8.42468E-023 -1.00000E+001 -6.00000E+000 -8.60000E+001 4.29594E-009 2.00000E-001 9.31481E-009 1.15000E-001 3.62618E+000 8 | -------------------------------------------------------------------------------- /Param/HELIOSPHERE/Arcade: -------------------------------------------------------------------------------- 1 | !- Arcade Instability Model based on Chip's 2D model 2 | ! October 12 3 | 4 | #PROBLEMTYPE 5 | 25 problem_type: 1=MHD Uniform Flow 6 | 2=Brio-Wu Initial Value 7 | 3=Solar Wind and Inner Heliosphere 8 | 4=Magnetotail Flow 9 | 5=Mass-Loaded Comet 10 | 11=Earth Magnetosphere 11 | 12=Saturn Magnetosphere 12 | 13=Jupiter Magnetosphere 13 | 14=Venus Ionosphere 14 | 21=Conducting Cylinder (2-D) 15 | 22=Conducting Cylinder (3-D) 16 | 25=Coronal Arcade 17 | 18 | #ARCADE 19 | 1.0E6 TArcDim [K] 20 | 1.0E-12 RhoArcDim [kg/m^3] 21 | 0.71814 BArcDim [Gauss] 22 | 0.0 ByArcDim [Gauss] 23 | 5.0E3 UzArcDim [5.0E3 m/s] 24 | 0.5 phi0Arc 25 | 1.3 muArc 26 | 3 expArc 27 | 0.5 widthArc 28 | 29 | #GAMMA 30 | 1.6666666667 g 31 | 32 | #GRID 33 | 2 proc_dims(1) 34 | 1 proc_dims(2) 35 | 6 proc_dims(3), 36 | -3.00 x1 2.58 37 | 3.00 x2 2.58 38 | 0.00 y1 -8.00 39 | 0.05 y2 8.00 40 | -3.00 z1 -3.45 41 | 13.00 z2 14.20 42 | 43 | #OUTERBOUNDARY 44 | reflect outerBCtype(east_) ('inflow','outflow','float') 45 | reflect outerBCtype(west_) 46 | periodic outerBCtype(south_) (periodic) 47 | periodic outerBCtype(north_) (periodic) 48 | arcadebot outerBCtype(bot_) (linetied) 49 | float outerBCtype(top_) (set_P) 50 | 51 | #AMRINIT 52 | all InitialRefineType "none", "all" and others 53 | 2 initial_refinement_levels 54 | 55 | #AMRCRITERIA 56 | 2 nRefineCrit (number of refinement criteria: 1,2 or 3) 57 | gradlogP RefineCrit(1) 58 | curlB RefineCrit(2) 59 | 60 | #AMRINITPHYSICS 61 | 12 nRefineLevelIC 62 | 63 | #AMR 64 | 100 dn_refine (frequency in terms of total steps n_step) 65 | .true. automatic_refinement 66 | 1.0 percentCoarsen 67 | 10.0 percentRefine 68 | 7200 maxTotalBlocks 7200 69 | 70 | #AMRLEVELS 71 | 0 min_block_level 72 | 99 max_block_level 73 | .false. fix_body_level 74 | 75 | #END 76 | -------------------------------------------------------------------------------- /Param/HELIOSPHERE/CME_fluxrope: -------------------------------------------------------------------------------- 1 | !- Arcade Instability Model based on BC Low's 3D model 2 | ! January 18 2001 3 | 4 | #PROBLEMTYPE 5 | 26 problem_type: 1=MHD Uniform Flow 6 | 2=Brio-Wu Initial Value 7 | 3=Solar Wind and Inner Heliosphere 8 | 4=Magnetotail Flow 9 | 5=Mass-Loaded Comet 10 | 11=Earth Magnetosphere 11 | 12=Saturn Magnetosphere 12 | 13=Jupiter Magnetosphere 13 | 14=Venus Ionosphere 14 | 21=Conducting Cylinder (2-D) 15 | 22=Conducting Cylinder (3-D) 16 | 25=Coronal Arcade Eruption (2-D) 17 | 26=CME driven by fluxrope (3-D) 18 | 19 | #BODY 20 | T 21 | 1. Rbody 22 | 4. Rcurrents 23 | 1.50E+14 NDsun [n.cm-3] 24 | 2.85E+06 Tsun [K] 25 | 26 | 27 | #CME 28 | Low 29 | 0.7 cme_a [scaled] 0.7 30 | 1.2 cme_r1 [scaled] 1.2 31 | 1.0 cme_r0 [scaled] 1.0 32 | 0.23 cme_a1 [Gauss] 0.23 33 | 0.0 alpha [scaled] 0.0 34 | 2.5E-12 cme_rho1[kg/m^3] 2.5E-12 35 | 2.0E-13 cme_rho2[kg/m^3] 2.0E-13 36 | 1.0 B1_dim [Gauss] 37 | 1.5E5 v_erupt [m/s] 38 | 39 | #DIPOLE 40 | -3.0E00 DipoleStrength (equatorial field strength in Gauss) 41 | .false. SetDipoleTilt 42 | -15.0 ThetaTilt 43 | 0.0 dt_UpdateB0 44 | 45 | #COROTATION 46 | .false. UseCorotation 47 | .false. AlignDipoleCorotation 48 | .false. SetCorotationTilt 49 | 0.0 CorotationTilt 50 | 51 | #GAMMA 52 | 1.3333333 g 53 | 54 | #GRID 55 | 1 proc_dims(1) 56 | 1 proc_dims(2) 57 | 1 proc_dims(3), 58 | -9.0 x1 59 | 9.0 x2 60 | -9.0 y1 61 | 9.0 y2 62 | -6.0 z1 63 | 12.0 z2 64 | 65 | #OUTERBOUNDARY 66 | float outerBCtype(east_) ('inflow','outflow','float') 67 | float outerBCtype(west_) 68 | float outerBCtype(south_) 69 | float outerBCtype(north_) 70 | float outerBCtype(bot_) 71 | float outerBCtype(top_) 72 | 73 | #BODY 74 | true body1 !Others read only if Body1=.true. 75 | 0.4 Rbody 76 | 0.0 Rcurrents 77 | 0.0 Magneto_rho_dim (/ccm) density fixed BC for rho_BLK 78 | 0.0 Magneto_T_dim (K) temperature for fixed BC for P_BLK 79 | 80 | #INNERBOUNDARY 81 | linetied innerBCtype 82 | 83 | #END 84 | -------------------------------------------------------------------------------- /Param/HELIOSPHERE/EARTH_TRAJ.in: -------------------------------------------------------------------------------- 1 | !A sample of a description file for the Earth trajectory 2 | #COOR 3 | HGR ! Coordinate system; only HAE, HGI, or HGR allowed 4 | 5 | #DATE 6 | 1998 ! Year 7 | 5 ! Month 8 | 1 ! Day 9 | 0 ! Hour 10 | 0 ! Minute 11 | 0 ! Second 12 | 0 ! Millisecond 13 | 14 | #TIMELOOP 15 | -86400.000 86400.000 1800.000 ! Starttime, Endtime, Timestep 16 | 17 | #END 18 | -------------------------------------------------------------------------------- /Param/HELIOSPHERE/EARTH_TRAJ.in.CR: -------------------------------------------------------------------------------- 1 | !A description file for the Earth trajectory for the whole Carrington rotation 2 | #COOR 3 | HGR ! Coordinate system; only HAE, HGI, or HGR allowed 4 | 5 | !!NOTE, this should be a START TIME for the CARRINGTON ROTATION 6 | 7 | #INCLUDE 8 | TIME.in 9 | 10 | 11 | #TIMELOOP 12 | -0.1 2419200.100 1800.000 ! Starttime, Endtime, Timestep 13 | 14 | #END 15 | -------------------------------------------------------------------------------- /Param/HELIOSPHERE/Heliosphere: -------------------------------------------------------------------------------- 1 | !- Arcade Instability Model based on Chip's 2D model 2 | ! October 12 3 | 4 | #PROBLEMTYPE 5 | 3 problem_type: 1=MHD Uniform Flow 6 | 2=Brio-Wu Initial Value 7 | 3=Solar Wind and Inner Heliosphere 8 | 4=Magnetotail Flow 9 | 5=Mass-Loaded Comet 10 | 11=Earth Magnetosphere 11 | 12=Saturn Magnetosphere 12 | 13=Jupiter Magnetosphere 13 | 14=Venus Ionosphere 14 | 21=Conducting Cylinder (2-D) 15 | 22=Conducting Cylinder (3-D) 16 | 25=Coronal Arcade 17 | 18 | #HELIOSPHERE 19 | 2.85E+06 Tsun [K] 1.0E6 20 | 1.50E+14 NDsun [n/m-3] 21 | 25.00 Qsun 22 | 1.75 Theat 23 | 1.00 Rheat 24 | 4.50 SIGMAheat 25 | 26 | #INERTIAL 27 | .true. UseInertial, if false then corotating frame 28 | 29 | #DIPOLE 30 | -3.0E00 DipoleStrength (equatorial field strength in Gauss) 31 | .true. SetDipoleTilt 32 | 0.0 ThetaTilt 33 | 0.0 dt_UpdateB0 34 | 35 | #COROTATION 36 | .true. UseCorotation 37 | .false. AlignDipoleCorotation 38 | .false. SetCorotationTilt 39 | 0.0 CorotationTilt 40 | 41 | #GAMMA 42 | 1.666666666666 g 43 | 44 | #SOLARWIND 45 | 5.0 SW_rho_dim [n/cc] 46 | 181712.175 SW_T_dim [K] 47 | -400.0 SW_Ux_dim [km/s] 48 | 0.0 SW_Uy_dim [km/s] 49 | 0.0 SW_Uz_dim [km/s] 50 | 0.0 SW_Bx_dim [nT] 51 | 0.0 SW_By_dim [nT] 52 | 5.0 SW_Bz_dim [nT] 53 | 54 | 55 | #GRID 56 | 1 proc_dims(1) 57 | 1 proc_dims(2) 58 | 1 proc_dims(3), 59 | -32.0 x1 60 | 224.0 x2 61 | -128.0 y1 62 | 128.0 y2 63 | -128.0 z1 64 | 128.0 z2 65 | 66 | #OUTERBOUNDARY 67 | float outerBCtype(east_) ('inflow','outflow','float') 68 | float outerBCtype(west_) 69 | float outerBCtype(south_) 70 | float outerBCtype(north_) 71 | float outerBCtype(bot_) 72 | float outerBCtype(top_) 73 | 74 | #INNERBOUNDARY 75 | inner_helio innerBCtype 76 | 77 | #END 78 | -------------------------------------------------------------------------------- /Param/HELIOSPHERE/SATELLITE.in: -------------------------------------------------------------------------------- 1 | #SATELLITE 2 | 3 nSatelliteFile 3 | MHD file StringSatellite 4 | -1 DnOutput ! controls frequency with which the whole 5 | 1800.0 DtOutput ! trajectory solution is provided 6 | IH/earthtraj.dat 7 | MHD file StringSatellite 8 | -1 DnOutput ! controls frequency with which the whole 9 | 1800.0 DtOutput ! trajectory solution is provided 10 | IH/Param/stereoa.dat 11 | MHD file StringSatellite 12 | -1 DnOutput ! controls frequency with which the whole 13 | 1800.0 DtOutput ! trajectory solution is provided 14 | IH/Param/stereob.dat 15 | 16 | -------------------------------------------------------------------------------- /Param/JUPITER/Cassini_nov18_IMF: -------------------------------------------------------------------------------- 1 | ! solarwind input data for Nov. 18 at jupiter 2 | 3 | #SOLARWIND 4 | 0.05 SW_rho_dim [n/cc] 5 | 25000. SW_T_dim [K] 6 | -405.0 SW_Ux_dim [km/s] 7 | 0.0 SW_Uy_dim [km/s] 8 | 0.0 SW_Uz_dim [km/s] 9 | 0.1 SW_Bx_dim [nT] 10 | -0.3 SW_By_dim [nT] 11 | 0.4 SW_Bz_dim [nT] 12 | 13 | #END 14 | -------------------------------------------------------------------------------- /Param/JUPITER/Ideal_Corotation: -------------------------------------------------------------------------------- 1 | 2 | #AXES 3 | GSM CoordinateSystemString 4 | ideal ideal AxesType 5 | 6 | #END -------------------------------------------------------------------------------- /Param/JUPITER/Mag_SSSS_0.25Rj: -------------------------------------------------------------------------------- 1 | 2 | #INCLUDE 3 | Param/JUPITER/Mag_SSSS_0.5Rj 4 | 5 | #AMR 6 | 1501 7 | F automatic_refinement 8 | 9 | #STOP 10 | 1750 11 | 12 | #RUN ###################################################### 13 | 14 | #AMR 15 | -1 16 | 17 | #STOP 18 | 2000 19 | 20 | #RUN ###################################################### 21 | 22 | #END -------------------------------------------------------------------------------- /Param/JUPITER/Mag_SSSS_0.2Rj: -------------------------------------------------------------------------------- 1 | 2 | #INCLUDE 3 | Param/JUPITER/Mag_SSSS_0.4Rj 4 | 5 | #AMR 6 | 1501 7 | F automatic_refinement 8 | 9 | #STOP 10 | 1750 11 | 12 | #RUN ###################################################### 13 | 14 | #AMR 15 | -1 16 | 17 | #STOP 18 | 2000 19 | 20 | #RUN ###################################################### 21 | 22 | #END -------------------------------------------------------------------------------- /Param/JUPITER/Mag_SSSS_0.4Rj: -------------------------------------------------------------------------------- 1 | 2 | #INCLUDE 3 | Param/JUPITER/Mag_SSSS_0.8Rj 4 | 5 | #AMR 6 | 1001 7 | F automatic_refinement 8 | 9 | #STOP 10 | 1250 11 | 12 | #RUN ###################################################### 13 | 14 | #AMR 15 | -1 16 | 17 | #STOP 18 | 1500 19 | 20 | #RUN ###################################################### 21 | 22 | #END -------------------------------------------------------------------------------- /Param/JUPITER/Mag_SSSS_0.5Rj: -------------------------------------------------------------------------------- 1 | 2 | #INCLUDE 3 | Param/JUPITER/Mag_SSSS_1.0Rj 4 | 5 | #AMR 6 | 1001 7 | F automatic_refinement 8 | 9 | #STOP 10 | 1250 11 | 12 | #RUN ###################################################### 13 | 14 | #AMR 15 | -1 16 | 17 | #STOP 18 | 1500 19 | 20 | #RUN ###################################################### 21 | 22 | #END -------------------------------------------------------------------------------- /Param/JUPITER/Mag_SSSS_0.8Rj: -------------------------------------------------------------------------------- 1 | #MASSLOADING 2 | F UseMassLoading 3 | F AccelerateMassLoading 4 | 5 | #TIMESTEPPING 6 | F time_accurate 7 | 1 nSTAGE 8 | 0.60 cfl 9 | 10 | #AMRINIT 11 | magnetojupiter InitialRefineType 12 | 7 initial_refine_levels 13 | 14 | #SCHEME 15 | 1 nORDER 16 | 2 fluxfcn_type 1=Roe, 2=Rusanov, 3=Linde 17 | minmod 18 | 19 | #BORIS 20 | T boris_correction 21 | 1.0 boris_cLIGHT_factor 22 | 23 | #STOP 24 | 200 n_step 25 | 26 | #RUN ###################################################### 27 | 28 | #BORIS 29 | T boris_correction 30 | 0.1 boris_cLIGHT_factor 31 | 32 | #AMR 33 | 300 nrefine1 34 | F automatic_refinement 35 | 36 | #STOP 37 | 700 38 | 39 | #RUN ###################################################### 40 | 41 | #AMR 42 | -1 nrefine1 43 | 44 | #SCHEME 45 | 2 nORDER 46 | 2 fluxfcn_type 1=Roe, 2=Rusanov, 3=Linde 47 | minmod limiter_type ('minmod' or 'mc') 48 | 49 | #STOP 50 | 1000 nITER 51 | 52 | #RUN ###################################################### 53 | 54 | #END 55 | -------------------------------------------------------------------------------- /Param/JUPITER/Mag_SSSS_1.0Rj: -------------------------------------------------------------------------------- 1 | #MASSLOADING 2 | F UseMassLoading 3 | F AccelerateMassLoading 4 | 5 | #TIMESTEPPING 6 | F time_accurate 7 | 1 nSTAGE 8 | 0.60 cfl 9 | 10 | #AMRINIT 11 | magnetojupiter InitialRefineType 12 | 7 initial_refine_levels 13 | 14 | #SCHEME 15 | 1 nORDER 16 | 2 fluxfcn_type 1=Roe, 2=Rusanov, 3=Linde 17 | minmod 18 | 19 | #BORIS 20 | T boris_correction 21 | 1.0 boris_cLIGHT_factor 22 | 23 | #STOP 24 | 200 n_step 25 | 26 | #RUN ###################################################### 27 | 28 | #BORIS 29 | T boris_correction 30 | 0.1 boris_cLIGHT_factor 31 | 32 | #AMR 33 | 300 nrefine1 34 | F automatic_refinement 35 | 36 | #STOP 37 | 700 38 | 39 | #RUN ###################################################### 40 | 41 | #AMR 42 | -1 nrefine1 43 | 44 | #SCHEME 45 | 2 nORDER 46 | 2 fluxfcn_type 1=Roe, 2=Rusanov, 3=Linde 47 | minmod limiter_type ('minmod' or 'mc') 48 | 49 | #STOP 50 | 1000 nITER 51 | 52 | #RUN ###################################################### 53 | 54 | #END 55 | -------------------------------------------------------------------------------- /Param/JUPITER/Magnetosphere: -------------------------------------------------------------------------------- 1 | 2 | #PROBLEMTYPE 3 | 13 problem_type: 13=Jupiter Magnetosphere 4 | 5 | #GRID 6 | 2 proc_dims(1) 7 | 1 proc_dims(2) 8 | 1 proc_dims(3), 9 | -1792. x1 10 | 256. x2 11 | -512. y1 12 | 512. y2 13 | -512. z1 14 | 512. z2 15 | 16 | #MAGNETOSPHERE 17 | T body1 18 | 3.0 Rbody 19 | 4.0 Rcurrents 20 | 1.0 Magneto_rho_dim (/ccm) density for fixed BC for rho_BLK 21 | 35000.0 Magneto_T_dim (K) temperature for fixed BC for P_BLK 22 | 23 | #OUTERBOUNDARY 24 | outflow outerBCtype(east_) ('inflow','outflow','float') 25 | inflow outerBCtype(west_) 26 | float outerBCtype(south_) 27 | float outerBCtype(north_) 28 | float outerBCtype(bot_) 29 | float outerBCtype(top_) 30 | 31 | #INNERBOUNDARY 32 | ionosphere innerBCtype 33 | 34 | #END -------------------------------------------------------------------------------- /Param/JUPITER/North_IMF: -------------------------------------------------------------------------------- 1 | 2 | #SOLARWIND 3 | 0.4 SW_rho_dim [n/cc] 4 | 181712.175 SW_T_dim [K] 5 | -400.0 SW_Ux_dim [km/s] 6 | 0.0 SW_Uy_dim [km/s] 7 | 0.0 SW_Uz_dim [km/s] 8 | 0.0 SW_Bx_dim [nT] 9 | 0.0 SW_By_dim [nT] 10 | 1.0 SW_Bz_dim [nT] 11 | 12 | #END 13 | -------------------------------------------------------------------------------- /Param/JUPITER/Parker_IMF: -------------------------------------------------------------------------------- 1 | 2 | #SOLARWIND 3 | 0.4 SW_rho_dim [n/cc] 4 | 181712.175 SW_T_dim [K] 5 | -400.0 SW_Ux_dim [km/s] 6 | 0.0 SW_Uy_dim [km/s] 7 | 0.0 SW_Uz_dim [km/s] 8 | 0.173648 SW_Bx_dim [nT] 9 | 0.984808 SW_By_dim [nT] 10 | 0.0 SW_Bz_dim [nT] 11 | 12 | #END 13 | -------------------------------------------------------------------------------- /Param/JUPITER/Real_Corotation: -------------------------------------------------------------------------------- 1 | 2 | #AXES 3 | GSM CoordinateSystemString 4 | planet planet AxesType 5 | 6 | #END 7 | 8 | -------------------------------------------------------------------------------- /Param/JUPITER/South_IMF: -------------------------------------------------------------------------------- 1 | 2 | #SOLARWIND 3 | 0.4 SW_rho_dim [n/cc] 4 | 181712.175 SW_T_dim [K] 5 | -400.0 SW_Ux_dim [km/s] 6 | 0.0 SW_Uy_dim [km/s] 7 | 0.0 SW_Uz_dim [km/s] 8 | 0.0 SW_Bx_dim [nT] 9 | 0.0 SW_By_dim [nT] 10 | -1.0 SW_Bz_dim [nT] 11 | 12 | #END 13 | -------------------------------------------------------------------------------- /Param/MAGNETOMETER/magin.dat: -------------------------------------------------------------------------------- 1 | #COORD 2 | MAG 3 | 4 | #START 5 | st1 70.0 80.0 6 | st2 60.0 140.0 7 | st3 -55.0 200.0 8 | st4 20.0 310.0 9 | -------------------------------------------------------------------------------- /Param/MARS/Grid: -------------------------------------------------------------------------------- 1 | #GRIDGEOMETRY 2 | spherical_lnr TypeGeometry 3 | 4 | #GRID 5 | 4 nRootBlock1 6 | 2 nRootBlock2 7 | 1 nRootBlock3 8 | -24. xMin 9 | 8. xMax 10 | -16. yMin 11 | 16. yMax 12 | -16. zMin 13 | 16. zMax 14 | 15 | #GRIDBLOCKALL 16 | 1000 MaxBlock 17 | 18 | #LIMITRADIUS 19 | 1.029884232197 rMin 20 | 32.9849 rMax = sqrt(24^2+16^2+16^2) 21 | 22 | #GRIDLEVEL 23 | 2 nLevel 24 | initial StringShape 25 | 26 | #GRIDLEVEL 27 | 1 nLevel 28 | all StringShape 29 | 30 | #GRIDLEVEL 31 | 10 nLevel 32 | sphere0 StringShape 33 | 1.544 Radius = 1.5*rBody 34 | -------------------------------------------------------------------------------- /Param/MARS/PARAM.in.restartread: -------------------------------------------------------------------------------- 1 | #INCLUDE 2 | GM/restartIN/restart.H 3 | 4 | #DESCRIPTION 5 | Param/MARS/PARAM.in.restartread: Mars test StringDescription 6 | 7 | #ROTATION 8 | F UseRotation 9 | 10 | #GRAVITY 11 | T UseGravity 12 | 0 iDirGravity 13 | 14 | #TESTXYZ 15 | 1.2 xTest 16 | 0.1 yTest 17 | 0.1 zTest 18 | 19 | #TIMEACCURATE 20 | F IsTimeAccurate 21 | 22 | #TIMESTEPPING 23 | 2 nStage 24 | 0.8 CflExpl 25 | 26 | #SCHEME 27 | 2 nOrder 28 | Linde TypeFlux 29 | minmod TypeLimiter 30 | 31 | #SAVELOGFILE 32 | T DoSaveLogfile 33 | VAR step time StringLog 34 | 1 DnSaveLogfile 35 | -1. DtSaveLogfile 36 | rho p rhopnt ppnt bxpnt pmin pmax Hpflx Opflx O2pflx CO2pflx NameLogVars 37 | 3.0 StringLogRadii 38 | 39 | #INNERBOUNDARY 40 | user TypeBcBody 41 | 42 | #BOXBOUNDARY 43 | outflow TypeBc1 44 | inflow TypeBc2 45 | fixed TypeBc3 46 | fixed TypeBc4 47 | fixed TypeBc5 48 | fixed TypeBc6 49 | 50 | #USERSWITCH 51 | +init +ic +b0 +source StringSwitch 52 | 53 | #USERINPUTBEGIN ----------------------------------- 54 | 55 | #UseHotO 56 | T 57 | 58 | #UseTempCont 59 | F 60 | 61 | #UseMarsB0 62 | T UseMarsB0 63 | 60 NNm 64 | -180 Rot 65 | 0.0 Thetilt 66 | marsmgsp.txt NameFileB0 67 | 68 | -99.5 69 | 25.26 70 | 71 | #SOLARCON 72 | solarmax 73 | 74 | #USERINPUTEND ------------------------------------- 75 | 76 | #MULTISPECIES 77 | T DoReplaceDensity 78 | 5.0 SpeciesPercentCheck 79 | 80 | #SAVERESTART 81 | F DoSaveRestart 82 | 83 | #SAVEPLOT 84 | 6 nPlotFile 85 | x=0 MHD tec StringPlot 86 | 5000 DnSavePlot 87 | -1.0 DtSavePlot 88 | y=0 MHD tec StringPlot 89 | 5000 DnSavePlot 90 | -1.0 DtSavePlot 91 | z=0 MHD tec StringPlot 92 | 5000 DnSavePlot 93 | -1.0 DtSavePlot 94 | x=0 MHD idl StringPlot 95 | 2000 DnSavePlot 96 | -1.0 DtSavePlot 97 | -1.0 DxSavePlot 98 | y=0 MHD idl StringPlot 99 | 2000 DnSavePlot 100 | -1.0 DtSavePlot 101 | -1.0 DxSavePlot 102 | z=0 MHD idl StringPlot 103 | 2000 DnSavePlot 104 | -1.0 DtSavePlot 105 | -1.0 DxSavePlot 106 | 107 | cut MHD idl StringPlot 108 | 1000 DnSavePlot 109 | -1.0 DtSavePlot 110 | -10.0 xMin 111 | 20.0 xMax 112 | 0.0001 yMin 113 | 0.00011 yMax 114 | 0.0001 zMin 115 | 0.00011 zmax 116 | 0. DxSavePlot 117 | 118 | #STOP 119 | 25 MaxIteration 120 | -1.0 tSimulationMax 121 | -------------------------------------------------------------------------------- /Param/MARS/TestOutput/log_ta.log: -------------------------------------------------------------------------------- 1 | test point (X,Y,Z): 1.20000E+00 1.00000E-01 1.00000E-01 2 | it t X Y Z dt rho p rhopnt ppnt bxpnt Hpflx_R=6.0 Hpflx_R=2.0 Hpflx_R=1.117786 Opflx_R=6.0 Opflx_R=2.0 Opflx_R=1.117786 O2pflx_R=6.0 O2pflx_R=2.0 O2pflx_R=1.117786 CO2pflx_R=6.0 CO2pflx_R=2.0 CO2pflx_R=1.117786 3 | 163 2.00000E+00 1.20142E+00 1.18329E-01 1.18902E-01 2.60632E-003 1.71516E+001 1.71534E-001 8.83363E+000 1.03328E-001 -3.73049E+000 -7.48265E+025 -2.06051E+022 -2.54419E+023 -9.85751E+014 -4.54022E+018 -8.93121E+021 -7.54633E+014 -3.51759E+013 -1.10430E+022 -6.91604E+014 -4.67884E+010 -9.80843E+019 4 | -------------------------------------------------------------------------------- /Param/MARS/imf.dat: -------------------------------------------------------------------------------- 1 | Units: nT, km/s, amu/cc, K 2 | 3 | #COORD 4 | GSE 5 | 6 | year mo dy hr mn sc msc bx by bz ux uy uz rho T 7 | #START 8 | 2017 09 12 18 0 0 0 0.0 -2.81 1.85 -394.0 0.0 0.0 4.99 1.5E6 9 | 2017 09 12 18 0 0 2 0.0 -2.81 -1.85 -394.0 0.0 0.0 4.99 1.5E6 10 | 2017 09 15 0 0 0 0 0.0 -2.81 1.85 -394.0 0.0 0.0 4.99 1.5E6 11 | -------------------------------------------------------------------------------- /Param/MARSFLUIDS/PARAM.in.restartread: -------------------------------------------------------------------------------- 1 | #UPDATECHECK 2 | T UseUpdateCheck 3 | 40. RhoMinPercent 4 | 400. RhoMaxPercent 5 | 40. pMinPercent 6 | 400. pMaxPercent 7 | 8 | #INCLUDE 9 | GM/restartIN/restart.H 10 | 11 | #DESCRIPTION 12 | Param/MARSFLUIDS/PARAM.in.restartread: Mars test read restart StringDescription 13 | 14 | #ROTATION 15 | F UseRotation 16 | 17 | #GRAVITY 18 | T UseGravity 19 | 0 iDirGravity 20 | 21 | #POINTIMPLICIT 22 | T UsePointImplicit 23 | 1.0 BetaPointImplicit 24 | F IsAsymmetric 25 | T DoNormalizeCell 26 | 27 | #TIMEACCURATE 28 | F IsTimeAccurate 29 | 30 | #SCHEME 31 | 2 nOrder 32 | Linde TypeFlux 33 | minmod TypeLimiter 34 | 35 | #TIMESTEPPING 36 | 2 nStage 37 | 0.2 CflExpl 38 | 39 | #NONCONSERVATIVE 40 | T UseNonConservative 41 | 42 | #SAVELOGFILE 43 | T DoSaveLogfile 44 | RAW step time StringLog 45 | 1 DnSaveLogfile 46 | -1. DtSaveLogfile 47 | 48 | #INNERBOUNDARY 49 | user TypeBcBody 50 | 51 | #BOXBOUNDARY 52 | outflow TypeBc1 MinX 53 | inflow TypeBc2 MaxX 54 | fixed TypeBc3 MinY 55 | fixed TypeBc4 MaxY 56 | fixed TypeBc5 MinZ 57 | fixed TypeBc6 MaxZ 58 | 59 | #USERSWITCH 60 | +init +ic +source StringSwitch 61 | 62 | #USERINPUTBEGIN ----------------------------------- 63 | 64 | #UseHotO 65 | T 66 | 1.0 67 | 0.0 68 | 1.0 69 | 70 | #UseTempCont 71 | F 72 | 73 | UseMarsB0 74 | T 75 | 60 76 | -180 77 | 0.0 78 | -99.5 79 | 25.26 80 | 81 | UseMarsAtm 82 | T 83 | Mars_TGCM_max.dat TGCMFilename 84 | 19 NAlt 85 | 86 | #SOLARCON 87 | issiC 88 | 89 | #REACTIONS 90 | F 91 | T 92 | 93 | 94 | #USERINPUTEND ------------------------------------- 95 | 96 | #PLASMA 97 | 1.0 FluidMass H [amu] 98 | 32.0 FluidMass O2 [amu] 99 | 16.0 FluidMass O [amu] 100 | 44.0 FluidMass CO2 [amu] 101 | 1.0 IonCharge H [e] 102 | 1.0 IonCharge O2 [e] 103 | 1.0 IonCharge O [e] 104 | 1.0 IonCharge CO2 [e] 105 | 1.0 ElectronTemperatureRatio 106 | 107 | #MULTIION 108 | 1e-9 LowDensityRatio 109 | 1e-11 LowPressureRatio 110 | F DoRestrictMultiIon 111 | 112 | #SAVERESTART 113 | F DoSaveRestart 114 | 115 | #SAVEPLOT 116 | 3 nPlotFile 117 | x=0 MHD idl StringPlot 118 | 2000 DnSavePlot 119 | -1.0 DtSavePlot 120 | -1.0 DxSavePlot 121 | y=0 MHD idl StringPlot 122 | 2000 DnSavePlot 123 | -1.0 DtSavePlot 124 | -1.0 DxSavePlot 125 | z=0 MHD idl StringPlot 126 | 2000 DnSavePlot 127 | -1.0 DtSavePlot 128 | -1.0 DxSavePlot 129 | 130 | 1d MHD idl StringPlot 131 | 1000 DnSavePlot 132 | -1.0 DtSavePlot 133 | -1.0 DxSavePlot 134 | 135 | #STOP 136 | 10 MaxIteration 137 | -1.0 tSimulationMax 138 | -------------------------------------------------------------------------------- /Param/MHDNONCONS/PARAM.in: -------------------------------------------------------------------------------- 1 | #PLANET 2 | NONE NamePlanet 3 | 4 | ! Use planetary units so that the multi-ion effects are in reasonable range 5 | #IOUNITS 6 | PLANETARY TypeIoUnit 7 | 8 | #NORMALIZATION 9 | PLANETARY TypeNormalization 10 | 11 | A single block 256x4x4 grid 12 | #CHECKGRIDSIZE 13 | 64 nI 14 | 2 nJ 15 | 2 nK 16 | 4 MinBlockAll 17 | 18 | Distance units are in km 19 | #GRID 20 | 4 nRootBlock1 21 | 1 nRootBlock2 22 | 1 nRootBlock3 23 | -128e3 xMin 24 | 128e3 xMax 25 | -1e3 yMin 26 | 1e3 yMax 27 | -1e3 zMin 28 | 1e3 zMax 29 | 30 | ! 90 per cent H+ and 10 percent O+ 31 | #SHOCKTUBE 32 | 5.0 LeftState Rho 33 | 400. LeftState Ux 34 | 0. LeftState Uy 35 | 0. LeftState Uz 36 | 0.0 LeftState Bx 37 | 0.0 LeftState By 38 | 5.0 LeftState Bz 39 | 0.01 LeftState p 40 | 18.0 RightState Rho 41 | 100.0 RightState Ux 42 | 0. RightState Uy 43 | 0. RightState Uz 44 | 0.0 RightState Bx 45 | 0.0 RightState By 46 | 18.0 RightState Bz 47 | 1.0 RightState p 48 | 49 | Sheared boundary condition is used for non-zero shock slope 50 | #OUTERBOUNDARY 51 | float TypeBc1 52 | float TypeBc2 53 | float TypeBc3 54 | float TypeBc4 55 | float TypeBc5 56 | float TypeBc6 57 | 58 | #TIMESTEPPING 59 | 2 nStage 60 | 0.8 CflExpl 61 | 62 | #SCHEME 63 | 2 nOrder 64 | Rusanov TypeFlux 65 | mc TypeLimiter 66 | 1.5 LimiterBeta 67 | 68 | #SAVELOGFILE 69 | T DoSaveLogfile 70 | raw StringLog 71 | 1 DnSaveLogfile 72 | -1. DtSaveLogfile 73 | 74 | #SAVERESTART 75 | F DoSaveRestart 76 | 77 | #SAVEINITIAL 78 | T DoSaveInitial 79 | 80 | #SAVEPLOT 81 | 1 nPlotFile 82 | cut MHD idl_ascii StringPlot 83 | -1 DnSavePlot 84 | 60.0 DtSavePlot 85 | -128e3 Coord1MinCut 86 | +128e3 Coord1MaxCut 87 | -0.1 Coord2MinCut 88 | 0.1 Coord2MaxCut 89 | -0.1 Coord3MinCut 90 | 0.1 Coord3MaxCut 91 | 0.0 DxSavePlot 92 | 93 | #STOP 94 | -1 MaxIteration 95 | 600.0 tSimulationMax 96 | -------------------------------------------------------------------------------- /Param/MULTIFLUID/PARAM.in: -------------------------------------------------------------------------------- 1 | #PLANET 2 | NONE NamePlanet 3 | 4 | #IOUNITS 5 | NONE TypeIoUnit 6 | 7 | #NORMALIZATION 8 | NONE TypeNormalization 9 | 10 | #DESCRIPTION 11 | Param/MULTIFLUID/PARAM.in: Brio-Wu shock tube for ions, Sod shock tube for neutrals StringDescription 12 | 13 | 4 blocks with 64x2x2 cells, -128 scaling.diff 30 | cat runlog_128 runlog_256 > scaling.log 31 | setenv COPY Sites/CRASH_SCALING/`date +%Y%m%d` 32 | if(-s scaling.diff)then 33 | scp scaling.diff ${HEROT}:${COPY}.error 34 | endif 35 | scp scaling.log ${HEROT}:${COPY}.log 36 | ssh ${HEROT} chmod go+r ${COPY}.\* 37 | endif 38 | 39 | exit 40 | -------------------------------------------------------------------------------- /Scripts/Scaling.job.jaguarpf-ext: -------------------------------------------------------------------------------- 1 | #!/bin/csh 2 | 3 | # This job script should be copied and modified by Scaling.pl 4 | 5 | # specify project number, if necessary 6 | ###PBS -A ... 7 | 8 | #PBS -V 9 | #PBS -N test 10 | #PBS -j oe 11 | 12 | # The number of cores (=size) should be a multiple of 16 13 | #PBS -l walltime=0:30:00 14 | #PBS -l size=16 15 | #PBS -l gres=widow2%widow3 16 | 17 | cd $PBS_O_WORKDIR 18 | 19 | # run CRASH with aprun -n ... or aprun -mode VN 20 | aprun -n $PBS_NNODES ./CRASH.exe 21 | 22 | # load opacity tables in next run(s) 23 | perl -pi -e 's/^LOOKUPTABLE/\#LOOKUPTABLE/; s/save/load/' PARAM.in 24 | -------------------------------------------------------------------------------- /Scripts/Scaling.job.pfe: -------------------------------------------------------------------------------- 1 | #!/bin/csh 2 | 3 | # This job script should be copied and modified by Scaling.pl 4 | 5 | #PBS -S /bin/csh 6 | #PBS -N CRASH 7 | 8 | # Use -q wide for 1024 or more cores. 9 | #PBS -q normal 10 | 11 | # To run on the 8-core Harpertown nodes (default now) 12 | #PBS -l select=16:ncpus=8:model=har 13 | 14 | # To run on the 8-core Nehalem nodes (2.4 times faster than Hapertown) 15 | ### PBS -l select=16:ncpus=8:model=neh 16 | 17 | # To run on the 12-core Westmere nodes use this (2 times faster than Hapertown) 18 | ### PBS -l select=16:ncpus=12:model=wes 19 | #PBS -l walltime=1:00:00 20 | #PBS -j oe 21 | #PBS -m e 22 | 23 | # Specify group if necessary 24 | ### PBS -W group_list=... 25 | 26 | # cd into the run directory 27 | cd $PBS_O_WORKDIR 28 | 29 | # run CRASH (the number of processors is already specified above) 30 | # switch on BATL and set plot directory just before the run starts 31 | perl -pi -e 's/^BATL/\#BATL/' PARAM.in 32 | perl -pi -e 's/.*NamePlotDir/plot\t\tNamePlotDir/' PARAM.in 33 | mpiexec CRASH.exe > runlog 34 | 35 | # load opacity tables in next run(s) 36 | perl -pi -e 's/^LOOKUPTABLE/\#LOOKUPTABLE/; s/save/load/' PARAM.in 37 | 38 | # exit here unless we want to do more runs with modified PARAM.in 39 | exit 40 | 41 | # switch off BATL and use message_pass_cell 42 | perl -pi -e 's/^\#BATL/BATL/' PARAM.in 43 | perl -pi -e 's/^\#MESSAGEPASS/MESSAGEPASS/' PARAM.in 44 | perl -pi -e 's/.*NamePlotDir/plot\t\tNamePlotDir/' PARAM.in 45 | mpiexec CRASH.exe > runlog_bats_cell 46 | 47 | # switch off BATL and use message_pass_dir 48 | perl -pi -e 's/^\#BATL/BATL/' PARAM.in 49 | perl -pi -e 's/^MESSAGEPASS/\#MESSAGEPASS/' PARAM.in 50 | perl -pi -e 's/.*NamePlotDir/plot\t\tNamePlotDir/' PARAM.in 51 | mpiexec CRASH.exe > runlog_bats_dir 52 | -------------------------------------------------------------------------------- /Scripts/Scaling.job.ubgl: -------------------------------------------------------------------------------- 1 | #!/bin/csh 2 | 3 | # news job.lim.ubgl 4 | # mjstat 5 | # This job script should be copied and modified by Scaling.pl 6 | 7 | #MSUB -A michigan 8 | #MSUB -l partition=ubgl 9 | #MSUB -l nodes=128 10 | #MSUB -l walltime=0:30:00 11 | #MSUB -o runlog 12 | 13 | # Use pdebug for 128 or 512 nodes, pshort from 1k to 4k nodes, pbatch for more than 4k node 14 | #MSUB -q pbatch 15 | 16 | # cd into the run directory 17 | cd /g/g99/toth2/crash/CRASH_ubgl/SCALING/run 18 | 19 | # set plot directory just before the run starts 20 | ### perl -pi -e 's/.*NamePlotDir/plot\t\tNamePlotDir/' PARAM.in 21 | 22 | # run CRASH with mpirun -np ... or mpirun -mode VN 23 | mpirun -np 128 -exe `pwd`/CRASH.exe -cwd `pwd` 24 | 25 | # load opacity tables in next run(s) 26 | perl -pi -e 's/^LOOKUPTABLE/\#LOOKUPTABLE/; s/save/load/' PARAM.in 27 | -------------------------------------------------------------------------------- /Scripts/Scaling.job.yslogin: -------------------------------------------------------------------------------- 1 | #!/bin/csh 2 | 3 | # This job script should be copied and modified by Scaling.pl 4 | 5 | #BSUB -P P35201105 # project code 6 | #BSUB -W 0:10 # wall-clock time (hrs:mins) 7 | #BSUB -n 32 # number of tasks in job 8 | #BSUB -R "span[ptile=16]" # run 16 MPI tasks per node 9 | #BSUB -J CRASH_hi # job name 10 | #BSUB -o CRASH.%J.out # output file name: %J is replaced by the job ID 11 | #BSUB -e CRASH.%J.err # error file name: %J is replaced by the job ID 12 | #BSUB -q regular # queue 13 | #BSUB -N # notify me when job is finished 14 | 15 | # run the executable 16 | mpirun.lsf ./CRASH.exe > runlog_`date +%y%m%d%H%M` 17 | -------------------------------------------------------------------------------- /Scripts/TEC/ppA: -------------------------------------------------------------------------------- 1 | #!/bin/sh 2 | # Copyright (C) 2002 Regents of the University of Michigan, portions used with permission 3 | # For more information, see http://csem.engin.umich.edu/tools/swmf 4 | 5 | echo ' ' 6 | echo 'Proccessing *.dat files' 7 | echo ' ' 8 | 9 | dir='.' 10 | for i in `ls $dir/*.dat` 11 | do 12 | echo ' processing' $i 13 | /usr/contrib/bin/preplot ${i} 14 | done 15 | 16 | echo 'ls -l' 17 | echo ' ' 18 | ls -l 19 | 20 | -------------------------------------------------------------------------------- /Scripts/TEC/ppAr: -------------------------------------------------------------------------------- 1 | #!/bin/sh 2 | # Copyright (C) 2002 Regents of the University of Michigan, portions used with permission 3 | # For more information, see http://csem.engin.umich.edu/tools/swmf 4 | 5 | echo ' ' 6 | echo 'Proccessing *.dat files' 7 | echo ' ' 8 | 9 | dir='.' 10 | for i in `ls $dir/*.dat` 11 | do 12 | echo ' processing' $i 13 | /usr/contrib/bin/preplot ${i} 14 | /bin/rm ${i} 15 | done 16 | 17 | echo 'ls -l' 18 | echo ' ' 19 | ls -l 20 | 21 | -------------------------------------------------------------------------------- /Scripts/TEC/ppgz: -------------------------------------------------------------------------------- 1 | #!/bin/csh 2 | # Copyright (C) 2002 Regents of the University of Michigan, portions used with permission 3 | # For more information, see http://csem.engin.umich.edu/tools/swmf 4 | # tec2bin: *.dat files into *.plt files using preplot 5 | # 6 | 7 | if($1 == 'help' || $1 == '-h' || $1 == '-help')then 8 | echo "Usage:" 9 | echo "" 10 | echo "ppgz [KEEP]" 11 | echo "" 12 | echo "This routine is intended for use only with tecplot files, Which" 13 | echo "it converts from ascii(.dat) to binary(.plt). The script only 14 | echo "runs on platforms which have the preplot program installted in 15 | echo "/usr/contrib/bin/preplot." 16 | echo "" 17 | echo "The script looks for .dat.gz files to process. It unzips them," 18 | echo "runs preplot on them. By default it then deletes them." 19 | echo "" 20 | echo "If the 2nd argument, KEEP, exists the files are not deleted but" 21 | echo "are gzipped again for storage." 22 | echo "" 23 | echo "Examples:" 24 | echo "" 25 | echo "ppgz =>process all .dat.gz files in current dir." 26 | echo "ppgz KEEP =>process all .dat.gz files in IO2/ dir. Do not delete." 27 | echo "" 28 | exit 29 | endif 30 | 31 | set dir='./' 32 | 33 | echo 'Seaching for files in '${dir}' --- preplot logfile: '${dir}'preplot.log' 34 | echo '' 35 | 36 | foreach i (${dir}*.dat.gz) 37 | set f = `dirname $i`/`basename $i .dat.gz` 38 | echo '-----------------------------------------------------------------' 39 | echo 'processing file: '${i} 40 | echo ' gunzip ... ('${i}')' 41 | gunzip $i 42 | echo ' preplot ... ('${f}'.dat => '${f}'.plt)' 43 | if (`hostname` == rhodot.engin.umich.edu) \ 44 | /usr/local/bin/preplot ${f}.dat >> ${dir}preplot.log 45 | if (`hostname` != rhodot.engin.umich.edu) \ 46 | /usr/contrib/bin/preplot ${f}.dat >> ${dir}preplot.log 47 | if($1 == '')then 48 | echo ' deleting ... ('${f}'.dat)' 49 | rm -f ${f}.dat 50 | else 51 | echo ' gzip ... ('${f}'.dat => '${f}'.dat.gz)' 52 | gzip ${f}.dat 53 | endif 54 | end 55 | -------------------------------------------------------------------------------- /Scripts/cat_sat_log.pl: -------------------------------------------------------------------------------- 1 | #!/usr/bin/perl -p 2 | # Copyright (C) 2002 Regents of the University of Michigan, portions used with permission 3 | # For more information, see http://csem.engin.umich.edu/tools/swmf 4 | # Usage: cat_sat_log.pl logfile1 logfile2 .... > log_all.log 5 | # Example: cat_sat_log.pl log_*.log > log_all.log 6 | $_ = '' unless ($. == 1 or $. == 2 or /^[\s\d\.eEdD\+\-]+$/); 7 | -------------------------------------------------------------------------------- /Scripts/shift_satellite.pl: -------------------------------------------------------------------------------- 1 | #!/usr/bin/perl -p -s 2 | # Copyright (C) 2002 Regents of the University of Michigan, portions used with permission 3 | # For more information, see http://csem.engin.umich.edu/tools/swmf 4 | # Shift the satellite trajectory by a constant vector (dx, dy, dz) 5 | # Usage: shift_satellite.pl -dx=1.3 -dy=-1.2 -dz=0.1 ORIGINAL.sat > SHIFTED.sat 6 | s/([\-\d\.]+)\s+([\-\d\.]+)\s+([\-\d\.]+)$/sprintf("%10.5f %10.5f %10.5f",$1+$dx,$2+$dy,$3+$dz)/e; 7 | 8 | 9 | -------------------------------------------------------------------------------- /src/Makefile.RULES.all: -------------------------------------------------------------------------------- 1 | # Syntax: 2 | # Conditional in Perl syntax using the variables known by Config.pl 3 | # Compilation rules preceeded with a single TAB character. 4 | # Empty line before next conditional. 5 | 6 | # Special rule to compile HDF5 code 7 | $Hdf5 eq "yes" and $Compiler !~ /(cray|ifort)ftn/ 8 | h5pfc ${Cflag2} ModHdf5.f90 9 | 10 | # Particle field line module works with 0 optimization only 11 | $Compiler eq "crayftn" 12 | ${COMPILE.f90} ${Cflag0} ModParticleFieldLine.f90 13 | 14 | # There is no point optimizing this code, speed up compilation 15 | 1 16 | ${COMPILE.f90} ${Cflag0} ModSetParameters.f90 17 | 18 | # compiling ModWritePlot takes a long time with ifort with -O2 19 | $Compiler eq "ifort" 20 | ${COMPILE.f90} ${Cflag1} ModWritePlot.f90 21 | 22 | # compiling ModWritePlot with -O2 and nagfor + GNU gcc fails on M1 architecture 23 | $Compiler eq "nagfor" 24 | ${COMPILE.f90} ${Cflag0} ModWritePlot.f90 25 | 26 | # who knows why 27 | $Compiler eq "nagfor" 28 | ${COMPILE.f90} ${Cflag0} ModRadioWaveImage.f90 29 | -------------------------------------------------------------------------------- /src/Makefile.def: -------------------------------------------------------------------------------- 1 | # Local Makefile.def 2 | 3 | localdefault: help # default target 4 | 5 | DIR = ${GMDIR} 6 | SCRIPTDIR = ${GMDIR}/share/Scripts 7 | LIBDIR = ${GMDIR}/src 8 | RUNDIR = ${GMDIR}/run 9 | BINDIR = ${GMDIR}/src 10 | INCLDIR = ${DIR}/share/include 11 | SHAREDIR = ${GMDIR}/share/Library/src 12 | UTILDIR = ${GMDIR}/util 13 | 14 | TIMINGDIR = ${UTILDIR}/TIMING/src 15 | #TIMINGDIR = ${UTILDIR}/TIMING/srcEmpty 16 | NOMPIDIR = ${UTILDIR}/NOMPI/src 17 | MAGNETOGRAMDIR = ${UTILDIR}/DATAREAD/srcMagnetogram 18 | DEMTDIR = ${UTILDIR}/DATAREAD/srcDemt 19 | EMPIRICALEEDIR = ${UTILDIR}/EMPIRICAL/srcEE 20 | EMPIRICALSCDIR = ${UTILDIR}/EMPIRICAL/srcSC 21 | CRASHDIR = ${UTILDIR}/CRASH/src 22 | -------------------------------------------------------------------------------- /src/ModExtraVariables.f90: -------------------------------------------------------------------------------- 1 | ! Copyright (C) 2002 Regents of the University of Michigan, 2 | ! portions used with permission 3 | ! For more information, see http://csem.engin.umich.edu/tools/swmf 4 | module ModExtraVariables 5 | 6 | ! Define indexes that are not used in most equation modules 7 | 8 | ! Set impossible value for indexes 9 | integer, parameter :: Ppar_ = 1, iPparIon_I(1) = 1 10 | integer, parameter :: Pe_ = 1, Pepar_ = 1 11 | integer, parameter :: Hyp_ = 1 12 | integer, parameter :: Erad_ = 1, ExtraEint_ = 1 13 | integer, parameter :: Ew_ = 1 14 | integer, parameter :: SignB_ = 1, BperU_ = 1 15 | integer, parameter :: Te0_ = 1 16 | integer, parameter :: Ehot_ = 1 17 | 18 | ! Electric field and hyperbolic scalar field 19 | integer, parameter :: Ex_ = 1, Ey_ = 2, Ez_ = 3 20 | integer, parameter :: HypE_ = 1 21 | 22 | ! Number of electron fluids 23 | integer, parameter :: nElectronFluid = 0 24 | 25 | ! The named index range for frequency bins in multi-group 26 | integer, parameter :: nWave = 1 27 | integer, parameter :: WaveFirst_ = 1, WaveLast_ = 1 28 | ! For Alfven wave energy difference, \rho\delta u^2 - \delta B^2/\mu_0 29 | integer, parameter :: WDiff_ = 1 30 | 31 | ! Perpendicular correlation length (Assumed to be the same for major 32 | ! and minor turbulence energy density) 33 | integer, parameter :: Lperp_ = 1 34 | 35 | ! Correlation lengths (may be different for major 36 | ! and minor turbulence energy density) 37 | integer, parameter :: LcorrFirst_ = 1, LcorrLast_ = 1, nLcorr = 0 38 | 39 | ! The named index range for velocity bins in pickup ions 40 | integer, parameter :: nPui = 1 41 | integer, parameter :: PuiFirst_ = 1, PuiLast_ = 1 42 | 43 | ! The named index range for material levels 44 | integer, parameter :: nMaterial = 1 45 | integer, parameter :: MaterialFirst_ = 1, MaterialLast_ = 1 46 | 47 | ! The named index range for charge states 48 | integer, parameter :: nChargeStateAll = 1 49 | integer, parameter :: ChargeStateFirst_ = 1, ChargeStateLast_ = 1 50 | integer, parameter :: nElement = 1, nChargeState_I(nElement) = [1] 51 | character(len=2), parameter :: NameElement_I(1:nElement) = [' '] 52 | 53 | ! The named index range for species and their masses 54 | integer, parameter :: SpeciesFirst_ = 1, SpeciesLast_ = 1 55 | real :: MassSpecies_V(SpeciesFirst_:SpeciesLast_) 56 | 57 | end module ModExtraVariables 58 | !============================================================================== 59 | -------------------------------------------------------------------------------- /src/ModNodes.f90: -------------------------------------------------------------------------------- 1 | ! Copyright (C) 2002 Regents of the University of Michigan, 2 | ! portions used with permission 3 | ! For more information, see http://csem.engin.umich.edu/tools/swmf 4 | module ModNodes 5 | 6 | use BATL_lib, ONLY: & 7 | test_start, test_stop, iProc 8 | 9 | use ModSize 10 | use ModIO, ONLY: iUnitOut, write_prefix 11 | 12 | implicit none 13 | save 14 | 15 | ! Block node-centered MHD numberings 16 | integer :: nNodeALL 17 | integer, allocatable :: iNodeLocal_NB(:,:,:,:) 18 | integer, allocatable :: iNodeGlobal_NB(:,:,:,:) 19 | logical, allocatable :: IsNodeUnique_NB(:,:,:,:) 20 | 21 | contains 22 | !============================================================================ 23 | subroutine init_mod_nodes 24 | 25 | logical:: DoTest 26 | character(len=*), parameter:: NameSub = 'init_mod_nodes' 27 | !-------------------------------------------------------------------------- 28 | call test_start(NameSub, DoTest) 29 | if(allocated(iNodeLocal_NB)) RETURN 30 | allocate(iNodeLocal_NB(1:nI+1,1:nJ+1,1:nK+1,MaxBlock)) 31 | allocate(iNodeGlobal_NB(1:nI+1,1:nJ+1,1:nK+1,MaxBlock)) 32 | allocate(IsNodeUnique_NB(1:nI+1,1:nJ+1,1:nK+1,MaxBlock)) 33 | if(iProc==0)then 34 | call write_prefix 35 | write(iUnitOut,'(a)') 'init_mod_nodes allocated arrays' 36 | end if 37 | 38 | call test_stop(NameSub, DoTest) 39 | end subroutine init_mod_nodes 40 | !============================================================================ 41 | subroutine clean_mod_nodes 42 | 43 | logical:: DoTest 44 | character(len=*), parameter:: NameSub = 'clean_mod_nodes' 45 | !-------------------------------------------------------------------------- 46 | call test_start(NameSub, DoTest) 47 | if(.not.allocated(iNodeLocal_NB)) RETURN 48 | deallocate(iNodeLocal_NB) 49 | deallocate(iNodeGlobal_NB) 50 | deallocate(IsNodeUnique_NB) 51 | 52 | if(iProc==0)then 53 | call write_prefix 54 | write(iUnitOut,'(a)') 'clean_mod_nodes deallocated arrays' 55 | end if 56 | 57 | call test_stop(NameSub, DoTest) 58 | end subroutine clean_mod_nodes 59 | !============================================================================ 60 | 61 | end module ModNodes 62 | !============================================================================== 63 | -------------------------------------------------------------------------------- /src/ModSingleFluid.f90: -------------------------------------------------------------------------------- 1 | ! Copyright (C) 2002 Regents of the University of Michigan, 2 | ! portions used with permission 3 | ! For more information, see http://csem.engin.umich.edu/tools/swmf 4 | module ModSingleFluid 5 | 6 | ! Define single fluid values for multifluid parameters 7 | 8 | implicit none 9 | 10 | integer, parameter :: nFluid = 1 11 | integer, parameter :: nIonFluid = 1 12 | logical, parameter :: IsMhd = .true. 13 | character(len=1), parameter :: NameFluid_I(nFluid) = [' '] 14 | 15 | ! Default is proton mass, but it can be changed 16 | real :: MassFluid_I(nFluid) = 1.0 17 | !$acc declare create(MassFluid_I) 18 | 19 | end module ModSingleFluid 20 | !============================================================================== 21 | 22 | -------------------------------------------------------------------------------- /src/ModSize.f90: -------------------------------------------------------------------------------- 1 | module ModSize 2 | 3 | use BATL_size 4 | 5 | ! Named indexes for coordinate directions 6 | use BATL_geometry, ONLY: & 7 | x_, y_, z_, r_, Phi_, Theta_, Lon_, Lat_ 8 | 9 | end module ModSize 10 | !============================================================================== 11 | -------------------------------------------------------------------------------- /src/ModUaCoupling.f90: -------------------------------------------------------------------------------- 1 | ! Copyright (C) 2002 Regents of the University of Michigan, 2 | ! portions used with permission 3 | ! For more information, see http://csem.engin.umich.edu/tools/swmf 4 | 5 | module ModUaCoupling 6 | 7 | use BATL_lib, ONLY: 8 | 9 | ! Routines related to the coupline with the Upper Atmosphere component 10 | 11 | implicit none 12 | 13 | SAVE 14 | 15 | private ! except 16 | 17 | public:: gm_init_ua_array 18 | 19 | real, public, allocatable:: NumDenNuSpecies_CBI(:,:,:,:,:) 20 | 21 | contains 22 | !============================================================================ 23 | subroutine gm_init_ua_array(nI,nJ,nK,MaxBlock,MaxNuSpecies) 24 | 25 | integer, intent(in):: nI,nJ,nK,MaxBlock,MaxNuSpecies 26 | !-------------------------------------------------------------------------- 27 | allocate(NumDenNuSpecies_CBI(nI,nJ,nK,MaxBlock,MaxNuSpecies)) 28 | 29 | end subroutine gm_init_ua_array 30 | !============================================================================ 31 | end module ModUaCoupling 32 | !============================================================================== 33 | -------------------------------------------------------------------------------- /src/ModUserDefault.f90: -------------------------------------------------------------------------------- 1 | ! Copyright (C) 2002 Regents of the University of Michigan, 2 | ! portions used with permission 3 | ! For more information, see http://csem.engin.umich.edu/tools/swmf 4 | module ModUser 5 | ! This is the default user module which contains empty methods defined 6 | ! in ModUserEmpty.f90 7 | 8 | use ModUserEmpty 9 | 10 | include 'user_module.h' ! list of public methods 11 | 12 | character (len=*), parameter :: NameUserFile = 'src/ModUserDefault.f90' 13 | character (len=*), parameter :: NameUserModule = 'DEFAULT EMPTY ROUTINES' 14 | 15 | end module ModUser 16 | !============================================================================== 17 | -------------------------------------------------------------------------------- /src/show_git_info.h.orig: -------------------------------------------------------------------------------- 1 | 2 | write(*,'(100(a/))') & 3 | 'Git repo=git@github.com:SWMFsoftware/BATSRUS', & 4 | ' GIT REFERENCES HAVE NOT BEEN OBTAINED!' 5 | -------------------------------------------------------------------------------- /src/user_module.h: -------------------------------------------------------------------------------- 1 | ! Copyright (C) 2002 Regents of the University of Michigan, 2 | ! portions used with permission 3 | ! For more information, see http://csem.engin.umich.edu/tools/swmf 4 | 5 | use ModUtilities, ONLY: CON_stop 6 | 7 | implicit none 8 | 9 | SAVE ! save all module variables 10 | 11 | private ! except 12 | public :: NameUserFile 13 | public :: NameUserModule 14 | public :: user_set_boundary_cells 15 | public :: user_set_face_boundary 16 | public :: user_set_cell_boundary 17 | public :: user_initial_perturbation 18 | public :: user_set_ics 19 | public :: user_init_session 20 | public :: user_specify_region 21 | public :: user_amr_criteria 22 | public :: user_read_inputs 23 | public :: user_get_log_var 24 | public :: user_set_plot_var 25 | public :: user_calc_sources_expl 26 | public :: user_calc_sources_impl 27 | public :: user_init_point_implicit 28 | public :: user_get_b0 29 | public :: user_update_states 30 | public :: user_calc_timestep 31 | public :: user_io_units 32 | public :: user_normalization 33 | public :: user_set_resistivity 34 | public :: user_material_properties 35 | public :: user_action 36 | public :: i_type_block_user 37 | -------------------------------------------------------------------------------- /srcEquation/ModEquationHd.f90: -------------------------------------------------------------------------------- 1 | ! Copyright (C) 2002 Regents of the University of Michigan 2 | ! portions used with permission 3 | ! For more information, see http://csem.engin.umich.edu/tools/swmf 4 | module ModVarIndexes 5 | 6 | use ModSingleFluid, Redefine => IsMhd 7 | use ModExtraVariables 8 | 9 | implicit none 10 | 11 | save 12 | 13 | character (len=*), parameter :: & 14 | NameEquationFile = "ModEquationHd.f90" 15 | 16 | ! This equation module contains the standard hydro equations. 17 | character (len=*), parameter :: & 18 | NameEquation = 'HD' 19 | 20 | integer, parameter :: nVar = 5 21 | 22 | logical, parameter :: IsMhd = .false. 23 | 24 | ! Named indexes for State_VGB and other variables 25 | ! These indexes should go subsequently, from 1 to nVar+nFluid. 26 | ! The energies are handled as an extra variable, so that we can use 27 | ! both conservative and non-conservative scheme and switch between them. 28 | integer, parameter :: & 29 | Rho_ = 1, & 30 | RhoUx_ = 2, Ux_ = 2, & 31 | RhoUy_ = 3, Uy_ = 3, & 32 | RhoUz_ = 4, Uz_ = 4, & 33 | p_ = 5, & 34 | Energy_ = nVar+1 35 | 36 | ! This allows to calculate RhoUx_ as RhoU_+x_ and so on. 37 | integer, parameter :: U_ = Ux_ - 1, RhoU_ = RhoUx_-1 38 | 39 | ! The default values for the state variables: 40 | ! Variables which are physically positive should be set to 1, 41 | ! variables that can be positive or negative should be set to 0: 42 | real, parameter :: DefaultState_V(nVar+nFluid) = [ & 43 | 1.0, & ! Rho_ 44 | 0.0, & ! RhoUx_ 45 | 0.0, & ! RhoUy_ 46 | 0.0, & ! RhoUz_ 47 | 1.0, & ! p_ 48 | 1.0 ] ! Energy_ 49 | 50 | ! The names of the variables used in i/o 51 | character(len=5) :: NameVar_V(nVar+nFluid) = [ & 52 | 'Rho ', & ! Rho_ 53 | 'Mx ', & ! RhoUx_ 54 | 'My ', & ! RhoUy_ 55 | 'Mz ', & ! RhoUz_ 56 | 'P ', & ! p_ 57 | 'E '] ! Energy_ 58 | 59 | ! Bx_, By_, Bz_ have to be defined so that the code compiles 60 | ! but the Bx_ = Ux_ choice indicates that B is not used (see UseB in ModMain) 61 | integer, parameter :: Bx_ = Ux_, By_ = Uy_, Bz_ = Uz_, B_ = U_ 62 | 63 | ! There are no extra scalars 64 | integer, parameter :: ScalarFirst_ = 2, ScalarLast_ = 1 65 | 66 | integer, parameter :: iRho_I(nFluid) = [Rho_] 67 | integer, parameter :: iRhoUx_I(nFluid) = [RhoUx_] 68 | integer, parameter :: iRhoUy_I(nFluid) = [RhoUy_] 69 | integer, parameter :: iRhoUz_I(nFluid) = [RhoUz_] 70 | integer, parameter :: iP_I(nFluid) = [p_] 71 | 72 | end module ModVarIndexes 73 | !============================================================================== 74 | -------------------------------------------------------------------------------- /srcEquation/ModEquationMhd.f90: -------------------------------------------------------------------------------- 1 | ! Copyright (C) 2002 Regents of the University of Michigan 2 | ! portions used with permission 3 | ! For more information, see http://csem.engin.umich.edu/tools/swmf 4 | module ModVarIndexes 5 | 6 | use ModSingleFluid 7 | use ModExtraVariables 8 | 9 | implicit none 10 | 11 | save 12 | 13 | character (len=*), parameter :: & 14 | NameEquationFile = "ModEquationMhd.f90" 15 | 16 | character (len=*), parameter :: & 17 | NameEquation = 'MHD' 18 | 19 | ! Number of variables without energy: 20 | integer, parameter :: nVar = 8 21 | 22 | ! Named indexes for State_VGB and other variables 23 | ! These indexes should go subsequently, from 1 to nVar+1. 24 | ! The energy is handled as an extra variable, so that we can use 25 | ! both conservative and non-conservative scheme and switch between them. 26 | integer, parameter :: & 27 | Rho_ = 1, & 28 | RhoUx_ = 2, & 29 | RhoUy_ = 3, & 30 | RhoUz_ = 4, & 31 | Bx_ = 5, & 32 | By_ = 6, & 33 | Bz_ = 7, & 34 | p_ = nVar, & 35 | Energy_= nVar+1 36 | 37 | ! This allows to calculate RhoUx_ as rhoU_+x_ and so on. 38 | integer, parameter :: RhoU_ = RhoUx_-1, B_ = Bx_-1 39 | 40 | ! These arrays are useful for multifluid 41 | integer, parameter :: iRho_I(nFluid) = [Rho_] 42 | integer, parameter :: iRhoUx_I(nFluid) = [RhoUx_] 43 | integer, parameter :: iRhoUy_I(nFluid) = [RhoUy_] 44 | integer, parameter :: iRhoUz_I(nFluid) = [RhoUz_] 45 | integer, parameter :: iP_I(nFluid) = [p_] 46 | 47 | ! The default values for the state variables: 48 | ! Variables which are physically positive should be set to 1, 49 | ! variables that can be positive or negative should be set to 0: 50 | real, parameter :: DefaultState_V(nVar+1) = [ & 51 | 1.0, & ! Rho_ 52 | 0.0, & ! RhoUx_ 53 | 0.0, & ! RhoUy_ 54 | 0.0, & ! RhoUz_ 55 | 0.0, & ! Bx_ 56 | 0.0, & ! By_ 57 | 0.0, & ! Bz_ 58 | 1.0, & ! p_ 59 | 1.0 ] ! Energy_ 60 | 61 | ! The names of the variables used in i/o 62 | character(len=3) :: NameVar_V(nVar+1) = [ & 63 | 'Rho', & ! Rho_ 64 | 'Mx ', & ! RhoUx_ 65 | 'My ', & ! RhoUy_ 66 | 'Mz ', & ! RhoUz_ 67 | 'Bx ', & ! Bx_ 68 | 'By ', & ! By_ 69 | 'Bz ', & ! Bz_ 70 | 'p ', & ! p_ 71 | 'e ' ] ! Energy_ 72 | 73 | !$acc declare create (NameVar_V) 74 | 75 | ! Primitive variable names 76 | integer, parameter :: U_ = RhoU_, Ux_ = RhoUx_, Uy_ = RhoUy_, Uz_ = RhoUz_ 77 | 78 | ! There are no extra scalars 79 | integer, parameter :: ScalarFirst_ = 2, ScalarLast_ = 1 80 | 81 | end module ModVarIndexes 82 | !============================================================================== 83 | 84 | -------------------------------------------------------------------------------- /srcEquation/ModEquationScalar.f90: -------------------------------------------------------------------------------- 1 | ! Copyright (C) 2002 Regents of the University of Michigan 2 | ! portions used with permission 3 | ! For more information, see http://csem.engin.umich.edu/tools/swmf 4 | module ModVarIndexes 5 | 6 | use ModSingleFluid, Redefine => IsMhd 7 | use ModExtraVariables 8 | 9 | implicit none 10 | 11 | save 12 | 13 | character (len=*), parameter :: & 14 | NameEquationFile = "ModEquationScalar.f90" 15 | 16 | ! This equation module contains the standard hydro equations. 17 | character (len=*), parameter :: & 18 | NameEquation = 'Scalar' 19 | 20 | integer, parameter :: nVar = 1 21 | 22 | logical, parameter :: IsMhd = .false. 23 | 24 | ! Named indexes for State_VGB and other variables 25 | ! These indexes should go subsequently, from 1 to nVar+nFluid. 26 | ! The energies are handled as an extra variable, so that we can use 27 | ! both conservative and non-conservative scheme and switch between them. 28 | integer, parameter :: & 29 | Rho_ = 1, & 30 | Energy_ = nVar+1 31 | 32 | ! These need to be defined 33 | integer, parameter :: RhoUx_=1, RhoUy_=1, RhoUz_=3 34 | integer, parameter :: Ux_=1, Uy_=2, Uz_=3, p_=1 35 | 36 | ! This allows to calculate RhoUx_ as RhoU_+x_ and so on. 37 | integer, parameter :: U_ = Ux_ - 1, RhoU_ = RhoUx_-1 38 | 39 | ! The default values for the state variables: 40 | ! Variables which are physically positive should be set to 1, 41 | ! variables that can be positive or negative should be set to 0: 42 | real, parameter :: DefaultState_V(nVar+nFluid) = [ & 43 | 1.0, & ! Rho_ 44 | 1.0 ] ! Energy_ 45 | 46 | ! The names of the variables used in i/o 47 | character(len=3) :: NameVar_V(nVar+nFluid) = [ & 48 | 'Rho', & ! Rho_ 49 | 'E '] ! Energy_ 50 | 51 | ! Bx_, By_, Bz_ have to be defined so that the code compiles 52 | ! but the Bx_ = Ux_ choice indicates that B is not used (see UseB in ModMain) 53 | integer, parameter :: Bx_ = Ux_, By_ = Uy_, Bz_ = Uz_, B_ = U_ 54 | 55 | ! There are no extra scalars 56 | integer, parameter :: ScalarFirst_ = 2, ScalarLast_ = 1 57 | 58 | integer, parameter :: iRho_I(nFluid) = [Rho_] 59 | integer, parameter :: iRhoUx_I(nFluid) = [RhoUx_] 60 | integer, parameter :: iRhoUy_I(nFluid) = [RhoUy_] 61 | integer, parameter :: iRhoUz_I(nFluid) = [RhoUz_] 62 | integer, parameter :: iP_I(nFluid) = [p_] 63 | 64 | end module ModVarIndexes 65 | !============================================================================== 66 | -------------------------------------------------------------------------------- /srcInterface/Makefile: -------------------------------------------------------------------------------- 1 | # Copyright (C) 2002 Regents of the University of Michigan, portions used with permission 2 | # For more information, see http://csem.engin.umich.edu/tools/swmf 3 | 4 | SHELL =/bin/sh 5 | 6 | include ../Makefile.def 7 | include ../Makefile.conf 8 | -include Makefile.DEPEND 9 | 10 | # Configurable objects 11 | 12 | OBJ_IH = GM_couple_ih.o read_ih_buffer.o #^CMP IF IH 13 | OBJ_IM = GM_couple_im.o #^CMP IF IM 14 | OBJ_PS = GM_couple_ps.o #^CMP IF PS 15 | OBJ_RB = GM_couple_rb.o #^CMP IF RB 16 | OBJ_IE = GM_couple_ie.o #^CMP IF IE 17 | OBJ_PT = GM_couple_pt.o #^CMP IF PT 18 | OBJ_PC = GM_couple_pc.o #^CMP IF PC 19 | OBJ_PW = GM_couple_pw.o read_pw_buffer.o #^CMP IF PW 20 | OBJ_UA = GM_couple_ua.o #^CMP IF UA 21 | 22 | OBJECTS = \ 23 | ModGridDescriptor.o\ 24 | GM_wrapper.o\ 25 | $(OBJ_IH) $(OBJ_IM) $(OBJ_RB) $(OBJ_IE) $(OBJ_PT) $(OBJ_PW) $(OBJ_PC) $(OBJ_PS) $(OBJ_UA) 26 | 27 | SEARCHDIR = -I${SHAREDIR} -I${COUPLERDIR} -I${LIBRARYDIR} -I${MAGNETOGRAMDIR} \ 28 | -I../src -I../srcBATL -I${CRASHDIR} 29 | 30 | DEPEND: 31 | @perl ${SCRIPTDIR}/depend.pl ${SEARCHDIR} ${OBJECTS} 32 | 33 | # 34 | # Making executables and libraries 35 | # 36 | MY_LIB = ${LIBDIR}/libGM.a 37 | 38 | LIB: DEPEND 39 | $(MAKE) ${MY_LIB} 40 | @echo 41 | @echo ${MY_LIB} has been brought up to date. 42 | @echo 43 | 44 | BATSRUS_LIB = ../src/libBATSRUS.a 45 | 46 | ${MY_LIB}: ${BATSRUS_LIB} ${OBJECTS} 47 | cp -f ${BATSRUS_LIB} ${MY_LIB} 48 | ${AR} ${MY_LIB} ${OBJECTS} 49 | 50 | clean: cleanfiles 51 | 52 | distclean: clean 53 | -------------------------------------------------------------------------------- /srcPostProc/Makefile.RULES.all: -------------------------------------------------------------------------------- 1 | $OS eq "Linux" and $Compiler eq "ifort" 2 | ${COMPILE.f90} ${Cflag1} PostIDL.f90 3 | -------------------------------------------------------------------------------- /srcPostProc/test-SPECTRUM.in: -------------------------------------------------------------------------------- 1 | Test imput file for spectrum code 2 | #ECHO 3 | F DoEcho 4 | 5 | #VERBOSE 6 | F IsVerbose 7 | 8 | #OUTFILE 9 | test-spectrum.out NameSpectrumFile 10 | 11 | DATAFILE 12 | box.out NameDataFile 13 | real4 TypeDataFile 14 | 15 | #TABLEFILE 16 | SPECTRUM_chianti_tbl.dat NameTableFile 17 | 1000 nMaxLine 18 | 19 | #WAVELENGTHINTERVAL 20 | 5 nWavelengthInterval 21 | 200. IntervalMin 22 | 210. IntervalMax 23 | 300. IntervalMin 24 | 310. IntervalMax 25 | 400. IntervalMin 26 | 410. IntervalMax 27 | 945. IntervalMin 28 | 955. IntervalMax 29 | 960. IntervalMin 30 | 974. IntervalMax 31 | 0.001 SizeWavelengthBin 32 | 33 | #DATABLOCK 34 | 1 n1 35 | 1 n2 36 | 1 n3 37 | 0.001 Dx 38 | 0.001 Dy 39 | 0.001 Dz 40 | 41 | INSTRUMENT 42 | EIS NameInstrument 43 | 44 | #UNIFORMDATA 45 | T IsPpar 46 | T IsPe 47 | 8.3e7 RhoUni 48 | 0 UxUni 49 | 0 UyUni 50 | 0 UzUni 51 | 0 BxUni 52 | 0 ByUni 53 | 0 BzUni 54 | 1e6 TparUni 55 | 1e6 TperpUni 56 | 1e6 TeUni 57 | 0 I01Uni 58 | 0 I02Uni 59 | 60 | #NOALFVEN 61 | T IsNoAlfven 62 | 63 | ISPERMUTEAXIS 64 | F IsPermuteAxis 65 | 1 iDimLOS 66 | 2 iDimVertical 67 | 3 iDimHorizontal 68 | 69 | #ISALLLINES 70 | F IsAllLines 71 | 72 | #ISDOPPLER 73 | F IsDoppler 74 | -------------------------------------------------------------------------------- /srcPostProc/test-spectrum.ref.gz: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/srcPostProc/test-spectrum.ref.gz -------------------------------------------------------------------------------- /srcPostProc/test-spectrum2-input.tree: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/SWMFsoftware/BATSRUS/a2fc03a3bcc988a7831e62bcea18943d33dec3db/srcPostProc/test-spectrum2-input.tree --------------------------------------------------------------------------------