Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

[svn:einsteintoolkit] GRHydro_InitData/trunk/src/ (Rev. 151)
by schnetter＠cct.lsu.edu 09 Aug '12

09 Aug '12

User: eschnett Date: 2012/08/09 03:02 PM Modified: /trunk/src/ GRHydro_BondiM_new.F90, GRHydro_Bondi_new.F90 Log: Define Pi with double precision File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_BondiM_new.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_BondiM_new.F90 2012-08-09 06:27:00 UTC (rev 150) +++ trunk/src/GRHydro_BondiM_new.F90 2012-08-09 20:02:21 UTC (rev 151) @@ -31,7 +31,7 @@ #include "cctk_Functions.h" #include "GRHydro_Macros.h" -# define M_PI 3.14159265358979323846 /* pi */ +# define M_PI 3.14159265358979323846d0 /* pi */ !!$Newton-Raphson parameters: File [modified]: GRHydro_Bondi_new.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_Bondi_new.F90 2012-08-09 06:27:00 UTC (rev 150) +++ trunk/src/GRHydro_Bondi_new.F90 2012-08-09 20:02:21 UTC (rev 151) @@ -31,7 +31,7 @@ #include "cctk_Functions.h" #include "GRHydro_Macros.h" -# define M_PI 3.14159265358979323846 /* pi */ +# define M_PI 3.14159265358979323846d0 /* pi */ !!$Newton-Raphson parameters:

1 0

[svn:einsteintoolkit] GRHydro_InitData/trunk/ (Rev. 150)
by rhaas＠tapir.caltech.edu 09 Aug '12

09 Aug '12

User: rhaas Date: 2012/08/09 01:27 AM Modified: /trunk/ param.ccl /trunk/src/ GRHydro_AdvectedLoopM.F90, GRHydro_AlfvenWaveM.F90 Log: Fix to oblique case for Advected loop, new implementation of Alfven wave test including used-defined pressure (to alter wavespeed) and oblique case From: Josh Faber <jfaberrit(a)jfaber3.local> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_AdvectedLoopM.F90 Delta lines: +3 -3 =================================================================== --- trunk/src/GRHydro_AdvectedLoopM.F90 2012-08-09 06:26:57 UTC (rev 149) +++ trunk/src/GRHydro_AdvectedLoopM.F90 2012-08-09 06:27:00 UTC (rev 150) @@ -133,9 +133,9 @@ dz = CCTK_DELTA_SPACE(3) !!$ Note that the 3D test wasn't deviced to be used with AMR! - range_x = (cctk_gsh(1)-1)*dx - range_y = (cctk_gsh(2)-1)*dy - range_z = (cctk_gsh(3)-1)*dz + range_x = (cctk_gsh(1)-2*cctk_nghostzones(1))*dx + range_y = (cctk_gsh(2)-2*cctk_nghostzones(2))*dy + range_z = (cctk_gsh(3)-2*cctk_nghostzones(3))*dz range_d = sqrt(range_z**2+range_x**2) diaglength = range_x*range_z/range_d File [modified]: GRHydro_AlfvenWaveM.F90 Delta lines: +53 -30 =================================================================== --- trunk/src/GRHydro_AlfvenWaveM.F90 2012-08-09 06:26:57 UTC (rev 149) +++ trunk/src/GRHydro_AlfvenWaveM.F90 2012-08-09 06:27:00 UTC (rev 150) @@ -65,6 +65,7 @@ CCTK_REAL :: dx,dy,dz CCTK_REAL :: range_x,range_y,range_z,range_d CCTK_REAL :: cos_theta, sin_theta + CCTK_REAL :: diaglength,xnew,vparallel,vperp,Bparallel,Bperp CCTK_REAL :: Bvecx_d, Bvecz_d CCTK_REAL :: velx_d, velz_d CCTK_REAL :: det @@ -77,10 +78,13 @@ !!$pressure, density, B^x, specific internal energy and enthalpy: rhoval = 1.0d0 - pressval = 1.0d0 + pressval = alfvenwave_pressure Bxval = 1.0d0 epsval = pressval/(gam-1.0d0)/rhoval hval = 1.0d0 + epsval + pressval/rhoval + +!!$ DZ: rho=P=Bxval=AA = 1 +!!$ Using DZ parameters, epsval=1.5, hval=3.5 !!$ Alfven Wave Amplitude: AA=1.0d0 @@ -91,6 +95,14 @@ t3 = 0.5d0*(1.0d0+sqrt(1.0d0-t2**2)) valf = sqrt(Bxval**2/t1/t3) + write(*,*)'Alfven velocity:',valf + +!!$ Using DZ parameters with P=1.0, we have: +!!$ t1=5.5 +!!$ t2=4/11 +!!$ t3=0.96577 +!!$ valf=0.43389 + !!$ Vx value: vxval=0.0d0 @@ -103,45 +115,56 @@ dz = CCTK_DELTA_SPACE(3) !!$ Note that the 3D test wasn't deviced to be used with AMR! - range_x = (cctk_gsh(1)-1)*dx - range_y = (cctk_gsh(2)-1)*dy - range_z = (cctk_gsh(3)-1)*dz + range_x = (cctk_gsh(1)-2*cctk_nghostzones(1))*dx + range_y = (cctk_gsh(2)-2*cctk_nghostzones(2))*dy + range_z = (cctk_gsh(3)-2*cctk_nghostzones(3))*dz - range_d = sqrt(range_z**2+range_x**2) - - cos_theta = range_z/range_d - sin_theta = range_x/range_d - !!$ Alfven wave number - wnbr = 2.0d0*pi/range_x do i=1,nx do j=1,ny do k=1,nz - - rho(i,j,k)=rhoval + + rho(i,j,k)=rhoval press(i,j,k)=pressval eps(i,j,k)=epsval - velx(i,j,k)=vxval - vely(i,j,k)=-valf*AA*cos(wnbr*x(i,j,k)) - velz(i,j,k)=-valf*AA*sin(wnbr*x(i,j,k)) - Bvecx(i,j,k)=Bxval - Bvecy(i,j,k)=AA*Bxval*cos(wnbr*x(i,j,k)) - Bvecz(i,j,k)=AA*Bxval*sin(wnbr*x(i,j,k)) + + if (CCTK_EQUALS(alfvenwave_type,"1D")) then + wnbr = 2.0d0*pi/range_x -!!$ if (CCTK_EQUALS(advectedloop_type,"3D")) then -!!$ Bvecx_d = cos_theta*Bvecx(i,j,k)+sin_theta*Bvecz(i,j,k) -!!$ Bvecz_d = cos_theta*Bvecz(i,j,k)-sin_theta*Bvecx(i,j,k) -!!$ velx_d = cos_theta*velx(i,j,k)+sin_theta*velz(i,j,k) -!!$ velz_d = cos_theta*velz(i,j,k)-sin_theta*velx(i,j,k) -!!$ -!!$ Bvecx(i,j,k) = Bvecx_d -!!$ Bvecz(i,j,k) = Bvecz_d -!!$ velx(i,j,k) = velx_d -!!$ velz(i,j,k) = velz_d -!!$ end if - + velx(i,j,k)=vxval + vely(i,j,k)=-valf*AA*cos(wnbr*x(i,j,k)) + velz(i,j,k)=-valf*AA*sin(wnbr*x(i,j,k)) + Bvecx(i,j,k)=Bxval + Bvecy(i,j,k)=AA*Bxval*cos(wnbr*x(i,j,k)) + Bvecz(i,j,k)=AA*Bxval*sin(wnbr*x(i,j,k)) + + else if (CCTK_EQUALS(alfvenwave_type,"2D")) then + + diaglength=range_x*range_y/range_d + range_d = sqrt(range_x**2+range_y**2) + cos_theta = range_y/range_d + sin_theta = range_x/range_d + wnbr = 2.0d0*pi/diaglength + + xnew = cos_theta*x(i,j,k)+sin_theta*y(i,j,k) + + vparallel=vxval + vperp=-valf*AA*cos(wnbr*xnew) + velx(i,j,k)=vparallel*cos_theta-vperp*sin_theta + vely(i,j,k)=vparallel*sin_theta+vperp*cos_theta + velz(i,j,k)=-valf*AA*sin(wnbr*xnew) + Bparallel=Bxval + Bperp=AA*Bxval*cos(wnbr*xnew) + Bvecx(i,j,k)=Bparallel*cos_theta-Bperp*sin_theta + Bvecy(i,j,k)=Bparallel*sin_theta+Bperp*cos_theta + Bvecz(i,j,k)=AA*Bxval*sin(wnbr*xnew) + + else + call CCTK_WARN(0,"Alfven wave case not recognized!") + end if + det=SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),gyy(i,j,k),gyz(i,j,k),gzz(i,j,k)) if (CCTK_EQUALS(GRHydro_eos_type,"Polytype")) then Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +14 -0 =================================================================== --- trunk/param.ccl 2012-08-09 06:26:57 UTC (rev 149) +++ trunk/param.ccl 2012-08-09 06:27:00 UTC (rev 150) @@ -191,6 +191,20 @@ "Numeric" :: "Finite difference approximation of the derivatives applied" } "Exact" +#################################3 +# Alfven Wave test +################################3 +KEYWORD alfvenwave_type "1-dimensional or 2-dimensional?" +{ + "1D" :: "1-dimensional" + "2D" :: "2-dimensional (in x-y plane)" +} "1D" + +CCTK_REAL alfvenwave_pressure "P_gas for the Alfven wave" +{ + (0:* :: "positive" +} 1.0 + ################################################################################## # BONDI PARAMETERS: (black hole mass specified by parameters from the "Exact" thorn) ##################################################################################

1 0

[svn:einsteintoolkit] GRHydro_InitData/trunk/ (Rev. 149)
by rhaas＠tapir.caltech.edu 09 Aug '12

09 Aug '12

User: rhaas Date: 2012/08/09 01:26 AM Modified: /trunk/ param.ccl /trunk/src/ GRHydro_MonopoleM.F90 Log: Slight tweaks to the monopole routine options (patch by Josh) From: Bruno Coutinho Mundim <bcmsma(a)astro.rit.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_MonopoleM.F90 Delta lines: +26 -5 =================================================================== --- trunk/src/GRHydro_MonopoleM.F90 2012-08-09 06:26:55 UTC (rev 148) +++ trunk/src/GRHydro_MonopoleM.F90 2012-08-09 06:26:57 UTC (rev 149) @@ -71,9 +71,14 @@ vely(i,j,k) = 0.0 velz(i,j,k) = 0.0 eps(i,j,k) = 0.1 + + Bvecx(i,j,k)=0.0 + Bvecy(i,j,k)=0.0 + Bvecz(i,j,k)=0.0 + if(CCTK_EQUALS(monopole_type,"Point")) then if(i.eq.nx/2+1.and.j.eq.ny/2+1.and.k.eq.nz/2+1) then - Bvecx(i,j,k)=0.1 + Bvecx(i,j,k)=Monopole_point_Bx else Bvecx(i,j,k)=0.0 endif @@ -93,7 +98,7 @@ else Bvecx(i,j,k) = 0.0 endif - if(mod(i+j,2).eq.0)Bvecx(i,j,k)=-1.0*Bvecx(i,j,k) + if(mod(i+j+k,2).eq.0)Bvecx(i,j,k)=-1.0*Bvecx(i,j,k) else if(CCTK_EQUALS(monopole_type,"2dalt")) then rr2=x(i,j,k)**2+y(i,j,k)**2+z(i,j,k)**2 rg2=R_Gauss*R_Gauss @@ -104,17 +109,33 @@ Bvecx(i,j,k) = 0.0 Bvecy(i,j,k) = 0.0 endif + if(mod(i+j+k,2).eq.0)then + Bvecx(i,j,k)=-1.0*Bvecx(i,j,k) +!!$ Only vary one component, for different character +!!$ Bvecy(i,j,k)=-1.0*Bvecy(i,j,k) + endif + else if(CCTK_EQUALS(monopole_type,"3dalt")) then + rr2=x(i,j,k)**2+y(i,j,k)**2+z(i,j,k)**2 + rg2=R_Gauss*R_Gauss + if(rr2.lt.rg2) then + Bvecx(i,j,k) = exp(-1.0*rr2/rg2)-1.0/exp(1.0) + Bvecy(i,j,k) = exp(-1.0*rr2/rg2)-1.0/exp(1.0) + Bvecz(i,j,k) = exp(-1.0*rr2/rg2)-1.0/exp(1.0) + else + Bvecx(i,j,k) = 0.0 + Bvecy(i,j,k) = 0.0 + Bvecz(i,j,k) = 0.0 + endif +!!$ Different spatial pattern! if(mod(i+j,2).eq.0)then Bvecx(i,j,k)=-1.0*Bvecx(i,j,k) Bvecy(i,j,k)=-1.0*Bvecy(i,j,k) + Bvecz(i,j,k)=-1.0*Bvecz(i,j,k) endif else call CCTK_WARN(0,"Unrecognized monopole type!!!") endif - Bvecy(i,j,k)=0.0 - Bvecz(i,j,k)=0.0 - det=SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),gyy(i,j,k),gyz(i,j,k),gzz(i,j,k)) if (CCTK_EQUALS(GRHydro_eos_type,"Polytype")) then Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +6 -0 =================================================================== --- trunk/param.ccl 2012-08-09 06:26:55 UTC (rev 148) +++ trunk/param.ccl 2012-08-09 06:26:57 UTC (rev 149) @@ -147,6 +147,7 @@ "Gauss" :: "Gaussian w/radius R_Gauss" "1dalt" :: "1-d alternating" "2dalt" :: "2-d alternating" + "3dalt" :: "3-d alternating" } "Point" CCTK_REAL R_Gauss "Radius for a Gaussian monopole" @@ -154,6 +155,11 @@ 0:* :: "Any positive number" } 1.0 +CCTK_REAL Monopole_point_Bx "Pointlike Monopole Bx value" +{ + *:* :: "Any number" +} 1.0 + # For cylindrical explosion test: CCTK_REAL cyl_r_inner "Inner Radius" {

