Commit c0cc5ef5 authored by Lulu Zhao's avatar Lulu Zhao
Browse files

update PARAM.in for SCIHOHSP

parent cbfe2171
......@@ -27,7 +27,7 @@ IH 0 -1 1 CompMap
#TIMEACCURATE
F IsTimeAccurate
#STARTTIME ^
#STARTTIME
2013 iYear
4 iMonth
11 iDay
......@@ -38,7 +38,7 @@ F IsTimeAccurate
#SAVERESTART
T DoSaveRestart
5000 DnSaveRestart
10000 DnSaveRestart
-1.0 DtSaveRestart
#TEST
......@@ -73,9 +73,6 @@ T UseFieldLineThreads
200 nPointThreadMax
2.0e-3 DsThreadMin
#TEST
radiative_cooling StringTest
#BODY
T UseBody
1.0 rBody
......@@ -91,37 +88,23 @@ load NameCommand
SC/fdips_bxyz.out NameFile
real4 TypeFile
#HARMONICSFILE ^
#HARMONICSFILE
SC/harmonics_adapt.dat NameHarmonicsFile
#HARMONICSGRID ^
#HARMONICSGRID
1.0 rMagnetogram
2.5 rSourceSurface^
F IsLogRadius^
25 rSourceSurface^
T IsLogRadius^
180 MaxOrder
100 nR
400 nR
180 nLon
90 nLat
#CURLB0
T UseCurlB0
2.5 rCurrentFreeB0
T UseB0MomentumFlux
#B0SOURCE
T UseB0Source
T UseDivFullBSource
#ALIGNBANDU
T UseChGL
2.5 RSourceChGL
5.0 RMinChGL
FACTORB0 ^
FACTORB0
1.0 FactorB0 BrFactor^
#POYNTINGFLUX
1e6 PoyntingFluxPerBSi
1.0e6 PoyntingFluxPerBSi^ [J/m^2/s/T]
#CORONALHEATING
turbulentcascade TypeCoronalHeating
......@@ -235,7 +218,7 @@ T UseSpeedMin
10 h TauSpeedMinDim
#NONCONSERVATIVE
T UseNonConservative
F UseNonConservative^
#RESCHANGE
T UseAccurateReschange
......@@ -267,19 +250,19 @@ F UsePlanarTriangles
#SAVEPLOT
4 nPlotFile
x=0 VAR idl_ascii StringPlot
1000 DnSavePlot
5000 DnSavePlot
-1.0 DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
y=0 VAR idl_ascii StringPlot
1000 DnSavePlot
5000 DnSavePlot
-1. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
z=0 VAR idl_ascii StringPlot
1000 DnSavePlot
5000 DnSavePlot
-1. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
......@@ -353,11 +336,6 @@ Linde TypeFlux (Roe, Rusanov, Linde, Sokolov
mc3 TypeLimiter
1.2 LimiterBeta
ALIGNBANDU
T UseChGL
8.0 RSourceChGL
8.0 RMinChGL
#END_COMP SC -----------------------------------------------------------------
#STOP
......@@ -390,7 +368,7 @@ T UseComp
#BEGIN_COMP IH ---------------------------------------------------------------
#GRIDBLOCKALL
160000 MaxBlock
3240000 MaxBlock
INCLUDE
IH/restartIN/restart.H
......@@ -416,12 +394,12 @@ T UseBody
1 nRootBlock1
1 nRootBlock2
1 nRootBlock3
-250.0 xMin
250.0 xMax
-250.0 yMin
250.0 yMax
-250.0 zMin
250.0 zMax
-600.0 xMin
600.0 xMax
-600.0 yMin
600.0 yMax
-600.0 zMin
600.0 zMax
#AMRREGION
OuterShell NameRegion
......@@ -429,18 +407,50 @@ 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
3 nRefineCrit
dx StringRefine
8. RefineTo
4. CoarsenFrom
7 nRefineCrit
dx +Outershell StringRefine
0.5 RefineTo
0.25 CoarsenFrom
currentsheet StringRefine
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
......@@ -503,7 +513,7 @@ It might be important to prevent code crashing due to CMEs with strong shocks.
1E-14 PeMinDim
#NONCONSERVATIVE
T UseNonConservative
F UseNonConservative^
#TIMESTEPPING
2 nStage
......@@ -524,22 +534,22 @@ RAW StringLog
#SAVEPLOT
3 nPlotFile
x=0 VAR idl_ascii StringPlot
1000 DnSavePlot
1 DnSavePlot
-1. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{MHD} b1x b1y b1z te ti dx NameVars
{default} NamePars
y=0 VAR idl_ascii StringPlot
1000 DnSavePlot
1 DnSavePlot
-1. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{MHD} b1x b1y b1z te ti dx NameVars
{default} NamePars
z=0 VAR idl_ascii StringPlot
1000 DnSavePlot
1 DnSavePlot
-1. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{MHD} b1x b1y b1z te ti dx NameVars
{default} NamePars
3d VAR tcp StringPlot
5000 DnSavePlot
......@@ -574,11 +584,6 @@ SC/TRAJECTORY/stb.dat NameTrajectoryFile
1 h DtTraj
rho ux uy uz bx by bz p pe NameSatelliteVars
#ALIGNBANDU
T UseChGL
0.0 RSourceChGL
0.0 RMinChGL
#END_COMP IH -----------------------------------------------------------------
#COUPLE1
......@@ -613,6 +618,32 @@ T DoAmr
-1.0 DtAmr
T IsStrictAmr
#SAVEPLOT
3 nPlotFile
x=0 VAR idl_ascii StringPlot
1000 DnSavePlot
-1. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
y=0 VAR idl_ascii StringPlot
1000 DnSavePlot
-1. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
z=0 VAR idl_ascii StringPlot
1000 DnSavePlot
-1. DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
3d VAR tcp StringPlot
5000 DnSavePlot
-1. DtSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
#END_COMP IH -----------------------------------------------------------------
#STOP
......
For a three-temperature solar wind simulation, configure SWMF with
./Config.pl -v=OH/BATSRUS,SP/MFLAMPA;
./Config.pl -o=SC:u=Awsom,e=AwsomChGL,nG=2,g=6,8,8;
./Config.pl -o=IH:u=Awsom,e=AwsomChGL,nG=2,g=8,8,8;
./Config.pl -o=OH:u=Awsom,e=AwsomChGL,nG=2,g=4,4,4;
Default inner boundary is 1.1 Rs
Default rMinWaveReflection is 0.0
Default UseSurfaceWaveRefl is F
ChGL SC setup: 1.1Rs, 3.5Rs
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.SCIHOHSP.CME: create steady state solar corona StringDescription
#INCLUDE
RESTART.in
#COMPONENTMAP
SC 0 -1 1 CompMap
IH 0 -1 1 CompMap
OH 0 -1 1 CompMap
#TIMEACCURATE
T IsTimeAccurate
#STARTTIME ^
2013 iYear
4 iMonth
11 iDay
7 iHour
24 iMinute
06 iSecond
0.0 FracSecond
#SAVERESTART
T DoSaveRestart
-1.0 DnSaveRestart
1 h DtSaveRestart
ROTATEHGR
-1. dLongitudeHgr
ROTATEHGI
-1. dLongitudeHgi
#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
TEST
radiative_cooling StringTest
#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
2.5 rSourceSurface^
F IsLogRadius^
90 MaxOrder
100 nR
180 nLon
90 nLat
#CURLB0
T UseCurlB0
2.5 rCurrentFreeB0
T UseB0MomentumFlux
#B0SOURCE
T UseB0Source
T UseDivFullBSource
ALIGNBANDU
T UseChGL
1.1 RSourceChGL
3.5 RMinChGL
FACTORB0 ^
1.0 FactorB0 BrFactor^
#POYNTINGFLUX
0.3e6 PoyntingFluxPerBSi
#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.1 rMin rMin_AWSoMR^
24.0 rMax
----------------------------------------------------
#AMRREGION
InnerShell NameRegion
shell0 StringShape
1.1 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
F 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
mc3 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
-1. DnSaveLogfile
120 DtSaveLogfile
#OUTERBOUNDARY
fieldlinethreads TypeBc1 rMin
float TypeBc2 rMax
periodic TypeBc3 LonMin
periodic TypeBc4 LonMax
periodic TypeBc5 LatMin
periodic TypeBc6 LatMax
add the perturb switch for CME
#USERSWITCH
+init +ic +perturb StringSwitch
flux-rope parameter
#CME
T UseCme
T DoAddFluxRope
71.50 LongitudeCme
9.50 LatitudeCme
39.76 OrientationCme
SPHEROMAK TypeCme
-12.04 BStrength
0.30 Radius
0.60 Stretch
0.59 ApexHeight
710 uCme Expected CME speed in km/s
#HELIOUPDATEB0
600 DtUodateB0 ! Time intervals to update the magnetic field in the threaded gap
#PLOTTHREADS
T DoPlotThreads
10 nGUniform
T UseTriangulation
T DoTRCorrection
F UsePlanarTriangles
#SAVETECPLOT
T DoSaveOneTecFile
#SAVEPLOT
4 nPlotFile
x=0 VAR idl StringPlot
-1.0 DnSavePlot
600 DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti dx NameVars
{default} NamePars
y=0 VAR idl StringPlot
-1.0 DnSavePlot
600 DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti dx NameVars
{default} NamePars
z=0 VAR idl StringPlot
-1.0 DnSavePlot
600 DtSavePlot
-1.0 DxSavePlot
{MHD} b1x b1y b1z te ti dx NameVars
{default} NamePars
3d VAR tcp StringPlot
-1 DnSavePlot
600. DtSavePlot
{MHD} b1x b1y b1z te ti NameVars
{default} NamePars
los ins idl StringPlot
-1.0 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
#COMPONENT
OH NameComp
F UseComp