Commit 931eeeea authored by Lulu Zhao's avatar Lulu Zhao
Browse files

updates

parent 8c134a67
......@@ -40,7 +40,7 @@ IH NameSource
240 RIhMax
OH NameSource
220.0 RIhMin
1500.0 RIhMax
1000.0 RIhMax
#TIMEACCURATE
T IsTimeAccurate
......@@ -301,13 +301,12 @@ SPHEROMAK TypeCme
#PARTICLELINE
T UseParticles
5000 nFieldLineMax
2000 nFieldLineMax
3000 nParticlePerLine
-1 SpaceStepMin
-1 SpaceStepMax
import InitMode
F UseBRAlignment
0.7 ConBRAngleMax
F UseBUAlignment
0.99 ConBUAngleMax
......@@ -322,7 +321,7 @@ F UsePlanarTriangles
T DoSaveOneTecFile
#SAVEPLOT
4 nPlotFile
3 nPlotFile
x=0 VAR idl StringPlot
-1.0 DnSavePlot
600 DtSavePlot
......@@ -542,7 +541,7 @@ RAW StringLog
100. DtSaveLogfile
#SAVEPLOT
4 nPlotFile
3 nPlotFile
x=0 VAR idl StringPlot
-1.0 DnSavePlot
3600 DtSavePlot
......@@ -601,13 +600,12 @@ F UseChGL
#PARTICLELINE
T UseParticles
5000 nFieldLineMax
2000 nFieldLineMax
3000 nParticlePerLine
-1 SpaceStepMin
-1 SpaceStepMax
import InitMode
F UseBRAlignment
0.1 ConBRAngleMax
F UseBUAlignment
0.99 ConBUAngleMax
......@@ -768,13 +766,12 @@ F UseChGL
#PARTICLELINE
T UseParticles
5000 nFieldLineMax
2000 nFieldLineMax
3000 nParticlePerLine
-1 SpaceStepMin
-1 SpaceStepMax
import InitMode
F UseBRAlignment
0.1 ConBRAngleMax
F UseBUAlignment
0.99 ConBUAngleMax
......@@ -793,8 +790,8 @@ T DoTraceShock
T UseDiffusion
#GRIDNODE
4 nLat
8 nLon
4 nLat
8 nLon
GRID
1.15 RMin
......@@ -805,7 +802,7 @@ GRID
#ORIGIN
2.5 ROrigin
0 LonMin
-80 LatMin
-80 LatMin
360 LonMax
80 LatMax
......@@ -818,7 +815,7 @@ T UseFixedMFPUpstream
0.3 MeanFreePath0 [AU]
#SAVEPLOT
1
2
mh1d flux tec
mh2d flux spread ascii
215.0
......@@ -883,6 +880,33 @@ Begin session: 3
#USERSWITCH
+init +ic -perturb StringSwitch
Saving every hour
#SAVEPLOT
3 nPlotFile
x=0 VAR idl StringPlot
-1 DnSavePlot
3600.0 DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
y=0 VAR idl StringPlot
-1 DnSavePlot
3600. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
z=0 VAR idl StringPlot
-1 DnSavePlot
3600. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
3d VAR tcp StringPlot
-1 DnSavePlot
7200. DtSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
#END_COMP SC -----------------------------------------------------------------
#COUPLE1
......@@ -899,6 +923,36 @@ End session: 3
#RUN ##########################################################################
Begin session: 4
#BEGIN_COMP IH ---------------------------------------------------------------
#SAVEPLOT
3 nPlotFile
x=0 VAR idl StringPlot
-1 DnSavePlot
3600. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
y=0 VAR idl StringPlot
-1 DnSavePlot
3600. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
z=0 VAR idl StringPlot
-1 DnSavePlot
3600. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
3d VAR tcp StringPlot
-1 DnSavePlot
7200. DtSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
#END_COMP IH -----------------------------------------------------------------
#COUPLE1
SC NameSource
IH NameTarget
......@@ -939,4 +993,3 @@ OH NameTarget
End session: 6
#END #####################################################################
For a three-temperature solar wind simulation, configure SWMF with
./Config.pl -v=Empty,SC/BATSRUS,IH/BATSRUS
./Config.pl -o=SC:u=Awsom,e=Awsom,g=6,8,8,ng=3
./Config.pl -o=IH:u=Awsom,e=Awsom,g=8,8,8,ng=3
Default inner boundary is 1.05 Rs
Default rMinWaveReflection is 0.0
Default UseSurfaceWaveRefl is F
For changing the distance of the source surface, rMax in FDIPS.in
needs to change if using FDIPS and rSourceSurface in this file
needs to change if using harmonics.
Begin session: 1
#DESCRIPTION
SWMFSOLAR/Param/PARAM.in.awsomr: create steady state solar corona StringDescription
INCLUDE
RESTART.in
#COMPONENTMAP
SC 0 -1 1 CompMap
IH 0 -1 1 CompMap
#TIMEACCURATE
F IsTimeAccurate
#STARTTIME
2013 iYear
4 iMonth
11 iDay
6 iHour
4 iMinute
0 iSecond
0.0 FracSecond
#SAVERESTART
T DoSaveRestart
10000 DnSaveRestart
-1.0 DtSaveRestart
#TEST
init_axes StringTest
#BEGIN_COMP SC ---------------------------------------------------------------
#GRIDBLOCKALL
120000 MaxBlock
INCLUDE
SC/restartIN/restart.H
#RESTARTOUTFILE
one StringRestartOutFile
#COORDSYSTEM
HGR TypeCoordSystem
#PLASMA
1.0 FluidMass [amu]
1.0 IonCharge [e]
1.0 ElectronTemperatureRatio
THREADEDBC
T UseAlignedVelocity
T DoConvergenceCheck
limited TypeBc
#FIELDLINETHREAD
T UseFieldLineThreads
200 nPointThreadMax
2.0e-3 DsThreadMin
#BODY
T UseBody
1.0 rBody
1.5e8 BodyNDim
1.5e6 BodyTDim
#BODY
F UseBody
LOOKUPTABLE FDIPS^
B0 NameTable
load NameCommand
SC/fdips_bxyz.out NameFile
real4 TypeFile
#HARMONICSFILE
SC/harmonics_adapt.dat NameHarmonicsFile
#HARMONICSGRID
1.0 rMagnetogram
25 rSourceSurface^
T IsLogRadius^
180 MaxOrder
400 nR
180 nLon
90 nLat
FACTORB0
1.0 FactorB0 BrFactor^
#POYNTINGFLUX
1.0e6 PoyntingFluxPerBSi^ [J/m^2/s/T]
#CORONALHEATING
turbulentcascade TypeCoronalHeating
T UseWaveReflection
1.5e5 LperpTimesSqrtBSi^
0.0 rMinWaveReflection^
F UseSurfaceWaveRefl^
#HEATPARTITIONING
stochasticheating TypeHeatPartitioning
0.21 StochasticExponent^
0.18 StochasticAmplitude^
#RADIATIVECOOLING
T UseRadCooling
#LOOKUPTABLE
radcool NameTable
load NameCommand
SC/Param/RadCoolCorona_8.0.dat NameFile
ascii TypeFile
#LOOKUPTABLE
TR NameTable
load NameCommand
SC/Param/TR.dat NameFile
ascii TypeFile
#COARSEAXIS
T UseCoarseAxis
2 nCoarseLayer
turn these off when restarting
---------------------------------------------------
#GRIDGEOMETRY
spherical_lnr TypeGeometry
! Use a true sphere, so box is made much larger
#GRID
2 nRootBlock1
2 nRootBlock2
1 nRootBlock3
-100.0 xMin
100.0 xMax
-100.0 yMin
100.0 yMax
-100.0 zMin
100.0 zMax
#LIMITRADIUS
1.05 rMin rMin_AWSoMR^
24.0 rMax
----------------------------------------------------
#AMRREGION
InnerShell NameRegion
shell0 StringShape
1.05 RadiusInner rMin_AWSoMR^
1.7 Radius
#AMRCRITERIARESOLUTION
3 nRefineCrit
dphi StringRefine
3.0 RefineTo
1.5 CoarsenFrom
dphi Innershell StringRefine
1.5 RefineTo
0.75 CoarsenFrom
currentsheet StringRefine
0.5 CoarsenLimit
0.5 RefineLimit
1.5 MaxResolution
#GRIDRESOLUTION
1.5 Resolution
initial StringShape
#DOAMR
T DoAmr
30001 DnAmr
-1.0 DtAmr
T IsStrictAmr
#TIMESTEPPING
2 nStage
0.8 CflExpl
#SCHEME
2 nOrder (1 or 2)
Linde TypeFlux (Roe, Rusanov, Linde, Sokolov
minmod TypeLimiter
1.2 LimiterBeta
#LIMITER
T UseLogRhoLimiter
T UseLogPLimiter
F UseRhoRatioLimiter
#MINIMUMTEMPERATURE
5e4 TminDim [K] for fluid 1
5e4 TeMinDim [K]
#MINIMUMPRESSURE
1E-9 pMinDim for fluid 1
1E-9 PeMinDim
#MINIMUMRADIALSPEED
T UseSpeedMin
10 rSpeedMin
250 SpeedMinDim
10 h TauSpeedMinDim
#NONCONSERVATIVE
F UseNonConservative^
#RESCHANGE
T UseAccurateReschange
#SAVELOGFILE
T DoSaveLogfile
RAW StringLog
100 DnSaveLogfile
-1. DtSaveLogfile
#OUTERBOUNDARY
fieldlinethreads TypeBc1 rMin
float TypeBc2 rMax
periodic TypeBc3 LonMin
periodic TypeBc4 LonMax
periodic TypeBc5 LatMin
periodic TypeBc6 LatMax
#USERSWITCH
+init +ic StringSwitch
#PLOTTHREADS
T DoPlotThreads
10 nGUniform
T UseTriangulation
T DoTRCorrection
F UsePlanarTriangles
#SAVEPLOT
4 nPlotFile
x=0 VAR idl_ascii StringPlot
5000 DnSavePlot
-1.0 DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
y=0 VAR idl_ascii StringPlot
5000 DnSavePlot
-1. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
z=0 VAR idl_ascii StringPlot
5000 DnSavePlot
-1. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
los ins idl_ascii StringPlot
80000 DnSavePlot
-1.0 DtSavePlot
sta:euvi stb:euvi sdo:aia StringsInstrument
#SATELLITE
3 nSatellite
MHD StringSatellite
-1 DnOutput
-1.0 DtOutput [sec]
SC/TRAJECTORY/earth.dat NameTrajectoryFile
MHD StringSatellite
-1 DnOutput
-1.0 DtOutput [sec]
SC/TRAJECTORY/sta.dat NameTrajectoryFile
MHD StringSatellite
-1 DnOutput
-1.0 DtOutput [sec]
SC/TRAJECTORY/stb.dat NameTrajectoryFile
#LOOKUPTABLE
AiaXrt NameTable
load NameCommand
SC/Param/los_tbl.dat NameFile
ascii TypeFile
#LOOKUPTABLE
euv NameTable
load NameCommand
SC/Param/los_Eit_cor.dat NameFile
ascii TypeFile
#LOOKUPTABLE
EuviA NameTable
load NameCommand
SC/Param/los_EuviA.dat NameFile
ascii TypeFile
#LOOKUPTABLE
EuviB NameTable
load NameCommand
SC/Param/los_EuviB.dat NameFile
ascii TypeFile
#SAVEINITIAL
T DoSaveInitial
#END_COMP SC -----------------------------------------------------------------
#COMPONENT
IH NameComp
F UseComp
#STOP
1000 MaxIter
-1.0 TimeMax
End session: 1
#RUN ##########################################################################
Begin session: 2
#BEGIN_COMP SC ---------------------------------------------------------------
#SCHEME
2 nOrder (1 or 2)
Linde TypeFlux (Roe, Rusanov, Linde, Sokolov
mc3 TypeLimiter
1.2 LimiterBeta
#END_COMP SC -----------------------------------------------------------------
#STOP
70000 MaxIter
-1.0 TimeMax
End session: 2
#RUN ##########################################################################
Begin session: 3
#BEGIN_COMP SC ---------------------------------------------------------------
#MINIMUMRADIALSPEED
F UseSpeedMin
#END_COMP SC -----------------------------------------------------------------
#STOP
80000 MaxIter
-1.0 TimeMax
End session: 3
#RUN ##########################################################################
Begin session: 4
#COMPONENT
IH NameComp
T UseComp
#BEGIN_COMP IH ---------------------------------------------------------------
#GRIDBLOCKALL
3240000 MaxBlock
INCLUDE
IH/restartIN/restart.H
#RESTARTOUTFILE
one StringRestartOutFile
#COORDSYSTEM
HGC TypeCoordSystem
#BODY
T UseBody
19. rBody
3.0E+8 BodyNDim
3.5E+6 BodyTDim
#PLASMA
1.0 FluidMass [amu]
1.0 IonCharge [e]
1.0 ElectronTemperatureRatio
#GRID
1 nRootBlock1
1 nRootBlock2
1 nRootBlock3
-600.0 xMin
600.0 xMax
-600.0 yMin
600.0 yMax
-600.0 zMin
600.0 zMax
#AMRREGION
OuterShell NameRegion
shell0 StringShape
18.0 RadiusInner
21.0 Radius
#AMRREGION
Shell1 NameRegion
shell0 StringShape
21.0 RadiusInner
250.0 Radius
#AMRREGION
Shell2 NameRegion
shell0 StringShape
250.0 RadiusInner
400.0 Radius
#AMRREGION
Shell3 NameRegion
shell0 StringShape
400.0 RadiusInner
600.0 Radius
#AMRCRITERIARESOLUTION
7 nRefineCrit
dx +Outershell StringRefine
0.5 RefineTo
0.25 CoarsenFrom
dx +Shell1 StringRefine
2. RefineTo
1. CoarsenFrom
dx +Shell2 StringRefine
4. RefineTo
2. CoarsenFrom
dx +Shell3 StringRefine
8. RefineTo
4. CoarsenFrom
currentsheet +Shell1 StringRefine
0.5 CoarsenLimit
0.5 RefineLimit
1.0 MaxResolution
currentsheet +Shell2 StringRefine
0.5 CoarsenLimit
0.5 RefineLimit
2.0 MaxResolution
currentsheet +Shell3 StringRefine
0.5 CoarsenLimit
0.5 RefineLimit
4.0 MaxResolution
#GRIDRESOLUTION
0.5 Resolution
initial StringShape
#SAVEINITIAL
F DoSaveInitial
#INNERBOUNDARY
buffergrid TypeBcBody
#BUFFERGRID
5 nRBuff
240 nPhiBuff
120 nThetaBuff
18.0 rBuffMin
20.0 rBuffMax
0.0 PhiBuffMin
360. PhiBuffMax
-90.0 LatBuffMin
90.0 LatBuffMax
#OUTERBOUNDARY
float TypeBc1
float TypeBc2
float TypeBc3
float TypeBc4
float TypeBc5