Commits August 2012

commits@lists.einsteintoolkit.org

7 participants
71 discussions

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

08 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 08 Aug '12

08 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 08 Aug '12

08 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

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

08 Aug '12

User: rhaas Date: 2012/08/09 01:26 AM Modified: /trunk/par/ rotor.par Log: GRHydro_InitData: adapt rotor parfile to current code From: Roland Haas <roland.haas(a)physics.gatech.edu> File Changes: Directory: /trunk/par/ ====================== File [modified]: rotor.par Delta lines: +60 -27 =================================================================== --- trunk/par/rotor.par 2012-08-09 06:26:40 UTC (rev 140) +++ trunk/par/rotor.par 2012-08-09 06:26:42 UTC (rev 141) @@ -1,15 +1,16 @@ ## -ActiveThorns = "time MoL cartgrid3d carpetioascii ioutil Fortran boundary hydrobase grhydro grhydro_initdata nanchecker eos_omni CoordBase LoopControl CarpetLib Carpet CarpetReduce spacemask ADMCoupling ADMMacros ADMBASE STATICCONFORMAL TMUNUBASE SymBase localreduce CarpetIOBasic initbase " +ActiveThorns = "time MoL cartgrid3d carpetioascii ioutil boundary hydrobase grhydro grhydro_initdata nanchecker eos_omni CoordBase LoopControl CarpetLib Carpet CarpetReduce spacemask ADMCoupling ADMMacros ADMBASE STATICCONFORMAL TMUNUBASE SymBase localreduce CarpetIOBasic initbase " ActiveThorns = "CarpetIOHDF5" -#This is a hack. ghost_size should be the same as grhydro::grhydro_stencil. -driver::ghost_size=3 -grhydro::grhydro_stencil=3 +driver::ghost_size=2 # for TVD +grhydro::grhydro_stencil=2 -time::dtfac = 0.25 +Time::dtfac = 0.25 # 0.5 see p 409 (0.25 from David Neilsen, Eric WHirschmann, and R Steven Millward arXiv:gr-qc/0512147v1) -methodoflines::ODE_Method = "rk2" -methodoflines::MoL_Intermediate_Steps=2 +# rk2 dies quite quickly +methodoflines::ODE_Method = "rk3" +methodoflines::MoL_Intermediate_Steps=3 +#methodoflines::MoL_Num_Scratch_Levels=1 ADMBase::initial_shift = "zero" @@ -19,11 +20,16 @@ HydroBase::evolution_method="grhydro" HydroBase::Bvec_evolution_method="grhydro" +# rotor problem from DelZanna +# needs to be before EOS_Omni +grhydro_initdata::rotor_use_smoothing = no + SpaceMask::use_mask="yes" -EOS_Omni::gl_gamma=2.0 +EOS_Omni::gl_gamma=5./3. # relativistic rotor problem of DelZanna & Buccianti & Londrille AaA_400_397_2003 +# try to mimic the core fluid +EOS_Omni::poly_k = "grhydro_initdata::rotor_pressin/grhydro_initdata::rotor_rhoin**EOS_Omni::gl_gamma" +EOS_Omni::poly_gamma="EOS_Omni::gl_gamma" grhydro::recon_method="tvd" -grhydro::method_type="RSA FV" -grhydro::recon_vars="primitive" grhydro::riemann_solver="HLLE" grhydro::tvd_limiter="vanleerMC" grhydro::bound="flat" @@ -31,29 +37,60 @@ grhydro::grhydro_eos_table="Ideal_Fluid" grhydro::Grhydro_MaxNumConstrainedVars=13 grhydro::GRHydro_MaxNumEvolvedVars=10 +grhydro::GRHydro_MaxNumSandRVars=0 # ADMBase is not evolved +grhydro::rho_abs_min = 1e-14 # need some reasonable value here, even though no density ever gets this low. This is about the same as GRHydro's default value from rho_central and rho_rel_min grhydro::clean_divergence = "yes" -grhydro::kap_dc = 1.0 +grhydro::transport_constraints = "no" +grhydro::decouple_normal_Bfield = "yes" +grhydro::kap_dc = 5.0 grhydro::track_divB = "yes" -grid::type = "BySpacing" +grid::type = "coordbase" grid::domain = "full" -grid::dxyz = 0.004 -driver::global_nx = 250 -driver::global_ny = 250 -driver::global_nz = 10 +CoordBase::domainsize = "minmax" +CoordBase::xmin = 0. +CoordBase::ymin = 0. +CoordBase::zmin = -0.0125 +CoordBase::xmax = 1.0 +CoordBase::ymax = 1.0 +CoordBase::zmax = 0.0125 +CoordBase::boundary_size_x_lower = carpet::ghost_size +CoordBase::boundary_size_y_lower = carpet::ghost_size +CoordBase::boundary_size_z_lower = carpet::ghost_size +CoordBase::boundary_size_x_upper = carpet::ghost_size +CoordBase::boundary_size_y_upper = carpet::ghost_size +CoordBase::boundary_size_z_upper = carpet::ghost_size + +CoordBase::spacing = "numcells" +CoordBase::ncells_x = 400 +CoordBase::ncells_y = 400 +CoordBase::ncells_z = 10 + +carpet::domain_from_coordbase = "yes" + Cactus::terminate = "time" -Cactus::cctk_final_time = 0.4 -#Cactus::terminate = "iteration" -#cactus::cctk_itlast = 2 +Cactus::cctk_final_time = 0.4 # see DelZanna p 410 IO::out_dir = $parfile -CarpetIOBasic::outInfo_vars="hydrobase::rho Hydrobase::bvec GRHydro::bcons" +CarpetIOBasic::outInfo_vars=" + hydrobase::rho + hydrobase::press + hydrobase::eps + hydrobase::vel[0] hydrobase::vel[1] + hydrobase::w_lorentz + GRHydro::bcons[0] GRHydro::bcons[1] + GRHydro::divB GRHydro::psidc +" CarpetIOBasic::outInfo_every=1 -CarpetIOASCII::out1D_every = 20 -CarpetIOASCII::out1D_vars = "HydroBase::rho HydroBase::press HydroBase::eps HydroBase::vel grhydro::dens grhydro::tau grhydro::scon HydroBase::Bvec grhydro::Bcons" +CarpetIOASCII::out1D_every = 10 +CarpetIOASCII::out1D_z = "no" +CarpetIOASCII::out1D_d = "no" +CarpetIOASCII::out1D_xline_y = 0.5 +CarpetIOASCII::out1D_yline_x = 0.5 +CarpetIOASCII::out1D_vars = "HydroBase::rho HydroBase::press HydroBase::eps HydroBase::vel grhydro::dens grhydro::tau grhydro::scon HydroBase::Bvec grhydro::Bcons grhydro::divB grhydro::psidc" CarpetIOHDF5::one_file_per_group = "yes" CarpetIOHDF5::open_one_input_file_at_a_time = "yes" @@ -63,9 +100,6 @@ CarpetIOHDF5::out2D_yz = "no" CarpetIOHDF5::out2D_xyplane_z = 0.0 CarpetIOHDF5::out2D_vars = " -# CarpetReduce::weight -# GRHydro::GRHydro_tracers -# Grid::coordinates GRHydro::Bcons GrHydro::dens GrHydro::divB @@ -77,6 +111,5 @@ HydroBase::rho HydroBase::vel HydroBase::w_lorentz + GRHydro::psidc " - -

