[Commits] [svn:einsteintoolkit] GRHydro_InitData/trunk/src/ (Rev. 189)

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/01/11 09:04 AM

Modified:
 /trunk/src/
  GRHydro_BondiM_new.F90

Log:
 GRHydro_InitData: output bsonic and psonicmag as well.
 
 From: Bruno Coutinho Mundim <bcmsma at astro.rit.edu>

File Changes:

Directory: /trunk/src/
======================

File [modified]: GRHydro_BondiM_new.F90
Delta lines: +9 -3
===================================================================
--- trunk/src/GRHydro_BondiM_new.F90	2012-11-21 22:57:45 UTC (rev 188)
+++ trunk/src/GRHydro_BondiM_new.F90	2013-01-11 15:04:01 UTC (rev 189)
@@ -66,7 +66,7 @@
   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 @@
  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