1 0

[svn:einsteintoolkit] GRHydro_InitData/trunk/ (Rev. 148)
by rhaas＠tapir.caltech.edu 09 Aug '12

09 Aug '12

User: rhaas Date: 2012/08/09 01:26 AM Modified: /trunk/ param.ccl /trunk/src/ GRHydro_Bondi_new.F90 Log: Added option to smooth metric inside horizon Signed-off-by: Bruno Coutinho Mundim <bcmsma(a)astro.rit.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Bondi_new.F90 Delta lines: +11 -2 =================================================================== --- trunk/src/GRHydro_Bondi_new.F90 2012-08-09 06:26:54 UTC (rev 147) +++ trunk/src/GRHydro_Bondi_new.F90 2012-08-09 06:26:55 UTC (rev 148) @@ -55,7 +55,7 @@ CCTK_REAL :: M, Msq, Mdot, rs, gam, rmin_bondi, rmax_bondi, cs_sq,cs,vs_sq,vs,rhos,gtemp,hs, Kval, Qdot CCTK_REAL :: logrmin,dlogr,rhotmp,utmp,vtmp,rspher CCTK_REAL :: r_bondi(N_points), logr_bondi(N_points), rho_bondi(N_points), u_bondi(N_points), v_bondi(N_points) - CCTK_REAL :: drhodr, det, rhocheck, rhocheck2, riso, rnew, rsch, ucheck + CCTK_REAL :: drhodr, det, rhocheck, rhocheck2, riso, rnew, rsch, ucheck, psinew CCTK_REAL :: uiso, uisocheck, vcheck, ucheck2, vcheck2, xhat,yhat, zhat, xp, yp, zp !!$set_bondi_parameters @@ -169,7 +169,16 @@ yhat = yp/riso zhat = zp/riso - if(riso < 1.0e-7) then + if(use_smooth_puncture_data.ne.0 .and. riso.lt.M) then + psinew=2.875 - 5*(riso/M)**2 + 6*(riso/M)**4 + gxx(i,j,k) = psinew**4 + gyy(i,j,k) = psinew**4 + gzz(i,j,k) = psinew**4 + gxy(i,j,k) = 0.0 + gxz(i,j,k) = 0.0 + gyz(i,j,k) = 0.0 + + else if(riso < 1.0e-7) then gxx(i,j,k) = 1.0e4 gyy(i,j,k) = 1.0e4 gzz(i,j,k) = 1.0e4 Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +4 -0 =================================================================== --- trunk/param.ccl 2012-08-09 06:26:54 UTC (rev 147) +++ trunk/param.ccl 2012-08-09 06:26:55 UTC (rev 148) @@ -97,6 +97,10 @@ 0.0:* :: "Anything positive" } 1.0 +BOOLEAN use_smooth_puncture_data "Smooth puncture solution inside R=M using 4th order matching polynomial" +{ +} "no" + BOOLEAN change_shock_direction "Change the shock direction" { } "no"

1 0

[svn:einsteintoolkit] GRHydro_InitData/trunk/par/ (Rev. 147)
by rhaas＠tapir.caltech.edu 09 Aug '12

09 Aug '12