1 0

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

08 Aug '12

User: rhaas Date: 2012/08/09 01:26 AM Modified: /trunk/src/ GRHydro_Bondi.c, GRHydro_BondiM.c Log: GRHydro_InitData: change 2d array access to avoid compiler warnings * don't access 2d array as 1d array through second index * always assign something to rho_guess in Bondi solution From: Roland Haas <roland.haas(a)physics.gatech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Bondi.c Delta lines: +29 -12 =================================================================== --- trunk/src/GRHydro_Bondi.c 2012-08-09 06:26:39 UTC (rev 139) +++ trunk/src/GRHydro_Bondi.c 2012-08-09 06:26:40 UTC (rev 140) @@ -123,13 +123,17 @@ for Kerr-Schild coordinates. ******************************************************************************/ -static void dxc_dxs_ks_calc(CCTK_REAL *x_cart, CCTK_REAL *x_spher, CCTK_REAL dxc_dxs[][NDIM] ) +static void dxc_dxs_ks_calc(CCTK_REAL *x_cart, CCTK_REAL *x_spher, CCTK_REAL dxc_dxs[NDIM][NDIM] ) { - int i; + int i, j; CCTK_REAL r, th, ph; CCTK_REAL sth,cth,sph,cph; - for( i = 0 ; i < NDIM*NDIM ; i++ ) { dxc_dxs[0][i] = 0. ; } + for( i = 0 ; i < NDIM ; i++ ) { + for( j = 0 ; j < NDIM ; j++ ) { + dxc_dxs[0][i] = 0. ; + } + } r = x_spher[RR]; th = x_spher[TH]; @@ -165,13 +169,17 @@ for Boyer-Lindquist coordinates. ******************************************************************************/ -static void dxc_dxs_bl_calc(CCTK_REAL *x_cart, CCTK_REAL *x_spher, CCTK_REAL dxc_dxs[][NDIM] ) +static void dxc_dxs_bl_calc(CCTK_REAL *x_cart, CCTK_REAL *x_spher, CCTK_REAL dxc_dxs[NDIM][NDIM] ) { - int i; + int i, j; CCTK_REAL r, th, ph, rterm, dr; CCTK_REAL sth,cth,sph,cph; - for( i = 0 ; i < NDIM*NDIM ; i++ ) { dxc_dxs[0][i] = 0. ; } + for( i = 0 ; i < NDIM ; i++ ) { + for( j = 0 ; j < NDIM ; j++ ) { + dxc_dxs[j][i] = 0. ; + } + } r = x_spher[RR]; th = x_spher[TH]; @@ -210,13 +218,17 @@ for "Isotropic" coordinates. ******************************************************************************/ -static void dxc_dxs_iso_calc(CCTK_REAL *x_cart, CCTK_REAL *x_spher, CCTK_REAL dxc_dxs[][NDIM] ) +static void dxc_dxs_iso_calc(CCTK_REAL *x_cart, CCTK_REAL *x_spher, CCTK_REAL dxc_dxs[NDIM][NDIM] ) { - int i; + int i, j; CCTK_REAL th, ph,r_iso; CCTK_REAL sth,cth,sph,cph; - for( i = 0 ; i < NDIM*NDIM ; i++ ) { dxc_dxs[0][i] = 0. ; } + for( i = 0 ; i < NDIM ; i++ ) { + for( j = 0 ; j < NDIM ; j++ ) { + dxc_dxs[j][i] = 0. ; + } + } /* BL spherical coordinates : */ th = x_spher[TH]; @@ -355,13 +367,17 @@ --------------- -- Covariant Kerr metric in Boyer-Lindquist coordinates. ****************************************************************************/ -static void bl_gcov_func( CCTK_REAL *x, CCTK_REAL gcov[][NDIM]) +static void bl_gcov_func( CCTK_REAL *x, CCTK_REAL gcov[NDIM][NDIM]) { int i,j,k ; CCTK_REAL sth,cth,s2,r2,DD,mu ; CCTK_REAL r,th; - for( i = 0 ; i < NDIM*NDIM ; i++ ) { gcov[0][i] = 0. ; } + for( i = 0 ; i < NDIM ; i++ ) { + for( j = 0 ; j < NDIM ; j++ ) { + gcov[i][j] = 0. ; + } + } r = x[RR]; th = x[TH]; @@ -648,7 +664,8 @@ if( *rho < 0. ) { if( r > 0.9*rs && r < 1.1*rs ) { - *rho = rhos; + rho_guess = rhos; + *rho = rho_guess; } else { // rhotmp = (sqrt(Qdot) - 1.) * (gamma_eos - 1.) / ( gamma_eos * K ); File [modified]: GRHydro_BondiM.c Delta lines: +28 -11 =================================================================== --- trunk/src/GRHydro_BondiM.c 2012-08-09 06:26:39 UTC (rev 139) +++ trunk/src/GRHydro_BondiM.c 2012-08-09 06:26:40 UTC (rev 140) @@ -125,11 +125,15 @@ ******************************************************************************/ static void dxc_dxs_ks_calc(CCTK_REAL *x_cart, CCTK_REAL *x_spher, CCTK_REAL dxc_dxs[][NDIM] ) { - int i; + int i, j; CCTK_REAL r, th, ph; CCTK_REAL sth,cth,sph,cph; - for( i = 0 ; i < NDIM*NDIM ; i++ ) { dxc_dxs[0][i] = 0. ; } + for( i = 0 ; i < NDIM ; i++ ) { + for( j = 0 ; j < NDIM ; j++ ) { + dxc_dxs[j][i] = 0. ; + } + } r = x_spher[RR]; th = x_spher[TH]; @@ -165,13 +169,17 @@ for Boyer-Lindquist coordinates. ******************************************************************************/ -static void dxc_dxs_bl_calc(CCTK_REAL *x_cart, CCTK_REAL *x_spher, CCTK_REAL dxc_dxs[][NDIM] ) +static void dxc_dxs_bl_calc(CCTK_REAL *x_cart, CCTK_REAL *x_spher, CCTK_REAL dxc_dxs[NDIM][NDIM] ) { - int i; + int i, j; CCTK_REAL r, th, ph, rterm, dr; CCTK_REAL sth,cth,sph,cph; - for( i = 0 ; i < NDIM*NDIM ; i++ ) { dxc_dxs[0][i] = 0. ; } + for( i = 0 ; i < NDIM ; i++ ) { + for( j = 0 ; j < NDIM ; j++ ) { + dxc_dxs[j][i] = 0. ; + } + } r = x_spher[RR]; th = x_spher[TH]; @@ -210,13 +218,17 @@ for "Isotropic" coordinates. ******************************************************************************/ -static void dxc_dxs_iso_calc(CCTK_REAL *x_cart, CCTK_REAL *x_spher, CCTK_REAL dxc_dxs[][NDIM] ) +static void dxc_dxs_iso_calc(CCTK_REAL *x_cart, CCTK_REAL *x_spher, CCTK_REAL dxc_dxs[NDIM][NDIM] ) { - int i; + int i, j; CCTK_REAL th, ph,r_iso; CCTK_REAL sth,cth,sph,cph; - for( i = 0 ; i < NDIM*NDIM ; i++ ) { dxc_dxs[0][i] = 0. ; } + for( i = 0 ; i < NDIM ; i++ ) { + for( j = 0 ; j < NDIM ; j++ ) { + dxc_dxs[j][i] = 0. ; + } + } /* BL spherical coordinates : */ th = x_spher[TH]; @@ -355,13 +367,17 @@ --------------- -- Covariant Kerr metric in Boyer-Lindquist coordinates. ****************************************************************************/ -static void bl_gcov_func( CCTK_REAL *x, CCTK_REAL gcov[][NDIM]) +static void bl_gcov_func( CCTK_REAL *x, CCTK_REAL gcov[NDIM][NDIM]) { int i,j,k ; CCTK_REAL sth,cth,s2,r2,DD,mu ; CCTK_REAL r,th; - for( i = 0 ; i < NDIM*NDIM ; i++ ) { gcov[0][i] = 0. ; } + for( i = 0 ; i < NDIM ; i++ ) { + for( j = 0 ; j < NDIM ; j++ ) { + gcov[i][j] = 0. ; + } + } r = x[RR]; th = x[TH]; @@ -648,7 +664,8 @@ if( *rho < 0. ) { if( r > 0.9*rs && r < 1.1*rs ) { - *rho = rhos; + rho_guess = rhos; + *rho = rho_guess; } else { // rhotmp = (sqrt(Qdot) - 1.) * (gamma_eos - 1.) / ( gamma_eos * K );

