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