User: rhaas Date: 2012/08/09 01:26 AM Added: /trunk/par/ advectedloop2D_hi.par, advectedloop2D_vz_hi.par Log: high resolution advectedloop2D par files. From: Bruno Coutinho Mundim <bcmsma(a)astro.rit.edu> File Changes: Directory: /trunk/par/ ====================== File [added]: advectedloop2D_hi.par Delta lines: +398 -0 =================================================================== --- trunk/par/advectedloop2D_hi.par (rev 0) +++ trunk/par/advectedloop2D_hi.par 2012-08-09 06:26:54 UTC (rev 147) @@ -0,0 +1,398 @@ +#============================================================================== +# Advected magnetic field loop (2d version): GRHydro_InitData +#============================================================================== + +#------------------------------------------------------------------------------ +# Cactus parameters: +#------------------------------------------------------------------------------ +Cactus::cctk_run_title = "Advected magnetic field loop (2d version)" +Cactus::cctk_full_warnings = "yes" + +#Cactus::terminate = "never" +Cactus::terminate = "time" +Cactus::cctk_final_time = 25.0 +#Cactus::terminate = "iteration" +#Cactus::cctk_itlast = 0 +#Cactus::cctk_itlast = 256 +#Cactus::terminate="runtime" +#Cactus::max_runtime=28800 # 46 days +#Cactus::max_runtime=1410 # 23.5 hours +#Cactus::max_runtime=1380 # 23 hours +#Cactus::max_runtime=690 # 11.5 hours +#Cactus::max_runtime=450 # 7.5 hours +#Cactus::max_runtime=210 # 3.5 hours + +#------------------------------------------------------------------------------ +# Activate all necessary thorns: +#------------------------------------------------------------------------------ + +ActiveThorns = "Boundary CartGrid3D CoordBase Fortran InitBase IOUtil LocalReduce SymBase Time" +ActiveThorns = "AEILocalInterp LocalInterp" +ActiveThorns = "MoL Periodic Slab SpaceMask" +ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetRegrid2 CarpetSlab LoopControl" +#ActiveThorns = "Formaline NaNChecker TerminationTrigger TimerReport" +ActiveThorns = "NaNChecker TerminationTrigger TimerReport" +ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal" +ActiveThorns = "Constants TmunuBase HydroBase " +ActiveThorns = "EOS_Omni" +ActiveThorns = "GRHydro GRHydro_InitData" +ActiveThorns = "HDF5" +ActiveThorns = "SystemStatistics" + +#------------------------------------------------------------------------------ +# Diagnostic parameters: +#------------------------------------------------------------------------------ +#AHFinderDirect::verbose_level = "physics details" + +Carpet::output_timers_every = 1024 +Carpet::output_timer_tree_every = 256 +Carpet::output_initialise_timer_tree = "yes" +Carpet::storage_verbose = "no" +#Carpet::timers_verbose = "no" +#Carpet::verbose = "yes" # only for Initial Data +Carpet::verbose = "no" +Carpet::veryverbose = "no" +Carpet::grid_structure_filename = "carpet-grid-structure.asc" +Carpet::grid_coordinates_filename = "carpet-grid-coordinates.asc" + +CarpetLib::output_bboxes = "no" +CarpetLib::print_memstats_every = 1024 +CarpetLib::memstat_file = "carpetlib-memory-statistics.asc" +#CarpetLib::print_timestats_every = 1024 +#CarpetLib::timestat_file = "carpetlib-timing-statistics.asc" + +CarpetReduce::verbose = "no" +CarpetRegrid2::verbose = "yes" +#CarpetRegrid2::veryverbose = "yes" + +#NaNChecker::verbose = "all" + +TimerReport::out_every = 1024 +TimerReport::out_filename = "TimerReport" +#TimerReport::output_all_timers = "yes" +#TimerReport::output_all_timers_together = "yes" +#TimerReport::output_all_timers_readable = "yes" +#TimerReport::before_checkpoint = "yes" +TimerReport::n_top_timers = 40 + +#------------------------------------------------------------------------------ +# Utility parameters: +#------------------------------------------------------------------------------ + +#NaNChecker::check_every = 16 +NaNChecker::check_every = 1 +NaNChecker::check_vars = " + HydroBase::Bvec + HydroBase::rho + HydroBase::eps + HydroBase::press + HydroBase::vel + HydroBase::w_lorentz +" +NaNChecker::action_if_found = "terminate" +#NaNChecker::action_if_found = "abort" +#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort" +#NaNChecker::check_after=0 + +#TerminationTrigger::max_walltime = 0.0 # hours +TerminationTrigger::max_walltime = 12.0 # hours +TerminationTrigger::on_remaining_walltime = 10.0 # minutes +TerminationTrigger::termination_from_file = "yes" +TerminationTrigger::create_termination_file = "yes" +TerminationTrigger::termination_file = "../TERMINATE" + + +#------------------------------------------------------------------------------ +# Run parameters: +#------------------------------------------------------------------------------ + +#------ +# Grid: +#------ + +Time::dtfac = 0.4 +MoL::ODE_Method = "rk2" +MoL::MoL_Intermediate_Steps = 2 +#MoL::ODE_Method = "rk4" +#MoL::MoL_Intermediate_Steps = 4 +#MoL::MoL_Num_Scratch_Levels = 1 + +CartGrid3D::type = "coordbase" +CartGrid3D::domain = "full" +CartGrid3D::avoid_origin = "no" + +CoordBase::domainsize = "minmax" +CoordBase::spacing = "gridspacing" # "gridspacing" or "numcells" + +CoordBase::xmin = -1.0 +CoordBase::ymin = -0.5 +#CoordBase::zmin = -0.031250 +#CoordBase::zmin = -0.015625 +CoordBase::zmin = -0.0078125 +CoordBase::xmax = 1.0 +CoordBase::ymax = 0.5 +#CoordBase::zmax = 0.031250 +#CoordBase::zmax = 0.015625 +CoordBase::zmax = 0.0078125 +#CoordBase::dx = 0.0156250 # dx or ncells_x +#CoordBase::dy = 0.0156250 # dy or ncells_y +#CoordBase::dz = 0.0156250 # dz or ncells_z +CoordBase::dx = 0.00390625 # dx or ncells_x +CoordBase::dy = 0.00390625 # dy or ncells_y +CoordBase::dz = 0.00390625 # dz or ncells_z + +CoordBase::boundary_size_x_lower = 3 +CoordBase::boundary_size_y_lower = 3 +CoordBase::boundary_size_z_lower = 3 +CoordBase::boundary_size_x_upper = 3 +CoordBase::boundary_size_y_upper = 3 +CoordBase::boundary_size_z_upper = 3 +CoordBase::boundary_shiftout_x_lower = 0 +CoordBase::boundary_shiftout_y_lower = 0 +CoordBase::boundary_shiftout_z_lower = 0 +CoordBase::boundary_shiftout_x_upper = 0 +CoordBase::boundary_shiftout_y_upper = 0 +CoordBase::boundary_shiftout_z_upper = 0 +CoordBase::boundary_staggered_x_lower = "yes" +CoordBase::boundary_staggered_y_lower = "yes" +CoordBase::boundary_staggered_z_lower = "yes" +CoordBase::boundary_staggered_x_upper = "yes" +CoordBase::boundary_staggered_y_upper = "yes" +CoordBase::boundary_staggered_z_upper = "yes" + +Driver::ghost_size = 3 +GRHydro::GRHydro_stencil = 3 + +Periodic::periodic = "yes" +Periodic::periodic_x = "yes" +Periodic::periodic_y = "yes" +Periodic::periodic_z = "yes" + +# General Carpet parameters: +Carpet::domain_from_coordbase = "yes" +Carpet::enable_all_storage = "no" +Carpet::schedule_barriers = "no" + +CarpetLib::poison_new_memory = "yes" +CarpetLib::poison_value = 114 +CarpetLib::check_bboxes = "no" +CarpetLib::interleave_communications = "yes" +#CarpetLib::vary_tags = "yes" # git version only! +CarpetLib::combine_sends = "yes" +#CarpetLib::reduce_mpi_waitall = "yes" # git version only! + +CarpetInterp::tree_search = "yes" +CarpetInterp::check_tree_search = "no" + +#------ +# MODEL: +#------ + +HydroBase::initial_hydro = "advectedloop" + +GRHydro_InitData::advectedloop_type = "2D" +#GRHydro_InitData::advectedloop_case = "V^z/=0" +GRHydro_InitData::advectedloop_case = "V^z=0" +GRHydro_InitData::advectedloop_delA = "Exact" +#GRHydro_InitData::advectedloop_delA = "Numeric" + +EOS_Omni::gl_gamma = 1.666666666666667 +GRHydro::initial_Gamma = 1.666666666666667 + +#---------- +# Numerics: +#---------- + +#TmunuBase::stress_energy_storage = "yes" +#TmunuBase::stress_energy_at_RHS = "no" +#TmunuBase::timelevels = 1 +#TmunuBase::prolongation_type = "none" +#TmunuBase::support_old_CalcTmunu_mechanism = "no" + +HydroBase::timelevels = 2 +#HydroBase::timelevels = 3 + +SpaceMask::use_mask = "yes" + +#----------- +# Evolution: +#----------- + +HydroBase::evolution_method = "GRHydro" +HydroBase::Bvec_evolution_method = "GRHydro" + +GRHydro::method_type = "RSA FV" +GRHydro::recon_method = "tvd" +GRHydro::tvd_limiter = "vanleerMC2" +GRHydro::recon_vars = "primitive" +GRHydro::riemann_solver = "HLLE" +GRHydro::bound = "none" +GRHydro::GRHydro_eos_type = "General" +GRHydro::GRHydro_eos_table = "Ideal_Fluid" + +GRHydro::Grhydro_MaxNumConstrainedVars = 13 +GRHydro::GRHydro_MaxNumEvolvedVars = 10 + +GRHydro::clean_divergence = "yes" +GRHydro::kap_dc = 10.0 +GRHydro::track_divB = "yes" +GRHydro::calculate_bcom = "yes" + +GRHydro::psidcspeed = "light speed" +GRHydro::decouple_normal_Bfield = "yes" + +#------------------------------------------------------------------------------ +# Output: +#------------------------------------------------------------------------------ + +IO::out_dir = $parfile +IO::out_fileinfo = "none" + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_reductions = "maximum" +IOBasic::outInfo_vars = " + Carpet::physical_time_per_hour + GRHydro::bcom_sq + HydroBase::press + HydroBase::rho + SystemStatistics::maxrss_mb +" + +IOScalar::outScalar_every = 640 #t=1 +IOScalar::one_file_per_group = "yes" +IOScalar::outScalar_reductions = "minimum maximum average norm1 norm2" +IOScalar::outScalar_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GrHydro::divB + GRHydro::scon + GRHydro::tau + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz + SystemStatistics::process_memory_mb +" + +IOASCII::one_file_per_group = "yes" +IOASCII::output_symmetry_points = "no" +IOASCII::output_ghost_points = "no" #hg version only! +IOASCII::output_boundary_points = "no" #hg version only! +#IOASCII::out3D_ghosts = "no" #hg (DEPRECATED) +#IOASCII::out3D_outer_ghosts = "no" #hg (DEPRECATED) + +IOASCII::out0D_every = 640 +IOASCII::out0D_vars = " + Carpet::timing +" + +#Set these IOASCII options for initial data only: +IOASCII::out1D_every = 15360 +#IOASCII::out1D_dt = 10 +IOASCII::out1D_d = "yes" +IOASCII::out1D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GrHydro::divB + GrHydro::psidc + GRHydro::scon + GRHydro::tau + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +#IOASCII::out2D_every = 256000000 +#IOASCII::out2D_vars = " +# GRHydro::bcom +# GRHydro::bcom0 +# GRHydro::bcom_sq +# GRHydro::Bcons +# GRHydro::dens +# GrHydro::divB +# GRHydro::scon +# GRHydro::tau +# HydroBase::Bvec +# HydroBase::eps +# HydroBase::press +# HydroBase::rho +# HydroBase::vel +# HydroBase::w_lorentz +#" + +CarpetIOHDF5::one_file_per_group = "yes" +CarpetIOHDF5::open_one_input_file_at_a_time = "yes" + +CarpetIOHDF5::output_symmetry_points = "yes" +CarpetIOHDF5::output_ghost_points = "yes" +CarpetIOHDF5::output_boundary_points = "yes" +CarpetIOHDF5::output_buffer_points = "yes" + +CarpetIOHDF5::out2D_every = 320 +CarpetIOHDF5::out2D_xy = "yes" +CarpetIOHDF5::out2D_xz = "no" +CarpetIOHDF5::out2D_yz = "no" +CarpetIOHDF5::out2D_xyplane_z = 0.0 +CarpetIOHDF5::out2D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GRHydro::divB + GrHydro::psidc + GRHydro::scon + GRHydro::tau + Grid::coordinates + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +CarpetIOHDF5::out3D_every = 3840 +CarpetIOHDF5::out3D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::divB + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +#------------------------------------------------------------------------------ +# Analysis: +#------------------------------------------------------------------------------ + +#------------------------------------------------------------------------------ +# Checkpoint/Recovery: +#------------------------------------------------------------------------------ +#IOHDF5::checkpoint = "yes" +#IO::checkpoint_dir = $parfile +#IO::checkpoint_ID = "yes" +#IO::checkpoint_every = 6144 +##IO::checkpoint_every_walltime_hours = 6.0 +#IO::checkpoint_keep=3 +#IO::checkpoint_on_terminate = "yes" + +#IO::recover = "autoprobe" +#IO::recover_dir = $parfile + + File [added]: advectedloop2D_vz_hi.par Delta lines: +398 -0 =================================================================== --- trunk/par/advectedloop2D_vz_hi.par (rev 0) +++ trunk/par/advectedloop2D_vz_hi.par 2012-08-09 06:26:54 UTC (rev 147) @@ -0,0 +1,398 @@ +#============================================================================== +# Advected magnetic field loop (2d version): GRHydro_InitData +#============================================================================== + +#------------------------------------------------------------------------------ +# Cactus parameters: +#------------------------------------------------------------------------------ +Cactus::cctk_run_title = "Advected magnetic field loop (2d version)" +Cactus::cctk_full_warnings = "yes" + +#Cactus::terminate = "never" +Cactus::terminate = "time" +Cactus::cctk_final_time = 25.0 +#Cactus::terminate = "iteration" +#Cactus::cctk_itlast = 0 +#Cactus::cctk_itlast = 256 +#Cactus::terminate="runtime" +#Cactus::max_runtime=28800 # 46 days +#Cactus::max_runtime=1410 # 23.5 hours +#Cactus::max_runtime=1380 # 23 hours +#Cactus::max_runtime=690 # 11.5 hours +#Cactus::max_runtime=450 # 7.5 hours +#Cactus::max_runtime=210 # 3.5 hours + +#------------------------------------------------------------------------------ +# Activate all necessary thorns: +#------------------------------------------------------------------------------ + +ActiveThorns = "Boundary CartGrid3D CoordBase Fortran InitBase IOUtil LocalReduce SymBase Time" +ActiveThorns = "AEILocalInterp LocalInterp" +ActiveThorns = "MoL Periodic Slab SpaceMask" +ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetRegrid2 CarpetSlab LoopControl" +#ActiveThorns = "Formaline NaNChecker TerminationTrigger TimerReport" +ActiveThorns = "NaNChecker TerminationTrigger TimerReport" +ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal" +ActiveThorns = "Constants TmunuBase HydroBase " +ActiveThorns = "EOS_Omni" +ActiveThorns = "GRHydro GRHydro_InitData" +ActiveThorns = "HDF5" +ActiveThorns = "SystemStatistics" + +#------------------------------------------------------------------------------ +# Diagnostic parameters: +#------------------------------------------------------------------------------ +#AHFinderDirect::verbose_level = "physics details" + +Carpet::output_timers_every = 1024 +Carpet::output_timer_tree_every = 256 +Carpet::output_initialise_timer_tree = "yes" +Carpet::storage_verbose = "no" +#Carpet::timers_verbose = "no" +#Carpet::verbose = "yes" # only for Initial Data +Carpet::verbose = "no" +Carpet::veryverbose = "no" +Carpet::grid_structure_filename = "carpet-grid-structure.asc" +Carpet::grid_coordinates_filename = "carpet-grid-coordinates.asc" + +CarpetLib::output_bboxes = "no" +CarpetLib::print_memstats_every = 1024 +CarpetLib::memstat_file = "carpetlib-memory-statistics.asc" +#CarpetLib::print_timestats_every = 1024 +#CarpetLib::timestat_file = "carpetlib-timing-statistics.asc" + +CarpetReduce::verbose = "no" +CarpetRegrid2::verbose = "yes" +#CarpetRegrid2::veryverbose = "yes" + +#NaNChecker::verbose = "all" + +TimerReport::out_every = 1024 +TimerReport::out_filename = "TimerReport" +#TimerReport::output_all_timers = "yes" +#TimerReport::output_all_timers_together = "yes" +#TimerReport::output_all_timers_readable = "yes" +#TimerReport::before_checkpoint = "yes" +TimerReport::n_top_timers = 40 + +#------------------------------------------------------------------------------ +# Utility parameters: +#------------------------------------------------------------------------------ + +#NaNChecker::check_every = 16 +NaNChecker::check_every = 1 +NaNChecker::check_vars = " + HydroBase::Bvec + HydroBase::rho + HydroBase::eps + HydroBase::press + HydroBase::vel + HydroBase::w_lorentz +" +NaNChecker::action_if_found = "terminate" +#NaNChecker::action_if_found = "abort" +#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort" +#NaNChecker::check_after=0 + +#TerminationTrigger::max_walltime = 0.0 # hours +TerminationTrigger::max_walltime = 12.0 # hours +TerminationTrigger::on_remaining_walltime = 10.0 # minutes +TerminationTrigger::termination_from_file = "yes" +TerminationTrigger::create_termination_file = "yes" +TerminationTrigger::termination_file = "../TERMINATE" + + +#------------------------------------------------------------------------------ +# Run parameters: +#------------------------------------------------------------------------------ + +#------ +# Grid: +#------ + +Time::dtfac = 0.4 +MoL::ODE_Method = "rk2" +MoL::MoL_Intermediate_Steps = 2 +#MoL::ODE_Method = "rk4" +#MoL::MoL_Intermediate_Steps = 4 +#MoL::MoL_Num_Scratch_Levels = 1 + +CartGrid3D::type = "coordbase" +CartGrid3D::domain = "full" +CartGrid3D::avoid_origin = "no" + +CoordBase::domainsize = "minmax" +CoordBase::spacing = "gridspacing" # "gridspacing" or "numcells" + +CoordBase::xmin = -1.0 +CoordBase::ymin = -0.5 +#CoordBase::zmin = -0.031250 +#CoordBase::zmin = -0.015625 +CoordBase::zmin = -0.0078125 +CoordBase::xmax = 1.0 +CoordBase::ymax = 0.5 +#CoordBase::zmax = 0.031250 +#CoordBase::zmax = 0.015625 +CoordBase::zmax = 0.0078125 +#CoordBase::dx = 0.0156250 # dx or ncells_x +#CoordBase::dy = 0.0156250 # dy or ncells_y +#CoordBase::dz = 0.0156250 # dz or ncells_z +CoordBase::dx = 0.00390625 # dx or ncells_x +CoordBase::dy = 0.00390625 # dy or ncells_y +CoordBase::dz = 0.00390625 # dz or ncells_z + +CoordBase::boundary_size_x_lower = 3 +CoordBase::boundary_size_y_lower = 3 +CoordBase::boundary_size_z_lower = 3 +CoordBase::boundary_size_x_upper = 3 +CoordBase::boundary_size_y_upper = 3 +CoordBase::boundary_size_z_upper = 3 +CoordBase::boundary_shiftout_x_lower = 0 +CoordBase::boundary_shiftout_y_lower = 0 +CoordBase::boundary_shiftout_z_lower = 0 +CoordBase::boundary_shiftout_x_upper = 0 +CoordBase::boundary_shiftout_y_upper = 0 +CoordBase::boundary_shiftout_z_upper = 0 +CoordBase::boundary_staggered_x_lower = "yes" +CoordBase::boundary_staggered_y_lower = "yes" +CoordBase::boundary_staggered_z_lower = "yes" +CoordBase::boundary_staggered_x_upper = "yes" +CoordBase::boundary_staggered_y_upper = "yes" +CoordBase::boundary_staggered_z_upper = "yes" + +Driver::ghost_size = 3 +GRHydro::GRHydro_stencil = 3 + +Periodic::periodic = "yes" +Periodic::periodic_x = "yes" +Periodic::periodic_y = "yes" +Periodic::periodic_z = "yes" + +# General Carpet parameters: +Carpet::domain_from_coordbase = "yes" +Carpet::enable_all_storage = "no" +Carpet::schedule_barriers = "no" + +CarpetLib::poison_new_memory = "yes" +CarpetLib::poison_value = 114 +CarpetLib::check_bboxes = "no" +CarpetLib::interleave_communications = "yes" +#CarpetLib::vary_tags = "yes" # git version only! +CarpetLib::combine_sends = "yes" +#CarpetLib::reduce_mpi_waitall = "yes" # git version only! + +CarpetInterp::tree_search = "yes" +CarpetInterp::check_tree_search = "no" + +#------ +# MODEL: +#------ + +HydroBase::initial_hydro = "advectedloop" + +GRHydro_InitData::advectedloop_type = "2D" +GRHydro_InitData::advectedloop_case = "V^z/=0" +#GRHydro_InitData::advectedloop_case = "V^z=0" +GRHydro_InitData::advectedloop_delA = "Exact" +#GRHydro_InitData::advectedloop_delA = "Numeric" + +EOS_Omni::gl_gamma = 1.666666666666667 +GRHydro::initial_Gamma = 1.666666666666667 + +#---------- +# Numerics: +#---------- + +#TmunuBase::stress_energy_storage = "yes" +#TmunuBase::stress_energy_at_RHS = "no" +#TmunuBase::timelevels = 1 +#TmunuBase::prolongation_type = "none" +#TmunuBase::support_old_CalcTmunu_mechanism = "no" + +HydroBase::timelevels = 2 +#HydroBase::timelevels = 3 + +SpaceMask::use_mask = "yes" + +#----------- +# Evolution: +#----------- + +HydroBase::evolution_method = "GRHydro" +HydroBase::Bvec_evolution_method = "GRHydro" + +GRHydro::method_type = "RSA FV" +GRHydro::recon_method = "tvd" +GRHydro::tvd_limiter = "vanleerMC2" +GRHydro::recon_vars = "primitive" +GRHydro::riemann_solver = "HLLE" +GRHydro::bound = "none" +GRHydro::GRHydro_eos_type = "General" +GRHydro::GRHydro_eos_table = "Ideal_Fluid" + +GRHydro::Grhydro_MaxNumConstrainedVars = 13 +GRHydro::GRHydro_MaxNumEvolvedVars = 10 + +GRHydro::clean_divergence = "yes" +GRHydro::kap_dc = 10.0 +GRHydro::track_divB = "yes" +GRHydro::calculate_bcom = "yes" + +GRHydro::psidcspeed = "light speed" +GRHydro::decouple_normal_Bfield = "yes" + +#------------------------------------------------------------------------------ +# Output: +#------------------------------------------------------------------------------ + +IO::out_dir = $parfile +IO::out_fileinfo = "none" + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_reductions = "maximum" +IOBasic::outInfo_vars = " + Carpet::physical_time_per_hour + GRHydro::bcom_sq + HydroBase::press + HydroBase::rho + SystemStatistics::maxrss_mb +" + +IOScalar::outScalar_every = 640 #t=1 +IOScalar::one_file_per_group = "yes" +IOScalar::outScalar_reductions = "minimum maximum average norm1 norm2" +IOScalar::outScalar_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GrHydro::divB + GRHydro::scon + GRHydro::tau + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz + SystemStatistics::process_memory_mb +" + +IOASCII::one_file_per_group = "yes" +IOASCII::output_symmetry_points = "no" +IOASCII::output_ghost_points = "no" #hg version only! +IOASCII::output_boundary_points = "no" #hg version only! +#IOASCII::out3D_ghosts = "no" #hg (DEPRECATED) +#IOASCII::out3D_outer_ghosts = "no" #hg (DEPRECATED) + +IOASCII::out0D_every = 640 +IOASCII::out0D_vars = " + Carpet::timing +" + +#Set these IOASCII options for initial data only: +IOASCII::out1D_every = 15360 +#IOASCII::out1D_dt = 10 +IOASCII::out1D_d = "yes" +IOASCII::out1D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GrHydro::divB + GrHydro::psidc + GRHydro::scon + GRHydro::tau + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +#IOASCII::out2D_every = 256000000 +#IOASCII::out2D_vars = " +# GRHydro::bcom +# GRHydro::bcom0 +# GRHydro::bcom_sq +# GRHydro::Bcons +# GRHydro::dens +# GrHydro::divB +# GRHydro::scon +# GRHydro::tau +# HydroBase::Bvec +# HydroBase::eps +# HydroBase::press +# HydroBase::rho +# HydroBase::vel +# HydroBase::w_lorentz +#" + +CarpetIOHDF5::one_file_per_group = "yes" +CarpetIOHDF5::open_one_input_file_at_a_time = "yes" + +CarpetIOHDF5::output_symmetry_points = "yes" +CarpetIOHDF5::output_ghost_points = "yes" +CarpetIOHDF5::output_boundary_points = "yes" +CarpetIOHDF5::output_buffer_points = "yes" + +CarpetIOHDF5::out2D_every = 320 +CarpetIOHDF5::out2D_xy = "yes" +CarpetIOHDF5::out2D_xz = "no" +CarpetIOHDF5::out2D_yz = "no" +CarpetIOHDF5::out2D_xyplane_z = 0.0 +CarpetIOHDF5::out2D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GRHydro::divB + GrHydro::psidc + GRHydro::scon + GRHydro::tau + Grid::coordinates + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +CarpetIOHDF5::out3D_every = 3840 +CarpetIOHDF5::out3D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::divB + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +#------------------------------------------------------------------------------ +# Analysis: +#------------------------------------------------------------------------------ + +#------------------------------------------------------------------------------ +# Checkpoint/Recovery: +#------------------------------------------------------------------------------ +#IOHDF5::checkpoint = "yes" +#IO::checkpoint_dir = $parfile +#IO::checkpoint_ID = "yes" +#IO::checkpoint_every = 6144 +##IO::checkpoint_every_walltime_hours = 6.0 +#IO::checkpoint_keep=3 +#IO::checkpoint_on_terminate = "yes" + +#IO::recover = "autoprobe" +#IO::recover_dir = $parfile + +