1 0

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

08 Aug '12

User: rhaas Date: 2012/08/09 01:26 AM Modified: /trunk/ param.ccl /trunk/src/ GRHydro_RotorM.F90 Log: GRHydro_InitData: make rotor parameters Cactus parameters * cite paper with description of system * reverse sense of rotation * add openmp statements From: Roland Haas <roland.haas(a)physics.gatech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_RotorM.F90 Delta lines: +56 -74 =================================================================== --- trunk/src/GRHydro_RotorM.F90 2012-08-09 06:26:37 UTC (rev 138) +++ trunk/src/GRHydro_RotorM.F90 2012-08-09 06:26:39 UTC (rev 139) @@ -44,93 +44,82 @@ subroutine GRHydro_RotorM(CCTK_ARGUMENTS) + use cctk + implicit none - + DECLARE_CCTK_ARGUMENTS DECLARE_CCTK_PARAMETERS DECLARE_CCTK_FUNCTIONS - + CCTK_INT :: i, j, k, nx, ny, nz - CCTK_REAL :: radius, det, v_max - CCTK_REAL :: rhoin, rhoout, pressin, pressout - CCTK_REAL :: bvcxl,bvcyl,bvczl - CCTK_REAL :: tmp,tmp2,gam,r_rot + CCTK_REAL :: radius, det, rfact - CCTK_INT :: use_smoothing - CCTK_REAL :: rsmooth_rel,rfact + !begin EOS_Omni stuff + CCTK_REAL :: tempEOS(1), yeEOS(1), xepsEOS(1) + CCTK_INT :: keyerr(1), anyerr + CCTK_INT, parameter :: have_temp = 0, n = 1 + !end EOS_Omni -!!$Hardwired bfield and gamma values - bvcxl = 1.0 - bvcyl = 0.0 - bvczl = 0.0 - -!!$Adiabatic index for test - gam = (5.d0/3.d0) - -!!$radius of rotor is 0.1 - r_rot = 0.1 - -!!$Maximum velocity is 0.995 - v_max = 0.995 - -!!$Inner values - rhoin = 10.d0 - pressin = 1.d0 - -!!$Outer values - rhoout = 1.d0 - pressout = 1.d0 - - -!!$ Smooth the edge? Define the radius in relative terms if so - use_smoothing = 0 - rsmooth_rel = 0.05 - nx = cctk_lsh(1) ny = cctk_lsh(2) nz = cctk_lsh(3) - - if(v_max>1.0) then - write(*,*)"No superluminal speeds!!!! v_max=",v_max," reset to 0.995" - v_max=0.995 - endif + !$OMP PARALLEL DO PRIVATE(i,j,k,radius,det,rfact,tempEOS,yeEOS,xepsEOS,keyerr,anyerr) & + !$OMP default(none) & + !$OMP firstprivate(nx,ny,nz,GRHydro_eos_handle, GRHydro_eos_rf_prec, grhydro_eos_type, & + !$OMP rotor_xc,rotor_yc, rotor_rsmooth_rel, rotor_use_smoothing, & + !$OMP rotor_r_rot, rotor_rhoin,rotor_pressin,rotor_rhoout,rotor_pressout, & + !$OMP rotor_v_max, rotor_bvcxl, rotor_bvcyl, rotor_bvczl) & + !$OMP shared(x,y,z,rho,press,eps,vel,Bvec,w_lorentz,dens,Scon,tau,Bcons, & + !$OMP gxx,gxy,gxz,gyy,gyz,gzz) do i=1,nx do j=1,ny do k=1,nz - - radius = sqrt(x(i,j,k)**2+y(i,j,k)**2) - - if(radius.le.r_rot) then - rho(i,j,k) = rhoin - press(i,j,k) = pressin - velx(i,j,k) = v_max/r_rot*y(i,j,k) - vely(i,j,k) = -1.d0*v_max/r_rot*x(i,j,k) + + radius = sqrt((x(i,j,k)-rotor_xc)**2+(y(i,j,k)-rotor_yc)**2) + + if(radius.le.rotor_r_rot) then + rho(i,j,k) = rotor_rhoin + press(i,j,k) = rotor_pressin + velx(i,j,k) = -1.d0*rotor_v_max/rotor_r_rot*(y(i,j,k)-rotor_yc) + vely(i,j,k) = +1.d0*rotor_v_max/rotor_r_rot*(x(i,j,k)-rotor_xc) velz(i,j,k) = 0.d0 - else if((use_smoothing.eq.1) .and. & - ((radius.gt.r_rot) .and. & - (radius.le.((1.0+rsmooth_rel)*r_rot)))) then - rfact = (radius/r_rot - 1.d0) / rsmooth_rel - rho(i,j,k) = rfact*rhoout + (1.d0-rfact)*rhoin - press(i,j,k) = rfact*pressout + (1.d0-rfact)*pressin - velx(i,j,k) = (1.d0-rfact)*v_max * y(i,j,k) / radius - velx(i,j,k) = -1.d0*(1.d0-rfact)*v_max * x(i,j,k) / radius + else if((rotor_use_smoothing.eq.1) .and. & + ((radius.gt.rotor_r_rot) .and. & + (radius.le.((1.0+rotor_rsmooth_rel)*rotor_r_rot)))) then + rfact = (radius/rotor_r_rot - 1.d0) / rotor_rsmooth_rel + rho(i,j,k) = rfact*rotor_rhoout + (1.d0-rfact)*rotor_rhoin + press(i,j,k) = rfact*rotor_pressout + (1.d0-rfact)*rotor_pressin + velx(i,j,k) = -1.d0*(1.d0-rfact)*rotor_v_max * (y(i,j,k)-rotor_yc) / radius + velx(i,j,k) = +1.d0*(1.d0-rfact)*rotor_v_max * (x(i,j,k)-rotor_xc) / radius velz(i,j,k) = 0.d0 else - rho(i,j,k) = rhoout - press(i,j,k) = pressout + rho(i,j,k) = rotor_rhoout + press(i,j,k) = rotor_pressout velx(i,j,k) = 0.d0 vely(i,j,k) = 0.d0 velz(i,j,k) = 0.d0 endif - eps(i,j,k)=press(i,j,k)/(gam-1.d0)/rho(i,j,k) - - Bvecx(i,j,k)=bvcxl - Bvecy(i,j,k)=bvcyl - Bvecz(i,j,k)=bvczl - + keyerr = 0; anyerr = 0 + call EOS_Omni_EpsFromPress(GRHydro_eos_handle, have_temp, GRHydro_eos_rf_prec, & + n, rho(i:i,j:j,k:k), eps(i:i,j:j,k:k), & + tempEOS, yeEOS, press(i:i,j:j,k:k), & + eps(i:i,j:j,k:k), & ! xeps argument + keyerr, anyerr) + if(anyerr .ne. 0) then + call CCTK_WARN(0, "Error when calling EOS. After stopping swearing at me, add a decent output text.") + end if + + Bvecx(i,j,k)=rotor_bvcxl + Bvecy(i,j,k)=rotor_bvcyl + Bvecz(i,j,k)=rotor_bvczl + 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(abs(det - 1d0) .gt. 1e-8) then + call CCTK_WARN(0, "Rotor initial data only supports flat spacetime right now") + end if + if (CCTK_EQUALS(GRHydro_eos_type,"Polytype")) then call Prim2ConPolyM(GRHydro_eos_handle,gxx(i,j,k),gxy(i,j,k),& gxz(i,j,k),gyy(i,j,k),gyz(i,j,k),gzz(i,j,k),& @@ -148,19 +137,12 @@ eps(i,j,k),press(i,j,k),Bvecx(i,j,k),Bvecy(i,j,k),Bvecz(i,j,k),& w_lorentz(i,j,k)) end if - + enddo enddo enddo - - densrhs = 0.d0 - srhs = 0.d0 - taurhs = 0.d0 - Bconsrhs = 0.d0 + !$OMP END PARALLEL DO - - return - end subroutine GRHydro_RotorM Directory: /trunk/ ================== File [modified]: param.ccl Delta lines: +70 -1 =================================================================== --- trunk/param.ccl 2012-08-09 06:26:37 UTC (rev 138) +++ trunk/param.ccl 2012-08-09 06:26:39 UTC (rev 139) @@ -14,7 +14,7 @@ "simple_wave" :: "Set initial data from Anile Miller Motta, Phys.Fluids. 26, 1450 (1983)" "monopole" :: "Monopole at the center" "cylexp" :: "Cylindrical Explosion" - "rotor" :: "Magnetic Rotor test" + "rotor" :: "Magnetic Rotor test from DelZanna,Bucciantini, and Londrillo A&A 400, 397â413 (2003)" "advectedloop":: "Magnetic advected loop test" "alfvenwave" :: "Circularly polarized Alfven wave" "hydro_bondi_solution" :: "Spherical single black hole Bondi solution" @@ -275,6 +275,75 @@ (0:* :: "Anything positive." } 1.0e-3 +# for the Magnetic Rotor test: +# default values are from DelZanna,Bucciantini, and Londrillo A&A 400, 397â413 (2003) though notation differs + +CCTK_REAL rotor_xc "center of rotation" +{ + *:* :: "Any location" +} 0.5 + +CCTK_REAL rotor_yc "center of rotation" +{ + *:* :: "Any location" +} 0.5 + +CCTK_REAL rotor_bvcxl "intial component of Bvec[0]" +{ + *:* :: "any real number" +} 1.0 + +CCTK_REAL rotor_bvcyl "intial component of Bvec[1]" +{ + *:* :: "any real number" +} 0.0 + +CCTK_REAL rotor_bvczl "intial component of Bvec[2]" +{ + *:* :: "any real number" +} 0.0 + +CCTK_REAL rotor_r_rot "radius of rotor" +{ + (0:* :: "any positive number" +} 0.1 + +CCTK_REAL rotor_v_max "Maximum velocity" +{ + (-1:1) :: "any subluminal speed (negative is clockwise)" +} 0.995 + + +CCTK_REAL rotor_rhoin "initial density inside rotor" +{ + (0:* :: "any positive number" +} 10.d0 + +CCTK_REAL rotor_pressin "initial pressure inside rotor" +{ + (0:* :: "any positive number" +} 1.d0 + +CCTK_REAL rotor_rhoout "initial density outside rotor" +{ + (0:* :: "any positive number" +} 1.d0 + +CCTK_REAL rotor_pressout "initial pressure outside rotor" +{ + (0:* :: "any positive number" +} 1.d0 + +CCTK_BOOLEAN rotor_use_smoothing "Smooth the edge?" +{ +} yes + +CCTK_REAL rotor_rsmooth_rel "Define the radius in relative terms if so" +{ + (0:* :: "any positive number" +} 0.05 + + shares:GRHydro USES real GRHydro_eos_rf_prec

