Commit 4ec5669a authored by Udo Ziegler's avatar Udo Ziegler
Browse files

version 4.1 update

parent cfa1e637
......@@ -48,10 +48,6 @@ for(lev=0; lev<=_C.level; lev++)
for(iy=0; iy<=g->ny; iy++)
for(ix=0; ix<=g->nx; ix++)
{
g->mx[iz][iy][ix]/=g->rho[iz][iy][ix];
g->my[iz][iy][ix]/=g->rho[iz][iy][ix];
g->mz[iz][iy][ix]/=g->rho[iz][iy][ix]*g->hyc[ix]*g->hzyc[iy];
x=xc=g->xc[ix];
y=yc=g->yc[iy];
z=zc=g->zc[iz];
......@@ -64,9 +60,11 @@ for(lev=0; lev<=_C.level; lev++)
dr=sqrt((xc-x0)*(xc-x0)+(yc-y0)*(yc-y0)+(zc-z0)*(zc-z0));
dv=1.-MIN(dt,1.)*(1.-tanh(4.*MAX(2.*dr-1.,0.)));
g->mx[iz][iy][ix]*=dv*g->rho[iz][iy][ix];
g->my[iz][iy][ix]*=dv*g->rho[iz][iy][ix];
g->mz[iz][iy][ix]*=dv*g->rho[iz][iy][ix]*g->hyc[ix]*g->hzyc[iy];
g->mx[iz][iy][ix]*=dv;
g->my[iz][iy][ix]*=dv;
g->mz[iz][iy][ix]*=dv;
}
return;
......
......@@ -6,7 +6,7 @@
*
* Category: Numerics.User
*
* Version: 4
* Version: 4.1
*
* Public functions: bcbUser()
*/
......@@ -30,11 +30,11 @@ void bcbUser(GRD *g, flag_t boundary)
* @note Testproblem: MHD/problem26: Alfven wings propagation
*/
{
int ix,iy,iz,iyp,izp,ib=0,iby,ibz=0;
double sgn,bx0,y;
int ix,iy,iz,iyp,izp;
double bx0,y;
bx0=_C.param[0];
bx0=B_FIELD_SPLITTING ? 0. : _C.param[0];
switch(boundary)
{
......@@ -57,12 +57,14 @@ switch(boundary)
}
for(iz=0; iz<=g->nz+1; iz++)
for(iy=0; iy<=g->ny; iy++) {
y=g->y[iy]+0.5*g->dy;
g->bz[iz][iy][ix]=-bx0*sin(g->z[iz])*(2.*sin(y)*sin(y)
+(sin(g->y[iy+1])*cos(g->y[iy+1])
-sin(g->y[iy])*cos(g->y[iy]))*g->dyi);
}
for(iy=0; iy<=g->ny; iy++)
{
y=g->y[iy]+0.5*g->dy;
g->bz[iz][iy][ix]=-bx0*sin(g->z[iz])*(2.*sin(y)*sin(y)
+(sin(g->y[iy+1])*cos(g->y[iy+1])
-sin(g->y[iy])*cos(g->y[iy]))*g->dyi);
}
}
break;
......
......@@ -31,9 +31,14 @@ void bceUser(GRD *g, flag_t boundary)
* @note Testproblem: MHD/problem26: Alfven wings propagation
*/
{
int ix,iy,iz,ib=0,iby,ibz=0;
double e0,bx,by,bz,dy,pih;
int ix,iy,iz;
pih=0.5/_C.permeability;
e0=_C.param[12]+0.5;
switch(boundary)
{
case X_LOWER: /* LOWER X-FACE --------------------------------------------- */
......@@ -41,22 +46,16 @@ switch(boundary)
for(iz=0; iz<=g->nz; iz++)
for(iy=0; iy<=g->ny; iy++)
{
iby=iy;
ibz=iz;
if(iy<g->iys)
{
iby=2*g->iys-iy-1;
if(_C.idz) ibz=g->izs+(iz-g->izs+_C.dim[2]/2)%_C.dim[2];
}
if(iy>g->iye)
{
iby=2*g->iye-iy+1;
if(_C.idz) ibz=g->izs+(iz-g->izs+_C.dim[2]/2)%_C.dim[2];
if(_C.mf)
{
dy=(g->yc[iy]-g->y[iy])*g->dyi;
bx=g->bx[iz][iy][ix]+(g->bx[iz][iy][ix+1]
-g->bx[iz][iy][ix])*(g->xc[ix]-g->x[ix])*g->dxi;
by=g->by[iz][iy][ix]+(g->by[iz][iy+1][ix]
-g->by[iz][iy][ix])*dy;
bz=0.5*(g->bz[iz][iy][ix]+g->bz[iz+_C.idz][iy][ix]);
}
if(_C.idz) ibz=g->izs+(ibz-g->izs+_C.dim[2]/2)%_C.dim[2];
iby=2*g->iys+_C.dim[1]-iby-1;
g->e[iz][iy][ix]=g->e[ibz][iby][2*g->ixs-1-ix];
g->e[iz][iy][ix]=e0+pih*SPQR(bx,by,bz);
}
break;
......@@ -65,6 +64,16 @@ switch(boundary)
for(iz=0; iz<=g->nz; iz++)
for(iy=0; iy<=g->ny; iy++)
{
if(_C.mf)
{
dy=(g->yc[iy]-g->y[iy])*g->dyi;
bx=g->bx[iz][iy][ix]+(g->bx[iz][iy][ix+1]
-g->bx[iz][iy][ix])*(g->xc[ix]-g->x[ix])*g->dxi;
by=g->by[iz][iy][ix]+(g->by[iz][iy+1][ix]
-g->by[iz][iy][ix])*dy;
bz=0.5*(g->bz[iz][iy][ix]+g->bz[iz+_C.idz][iy][ix]);
}
g->e[iz][iy][ix]=e0+pih*SPQR(bx,by,bz);
}
break;
......@@ -73,8 +82,16 @@ switch(boundary)
for(iz=0; iz<=g->nz; iz++)
for(ix=g->ixs; ix<=g->nx; ix++)
{
if(_C.idz) ib=g->izs+(iz-g->izs+_C.dim[2]/2)%_C.dim[2];
g->e[iz][iy][ix]=g->e[ib][2*g->iys-1-iy][ix];
if(_C.mf)
{
dy=(g->yc[iy]-g->y[iy])*g->dyi;
bx=g->bx[iz][iy][ix]+(g->bx[iz][iy][ix+1]
-g->bx[iz][iy][ix])*(g->xc[ix]-g->x[ix])*g->dxi;
by=g->by[iz][iy][ix]+(g->by[iz][iy+1][ix]
-g->by[iz][iy][ix])*dy;
bz=0.5*(g->bz[iz][iy][ix]+g->bz[iz+_C.idz][iy][ix]);
}
g->e[iz][iy][ix]=e0+pih*SPQR(bx,by,bz);
}
break;
......@@ -83,8 +100,16 @@ switch(boundary)
for(iz=0; iz<=g->nz; iz++)
for(ix=g->ixs; ix<=g->nx; ix++)
{
if(_C.idz) ib=g->izs+(iz-g->izs+_C.dim[2]/2)%_C.dim[2];
g->e[iz][iy][ix]=g->e[ib][2*g->iye+1-iy][ix];
if(_C.mf)
{
dy=(g->yc[iy]-g->y[iy])*g->dyi;
bx=g->bx[iz][iy][ix]+(g->bx[iz][iy][ix+1]
-g->bx[iz][iy][ix])*(g->xc[ix]-g->x[ix])*g->dxi;
by=g->by[iz][iy][ix]+(g->by[iz][iy+1][ix]
-g->by[iz][iy][ix])*dy;
bz=0.5*(g->bz[iz][iy][ix]+g->bz[iz+_C.idz][iy][ix]);
}
g->e[iz][iy][ix]=e0+pih*SPQR(bx,by,bz);
}
break;
......@@ -93,6 +118,16 @@ switch(boundary)
for(iy=0; iy<=g->ny; iy++)
for(ix=0; ix<=g->nx; ix++)
{
if(_C.mf)
{
dy=(g->yc[iy]-g->y[iy])*g->dyi;
bx=g->bx[iz][iy][ix]+(g->bx[iz][iy][ix+1]
-g->bx[iz][iy][ix])*(g->xc[ix]-g->x[ix])*g->dxi;
by=g->by[iz][iy][ix]+(g->by[iz][iy+1][ix]
-g->by[iz][iy][ix])*dy;
bz=0.5*(g->bz[iz][iy][ix]+g->bz[iz+_C.idz][iy][ix]);
}
g->e[iz][iy][ix]=e0+pih*SPQR(bx,by,bz);
}
break;
......@@ -101,6 +136,16 @@ switch(boundary)
for(iy=0; iy<=g->ny; iy++)
for(ix=0; ix<=g->nx; ix++)
{
if(_C.mf)
{
dy=(g->yc[iy]-g->y[iy])*g->dyi;
bx=g->bx[iz][iy][ix]+(g->bx[iz][iy][ix+1]
-g->bx[iz][iy][ix])*(g->xc[ix]-g->x[ix])*g->dxi;
by=g->by[iz][iy][ix]+(g->by[iz][iy+1][ix]
-g->by[iz][iy][ix])*dy;
bz=0.5*(g->bz[iz][iy][ix]+g->bz[iz+_C.idz][iy][ix]);
}
g->e[iz][iy][ix]=e0+pih*SPQR(bx,by,bz);
}
break;
}
......
......@@ -6,7 +6,7 @@
*
* Category: Numerics.User
*
* Version: 4
* Version: 4.1
*
* Public functions: bcetUser()
*/
......@@ -28,80 +28,57 @@ void bcetUser(GRD *g, flag_t boundary)
*
* @author Udo Ziegler, AIP
*
* @note Testproblem: MHD/problem26: Alfveb wings propagation
* @note Testproblem: MHD/problem26: Alfven wings propagation
*/
{
int ix,iy,iz,ib=0,iby,ibz=0;
int ix,iy,iz;
double e0;
e0=_C.param[12];
switch(boundary)
{
case X_LOWER: /* LOWER X-FACE --------------------------------------------- */
for(ix=0; ix<g->ixs; ix++)
for(iz=0; iz<=g->nz; iz++)
for(iy=0; iy<=g->ny; iy++)
{
iby=iy;
ibz=iz;
if(iy<g->iys)
{
iby=2*g->iys-iy-1;
if(_C.idz) ibz=g->izs+(iz-g->izs+_C.dim[2]/2)%_C.dim[2];
}
if(iy>g->iye)
{
iby=2*g->iye-iy+1;
if(_C.idz) ibz=g->izs+(iz-g->izs+_C.dim[2]/2)%_C.dim[2];
}
if(_C.idz) ibz=g->izs+(ibz-g->izs+_C.dim[2]/2)%_C.dim[2];
iby=2*g->iys+_C.dim[1]-iby-1;
g->et[iz][iy][ix]=g->et[ibz][iby][2*g->ixs-1-ix];
}
g->et[iz][iy][ix]=e0;
break;
case X_UPPER: /* UPPER X-FACE --------------------------------------------- */
for(ix=g->ixe+1; ix<=g->nx; ix++)
for(iz=0; iz<=g->nz; iz++)
for(iy=0; iy<=g->ny; iy++)
{
}
g->et[iz][iy][ix]=e0;
break;
case Y_LOWER: /* LOWER Y-FACE --------------------------------------------- */
for(iy=0; iy<g->iys; iy++)
for(iz=0; iz<=g->nz; iz++)
for(ix=g->ixs; ix<=g->nx; ix++)
{
if(_C.idz) ib=g->izs+(iz-g->izs+_C.dim[2]/2)%_C.dim[2];
g->et[iz][iy][ix]=g->et[ib][2*g->iys-1-iy][ix];
}
g->et[iz][iy][ix]=e0;
break;
case Y_UPPER: /* UPPER Y-FACE --------------------------------------------- */
for(iy=g->iye+1; iy<=g->ny; iy++)
for(iz=0; iz<=g->nz; iz++)
for(ix=g->ixs; ix<=g->nx; ix++)
{
if(_C.idz) ib=g->izs+(iz-g->izs+_C.dim[2]/2)%_C.dim[2];
g->et[iz][iy][ix]=g->et[ib][2*g->iye+1-iy][ix];
}
g->et[iz][iy][ix]=e0;
break;
case Z_LOWER: /* LOWER Z-FACE --------------------------------------------- */
for(iz=0; iz<g->izs; iz++)
for(iy=0; iy<=g->ny; iy++)
for(ix=0; ix<=g->nx; ix++)
{
}
g->et[iz][iy][ix]=e0;
break;
case Z_UPPER: /* UPPER Z-FACE --------------------------------------------- */
for(iz=g->ize+1; iz<=g->nz; iz++)
for(iy=0; iy<=g->ny; iy++)
for(ix=0; ix<=g->nx; ix++)
{
}
g->et[iz][iy][ix]=e0;
break;
}
......
......@@ -6,7 +6,7 @@
*
* Category: Numerics.User
*
* Version: 4
* Version: 4.1
*
* Public functions: configUser()
*/
......@@ -26,11 +26,12 @@ void configUser(GRD **gm)
* @author Udo Ziegler, AIP
*
* @note Testproblem: MHD/problem26: Alfven wings propagation
* (cf. Ziegler 12, SIAM J. SCI. COMPUT. 34, C102)
* (cf. Ziegler 12, SIAM J. Sci. Comput. 34, C102)
* - modified to test Boris corr and B-field splitting
*/
{
GRD *g;
double x0,y0,z0,x,y,z,xc,yc,zc,dy,r2,r0,e0,bx0,bx,by0,by,bz0,bz,g1,pih;
double x,y,z,dy,e0,bx0,bx,by0,by,bz0,bz,g1,pih;
int il,ix,iy,iz;
......@@ -40,14 +41,15 @@ g1=_C.gamma-1.;
pih=0.5/_C.permeability;
e0=10./g1;
_C.param[12]=e0;
bx0=_C.param[0];
by0=_C.param[1];
bz0=_C.param[2];
x0=_C.param[3];
y0=_C.param[4];
z0=_C.param[5];
r0=_C.param[6];
if(!B_FIELD_SPLITTING)
{
bx0=_C.param[0];
by0=_C.param[1];
bz0=_C.param[2];
}
else bx0=by0=bz0=0.;
for(il=_C.level; il>=0; il--)
for(g=gm[il]; (g); g=g->next)
......@@ -115,10 +117,6 @@ for(il=_C.level; il>=0; il--)
x=g->xc[ix];
y=g->yc[iy];
z=g->zc[iz];
xc=GRID_SPH ? x*sin(y)*cos(z) : (GRID_CYL ? y*cos(z) : x);
yc=GRID_SPH ? x*sin(y)*sin(z) : (GRID_CYL ? y*sin(z) : y);
zc=GRID_SPH ? x*cos(y) : (GRID_CYL ? x : z);
r2=(xc-x0)*(xc-x0)+(yc-y0)*(yc-y0)+(zc-z0)*(zc-z0);
g->rho[iz][iy][ix]=1.;
if(GRID_SPH)
......@@ -132,17 +130,22 @@ for(il=_C.level; il>=0; il--)
g->my[iz][iy][ix]=1.;
}
g->mz[iz][iy][ix]=0.;
if(_C.energy) g->e[iz][iy][ix]=e0+0.5;
if(_C.mf)
if(_C.energy)
{
dy=(g->yc[iy]-g->y[iy])*g->dyi;
bx=g->bx[iz][iy][ix]+(g->bx[iz][iy][ix+1]
-g->bx[iz][iy][ix])*(g->xc[ix]-g->x[ix])*g->dxi;
by=g->by[iz][iy][ix]+(g->by[iz][iy+1][ix]
-g->by[iz][iy][ix])*dy;
bz=0.5*(g->bz[iz][iy][ix]+g->bz[iz+_C.idz][iy][ix]);
if(_C.energy) g->e[iz][iy][ix]+=pih*SPQR(bx,by,bz);
g->e[iz][iy][ix]=e0+0.5;
if(_C.mf)
{
dy=(g->yc[iy]-g->y[iy])*g->dyi;
bx=g->bx[iz][iy][ix]+(g->bx[iz][iy][ix+1]
-g->bx[iz][iy][ix])*(g->xc[ix]-g->x[ix])*g->dxi;
by=g->by[iz][iy][ix]+(g->by[iz][iy+1][ix]
-g->by[iz][iy][ix])*dy;
bz=0.5*(g->bz[iz][iy][ix]+g->bz[iz+_C.idz][iy][ix]);
g->e[iz][iy][ix]+=pih*SPQR(bx,by,bz);
}
}
}
}
......
......@@ -16,14 +16,14 @@
*/
>SIMULATION I/O ----------------------------------------------------------------------------------
NEW ./NIRLAST >mode:{NEW,RES,MOD,ANA},fname
9000120 7.75e+00 >mod_max,time_max
0020 00500 00001 3.60e+03 >freq_{log,nir,ana},freq_walltime
1000050 7.75e+00 >mod_max,time_max
0020 00050 00001 3.60e+03 >freq_{log,nir,ana},freq_walltime
>GEOMETRY ----------------------------------------------------------------------------------------
CART 0.00e+00 0.00e+00 0.00e+00 >geometry:{CART,CYL,SPH},omega[0-2]
>DOMAIN SETTINGS ---------------------------------------------------------------------------------
-1.00000000e+01 1.00000000e+01 0060 >lo[0],up[0],dim[0]
-1.00000000e+01 1.00000000e+01 0060 >lo[1],up[1],dim[1]
-1.00000000e+01 1.00000000e+01 0060 >lo[2],up[2],dim[2]
-1.00000000e+01 1.00000000e+01 0064 >lo[0],up[0],dim[0]
-1.00000000e+01 1.00000000e+01 0064 >lo[1],up[1],dim[1]
-1.00000000e+01 1.00000000e+01 0064 >lo[2],up[2],dim[2]
SLOCK 1 1 0 >partitioning_type:{SFC,BLOCK # # #}
>BOUNDARY CONDITIONS -----------------------------------------------------------------------------
PPPPPP >bc[0-5]:{I,O,A,M,R,C,P,D,F,U}
......@@ -48,17 +48,17 @@ PPPPPP >bc[0-5]:{I,O,A,M,R,C,P,D,F,U}
STD 4.0e-01 >APdiffusion_{solver,courant}
0.1 >heatloss_max_change
0.1 0.1 >reactions_max_change{X,T}
1. >dt0_reduce:{<=1.}
0.1 >dt0_reduce:{<=1.}
>PHYSICS SPECIFICATIONS --------------------------------------------------------------------------
Y 7.94e+05 >mhd,permeability_rel
N 0.00e+00 >viscosity{,_coeff}
N 0.00e+00 >diffusion{,_coeff}
N 0.00e+00 0.00e+00 0.00e+00 >conduction{,_coeff,_coeff_perp,_sat}
N 0.00e+00 >APdiffusion{,_coeff}
N 0.01 >energy,energy_dual_sw:[0,1]
Y 0.005 >energy,energy_dual_sw:[0,1]
N >force
N >gravity
ISO 1.67e+00 1.00e+00 1.21e-02 >eos,gamma,poly_const,temperature
ADI 1.67e+00 1.00e+00 1.21e-04 >eos,gamma,poly_const,temperature
N >heatloss
0 >tracer
0 >fields
......
......@@ -16,14 +16,14 @@
*/
>SIMULATION I/O ----------------------------------------------------------------------------------
NEW ./NIRLAST >mode:{NEW,RES,MOD,ANA},fname
9000120 7.75e+00 >mod_max,time_max
0020 00500 00001 3.60e+03 >freq_{log,nir,ana},freq_walltime
1000050 7.75e+00 >mod_max,time_max
0020 00050 00001 3.60e+03 >freq_{log,nir,ana},freq_walltime
>GEOMETRY ----------------------------------------------------------------------------------------
CART 0.00e+00 0.00e+00 0.00e+00 >geometry:{CART,CYL,SPH},omega[0-2]
>DOMAIN SETTINGS ---------------------------------------------------------------------------------
-1.00000000e+01 1.00000000e+01 0060 >lo[0],up[0],dim[0]
-1.00000000e+01 1.00000000e+01 0060 >lo[1],up[1],dim[1]
-1.00000000e+01 1.00000000e+01 0060 >lo[2],up[2],dim[2]
-1.00000000e+01 1.00000000e+01 0064 >lo[0],up[0],dim[0]
-1.00000000e+01 1.00000000e+01 0064 >lo[1],up[1],dim[1]
-1.00000000e+01 1.00000000e+01 0064 >lo[2],up[2],dim[2]
SLOCK 1 1 0 >partitioning_type:{SFC,BLOCK # # #}
>BOUNDARY CONDITIONS -----------------------------------------------------------------------------
PPPPPP >bc[0-5]:{I,O,A,M,R,C,P,D,F,U}
......@@ -48,17 +48,17 @@ PPPPPP >bc[0-5]:{I,O,A,M,R,C,P,D,F,U}
STD 4.0e-01 >APdiffusion_{solver,courant}
0.1 >heatloss_max_change
0.1 0.1 >reactions_max_change{X,T}
1. >dt0_reduce:{<=1.}
0.1 >dt0_reduce:{<=1.}
>PHYSICS SPECIFICATIONS --------------------------------------------------------------------------
Y 7.94e+05 >mhd,permeability_rel
N 0.00e+00 >viscosity{,_coeff}
N 0.00e+00 >diffusion{,_coeff}
N 0.00e+00 0.00e+00 0.00e+00 >conduction{,_coeff,_coeff_perp,_sat}
N 0.00e+00 >APdiffusion{,_coeff}
N 0.01 >energy,energy_dual_sw:[0,1]
Y 0.005 >energy,energy_dual_sw:[0,1]
N >force
N >gravity
ISO 1.67e+00 1.00e+00 1.21e-02 >eos,gamma,poly_const,temperature
ADI 1.67e+00 1.00e+00 1.21e-04 >eos,gamma,poly_const,temperature
N >heatloss
0 >tracer
0 >fields
......
/** @file nirvana.par
*
* @brief NIRVANA parameter file
*
* Specification of basic code parameter for a problem.
*
* Copyright (C) 2020 Udo Ziegler
*
* Category: Infrastructure.Kernel
*
* Version: 4
*
* @author Udo Ziegler, AIP
*
* @note Testproblem: MHD/problem26: Alfven wings propagation
* - modified to test Boris correction and B-field splitting
*/
>SIMULATION I/O ----------------------------------------------------------------------------------
NEW ./NIRLAST >mode:{NEW,RES,MOD,ANA},fname
1000050 0.11e+00 >mod_max,time_max
0020 00050 00001 3.60e+03 >freq_{log,nir,ana},freq_walltime
>GEOMETRY ----------------------------------------------------------------------------------------
CART 0.00e+00 0.00e+00 0.00e+00 >geometry:{CART,CYL,SPH},omega[0-2]
>DOMAIN SETTINGS ---------------------------------------------------------------------------------
-1.00000000e+01 1.00000000e+01 0064 >lo[0],up[0],dim[0]
-1.00000000e+01 1.00000000e+01 0064 >lo[1],up[1],dim[1]
-1.00000000e+01 1.00000000e+01 0064 >lo[2],up[2],dim[2]
SLOCK 1 1 0 >partitioning_type:{SFC,BLOCK # # #}
>BOUNDARY CONDITIONS -----------------------------------------------------------------------------
PPPPPP >bc[0-5]:{I,O,A,M,R,C,P,D,F,U}
>MESH REFINEMENT ---------------------------------------------------------------------------------
0 2 >imr,amr
1.0e-01 3.0e-01 -5.0e-01 4.0e-01 2.0e-01 >amr_eps[0-4] / refinement thresholds
0.0e-00 1.0e-01 0.0e+00 3.0e-01 1.0e-02 >rhoref,vref,eref,Bref,Cref
0.6e+00 0.6e+00 >amr_{d1,exp}
2.0e-01 0.0e+00 -0.2e+00 >amr_{Jeans,dJeans,Field}
>USER-SPECIFIC PARAMETERS ------------------------------------------------------------------------
1.00000e+02 0.00000e+00 0.00000e+00 >param[0-2] / bx0,by0,bz0
0.00000e+00 0.00000e+00 0.00000e+00 >param[3-5] / x0,y0,z0
0.15000e+00 0.00000e+00 0.00000e+00 >param[6-8] / r0
0.00000e+00 0.00000e+00 0.00000e+00 >param[9-11] /
0.00000e+00 0.00000e+00 0.00000e+00 >param[12-14] /
0 0 0 0 0 0 0 0 >flag[0-7] /
>SOLVER SPECIFICATIONS ---------------------------------------------------------------------------
RK3 CU CCT 4.0e-01 >mhd_solver_{time,flux,ef},mhd_courant
STD 4.0e-01 >viscosity_{solver,courant}
STD 4.0e-01 >diffusion_{solver,courant}
STD 4.0e-01 >conduction_{solver,courant}
STD 4.0e-01 >APdiffusion_{solver,courant}
0.1 >heatloss_max_change
0.1 0.1 >reactions_max_change{X,T}
0.1 >dt0_reduce:{<=1.}
>PHYSICS SPECIFICATIONS --------------------------------------------------------------------------
Y 7.94e+05 >mhd,permeability_rel
N 0.00e+00 >viscosity{,_coeff}
N 0.00e+00 >diffusion{,_coeff}
N 0.00e+00 0.00e+00 0.00e+00 >conduction{,_coeff,_coeff_perp,_sat}
N 0.00e+00 >APdiffusion{,_coeff}
D 0.005 >energy,energy_dual_sw:[0,1]
N >force
N >gravity
ADI 1.67e+00 1.00e+00 1.21e-04 >eos,gamma,poly_const,temperature
N >heatloss
0 >tracer
0 >fields
N N N 1.0 0.0 >species,reac,energy_reac,mmw,mean_ionization
>>END INPUT --------------------------------------------------------------------------------------
......@@ -17,7 +17,7 @@
>SIMULATION I/O ----------------------------------------------------------------------------------
NEW ./NIRLAST >mode:{NEW,RES,MOD,ANA},fname
900020 7.75e+00 >mod_max,time_max
0050 00500 00001 3.60e+03 >freq_{log,nir,ana},freq_walltime
0010 00050 00001 3.60e+03 >freq_{log,nir,ana},freq_walltime
>GEOMETRY ----------------------------------------------------------------------------------------
SPH 0.00e+00 0.00e+00 0.00e+00 >geometry:{CART,CYL,SPH},omega[0-2]
>DOMAIN SETTINGS ---------------------------------------------------------------------------------
......@@ -29,7 +29,7 @@ SLOCK 4 4 4 >partitioning_type:{SFC,BLOCK # #
UUUUUU >bc[0-5]:{I,O,A,M,R,C,P,D,F,U}
>MESH REFINEMENT ---------------------------------------------------------------------------------
0 2 >imr,amr
8.0e-02 3.5e-01 -5.0e-01 3.0e-01 2.0e-01 >amr_eps[0-4] / refinement thresholds
8.0e-02 3.5e-01 -5.0e-01 -3.0e-01 2.0e-01 >amr_eps[0-4] / refinement thresholds
0.0e-00 1.0e-01 0.0e+00 1.0e-01 1.0e-02 >rhoref,vref,eref,Bref,Cref
0.6e+00 0.0e+00 >amr_{d1,exp}
2.0e-01 0.0e+00 -0.2e+00 >amr_{Jeans,dJeans,Field}
......@@ -41,7 +41,7 @@ UUUUUU >bc[0-5]:{I,O,A,M,R,C,P,D,F,U}
0.00000e+00 0.00000e+00 0.00000e+00 >param[12-14] /
0 0 0 0 0 0 0 0 >flag[0-7] /
>SOLVER SPECIFICATIONS ---------------------------------------------------------------------------
RK3 CU CCT 4.0e-01 >mhd_solver_{time,flux,ef},mhd_courant
RK3 CU CT 4.0e-01 >mhd_solver_{time,flux,ef},mhd_courant
STD 4.0e-01 >viscosity_{solver,courant}
STD 4.0e-01 >diffusion_{solver,courant}
STD 4.0e-01 >conduction_{solver,courant}
......@@ -55,10 +55,10 @@ UUUUUU >bc[0-5]:{I,O,A,M,R,C,P,D,F,U}