diff options
| author | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2013-01-11 15:04:01 +0000 |
|---|---|---|
| committer | rhaas <rhaas@ac85fae7-cede-4708-beff-ae01c7fa1c26> | 2013-01-11 15:04:01 +0000 |
| commit | 3d265ab3a727dc5b60eac81c17992cdaa0c6f221 (patch) | |
| tree | 7c52444663e1b3669ea9f161c71bb9843dd56a6d | |
| parent | fba24a1753d874047237908fa4b71ab70293a970 (diff) | |
GRHydro_InitData: output bsonic and psonicmag as well.
From: Bruno Coutinho Mundim <bcmsma@astro.rit.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/GRHydro_InitData/trunk@189 ac85fae7-cede-4708-beff-ae01c7fa1c26
| -rw-r--r-- | src/GRHydro_BondiM_new.F90 | 12 |
1 files 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 |