1 0

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

08 Aug '12

User: rhaas Date: 2012/08/09 01:26 AM Added: /trunk/ configuration.ccl Log: GRHydro_InitData: REQUIRE EOS_Omni From: Roland Haas <roland.haas(a)physics.gatech.edu> File Changes: Directory: /trunk/ ================== File [added]: configuration.ccl Delta lines: +2 -0 =================================================================== --- trunk/configuration.ccl (rev 0) +++ trunk/configuration.ccl 2012-08-09 06:26:37 UTC (rev 138) @@ -0,0 +1,2 @@ +REQUIRES EOS_Omni +

1 0

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

08 Aug '12

User: rhaas Date: 2012/08/09 01:21 AM Modified: /trunk/src/ GRHydro_PoloidalMagFieldM.F90 Log: A couple of changes: 1) Calculate divergence of B in MoL PseudoEvolution and set its Prolongation="Restrict". 2) Correct the source terms and fluxes for the Bfield and the divergence cleaning field when having a non-flat space-time. 3) Make sure alpha factors match between UpdateCalculation and fluxes definition. 4) Include 1/sqrt(detg) factor in calculation of \epsilon^{\muijk} in the cross product to obtain the Bfield form the vector potential. From: Philipp Moesta <pmoesta(a)tapir.caltech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_PoloidalMagFieldM.F90 Delta lines: +7 -4 =================================================================== --- trunk/src/GRHydro_PoloidalMagFieldM.F90 2012-08-09 06:11:31 UTC (rev 136) +++ trunk/src/GRHydro_PoloidalMagFieldM.F90 2012-08-09 06:21:04 UTC (rev 137) @@ -58,6 +58,7 @@ CCTK_INT :: i, j, k, nx, ny, nz CCTK_REAL :: det + CCTK_REAL :: sdet CCTK_REAL :: dx,dy,dz CCTK_REAL :: rhofac, delPcut, maxP_Pcut CCTK_REAL :: Aphi, Ax, Ay, Az @@ -112,12 +113,14 @@ Ay_dx = Aphi + x(i,j,k)*Aphi_dx Ay_dz = x(i,j,k)*Aphi_dz - Bvecx(i,j,k) = alp(i,j,k)*Ay_dz - Bvecy(i,j,k) = -alp(i,j,k)*Ax_dz - Bvecz(i,j,k) = alp(i,j,k)*(Ax_dy-Ay_dx) 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)) - + sdet = sqrt(det) + + Bvecx(i,j,k) = Ay_dz/sdet + Bvecy(i,j,k) = - Ax_dz/sdet + Bvecz(i,j,k) = (Ax_dy-Ay_dx)/sdet + if (CCTK_EQUALS(GRHydro_eos_type,"Polytype")) then call Prim2ConPolyM(GRHydro_eos_handle,gxx(i,j,k),gxy(i,j,k),& gxz(i,j,k),gyy(i,j,k),gyz(i,j,k),gzz(i,j,k),&

1 0

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

08 Aug '12

User: rhaas Date: 2012/08/09 01:11 AM Modified: /trunk/src/ GRHydro_PoloidalMagFieldM.F90 Log: Fix a sign error in the poloidal mag field idata From: Philipp Moesta <pmoesta(a)tapir.caltech.edu> File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_PoloidalMagFieldM.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_PoloidalMagFieldM.F90 2012-08-09 06:10:15 UTC (rev 135) +++ trunk/src/GRHydro_PoloidalMagFieldM.F90 2012-08-09 06:11:31 UTC (rev 136) @@ -108,7 +108,7 @@ ( rhofac*press_dz/delPcut - poloidal_n_p*rho_dz/poloidal_rho_max ) Ax_dy = -Aphi - y(i,j,k)*Aphi_dy - Ax_dz = y(i,j,k)*Aphi_dz + Ax_dz = -y(i,j,k)*Aphi_dz Ay_dx = Aphi + x(i,j,k)*Aphi_dx Ay_dz = x(i,j,k)*Aphi_dz

1 0

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

08 Aug '12

User: rhaas Date: 2012/08/09 01:10 AM Modified: /trunk/ schedule.ccl /trunk/par/ balsara1d.par /trunk/src/ GRHydro_ShockTubeM.F90 Log: Introduce a better treatment of the shocktube boundaries. * The 2d diagonal boundary conditions are correctly applied when using multiprocessors. * The ymin (ymax) surfaces have their points in the region indsum < minsum (indsum > maxsum) fixed after a first sync operation. A second sync guarantees the neighbour grid receives the correct value of the boundary on those regions. * Mirror symmetries applied at the faces intercepting the upper right corner assumed cubic local domain. However, rectangular local domains are common after domain decomposition. This commit fixes the indexes applying this symmetry for these cases. * Upper right corner/edge/surface condition should use sten instead of stenp1. Thanks for catching this Josh! * Attempt at a 3d diagonal boundary treatment (Josh's commit) File Changes: Directory: /trunk/par/ ====================== File [modified]: balsara1d.par Delta lines: <Unavailable> Directory: /trunk/src/ ====================== File [modified]: GRHydro_ShockTubeM.F90 Delta lines: <Unavailable> Directory: /trunk/ ================== File [modified]: schedule.ccl Delta lines: <Unavailable>

1 0

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits August 2012