For a two-temperature solar wind simulation, configure SWMF with ./Config.pl -v=Empty,IH/BATSRUS,SC/BATSRUS ./Config.pl -o=SC:u=Awsom,e=Awsom,ng=2,g=8,8,4 ./Config.pl -o=IH:u=Awsom,e=Awsom,ng=2,g=4,4,4 Begin session: 1 #COMPONENTMAP SC 0 -1 1 CompMap IH 0 -1 1 CompMap #ECHO T DoEcho #DESCRIPTION Creating a steady state SC and IH with Threaded-Field-Line Low Corona #TIMEACCURATE F DoTimeAccurate ! CCMC version with the included start time file #STARTTIME 2022 03 28 20 23 0 0.0 #SAVERESTART T DoSaveRestart 20000 DnSaveRestart -1 DtSaveRestart #BEGIN_COMP SC -------------------------------------------------------------- #RESTARTOUTFILE one #GRIDBLOCKALL 200000 MaxBlock #COORDSYSTEM HGR TypeCoordSystem #PLASMA 1.0 IonMass [amu] 1.0 AverageIonCharge [e] 1.0 ElectronTemperatureRatio #FIELDLINETHREAD T UseFieldLineThreads 400 nPointInThreadMax 1.0e-3 DsThreadMin #BODY T UseBody 1.0 rBody 1.5e8 BodyNDim 1.5e6 BodyTDim #BODY F UseBody #HARMONICSGRID 1.0 rMagnetogram 2.5 rSourceSurface F IsLogRadius 90 MaxOrder 30 nR 90 nLon 90 nLat #HARMONICSFILE SC/mf.dat NameHarmonicsFile #CURLB0 T UseCurlB0 2.5 rCurrentFreeB0 T UseB0MomentumFlux #B0SOURCE T UseB0Source T UseDivFullBSource #USERSWITCH +init +ic StringSwitch #POYNTINGFLUX 1.0e6 PoyntingFluxPerBSi [J/m^2/s/T] #CORONALHEATING turbulentcascade TypeCoronalHeating 1.5e5 LperpTimesSqrtBSi 0.0 rMinWaveReflection F UseSurfaceWaveRefl #HEATPARTITIONING stochasticheating TypeHeatPartitioning 0.17 StochasticExponent 0.18 StochasticAmplitude #RADIATIVECOOLING T UseRadCooling #LOOKUPTABLE radcool NameTable load NameCommand SC/Param/RadCoolCorona.dat ascii TypeFile #LOOKUPTABLE TR NameTable load NameCommand SC/Param/TR.dat ascii TypeFile #COARSEAXIS T UseCoarseAxis 3 nCoarseLayer --------------------------------------------------- #GRIDGEOMETRY spherical_lnr ! Use a true sphere, so box is made much larger #GRID 2 nRootBlockX 2 nRootBlockY 2 nRootBlockZ -100.0 xMin 100.0 xMax -100.0 yMin 100.0 yMax -100.0 zMin 100.0 zMax #LIMITRADIUS 1.1 rMin 24.0 rMax ---------------------------------------------------- #AMRREGION InnerShell NameRegion box_gen StringShape 1.1 xMinBox 0.0 yMinBox -70.0 zMinBox 1.7 xMaxBox 360.0 yMaxBox 70.0 zMaxBox #REGION lowlat NameRegion box_gen StringShape 1.0 xMinBox 0.0 yMinBox -70.0 zMinBox 32.0 xMaxBox 360.0 yMaxBox +70.0 zMaxBox #AMRREGION CMEbox box_gen 1.15 80 LongMin 10 LatMin 20 88 LongMax 14 LatMax #AMRCRITERIALEVEL 4 nCriteria Level StringRefine 3 RefineTo 3 CoarsenFrom Level Innershell StringRefine 4 RefineTo 4 CoarsenFrom currentsheet +lowlat StringRefine 0.5 CoarsenLimit 0.5 RefineLimit 4 MaxLevel Level CMEbox StringRefine 5 RefineTo 5 CoarsenFrom #GRIDLEVEL 4 nLevel initial StringShape #DOAMR T DoAmr 20001 DnAmr -1.0 DtAmr F IsStrictAmr #TIMESTEPPING 2 nStage 0.8 CflExpl #MINIMUMPRESSURE 1E-9 pMinDim 1E-9 peMinDim #MINIMUMTEMPERATURE 5e4 TMinDim 5e4 TeMinDim #MINIMUMDENSITY 1e-24 #MINIMUMRADIALSPEED T UseSpeedMin 10 rSpeedMin 250 SpeedMinDim 10 h TauSpeedMinDim #SCHEME 2 nOrder (1 or 2) Sokolov TypeFlux (Roe, Rusanov, Linde, Sokolov mc3 TypeLimiter 1.2 LimiterBeta #LIMITER T UseLogRhoLimiter T UseLogPLimiter F UseRhoRatioLimiter #NONCONSERVATIVE T UseNonConservative #MESSAGEPASS all TypeMessagePass #TVDRESCHANGE T UseTVDResChange #SAVELOGFILE T DoSaveLogfile RAW StringLogfile 1 DnSaveLogfile -1. DtSaveLogfile #OUTERBOUNDARY fieldlinethreads TypeBc1 user for spherical and user_outerbcs float TypeBc2 none TypeBc3 none TypeBc4 none TypeBc5 none TypeBc6 #SAVEPLOT 2 nPlotFiles 3d FUL idl 60000 DnSavePlot -1. DtSavePlot -1 DxSavePlot x=0 VAR idl StringPlot 1000 DnSavePlot -1. DtSavePlot -1.0 DxSavePlot rho ux uy uz ti te bx by bz I01 I02 p pe b1x b1y b1z g rbody #SAVEINITIAL F DoSaveInitial #PLOTTHREADS T DoPlotThreads 10 nGUniform T UseTriangulation T DoTRCorrection F UsePlanarTriangles #ALIGNBANDU T UseChGL 1.1 RSourceChGL 3.5 RMinChGL #END_COMP SC ---------------------------------------------------------------- #COMPONENT IH NameComp F UseComp #STOP 60000 MaxIteration -1.0 tSimulationMax End session: 1 #RUN ########################################################################## Begin session: 2 #BEGIN_COMP SC ---------------------------------------------------------------- #AMR -1 DnRefine #END_COMP SC ---------------------------------------------------------------- #COMPONENT IH NameComp T UseComp #BEGIN_COMP IH ------------------------------------------------------------- #RESTARTOUTFILE one #GRIDBLOCKALL 450000 MaxBlock #COORDSYSTEM HGR TypeCoordSystem #BODY T UseBody 19.0 rBody 3.0E+8 BodyNDim 3.5E+6 BodyTDim #PLASMA 1.0 IonMass [amu] 1.0 AverageIonCharge [e] 1.0 ElectronTemperatureRatio #GRID 1 nRootBlockX 1 nRootBlockY 1 nRootBlockZ -250.0 xMin 250.0 xMax -250.0 yMin 250.0 yMax -250.0 zMin 250.0 zMax #AMRREGION OuterShell NameRegion shell0 StringShape 19.0 RadiusInner 21.55 Radius #AMRREGION coneIH_CME NameRegion conex rotated StringShape 0.0 xPosition 0.0 yPosition 0.0 zPosition 220.0 Height 127.0 Radius 0.0 xrotate -12.0 yrotate 84.0 zrotate #AMRCRITERIALEVEL 4 nCriteria Level StringRefine 4 RefineTo 4 CoarsenFrom Level OuterShell StringRefine 8 RefineTo 8 CoarsenFrom currentsheet StringRefine 0.5 CoarsenLimit 0.5 RefineLimit 8 MaxLevel level +coneIH_CME StringRefine 7 7 #GRIDLEVEL 8 nLevel initial StringRefine #INNERBOUNDARY buffergrid TypeBcInner #HELIOBUFFERGRID 5 nRBuff 180 nLonBuff 90 nLatBuff 18.0 RBuffMin 21.55 RBuffMax #OUTERBOUNDARY float TypeBc1 float TypeBc2 float TypeBc3 float TypeBc4 float TypeBc5 float TypeBc6 #RESCHANGE T UseAccurateResChange #USERSWITCH +init StringSwitch #POYNTINGFLUX 1.0e6 PoyntingFluxPerBSi [J/m^2/s/T] #CORONALHEATING turbulentcascade TypeCoronalHeating 1.5e5 LperpTimesSqrtBSi 0.0 rMinWaveReflection F UseSurfaceWaveRefl #SCHEME 2 nOrder Sokolov TypeFlux mc3 TypeLimiter 1.2 LimiterBeta #NONCONSERVATIVE T UseNonConservative #TIMESTEPPING 2 nStage 0.8 CflExpl #DIVB T UseDivbSource T UseDivbDiffusion F UseProjection F UseConstrainB #SAVELOGFILE T DoSaveLogfile RAW StringLogfile 1 DnSaveLogfile -1. DtSaveLogfile #SAVEINITIAL T #SAVEPLOT 4 nPlotFiles 3d FUL idl StringPlot 5000 DnSavePlot -1. DtSavePlot -1 DxSavePlot x=0 VAR idl StringPlot 1000 DnSavePlot -1. DtSavePlot -1.0 DxSavePlot rho ux uy uz ti te bx by bz I01 I02 p pe b1x b1y b1z g rbody z=0 VAR idl StringPlot 1000 DnSavePlot -1. DtSavePlot -1.0 DxSavePlot rho ux uy uz ti te bx by bz I01 I02 p pe b1x b1y b1z g rbody buf MHD idl 5000 -1.0 #ALIGNBANDU T UseChGL 0 RSourceChGL 0 RMinChGL #END_COMP IH --------------------------------------------------------------- #COUPLE1 SC NameSource IH NameTarget 1 DnCouple -1.0 DtCouple #STOP 60001 MaxIter -1 TimeMax End session: 2 #RUN ########################################################################## Begin session: 3 #COUPLE1 SC NameSource IH NameTarget -1 DnCouple -1.0 DtCouple #COMPONENT SC NameComp F UseComp #BEGIN_COMP IH ---------------------------------------------------------------- #DOAMR T DoAmr 1001 DnAmr -1.0 DtAmr F IsStrictAmr #END_COMP IH ---------------------------------------------------------------- #STOP 64000 MaxIter -1.0 TimeMax End session: 3 #RUN ########################################################################## Begin session: 4 #BEGIN_COMP IH ---------------------------------------------------------------- #AMR -1 DnRefine #END_COMP IH ---------------------------------------------------------------- #STOP 65000 MaxIter -1.0 TimeMax #END ############################################################