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@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"