1 0

[svn:einsteintoolkit] GRHydro_InitData/trunk/par/ (Rev. 146)
by rhaas＠tapir.caltech.edu 09 Aug '12

09 Aug '12

User: rhaas Date: 2012/08/09 01:26 AM Added: /trunk/par/ balsara2_2d.par, balsara3_2d.par, balsara4_2d.par, balsara5_2d.par Log: update balsara par files From: Bruno Coutinho Mundim <bcmsma(a)astro.rit.edu> File Changes: Directory: /trunk/par/ ====================== File [added]: balsara2_2d.par Delta lines: <Unavailable> File [added]: balsara3_2d.par Delta lines: <Unavailable> File [added]: balsara4_2d.par Delta lines: <Unavailable> File [added]: balsara5_2d.par Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] GRHydro_InitData/trunk/par/ (Rev. 145)
by rhaas＠tapir.caltech.edu 09 Aug '12

09 Aug '12

User: rhaas Date: 2012/08/09 01:26 AM Added: /trunk/par/ balsara1_2d.par, balsara1_2d_low.par Log: GRHydro_InitData: update balsara1_2d par files. From: Bruno Coutinho Mundim <bcmsma(a)astro.rit.edu> File Changes: Directory: /trunk/par/ ====================== File [added]: balsara1_2d.par Delta lines: +384 -0 =================================================================== --- trunk/par/balsara1_2d.par (rev 0) +++ trunk/par/balsara1_2d.par 2012-08-09 06:26:49 UTC (rev 145) @@ -0,0 +1,384 @@ +#============================================================================== +# Balsara magnetized shock tube test #1: GRHydro_InitData +#============================================================================== + +#------------------------------------------------------------------------------ +# Cactus parameters: +#------------------------------------------------------------------------------ +Cactus::cctk_run_title = "Balsara magnetized shock tube test #1 along z=0 plane square diagonal" +Cactus::cctk_full_warnings = "yes" + +#Cactus::terminate = "never" +Cactus::terminate = "time" +Cactus::cctk_final_time = 0.4 +#Cactus::terminate = "iteration" +#Cactus::cctk_itlast = 0 +#Cactus::cctk_itlast = 100 +#Cactus::cctk_itlast = 256 +#Cactus::cctk_itlast = 1024 +#Cactus::terminate="runtime" +#Cactus::max_runtime=28800 # 46 days +#Cactus::max_runtime=1410 # 23.5 hours +#Cactus::max_runtime=1380 # 23 hours +#Cactus::max_runtime=690 # 11.5 hours +#Cactus::max_runtime=450 # 7.5 hours +#Cactus::max_runtime=210 # 3.5 hours + +#------------------------------------------------------------------------------ +# Activate all necessary thorns: +#------------------------------------------------------------------------------ + +ActiveThorns = "Boundary CartGrid3D CoordBase Fortran InitBase IOUtil LocalReduce SymBase Time" +ActiveThorns = "AEILocalInterp LocalInterp" +ActiveThorns = "MoL Slab SpaceMask" +ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetSlab LoopControl" +#ActiveThorns = "Formaline NaNChecker TerminationTrigger TimerReport" +ActiveThorns = "NaNChecker TerminationTrigger TimerReport" +ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal" +ActiveThorns = "Constants TmunuBase HydroBase " +#ActiveThorns = "EOS_Base EOS_IdealFluid EOS_Polytrope" +ActiveThorns = "EOS_Omni" +ActiveThorns = "GRHydro GRHydro_InitData" +ActiveThorns = "HDF5" +ActiveThorns = "SystemStatistics" + +#------------------------------------------------------------------------------ +# Diagnostic parameters: +#------------------------------------------------------------------------------ +Carpet::output_timers_every = 1024 +Carpet::output_timer_tree_every = 256 +Carpet::output_initialise_timer_tree = "yes" +Carpet::storage_verbose = "no" +#Carpet::timers_verbose = "no" +#Carpet::verbose = "yes" # only for Initial Data +Carpet::verbose = "no" +Carpet::veryverbose = "no" +Carpet::grid_structure_filename = "carpet-grid-structure.asc" +Carpet::grid_coordinates_filename = "carpet-grid-coordinates.asc" + +CarpetLib::output_bboxes = "no" +CarpetLib::print_memstats_every = 1024 +CarpetLib::memstat_file = "carpetlib-memory-statistics.asc" +#CarpetLib::print_timestats_every = 1024 +#CarpetLib::timestat_file = "carpetlib-timing-statistics.asc" + +CarpetReduce::verbose = "no" +#CarpetRegrid2::verbose = "yes" +#CarpetRegrid2::veryverbose = "yes" + +#NaNChecker::verbose = "all" + +TimerReport::out_every = 1024 +TimerReport::out_filename = "TimerReport" +#TimerReport::output_all_timers = "yes" +#TimerReport::output_all_timers_together = "yes" +#TimerReport::output_all_timers_readable = "yes" +#TimerReport::before_checkpoint = "yes" +TimerReport::n_top_timers = 40 + +#------------------------------------------------------------------------------ +# Utility parameters: +#------------------------------------------------------------------------------ + +#NaNChecker::check_every = 256 +NaNChecker::check_every = 1 +NaNChecker::check_vars = " + ADMBase::curv + ADMBase::metric + ADMBase::lapse + ADMBase::shift + HydroBase::Bvec + HydroBase::rho + HydroBase::eps + HydroBase::press + HydroBase::vel + HydroBase::w_lorentz +" +NaNChecker::action_if_found = "terminate" +#NaNChecker::action_if_found = "abort" +#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort" +#NaNChecker::check_after=0 + +#TerminationTrigger::max_walltime = 0.0 # hours +TerminationTrigger::max_walltime = 24.0 # hours +TerminationTrigger::on_remaining_walltime = 10.0 # minutes +TerminationTrigger::termination_from_file = "yes" +TerminationTrigger::create_termination_file = "yes" +TerminationTrigger::termination_file = "../TERMINATE" + + +#------------------------------------------------------------------------------ +# Run parameters: +#------------------------------------------------------------------------------ + +#------ +# Grid: +#------ + +#Time::dtfac = 0.8 +Time::dtfac = 0.4 +MoL::ODE_Method = "rk2" +MoL::MoL_Intermediate_Steps = 2 +#MoL::ODE_Method = "rk4" +#MoL::MoL_Intermediate_Steps = 4 +#MoL::MoL_Num_Scratch_Levels = 1 + +CartGrid3D::type = "coordbase" +CartGrid3D::domain = "full" +CartGrid3D::avoid_origin = "no" + +CoordBase::domainsize = "minmax" +CoordBase::spacing = "numcells" # "gridspacing" or "numcells" + +#CoordBase::xmin = -0.035355339059327376220 +#CoordBase::ymin = -0.035355339059327376220 +#CoordBase::zmin = -0.0017677669529663688110 +#CoordBase::xmax = 0.035355339059327376220 +#CoordBase::ymax = 0.035355339059327376220 +#CoordBase::zmax = 0.0017677669529663688110 +#CoordBase::ncells_x = 160 # dx or ncells_x +#CoordBase::ncells_y = 160 # dy or ncells_y +#CoordBase::ncells_z = 8 # dz or ncells_z + +# Approximate memory footprint: 1.3GB x 16 = 20.8GB + +CoordBase::xmin = -0.35355339059327376220 +CoordBase::ymin = -0.35355339059327376220 +CoordBase::zmin = -0.0017677669529663688110 +CoordBase::xmax = 0.35355339059327376220 +CoordBase::ymax = 0.35355339059327376220 +CoordBase::zmax = 0.0017677669529663688110 +CoordBase::ncells_x = 1600 # dx or ncells_x +CoordBase::ncells_y = 1600 # dy or ncells_y +CoordBase::ncells_z = 8 # dz or ncells_z + +Driver::ghost_size = 3 +GRHydro::GRHydro_stencil = 3 + +# General Carpet parameters: +Carpet::domain_from_coordbase = "yes" + +CarpetLib::poison_new_memory = "yes" +CarpetLib::poison_value = 114 + +#------ +# MODEL: +#------ + +HydroBase::initial_hydro = "shocktube" +GRHydro_InitData::shocktube_type = "diagshock2d" +GRHydro_InitData::shock_case = "Balsara1" + +ADMBase::metric_type = "physical" +ADMBase::initial_data = "Cartesian Minkowski" +ADMBase::initial_lapse = "one" +ADMBase::initial_shift = "zero" +ADMBase::initial_dtlapse = "none" +ADMBase::initial_dtshift = "none" + +EOS_Omni::poly_K = 100.0 +EOS_Omni::poly_gamma = 2.0 +EOS_Omni::gl_K = 100.0 +EOS_Omni::gl_gamma = 2.0 +GRHydro::initial_Gamma = 2.0 + +#---------- +# Numerics: +#---------- +#Boundary::radpower = 2 +#HydroBase::timelevels = 3 +HydroBase::timelevels = 2 +#ADMMacros::spatial_order = 4 +SpaceMask::use_mask = "yes" + +#----------- +# Evolution: +#----------- + +HydroBase::evolution_method = "GRHydro" +HydroBase::Bvec_evolution_method = "GRHydro" + +GRHydro::method_type = "RSA FV" +GRHydro::recon_method = "tvd" +#GRHydro::recon_method = "ppm" +#GRHydro::recon_method = "eno" +GRHydro::tvd_limiter = "vanleerMC2" +#GRHydro::tvd_limiter = "minmod" +GRHydro::recon_vars = "primitive" +GRHydro::riemann_solver = "HLLE" +GRHydro::decouple_normal_Bfield = "yes" +#GRHydro::riemann_solver = "LLF" +GRHydro::bound = "none" +GRHydro::GRHydro_eos_type = "General" +GRHydro::GRHydro_eos_table = "Ideal_Fluid" + +GRHydro::Grhydro_MaxNumConstrainedVars = 13 +GRHydro::GRHydro_MaxNumEvolvedVars = 10 + +GRHydro::clean_divergence = "yes" +GRHydro::psidcspeed = "light speed" +GRHydro::kap_dc = 10.0 +GRHydro::track_divB = "yes" +GRHydro::calculate_bcom = "yes" + +#------------------------------------------------------------------------------ +# Output: +#------------------------------------------------------------------------------ + +IO::out_dir = $parfile +IO::out_fileinfo = "none" + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_reductions = "maximum" +IOBasic::outInfo_vars = " + Carpet::physical_time_per_hour + GRHydro::bcom_sq + HydroBase::rho + HydroBase::press + SystemStatistics::maxrss_mb +" + +#IOScalar::outScalar_every = 1024 +#IOScalar::one_file_per_group = "yes" +#IOScalar::outScalar_reductions = "minimum maximum average norm1 norm2" +#IOScalar::outScalar_vars = " +# GRHydro::bcom +# GRHydro::bcom0 +# GRHydro::bcom_sq +# GRHydro::Bcons +# GRHydro::dens +# GrHydro::divB +# GRHydro::scon +# GRHydro::tau +# HydroBase::Bvec +# HydroBase::eps +# HydroBase::press +# HydroBase::rho +# HydroBase::vel +# HydroBase::w_lorentz +# SystemStatistics::process_memory_mb +#" + +IOASCII::one_file_per_group = "yes" +IOASCII::output_symmetry_points = "no" +IOASCII::output_ghost_points = "no" #hg version only! +IOASCII::output_boundary_points = "no" #hg version only! +#IOASCII::out3D_ghosts = "no" #hg (DEPRECATED) +#IOASCII::out3D_outer_ghosts = "no" #hg (DEPRECATED) + +IOASCII::out0D_every = 200 +IOASCII::out0D_vars = " + Carpet::timing +" + +#Set these IOASCII options for initial data only: +#IOASCII::out1D_every = 1 +IOASCII::out1D_every = 2262 +#IOASCII::out1D_dt = 0.01 +IOASCII::out1D_d = "yes" +IOASCII::out1D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GrHydro::divB + GrHydro::psidc + GRHydro::scon + GRHydro::tau + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +IOASCII::out2D_every = 2262 +IOASCII::out2D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GrHydro::divB + GrHydro::psidc + GRHydro::scon + GRHydro::tau + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +CarpetIOHDF5::one_file_per_group = "yes" +CarpetIOHDF5::open_one_input_file_at_a_time = "yes" + +#CarpetIOHDF5::output_symmetry_points = "no" +#CarpetIOHDF5::output_ghost_points = "no" +#CarpetIOHDF5::output_boundary_points = "no" +#CarpetIOHDF5::output_buffer_points = "no" + +#CarpetIOHDF5::out2D_every = 1 +CarpetIOHDF5::out2D_every = 100 +#CarpetIOHDF5::out2D_dt = 0.01 +CarpetIOHDF5::out2D_xy = "yes" +CarpetIOHDF5::out2D_xz = "no" +CarpetIOHDF5::out2D_yz = "no" +CarpetIOHDF5::out2D_xyplane_z = 0.0 +CarpetIOHDF5::out2D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GRHydro::divB + GrHydro::psidc + GRHydro::scon + GRHydro::tau + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +#CarpetIOHDF5::out3D_every = 1 +CarpetIOHDF5::out3D_every = 100 +#CarpetIOHDF5::out3D_dt = 0.01 +CarpetIOHDF5::out3D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::divB + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +#------------------------------------------------------------------------------ +# Analysis: +#------------------------------------------------------------------------------ + + +#------------------------------------------------------------------------------ +# Checkpoint/Recovery: +#------------------------------------------------------------------------------ +#IOHDF5::checkpoint = "yes" +#IO::checkpoint_dir = $parfile +#IO::checkpoint_ID = "yes" +#IO::checkpoint_every = 6144 +##IO::checkpoint_every_walltime_hours = 6.0 +#IO::checkpoint_keep=3 +#IO::checkpoint_on_terminate = "yes" +# +#IO::recover = "autoprobe" +#IO::recover_dir = $parfile + + File [added]: balsara1_2d_low.par Delta lines: +384 -0 =================================================================== --- trunk/par/balsara1_2d_low.par (rev 0) +++ trunk/par/balsara1_2d_low.par 2012-08-09 06:26:49 UTC (rev 145) @@ -0,0 +1,384 @@ +#============================================================================== +# Balsara magnetized shock tube test #1: GRHydro_InitData +#============================================================================== + +#------------------------------------------------------------------------------ +# Cactus parameters: +#------------------------------------------------------------------------------ +Cactus::cctk_run_title = "Balsara magnetized shock tube test #1 along z=0 plane square diagonal" +Cactus::cctk_full_warnings = "yes" + +#Cactus::terminate = "never" +Cactus::terminate = "time" +Cactus::cctk_final_time = 0.4 +#Cactus::terminate = "iteration" +#Cactus::cctk_itlast = 0 +#Cactus::cctk_itlast = 100 +#Cactus::cctk_itlast = 256 +#Cactus::cctk_itlast = 1024 +#Cactus::terminate="runtime" +#Cactus::max_runtime=28800 # 46 days +#Cactus::max_runtime=1410 # 23.5 hours +#Cactus::max_runtime=1380 # 23 hours +#Cactus::max_runtime=690 # 11.5 hours +#Cactus::max_runtime=450 # 7.5 hours +#Cactus::max_runtime=210 # 3.5 hours + +#------------------------------------------------------------------------------ +# Activate all necessary thorns: +#------------------------------------------------------------------------------ + +ActiveThorns = "Boundary CartGrid3D CoordBase Fortran InitBase IOUtil LocalReduce SymBase Time" +ActiveThorns = "AEILocalInterp LocalInterp" +ActiveThorns = "MoL Slab SpaceMask" +ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetSlab LoopControl" +#ActiveThorns = "Formaline NaNChecker TerminationTrigger TimerReport" +ActiveThorns = "NaNChecker TerminationTrigger TimerReport" +ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal" +ActiveThorns = "Constants TmunuBase HydroBase " +#ActiveThorns = "EOS_Base EOS_IdealFluid EOS_Polytrope" +ActiveThorns = "EOS_Omni" +ActiveThorns = "GRHydro GRHydro_InitData" +ActiveThorns = "HDF5" +ActiveThorns = "SystemStatistics" + +#------------------------------------------------------------------------------ +# Diagnostic parameters: +#------------------------------------------------------------------------------ +Carpet::output_timers_every = 1024 +Carpet::output_timer_tree_every = 256 +Carpet::output_initialise_timer_tree = "yes" +Carpet::storage_verbose = "no" +#Carpet::timers_verbose = "no" +#Carpet::verbose = "yes" # only for Initial Data +Carpet::verbose = "no" +Carpet::veryverbose = "no" +Carpet::grid_structure_filename = "carpet-grid-structure.asc" +Carpet::grid_coordinates_filename = "carpet-grid-coordinates.asc" + +CarpetLib::output_bboxes = "no" +CarpetLib::print_memstats_every = 1024 +CarpetLib::memstat_file = "carpetlib-memory-statistics.asc" +#CarpetLib::print_timestats_every = 1024 +#CarpetLib::timestat_file = "carpetlib-timing-statistics.asc" + +CarpetReduce::verbose = "no" +#CarpetRegrid2::verbose = "yes" +#CarpetRegrid2::veryverbose = "yes" + +#NaNChecker::verbose = "all" + +TimerReport::out_every = 1024 +TimerReport::out_filename = "TimerReport" +#TimerReport::output_all_timers = "yes" +#TimerReport::output_all_timers_together = "yes" +#TimerReport::output_all_timers_readable = "yes" +#TimerReport::before_checkpoint = "yes" +TimerReport::n_top_timers = 40 + +#------------------------------------------------------------------------------ +# Utility parameters: +#------------------------------------------------------------------------------ + +#NaNChecker::check_every = 256 +NaNChecker::check_every = 1 +NaNChecker::check_vars = " + ADMBase::curv + ADMBase::metric + ADMBase::lapse + ADMBase::shift + HydroBase::Bvec + HydroBase::rho + HydroBase::eps + HydroBase::press + HydroBase::vel + HydroBase::w_lorentz +" +NaNChecker::action_if_found = "terminate" +#NaNChecker::action_if_found = "abort" +#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort" +#NaNChecker::check_after=0 + +#TerminationTrigger::max_walltime = 0.0 # hours +TerminationTrigger::max_walltime = 24.0 # hours +TerminationTrigger::on_remaining_walltime = 10.0 # minutes +TerminationTrigger::termination_from_file = "yes" +TerminationTrigger::create_termination_file = "yes" +TerminationTrigger::termination_file = "../TERMINATE" + + +#------------------------------------------------------------------------------ +# Run parameters: +#------------------------------------------------------------------------------ + +#------ +# Grid: +#------ + +#Time::dtfac = 0.8 +Time::dtfac = 0.4 +MoL::ODE_Method = "rk2" +MoL::MoL_Intermediate_Steps = 2 +#MoL::ODE_Method = "rk4" +#MoL::MoL_Intermediate_Steps = 4 +#MoL::MoL_Num_Scratch_Levels = 1 + +CartGrid3D::type = "coordbase" +CartGrid3D::domain = "full" +CartGrid3D::avoid_origin = "no" + +CoordBase::domainsize = "minmax" +CoordBase::spacing = "numcells" # "gridspacing" or "numcells" + +#CoordBase::xmin = -0.035355339059327376220 +#CoordBase::ymin = -0.035355339059327376220 +#CoordBase::zmin = -0.0017677669529663688110 +#CoordBase::xmax = 0.035355339059327376220 +#CoordBase::ymax = 0.035355339059327376220 +#CoordBase::zmax = 0.0017677669529663688110 +#CoordBase::ncells_x = 160 # dx or ncells_x +#CoordBase::ncells_y = 160 # dy or ncells_y +#CoordBase::ncells_z = 8 # dz or ncells_z + +# Approximate memory footprint: 650MB x 8 = 5.2GB + +CoordBase::xmin = -0.35355339059327376220 +CoordBase::ymin = -0.35355339059327376220 +CoordBase::zmin = -0.003535533905932737622 +CoordBase::xmax = 0.35355339059327376220 +CoordBase::ymax = 0.35355339059327376220 +CoordBase::zmax = 0.003535533905932737622 +CoordBase::ncells_x = 800 # dx or ncells_x +CoordBase::ncells_y = 800 # dy or ncells_y +CoordBase::ncells_z = 8 # dz or ncells_z + +Driver::ghost_size = 3 +GRHydro::GRHydro_stencil = 3 + +# General Carpet parameters: +Carpet::domain_from_coordbase = "yes" + +CarpetLib::poison_new_memory = "yes" +CarpetLib::poison_value = 114 + +#------ +# MODEL: +#------ + +HydroBase::initial_hydro = "shocktube" +GRHydro_InitData::shocktube_type = "diagshock2d" +GRHydro_InitData::shock_case = "Balsara1" + +ADMBase::metric_type = "physical" +ADMBase::initial_data = "Cartesian Minkowski" +ADMBase::initial_lapse = "one" +ADMBase::initial_shift = "zero" +ADMBase::initial_dtlapse = "none" +ADMBase::initial_dtshift = "none" + +EOS_Omni::poly_K = 100.0 +EOS_Omni::poly_gamma = 2.0 +EOS_Omni::gl_K = 100.0 +EOS_Omni::gl_gamma = 2.0 +GRHydro::initial_Gamma = 2.0 + +#---------- +# Numerics: +#---------- +#Boundary::radpower = 2 +#HydroBase::timelevels = 3 +HydroBase::timelevels = 2 +#ADMMacros::spatial_order = 4 +SpaceMask::use_mask = "yes" + +#----------- +# Evolution: +#----------- + +HydroBase::evolution_method = "GRHydro" +HydroBase::Bvec_evolution_method = "GRHydro" + +GRHydro::method_type = "RSA FV" +GRHydro::recon_method = "tvd" +#GRHydro::recon_method = "ppm" +#GRHydro::recon_method = "eno" +GRHydro::tvd_limiter = "vanleerMC2" +#GRHydro::tvd_limiter = "minmod" +GRHydro::recon_vars = "primitive" +GRHydro::riemann_solver = "HLLE" +GRHydro::decouple_normal_Bfield = "yes" +#GRHydro::riemann_solver = "LLF" +GRHydro::bound = "none" +GRHydro::GRHydro_eos_type = "General" +GRHydro::GRHydro_eos_table = "Ideal_Fluid" + +GRHydro::Grhydro_MaxNumConstrainedVars = 13 +GRHydro::GRHydro_MaxNumEvolvedVars = 10 + +GRHydro::clean_divergence = "yes" +GRHydro::psidcspeed = "light speed" +GRHydro::kap_dc = 10.0 +GRHydro::track_divB = "yes" +GRHydro::calculate_bcom = "yes" + +#------------------------------------------------------------------------------ +# Output: +#------------------------------------------------------------------------------ + +IO::out_dir = $parfile +IO::out_fileinfo = "none" + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_reductions = "maximum" +IOBasic::outInfo_vars = " + Carpet::physical_time_per_hour + GRHydro::bcom_sq + HydroBase::rho + HydroBase::press + SystemStatistics::maxrss_mb +" + +#IOScalar::outScalar_every = 1024 +#IOScalar::one_file_per_group = "yes" +#IOScalar::outScalar_reductions = "minimum maximum average norm1 norm2" +#IOScalar::outScalar_vars = " +# GRHydro::bcom +# GRHydro::bcom0 +# GRHydro::bcom_sq +# GRHydro::Bcons +# GRHydro::dens +# GrHydro::divB +# GRHydro::scon +# GRHydro::tau +# HydroBase::Bvec +# HydroBase::eps +# HydroBase::press +# HydroBase::rho +# HydroBase::vel +# HydroBase::w_lorentz +# SystemStatistics::process_memory_mb +#" + +IOASCII::one_file_per_group = "yes" +IOASCII::output_symmetry_points = "no" +IOASCII::output_ghost_points = "no" #hg version only! +IOASCII::output_boundary_points = "no" #hg version only! +#IOASCII::out3D_ghosts = "no" #hg (DEPRECATED) +#IOASCII::out3D_outer_ghosts = "no" #hg (DEPRECATED) + +IOASCII::out0D_every = 200 +IOASCII::out0D_vars = " + Carpet::timing +" + +#Set these IOASCII options for initial data only: +#IOASCII::out1D_every = 1 +IOASCII::out1D_every = 1131 +#IOASCII::out1D_dt = 0.01 +IOASCII::out1D_d = "yes" +IOASCII::out1D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GrHydro::divB + GrHydro::psidc + GRHydro::scon + GRHydro::tau + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +IOASCII::out2D_every = 1131 +IOASCII::out2D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GrHydro::divB + GrHydro::psidc + GRHydro::scon + GRHydro::tau + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +CarpetIOHDF5::one_file_per_group = "yes" +CarpetIOHDF5::open_one_input_file_at_a_time = "yes" + +#CarpetIOHDF5::output_symmetry_points = "no" +#CarpetIOHDF5::output_ghost_points = "no" +#CarpetIOHDF5::output_boundary_points = "no" +#CarpetIOHDF5::output_buffer_points = "no" + +#CarpetIOHDF5::out2D_every = 1 +CarpetIOHDF5::out2D_every = 100 +#CarpetIOHDF5::out2D_dt = 0.01 +CarpetIOHDF5::out2D_xy = "yes" +CarpetIOHDF5::out2D_xz = "no" +CarpetIOHDF5::out2D_yz = "no" +CarpetIOHDF5::out2D_xyplane_z = 0.0 +CarpetIOHDF5::out2D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GRHydro::divB + GrHydro::psidc + GRHydro::scon + GRHydro::tau + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +#CarpetIOHDF5::out3D_every = 1 +CarpetIOHDF5::out3D_every = 100 +#CarpetIOHDF5::out3D_dt = 0.01 +CarpetIOHDF5::out3D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::divB + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +#------------------------------------------------------------------------------ +# Analysis: +#------------------------------------------------------------------------------ + + +#------------------------------------------------------------------------------ +# Checkpoint/Recovery: +#------------------------------------------------------------------------------ +#IOHDF5::checkpoint = "yes" +#IO::checkpoint_dir = $parfile +#IO::checkpoint_ID = "yes" +#IO::checkpoint_every = 6144 +##IO::checkpoint_every_walltime_hours = 6.0 +#IO::checkpoint_keep=3 +#IO::checkpoint_on_terminate = "yes" +# +#IO::recover = "autoprobe" +#IO::recover_dir = $parfile + +

