#-------------------------------------- # Inner Boundary at rScale (in AU) # Outer Boundary at rScale + (flowMag/C) * tSunDwell # *(NUM_SHELLS + (N_PROCS -1)*(NUM_SHELLS-1)) ; # dimension of distribution cm^-3 (u/c)^(-3) # kperxkpar = kappa_perp/kappa_par # BC0 = 0 implies gradient free BC on inner boundary # BC1 = 0 implies gradient free BC on outer boundary # ALPHA = 0.5, 2nd Order CN (fn+1 and fn); # ALPHA = 1, 1st Order Forward Differencing (fn+1); # ALPHA = 0, 1st order Reverse Differencing (fn). # rScale sets the scale of the inner boundary (now 1.0 AU) # all time quantities are normalized by TAU = AU/C # Sun's rotation rate 0.001399 = 2 * PI * TAU /(26 days) # 1 day / 8.3 min = 174 light crossing times # There are 4511 light crossing times in 26 days. # This will be the order of a standard run time # N_EP_STEPS defines the number of Energetic Particle Steps for each # move time-step #--------------------------------------------- rScale=0.9 # defines the simulation inner scale [AU] flowMag=450.0e5 # flow speed in cm/s mhdDensityAu=5.0 # density cm-3 at 1 AU mhdBAu=5.0e-5 # B at 1 AU [G] numEpSteps=10 # The number of substeps for EP computation aziSunStart=0.0 # Longitude of the starting point facing Earth omegaSun=0.001399 # Solar rotation rate lamo=0.4 # Scattering mean free path at 1 GV rig in AU kperxkpar=0.01 # Kper/Kpar useBoundaryRegion=180.0 # boundary region for running the source mass=1.0 # masses in nucleons charge=1.0 # charges in e- obsDtheta=1.0 # cell size in theta (degrees). obsDphi=1.0 # cell size in phi (degrees). d0Dpp=4.0e-6 # defines the stochastic accel rate (s-1) at 5 AU eta2Dpp=0.07 # define the (dB/B0)^2 from stochastic acceleration gamDpp=2.0 # Dpp propto (r)^-gamDpp v1Dpp=40000.0e5 # v [cm/s] for roll-over in Dpp alphaDpp=0.5 # used for implict (alp=1), crank-nicholson (alp=0.5), # or explicit (alp=0) inversion in difusive solver sepFilename=earthSepIn.txt # input filename for seps obsPosFilename=spiceObservers.pos # input filename for observer positions #mhdMapFilename=tim.0001.nc # netcdf filename #mhdMapsList=timList.txt # Loads a list of ENLIL nc files. debugCode=0 # Codes for debug printouts. Different codes correspond # to printouts from different parts of the code. # 5 - MHD (flow.c) interface functions. # 2 - EP transport (energeticParticles.c) functions. dumpFreq=50 # Number of timesteps before the grid is saved to a restart file. procMultiGrid0=1 # First processor in which to expand scalesize procMultiGrid1=5 # Last proc to expand scalesize flareSource=0 # A flag for a solar flare SEP source useCorotBoundary=0 # flag states corotating EP boundary conditions used useMixedBoundary=0 # flag to use corotating EP boundary for targets up to 45 deg from earth. oneauSource=1 # states that we are reading in data at 1 AU useEPBoundary=1 # flag that states EP boundary conditions will be used dumpStreamLines=0 # flag to read the positions of the grid streamlines into a file. dumpMhdAtObservers=0 # flag to read the solar wind parameters at the observers into a file. observerHistory=1 # moves observers in time and records data useParallelDiffusion = 1 # Switch to turn on/off parallel diffusion. useDrift = 1 # Switch to turn on/off particle drifts. useShellDiffusion = 1 # Switch to turn on/off perp. diffusion. useAdiabaticChange = 1 # Switch to turn on/off adiabatic change. useStochastic=1 # Switch to turn on/off stochastic acceleration dumpOnAbort=1 # Switch to turn on/off saving the grid information to a restart file when aborting a run. saveRestartFile=1 # Switch to turn on/off saving the restart file after an EPREM run has completed. simStartTime=2455569.50694444 simStopTime=2455570.49304398 eMin=2.5 eMax=433 tDel=5.0 tSunDwell=10.0 tCounterDwell=10.0 tOberverDwell=10.0 observers=earth,moon,mars # observer names for the output files. multiScale=1.0