From 3d265ab3a727dc5b60eac81c17992cdaa0c6f221 Mon Sep 17 00:00:00 2001 From: rhaas Date: Fri, 11 Jan 2013 15:04:01 +0000 Subject: GRHydro_InitData: output bsonic and psonicmag as well. From: Bruno Coutinho Mundim git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@189 ac85fae7-cede-4708-beff-ae01c7fa1c26 --- src/GRHydro_BondiM_new.F90 | 12 +++++++++--- 1 file changed, 9 insertions(+), 3 deletions(-) diff --git a/src/GRHydro_BondiM_new.F90 b/src/GRHydro_BondiM_new.F90 index 9483385..6991bbc 100644 --- a/src/GRHydro_BondiM_new.F90 +++ b/src/GRHydro_BondiM_new.F90 @@ -66,7 +66,7 @@ subroutine GRHydro_BondiM_Iso(CCTK_ARGUMENTS) CCTK_REAL :: drhodr, det, rhocheck, rhocheck2, riso, rnew, rsch, ucheck CCTK_REAL :: uiso, uisocheck, vcheck, ucheck2, vcheck2, xhat,yhat, zhat, xp, yp, zp CCTK_REAL :: f,df,ddf,a,b,c,rsm,roverm,dudr,uisocheck2,auiso,buiso - CCTK_REAL :: bondi_bsmooth, bmag + CCTK_REAL :: bondi_bsmooth, bmag, bsonic, psonicmag character(400) :: debug_message @@ -227,11 +227,17 @@ subroutine GRHydro_BondiM_Iso(CCTK_ARGUMENTS) else bmag = bondi_bmag end if - write(debug_message,'(a,2f22.14)')'Bondi pars: rs, psonic',rs, psonic - call CCTK_INFO(debug_message) + bsonic = bmag*(M/rs)**2 * sqrt(1.0d0-2.0d0*M/rs) + psonicmag = 0.5d0*bsonic**2 + write(debug_message,'(a,2f22.14)')'Bondi pars: bondi_bmag,bondi_beta_sonicpt',& bmag,bondi_beta_sonicpt call CCTK_INFO(debug_message) + write(debug_message,'(a,2f22.14)')'Bondi pars: rs,bsonic',rs,bsonic + call CCTK_INFO(debug_message) + write(debug_message,'(a,2f22.14)')'Bondi pars: psonic,psonicmag',& + psonic,psonicmag + call CCTK_INFO(debug_message) do i=1,nx do j=1,ny -- cgit v1.2.3