1 0

[svn:einsteintoolkit] GRHydro_InitData/trunk/par/ (Rev. 144)
by rhaas＠tapir.caltech.edu 09 Aug '12

09 Aug '12

User: rhaas Date: 2012/08/09 01:26 AM Added: /trunk/par/ advectedloop2D_vz.par Log: GRHydro_InitData: par file for advectedloop test with v^z \= 0. From: Bruno Coutinho Mundim <bcmsma(a)astro.rit.edu> File Changes: Directory: /trunk/par/ ====================== File [added]: advectedloop2D_vz.par Delta lines: +397 -0 =================================================================== --- trunk/par/advectedloop2D_vz.par (rev 0) +++ trunk/par/advectedloop2D_vz.par 2012-08-09 06:26:47 UTC (rev 144) @@ -0,0 +1,397 @@ +#============================================================================== +# Advected magnetic field loop (2d version): GRHydro_InitData +#============================================================================== + +#------------------------------------------------------------------------------ +# Cactus parameters: +#------------------------------------------------------------------------------ +Cactus::cctk_run_title = "Advected magnetic field loop (2d version)" +Cactus::cctk_full_warnings = "yes" + +#Cactus::terminate = "never" +Cactus::terminate = "time" +Cactus::cctk_final_time = 50.0 +#Cactus::terminate = "iteration" +#Cactus::cctk_itlast = 0 +#Cactus::cctk_itlast = 256 +#Cactus::terminate="runtime" +#Cactus::max_runtime=28800 # 46 days +#Cactus::max_runtime=1410 # 23.5 hours +#Cactus::max_runtime=1380 # 23 hours +#Cactus::max_runtime=690 # 11.5 hours +#Cactus::max_runtime=450 # 7.5 hours +#Cactus::max_runtime=210 # 3.5 hours + +#------------------------------------------------------------------------------ +# Activate all necessary thorns: +#------------------------------------------------------------------------------ + +ActiveThorns = "Boundary CartGrid3D CoordBase Fortran InitBase IOUtil LocalReduce SymBase Time" +ActiveThorns = "AEILocalInterp LocalInterp" +ActiveThorns = "MoL Periodic Slab SpaceMask" +ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetRegrid2 CarpetSlab LoopControl" +#ActiveThorns = "Formaline NaNChecker TerminationTrigger TimerReport" +ActiveThorns = "NaNChecker TerminationTrigger TimerReport" +ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal" +ActiveThorns = "Constants TmunuBase HydroBase " +ActiveThorns = "EOS_Omni" +ActiveThorns = "GRHydro GRHydro_InitData" +ActiveThorns = "HDF5" +ActiveThorns = "SystemStatistics" + +#------------------------------------------------------------------------------ +# Diagnostic parameters: +#------------------------------------------------------------------------------ +#AHFinderDirect::verbose_level = "physics details" + +Carpet::output_timers_every = 1024 +Carpet::output_timer_tree_every = 256 +Carpet::output_initialise_timer_tree = "yes" +Carpet::storage_verbose = "no" +#Carpet::timers_verbose = "no" +#Carpet::verbose = "yes" # only for Initial Data +Carpet::verbose = "no" +Carpet::veryverbose = "no" +Carpet::grid_structure_filename = "carpet-grid-structure.asc" +Carpet::grid_coordinates_filename = "carpet-grid-coordinates.asc" + +CarpetLib::output_bboxes = "no" +CarpetLib::print_memstats_every = 1024 +CarpetLib::memstat_file = "carpetlib-memory-statistics.asc" +#CarpetLib::print_timestats_every = 1024 +#CarpetLib::timestat_file = "carpetlib-timing-statistics.asc" + +CarpetReduce::verbose = "no" +CarpetRegrid2::verbose = "yes" +#CarpetRegrid2::veryverbose = "yes" + +#NaNChecker::verbose = "all" + +TimerReport::out_every = 1024 +TimerReport::out_filename = "TimerReport" +#TimerReport::output_all_timers = "yes" +#TimerReport::output_all_timers_together = "yes" +#TimerReport::output_all_timers_readable = "yes" +#TimerReport::before_checkpoint = "yes" +TimerReport::n_top_timers = 40 + +#------------------------------------------------------------------------------ +# Utility parameters: +#------------------------------------------------------------------------------ + +#NaNChecker::check_every = 16 +NaNChecker::check_every = 1 +NaNChecker::check_vars = " + HydroBase::Bvec + HydroBase::rho + HydroBase::eps + HydroBase::press + HydroBase::vel + HydroBase::w_lorentz +" +NaNChecker::action_if_found = "terminate" +#NaNChecker::action_if_found = "abort" +#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort" +#NaNChecker::check_after=0 + +#TerminationTrigger::max_walltime = 0.0 # hours +TerminationTrigger::max_walltime = 24.0 # hours +TerminationTrigger::on_remaining_walltime = 10.0 # minutes +TerminationTrigger::termination_from_file = "yes" +TerminationTrigger::create_termination_file = "yes" +TerminationTrigger::termination_file = "../TERMINATE" + + +#------------------------------------------------------------------------------ +# Run parameters: +#------------------------------------------------------------------------------ + +#------ +# Grid: +#------ + +Time::dtfac = 0.4 +MoL::ODE_Method = "rk2" +MoL::MoL_Intermediate_Steps = 2 +#MoL::ODE_Method = "rk4" +#MoL::MoL_Intermediate_Steps = 4 +#MoL::MoL_Num_Scratch_Levels = 1 + +CartGrid3D::type = "coordbase" +CartGrid3D::domain = "full" +CartGrid3D::avoid_origin = "no" + +CoordBase::domainsize = "minmax" +CoordBase::spacing = "gridspacing" # "gridspacing" or "numcells" + +CoordBase::xmin = -1.0 +CoordBase::ymin = -0.5 +#CoordBase::zmin = -0.031250 +CoordBase::zmin = -0.015625 +CoordBase::xmax = 1.0 +CoordBase::ymax = 0.5 +#CoordBase::zmax = 0.031250 +CoordBase::zmax = 0.015625 +#CoordBase::dx = 0.0156250 # dx or ncells_x +#CoordBase::dy = 0.0156250 # dy or ncells_y +#CoordBase::dz = 0.0156250 # dz or ncells_z +CoordBase::dx = 0.0078125 # dx or ncells_x +CoordBase::dy = 0.0078125 # dy or ncells_y +CoordBase::dz = 0.0078125 # dz or ncells_z + +CoordBase::boundary_size_x_lower = 3 +CoordBase::boundary_size_y_lower = 3 +CoordBase::boundary_size_z_lower = 3 +CoordBase::boundary_size_x_upper = 3 +CoordBase::boundary_size_y_upper = 3 +CoordBase::boundary_size_z_upper = 3 +CoordBase::boundary_shiftout_x_lower = 0 +CoordBase::boundary_shiftout_y_lower = 0 +CoordBase::boundary_shiftout_z_lower = 0 +CoordBase::boundary_shiftout_x_upper = 0 +CoordBase::boundary_shiftout_y_upper = 0 +CoordBase::boundary_shiftout_z_upper = 0 +CoordBase::boundary_staggered_x_lower = "yes" +CoordBase::boundary_staggered_y_lower = "yes" +CoordBase::boundary_staggered_z_lower = "yes" +CoordBase::boundary_staggered_x_upper = "yes" +CoordBase::boundary_staggered_y_upper = "yes" +CoordBase::boundary_staggered_z_upper = "yes" + +Driver::ghost_size = 3 +GRHydro::GRHydro_stencil = 3 + +Periodic::periodic = "yes" +Periodic::periodic_x = "yes" +Periodic::periodic_y = "yes" +Periodic::periodic_z = "yes" + +# General Carpet parameters: +Carpet::domain_from_coordbase = "yes" +Carpet::enable_all_storage = "no" +Carpet::schedule_barriers = "no" + +CarpetLib::poison_new_memory = "yes" +CarpetLib::poison_value = 114 +CarpetLib::check_bboxes = "no" +CarpetLib::interleave_communications = "yes" +#CarpetLib::vary_tags = "yes" # git version only! +CarpetLib::combine_sends = "yes" +#CarpetLib::reduce_mpi_waitall = "yes" # git version only! + +CarpetInterp::tree_search = "yes" +CarpetInterp::check_tree_search = "no" + +#------ +# MODEL: +#------ + +HydroBase::initial_hydro = "advectedloop" + +GRHydro_InitData::advectedloop_type = "2D" +GRHydro_InitData::advectedloop_case = "V^z/=0" +#GRHydro_InitData::advectedloop_case = "V^z=0" +GRHydro_InitData::advectedloop_delA = "Exact" +#GRHydro_InitData::advectedloop_delA = "Numeric" + +EOS_Omni::gl_gamma = 1.666666666666667 +GRHydro::initial_Gamma = 1.666666666666667 + +#---------- +# Numerics: +#---------- + +#TmunuBase::stress_energy_storage = "yes" +#TmunuBase::stress_energy_at_RHS = "no" +#TmunuBase::timelevels = 1 +#TmunuBase::prolongation_type = "none" +#TmunuBase::support_old_CalcTmunu_mechanism = "no" + +HydroBase::timelevels = 2 +#HydroBase::timelevels = 3 + +SpaceMask::use_mask = "yes" + +#----------- +# Evolution: +#----------- + +HydroBase::evolution_method = "GRHydro" +HydroBase::Bvec_evolution_method = "GRHydro" + +GRHydro::method_type = "RSA FV" +GRHydro::recon_method = "tvd" +GRHydro::tvd_limiter = "vanleerMC2" +GRHydro::recon_vars = "primitive" +GRHydro::riemann_solver = "HLLE" +GRHydro::bound = "none" +GRHydro::GRHydro_eos_type = "General" +GRHydro::GRHydro_eos_table = "Ideal_Fluid" + +GRHydro::Grhydro_MaxNumConstrainedVars = 13 +GRHydro::GRHydro_MaxNumEvolvedVars = 10 + +GRHydro::clean_divergence = "yes" +GRHydro::kap_dc = 10.0 +GRHydro::track_divB = "yes" +GRHydro::calculate_bcom = "yes" + +GRHydro::psidcspeed = "light speed" +GRHydro::decouple_normal_Bfield = "yes" + +#------------------------------------------------------------------------------ +# Output: +#------------------------------------------------------------------------------ + +IO::out_dir = $parfile +IO::out_fileinfo = "none" + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_reductions = "maximum" +IOBasic::outInfo_vars = " + Carpet::physical_time_per_hour + GRHydro::bcom_sq + HydroBase::press + HydroBase::rho + SystemStatistics::maxrss_mb +" + +#IOScalar::outScalar_every = 1024 +IOScalar::outScalar_every = 320 #t=1 +IOScalar::one_file_per_group = "yes" +IOScalar::outScalar_reductions = "minimum maximum average norm1 norm2" +IOScalar::outScalar_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GrHydro::divB + GRHydro::scon + GRHydro::tau + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz + SystemStatistics::process_memory_mb +" + +IOASCII::one_file_per_group = "yes" +IOASCII::output_symmetry_points = "no" +IOASCII::output_ghost_points = "no" #hg version only! +IOASCII::output_boundary_points = "no" #hg version only! +#IOASCII::out3D_ghosts = "no" #hg (DEPRECATED) +#IOASCII::out3D_outer_ghosts = "no" #hg (DEPRECATED) + +IOASCII::out0D_every = 320 +IOASCII::out0D_vars = " + Carpet::timing +" + +#Set these IOASCII options for initial data only: +IOASCII::out1D_every = 7680 +#IOASCII::out1D_dt = 10 +IOASCII::out1D_d = "yes" +IOASCII::out1D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GrHydro::divB + GrHydro::psidc + GRHydro::scon + GRHydro::tau + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +#IOASCII::out2D_every = 256000000 +#IOASCII::out2D_vars = " +# GRHydro::bcom +# GRHydro::bcom0 +# GRHydro::bcom_sq +# GRHydro::Bcons +# GRHydro::dens +# GrHydro::divB +# GRHydro::scon +# GRHydro::tau +# HydroBase::Bvec +# HydroBase::eps +# HydroBase::press +# HydroBase::rho +# HydroBase::vel +# HydroBase::w_lorentz +#" + +CarpetIOHDF5::one_file_per_group = "yes" +CarpetIOHDF5::open_one_input_file_at_a_time = "yes" + +CarpetIOHDF5::output_symmetry_points = "yes" +CarpetIOHDF5::output_ghost_points = "yes" +CarpetIOHDF5::output_boundary_points = "yes" +CarpetIOHDF5::output_buffer_points = "yes" + +CarpetIOHDF5::out2D_every = 160 +CarpetIOHDF5::out2D_xy = "yes" +CarpetIOHDF5::out2D_xz = "no" +CarpetIOHDF5::out2D_yz = "no" +CarpetIOHDF5::out2D_xyplane_z = 0.0 +CarpetIOHDF5::out2D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GRHydro::divB + GrHydro::psidc + GRHydro::scon + GRHydro::tau + Grid::coordinates + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +CarpetIOHDF5::out3D_every = 1920 +CarpetIOHDF5::out3D_vars = " + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::divB + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz +" + +#------------------------------------------------------------------------------ +# Analysis: +#------------------------------------------------------------------------------ + +#------------------------------------------------------------------------------ +# Checkpoint/Recovery: +#------------------------------------------------------------------------------ +#IOHDF5::checkpoint = "yes" +#IO::checkpoint_dir = $parfile +#IO::checkpoint_ID = "yes" +#IO::checkpoint_every = 6144 +##IO::checkpoint_every_walltime_hours = 6.0 +#IO::checkpoint_keep=3 +#IO::checkpoint_on_terminate = "yes" + +#IO::recover = "autoprobe" +#IO::recover_dir = $parfile + +

