parameters for hlc_graxi ######################################################################## # # # id_t sample initial data file to graxi (unigrid/no rotation) # # # # To construct otherwise identical ID files, run # # "mi id_t" which will create files id_td0, id_td1, etc # # which then differ only by the resolution. # # # ######################################################################## iter := 10000 output := 0-*/100 #====================================================================== # Matter Parameters #====================================================================== #====================================================================== # scalar field 1. 2 pieces in initial profile. # # The sinnX parameters multiply the distributions by sin(theta)^sinnX. # if ic_v1, or ic_v1_b !=0, then the flat-space scalar field solution # is Lorentz boosted along the z axis with the given velocity. #====================================================================== amp1 := 0.0 r01 := 2.0 delta1 := 1.0 # +1 --- in-going # 0 --- time-symmetric # -1 --- out-going signum1 := 0 rho_center1 := 0 z_center1 := 0.0 cyl1 := 0.75 sinn1 := 0 ic_v1 := 0 #====================================================================== # second part of scalar field profile #====================================================================== amp1_b := 0.0 r01_b := 0.0 delta1_b := 3.0 # +1 --- in-going # 0 --- time-symmetric # -1 --- out-going signum1_b := 0 rho_center1_b := 0 z_center1_b := 0.0 cyl1_b := 1 sinn1_b := 0 ic_v1_b := 0 #====================================================================== # complex scalar field. 2 pieces in initial profile. parameters # as with real s.f., except specify real and imaginary components # seperately, and can give a mass term (mass4) # and parameters in a "mexican hat" potential (lambda4, eta4). # # V = lambda4*(phi4*phi4_CC-eta4^2)^2+mass4^2*(phi4*phi4_CC) # # where phi4 = phi4_r + I*phi4_i # phi4_CC = phi4_r - I*phi4_i # #====================================================================== amp4_r := 0.03 r04_r := 0.0 delta4_r := 1.0 signum4_r := 0 rho_center4_r := 0 z_center4_r := 0.5 cyl4_r := 1.0 sinn4_r := 0 amp4_i := 0.03 r04_i := 2.0 delta4_i := 1.0 signum4_i := -1 rho_center4_i := 0 z_center4_i := -0.5 cyl4_i := 0.75 sinn4_i := 0 mass4 := 1 lambda4 := 0 eta4 := 0 W4 := 0.953910279 #====================================================================== # second part of complex scalar field profile #====================================================================== amp4_r_b := 0 r04_r_b := 0.0 delta4_r_b := 1.0 signum4_r_b := 1 rho_center4_r_b := 0 z_center4_r_b := 0.0 cyl4_r_b := 1 sinn4_r_b := 0 amp4_i_b := 0 r04_i_b := 0.0 delta4_i_b := 1.0 signum4_i_b := 0 rho_center4_i_b := 0 z_center4_i_b := 0.0 cyl4_i_b := 1 sinn4_i_b := 0 #====================================================================== # # Geometry Parameters # #====================================================================== #====================================================================== # sigma #====================================================================== ampsi := 0.06 r0si := 0.0 deltasi := 1.0 signumsi := 0 rho_centersi := 0.0 z_centersi := 0.0 cylsi := 1 sinnsi := 0 #====================================================================== # omega #====================================================================== ampo := 0.00000 r0o := 7.0 deltao := 3.0 signumo := 0 rho_centero := 0.0 z_centero := 0.0 cylo := 1 sinno := 0 #====================================================================== # Tracing level: # # > 0 ---- Canonical CN loop output # > 1 ---- Canonical MG tracing # > 2 ---- Output relativized errors in CN loop # > 3 ---- parameter dumping # > 4 ---- level output in MG #====================================================================== trace_level := 1 #====================================================================== # General Parameters #====================================================================== MAXITER := 50 maxstep := 500 Nrho0 := 32 Nz0 := 64 rhomin := 0 rhomax := 24 zmin := -24 zmax := 24 ser := 0 fout := 1 epsiter := 1.0e-3 ic_epsiter := 1.0e-5 tag := "t_b24_" lambda := .3 epsdis := 0.5 #====================================================================== # Specify form of Initial Data: # IDATA: # 0 ---- general Gaussians with parameters set below # 1 ---- specified ansatz, and computed source terms # 2 ---- load in scalar field profiles #====================================================================== IDATA := 0 #====================================================================== # Restrict model to Flat Space? # (if one, then metric restricted to flat metric) #====================================================================== flat_space := 0 #====================================================================== # force_cont_off == 1 do *not* use continuation methods in ID # != 1 do use them #====================================================================== force_cont_off := 0 #====================================================================== # Output initial time sdf's? (if non-zero) #====================================================================== output_init := 0 #======================================================================= # to evolve psi using the evolution equation #======================================================================= use_evo_for_psi := 1 #======================================================================= # An experimental 'over-relaxation' type parameter, scaling the # residual added to the rhs's on coarser grids in MG (though not to # try to speed things up as much as to get more stable evolution, # with w<1) #======================================================================= MG_w := 1 #======================================================================= # for excision: (currently not supported) #======================================================================= do_excise := 0 #======================================================================= #parameters for 'apparent horizon finder' #apph_search_dt --- how often to search (set <0 to disable) # # If (apph_find_a=1), then, if an AH is found, the parameter f_a1 # is modified (up to a maximum of apph_max_delta*apph_delta_f_a) # so that f_surf best fits the shape of the AH. # Similarly with apph_find_z, though the maximum change applied is # dz*apph_find_z*apph_max_delta. # In addition, apph_find_z=-1 is a special parameter for mergers, # and will cause the code to extrapolate the motion of the holes if # ever we loose track of the AH. # # mxstep is the maximum number of integration steps that can # be applied during a single integration (other parameters # controlling the integration are still hard-coded in apph2.f) #======================================================================= apph_search_dt := 2.0 apph_find_a := 0 apph_find_z := 0 apph_max_delta := 3 apph_delta_f_a := 0.025 apph_delta_f_z := 0.5 mxstep := 6000 #======================================================================= # Solve by relaxation ? (if equal to one) # nsweeps -- max number of relaxation sweeps # nout -- sweep interval for output (according to trace_level) #======================================================================= solve_by_relax := 0 nsweeps := 5000 nout := 1000 #======================================================================= # Multigrid parameters #======================================================================= nvcycle := 1 preswp := 3 pstswp := 3 #======================================================================= # For extrapolation before the first iteration at each time # step, except the initial time. (1=second order, 2=third order) # We should probably update RNPL to do this automatically for all # the variables. #======================================================================= extrapolate := 0 #======================================================================= # Alternative boundary condition flags # ------------------------------------ # # dbc_... : =1 for Dirichlet, t>0, else 'standard' b.c. # dbc_psi only has effect with the constrained solution for psi #======================================================================= dbc_alpha := 1 dbc_betarho := 1 dbc_betaz := 1 dbc_psi := 0 #======================================================================= # Alternative differencing of the MG outer b.c.'s #======================================================================= sym_diff := 1