#ECHO T #TIMEACCURATE T DoTimeAccurate #CHECKSTOP T DoCheckStop -1000 DnCheckStop 1.0 DtCheckStop #CPUTIMEMAX 84000.0 CpuTimeMax #INCLUDE RESTART.in NameRestartFile #SAVERESTART T DoSaveRestart -10000 DnSaveRestart 10.0 DtSaveRestart #COUPLE2 GM NameCompMaster PC NameCompSlave -1 DnCouple 0.05 DtCouple #COUPLETIME PC F #COUPLETIME GM F #BEGIN_COMP GM --------------------------------------------------------------- #CHECKTIMESTEP T DoCheckTimeStep 2 nCheckTimeStep 1e-4 TimeStepMin #INCLUDE GM/restartIN/restart.H #GRID 1 nRootBlock_D(x_) 1 nRootBlock_D(y_) 1 nRootBlock_D(z_) -224. x1 32. x2 -128. y1 128. y2 -128. z1 128. z2 #GRIDRESOLUTION 1.0 Resolution initial NameArea ! Minimum resolution in the computational domain #GRIDRESOLUTION 8.0 Resolution all NameArea #NONCONSERVATIVE T UseNonConservative #CONSERVATIVECRITERIA 1 nConservCrit parabola TypeConservCrit_I 5.0 xParabolaConserv 15.0 yParabolaConserv MAGNETOSPHERE T body1 2.5 Rbody 3.0 Rcurrents 8.0 Magneto_rho_dim (/ccm) density for fixed BC for rho_BLK 25000.0 Magneto_T_dim (K) temperature for fixed BC for P_BLK ! Pressures below pMinDim are replaced with pMinDim. #MINIMUMPRESSURE 1e-7 pMinDim [nPa] 1e-7 #MINIMUMDENSITY 1e-3 RhoMinDim [amu/cc] #OUTERBOUNDARY outflow TypeBC_I(east_) vary TypeBC_I(west_) float TypeBC_I(south_) float TypeBC_I(north_) float TypeBC_I(bot_) float TypeBC_I(top_) #INNERBOUNDARY ionosphere #PLASMA 1 Mass1 1 Charge1 1 ElectronTemperatureRatio #SCHEME 2 nORDER Rusanov TypeFlux mc3 TypeLimiter 1.5 #TIMESTEPPING 2 nStage 0.50 CflExlp #BORIS T 0.02 IMPLICIT F UsePointImplicit T UsePartImplicit F UseFullImplicit 100.0 CflImpl ImplCFL FIXEDTIMESTEP T UseDtFixed 0.02 DtFixedDim [sec] (read if UseDtFixed is true) #TEST krylov --------------Hall begin--------------- #HALLRESISTIVITY T UseHallResist (rest of parameters read only if true) 16.0 HallFactorMax 1.0 HallCmaxFactor 1.0 !!! #REGION box1 NameRegion box tapered NameHallRegion -20 xMinBox -20 yMinBox -20 zMinBox 15 xMaxBox 20 yMaxBox 20 zMaxBox 1 Taper #REGION sphere1 NameRegion sphere0 tapered StringShape 3 Radius 0.1 Taper REGION polars NameRegion doubleconez0 tapered StringShape 12 Height 2.0 Radius 0.2 Taper #HALLREGION +box1 -sphere1 #SEMIIMPLICIT T UseSemiImplicit resistivity TypeSemiImplicit #SEMIKRYLOV GMRES TypeKrylov (GMRES, BICGSTAB, CG) 0.001 ErrorMaxKrylov 100 MaxMatvecKrylov #SEMIPRECONDITIONER T DoPrecond (rest of parameters are read if true) MBILU TypePrecond (MBILU, BILU, DILU, GS, BLOCKJACOBI, JACOBI, HYPRE) 0.5 GustafssonPar (0 to 1, read for MBILU preconditioner only) --------------Hall end--------------- -------PC BEGIN--------------- PICGHOST 5 #PICUNIT 1.0 xUnitPicSi 2634.0e3 10000.0e3 uUnitPicSi 3000 km/s #PICREGIONROTATE 1 nPicRegion 5.5 xMinCut 7 xLen -8 yMinCut 16 yLen -3 zMinCut 12 zLen 1/25 DxPic 1/25 DyPic 1/25 DzPic 0 xRotate -15 yRotate 0 zRotate ------PC END--------------- #SAVEPLOT 3 nplotfiles y=0 VAR idl plot_strin StringPlot -10000 DnSavePlot 10 DtSavePlot -1. DxSavePlot {MHD} b1x b1y b1z eta jx jy jz dt dtblk cons impl dx hall pic absdivb NameVars {default} z=0 VAR idl plot_strin StringPlot -10000 DnSavePlot 100 DtSavePlot -1. DxSavePlot {MHD} b1x b1y b1z eta jx jy jz dt dtblk cons impl dx hall pic absdivb NameVars {default} cut VAR idl StringPlot ! 3d cell centered Tecplot with MHD data -1 DnSavePlot 100 DtSavePlot 3.5 Coord1MinCut 13 Coord1MaxCut -9 Coord2MinCut 9 Coord2MaxCut -4 Coord3MinCut 14. Coord3MaxCut 0.1 DxSavePlot {MHD} jx jy jz pic NameVars {default} NamePars box MHD idl StringPlot -1 DnSavePlot 5 DtSavePlot GSM TypeCoordPlot 10.0 x0 0.0 y0 0.0 z0 6.0 xSize .1 dX ! only read if xSize /= 0 12 ySize 0.1 dY ! only read if ySize /= 0 12.0 zSize 0.1 dZ ! only read if zSize /= 0 0.0 xAngle [deg] 0.0 yAngle [deg] 0.0 zAngle [deg] #SAVEINITIAL T DoSaveInitial #SAVELOGFILE T DoSaveLogfile RAW StringLogfile 1 DnSaveLogfile -1. DtSaveLogfile #END_COMP GM ----------------------------------------------------------------- #STOP -50000 MaxIter 14400 TimeMax #BEGIN_COMP PC --------------------------------------------------------------- #TIMESTEPPING F UseSWMFDt T UseFixedDt 0.04 CFL #PROCESSORS 4 XLEN 28 YLEN 25 ZLEN #RESTART T DoRrestart #CHECKSTOP 2.0 #ELECTRON -100.0 Charge/mass ratio (qom) #ENERGYCONSERVING T #DISCRETIZATION 0.6 th 0.0 gradRhoRatio 0.0 cDiff 0.1 divC2C #POISSON T doPoissonCorrection 1000000000 nCycle 0.1 20 #RESAMPLING T 0.8 3.0 #DIVE position_all divECleanType 1 nPower 1e-2 divECleanTol 20 divECleanIter 3 nIterNonLinear PARTICLEBC 4 #PARTICLES 5 Particles per cell in X region 1 4 Particles per cell in Y 5 Particles per cell in Z #SAVEIDL 3 3d var real4 planet -1 dn 20 dt 1 dx {fluid} qc divEc rgS0 rgS1 numS0 varName 3d var real4 planet -1 dn 2 dt 2 dx {fluid} qc divEc rgS0 rgS1 numS0 varName y=0 var real4 planet -100 10 0 {fluid} qc divEc rgS0 rgS1 numS0 varName #END_COMP PC -----------------------------------------------------------------