1 0

[svn:einsteintoolkit] GRHydro_InitData/trunk/par/ (Rev. 143)
by rhaas＠tapir.caltech.edu 09 Aug '12

09 Aug '12

User: rhaas Date: 2012/08/09 01:26 AM Modified: /trunk/par/ advectedloop2D.par Log: update advectedloop2D.par From: Bruno Coutinho Mundim <bcmsma(a)astro.rit.edu> File Changes: Directory: /trunk/par/ ====================== File [modified]: advectedloop2D.par Delta lines: +34 -65 =================================================================== --- trunk/par/advectedloop2D.par 2012-08-09 06:26:44 UTC (rev 142) +++ trunk/par/advectedloop2D.par 2012-08-09 06:26:46 UTC (rev 143) @@ -9,11 +9,11 @@ Cactus::cctk_full_warnings = "yes" #Cactus::terminate = "never" -#Cactus::terminate = "time" -#Cactus::cctk_final_time = 1500.0 -Cactus::terminate = "iteration" +Cactus::terminate = "time" +Cactus::cctk_final_time = 50.0 +#Cactus::terminate = "iteration" #Cactus::cctk_itlast = 0 -Cactus::cctk_itlast = 256 +#Cactus::cctk_itlast = 256 #Cactus::terminate="runtime" #Cactus::max_runtime=28800 # 46 days #Cactus::max_runtime=1410 # 23.5 hours @@ -28,33 +28,15 @@ ActiveThorns = "Boundary CartGrid3D CoordBase Fortran InitBase IOUtil LocalReduce SymBase Time" ActiveThorns = "AEILocalInterp LocalInterp" -#ActiveThorns = "MoL Periodic ReflectionSymmetry RotatingSymmetry90 Slab SpaceMask SphericalSurface" -#ActiveThorns = "MoL Periodic ReflectionSymmetry RotatingSymmetry180 Slab SpaceMask SphericalSurface" -#ActiveThorns = "MoL Periodic ReflectionSymmetry Slab SpaceMask SphericalSurface" -ActiveThorns = "MoL Periodic Slab SpaceMask SphericalSurface" -ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetRegrid2 CarpetSlab CarpetTracker CarpetMask LoopControl" +ActiveThorns = "MoL Periodic Slab SpaceMask" +ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetRegrid2 CarpetSlab LoopControl" #ActiveThorns = "Formaline NaNChecker TerminationTrigger TimerReport" ActiveThorns = "NaNChecker TerminationTrigger TimerReport" ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal" -#ActiveThorns = "PunctureTracker" -#ActiveThorns = "ADManalysis Constants TmunuBase HydroBase " ActiveThorns = "Constants TmunuBase HydroBase " -#ActiveThorns = "QuasiLocalMeasures" -#ActiveThorns = "ADMConstraints" -#ActiveThorns = "EOS_Base EOS_IdealFluid EOS_Polytrope" ActiveThorns = "EOS_Omni" ActiveThorns = "GRHydro GRHydro_InitData" -ActiveThorns = "BLAS LAPACK GSL HDF5" -#ActiveThorns = "TwoPunctures" -#ActiveThorns = "SummationByParts" -#ActiveThorns = "TOVSolver" -#ActiveThorns = "GenericFD NewRad" -#ActiveThorns = "ML_BSSN ML_BSSN_Helper ML_ADMConstraints" -#ActiveThorns = "Hydro_Analysis" -#ActiveThorns = "Dissipation" -#ActiveThorns = "AHFinderDirect" -#ActiveThorns = "WeylScal4 Multipole" -#ActiveThorns = "SetMask_SphericalSurface" +ActiveThorns = "HDF5" ActiveThorns = "SystemStatistics" #------------------------------------------------------------------------------ @@ -79,19 +61,12 @@ #CarpetLib::print_timestats_every = 1024 #CarpetLib::timestat_file = "carpetlib-timing-statistics.asc" -CarpetMask::verbose = "yes" CarpetReduce::verbose = "no" CarpetRegrid2::verbose = "yes" #CarpetRegrid2::veryverbose = "yes" -CarpetTracker::verbose = "yes" - -#Dissipation::verbose = "no" - #NaNChecker::verbose = "all" -#PunctureTracker::verbose = "yes" - TimerReport::out_every = 1024 TimerReport::out_filename = "TimerReport" #TimerReport::output_all_timers = "yes" @@ -100,11 +75,6 @@ #TimerReport::before_checkpoint = "yes" TimerReport::n_top_timers = 40 -#TwoPunctures::verbose = "yes" - -#QuasiLocalMeasures::verbose = "yes" -SphericalSurface::verbose = "yes" - #------------------------------------------------------------------------------ # Utility parameters: #------------------------------------------------------------------------------ @@ -196,13 +166,6 @@ Periodic::periodic_y = "yes" Periodic::periodic_z = "yes" -#ReflectionSymmetry::reflection_x = "no" -#ReflectionSymmetry::reflection_y = "no" -#ReflectionSymmetry::reflection_z = "yes" -#ReflectionSymmetry::avoid_origin_x = "no" -#ReflectionSymmetry::avoid_origin_y = "no" -#ReflectionSymmetry::avoid_origin_z = "no" - # General Carpet parameters: Carpet::domain_from_coordbase = "yes" Carpet::enable_all_storage = "no" @@ -226,8 +189,8 @@ HydroBase::initial_hydro = "advectedloop" GRHydro_InitData::advectedloop_type = "2D" -GRHydro_InitData::advectedloop_case = "V^z/=0" -#GRHydro_InitData::advectedloop_case = "V^z=0" +#GRHydro_InitData::advectedloop_case = "V^z/=0" +GRHydro_InitData::advectedloop_case = "V^z=0" GRHydro_InitData::advectedloop_delA = "Exact" #GRHydro_InitData::advectedloop_delA = "Numeric" @@ -238,11 +201,11 @@ # Numerics: #---------- -TmunuBase::stress_energy_storage = "yes" -TmunuBase::stress_energy_at_RHS = "no" -TmunuBase::timelevels = 1 -TmunuBase::prolongation_type = "none" -TmunuBase::support_old_CalcTmunu_mechanism = "no" +#TmunuBase::stress_energy_storage = "yes" +#TmunuBase::stress_energy_at_RHS = "no" +#TmunuBase::timelevels = 1 +#TmunuBase::prolongation_type = "none" +#TmunuBase::support_old_CalcTmunu_mechanism = "no" HydroBase::timelevels = 2 #HydroBase::timelevels = 3 @@ -258,7 +221,7 @@ GRHydro::method_type = "RSA FV" GRHydro::recon_method = "tvd" -GRHydro::tvd_limiter = "vanleerMC" +GRHydro::tvd_limiter = "vanleerMC2" GRHydro::recon_vars = "primitive" GRHydro::riemann_solver = "HLLE" GRHydro::bound = "none" @@ -269,10 +232,13 @@ GRHydro::GRHydro_MaxNumEvolvedVars = 10 GRHydro::clean_divergence = "yes" -GRHydro::kap_dc = 1.0 +GRHydro::kap_dc = 10.0 GRHydro::track_divB = "yes" GRHydro::calculate_bcom = "yes" +GRHydro::psidcspeed = "light speed" +GRHydro::decouple_normal_Bfield = "yes" + #------------------------------------------------------------------------------ # Output: #------------------------------------------------------------------------------ @@ -284,13 +250,14 @@ IOBasic::outInfo_reductions = "maximum" IOBasic::outInfo_vars = " Carpet::physical_time_per_hour - HydroBase::rho - HydroBase::press GRHydro::bcom_sq + HydroBase::press + HydroBase::rho SystemStatistics::maxrss_mb " -IOScalar::outScalar_every = 1024 +#IOScalar::outScalar_every = 1024 +IOScalar::outScalar_every = 320 #t=1 IOScalar::one_file_per_group = "yes" IOScalar::outScalar_reductions = "minimum maximum average norm1 norm2" IOScalar::outScalar_vars = " @@ -318,14 +285,14 @@ #IOASCII::out3D_ghosts = "no" #hg (DEPRECATED) #IOASCII::out3D_outer_ghosts = "no" #hg (DEPRECATED) -IOASCII::out0D_every = 256 +IOASCII::out0D_every = 320 IOASCII::out0D_vars = " Carpet::timing " #Set these IOASCII options for initial data only: -#IOASCII::out1D_every = 256000000 -IOASCII::out1D_dt = 10 +IOASCII::out1D_every = 7680 +#IOASCII::out1D_dt = 10 IOASCII::out1D_d = "yes" IOASCII::out1D_vars = " GRHydro::bcom @@ -334,6 +301,7 @@ GRHydro::Bcons GRHydro::dens GrHydro::divB + GrHydro::psidc GRHydro::scon GRHydro::tau HydroBase::Bvec @@ -365,12 +333,12 @@ CarpetIOHDF5::one_file_per_group = "yes" CarpetIOHDF5::open_one_input_file_at_a_time = "yes" -CarpetIOHDF5::output_symmetry_points = "no" -CarpetIOHDF5::output_ghost_points = "no" -CarpetIOHDF5::output_boundary_points = "no" -CarpetIOHDF5::output_buffer_points = "no" +CarpetIOHDF5::output_symmetry_points = "yes" +CarpetIOHDF5::output_ghost_points = "yes" +CarpetIOHDF5::output_boundary_points = "yes" +CarpetIOHDF5::output_buffer_points = "yes" -CarpetIOHDF5::out2D_every = 256 +CarpetIOHDF5::out2D_every = 160 CarpetIOHDF5::out2D_xy = "yes" CarpetIOHDF5::out2D_xz = "no" CarpetIOHDF5::out2D_yz = "no" @@ -382,6 +350,7 @@ GRHydro::Bcons GRHydro::dens GRHydro::divB + GrHydro::psidc GRHydro::scon GRHydro::tau Grid::coordinates @@ -393,7 +362,7 @@ HydroBase::w_lorentz " -CarpetIOHDF5::out3D_every = 4096 +CarpetIOHDF5::out3D_every = 1920 CarpetIOHDF5::out3D_vars = " GRHydro::bcom GRHydro::bcom0

1 0

[svn:einsteintoolkit] GRHydro_InitData/trunk/src/ (Rev. 142)
by rhaas＠tapir.caltech.edu 09 Aug '12

09 Aug '12

User: rhaas Date: 2012/08/09 01:26 AM Modified: /trunk/src/ GRHydro_AdvectedLoopM.F90 Log: Revised parameters for advected loop test From: Josh Faber <jfaberrit(a)jafsmamacbook.main.ad.rit.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_AdvectedLoopM.F90 Delta lines: +60 -12 =================================================================== --- trunk/src/GRHydro_AdvectedLoopM.F90 2012-08-09 06:26:42 UTC (rev 141) +++ trunk/src/GRHydro_AdvectedLoopM.F90 2012-08-09 06:26:44 UTC (rev 142) @@ -67,7 +67,7 @@ CCTK_REAL :: range_x,range_y,range_z,range_d CCTK_REAL :: cos_theta, sin_theta, tan_theta CCTK_REAL :: Bvecx_d, Bvecy_d, Bvecz_d - CCTK_REAL :: x_d, y_d, z_d + CCTK_REAL :: x_d, y_d, z_d,diaglength CCTK_REAL :: dx_d, dy_d, dz_d, dx_x, dz_x !!$Adiabatic index for test: @@ -83,19 +83,47 @@ rhoval = 1.0d0 pressval = 3.0d0 + if (CCTK_EQUALS(advectedloop_type,"2D")) then + + !!$ Vx, Vy and Vz values: - if (CCTK_EQUALS(advectedloop_case,"V^z/=0")) then - vxval=0.2d0/sqrt(6.0d0) - vyval=0.5d0*vxval - vzval=vyval - else if (CCTK_EQUALS(advectedloop_case,"V^z=0")) then - vxval=0.2d0/sqrt(6.0d0) - vyval=0.5d0*vxval - vzval=0.0d0 - else - call CCTK_WARN(0,"V^z component case not recognized!") - end if + if (CCTK_EQUALS(advectedloop_case,"V^z/=0")) then + +!!$vxval=0.2d0/sqrt(6.0d0) +!!$ This new choice yields a crossing time of exactly t=24 +!!$ assuming -1<x<1; -0.5<y<0.5 + vxval=1.d0/12.d0 + vyval=0.5d0*vxval + vzval=vyval + + else if (CCTK_EQUALS(advectedloop_case,"V^z=0")) then +!!$ vxval=0.2d0/sqrt(6.0d0) + vxval=1.d0/12.d0 + vyval=0.5d0*vxval + vzval=0.0d0 + else + call CCTK_WARN(0,"V^z component case not recognized!") + end if + else if (CCTK_EQUALS(advectedloop_type,"3D")) then + + vxval=0.2d0*sqrt(2.0d0) + vyval=0.2d0 + + if (CCTK_EQUALS(advectedloop_case,"V^z/=0")) then + +!!$vxval=0.2d0/sqrt(6.0d0) +!!$ This new choice yields a crossing time of exactly t=5 +!!$ assuming -1<x<1; -0.5<y<0.5 + vzval=0.1 + + else if (CCTK_EQUALS(advectedloop_case,"V^z=0")) then + vzval=0.0d0 + else + call CCTK_WARN(0,"V^z component case not recognized!") + end if + + endif nx = cctk_lsh(1) ny = cctk_lsh(2) nz = cctk_lsh(3) @@ -110,7 +138,10 @@ range_z = (cctk_gsh(3)-1)*dz range_d = sqrt(range_z**2+range_x**2) + diaglength = range_x*range_z/range_d +!!$ For 3-d case, the grid is going to be assumed to be a cube + cos_theta = range_z/range_d sin_theta = range_x/range_d tan_theta = sin_theta/cos_theta @@ -162,16 +193,33 @@ else if (CCTK_EQUALS(advectedloop_type,"3D")) then + +!!$ tangential velocity should be parallel to (1,1,1) plus a normal component (-1,0,1) + velx(i,j,k)=cos_theta*vxval-sin_theta*vzval vely(i,j,k)=vyval velz(i,j,k)=cos_theta*vzval+sin_theta*vxval Bvecz_d=0.0d0 +!!$ x_d = (x+z)/sqrt(2) => x_d=0 is equivalent to x+z=0 + x_d = cos_theta*x(i,j,k)+sin_theta*z(i,j,k) y_d = y(i,j,k) z_d = cos_theta*z(i,j,k)-sin_theta*x(i,j,k) +!!$ need to make x_d periodic! + + if(x_d.gt.1.5*diaglength) then + x_d=x_d-2.0*diaglength + else if (x_d.gt.0.5*diaglength .and. x_d.lt.1.5*diaglength) then + x_d=x_d-diaglength + else if(x_d.lt.-1.5*diaglength) then + x_d=x_d+2.0*diaglength + else if (x_d.lt.(-0.5*diaglength) .and. x_d.gt.(-1.5*diaglength)) then + x_d=x_d+diaglength + endif + radius = sqrt(x_d**2+y_d**2) if (CCTK_EQUALS(advectedloop_delA,"Exact")) then

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits