User: rhaas Date: 2012/08/04 02:41 AM
Modified: /trunk/src/ GRHydro_Eigenproblem.F90, GRHydro_PPMMReconstruct_drv.F90, GRHydro_UpdateMask.F90
Log: GRHydro: add forgotten label in error message, whitespace change
File Changes:
Directory: /trunk/src/ ======================
File [modified]: GRHydro_Eigenproblem.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_Eigenproblem.F90 2012-08-04 07:18:51 UTC (rev 406) +++ trunk/src/GRHydro_Eigenproblem.F90 2012-08-04 07:41:17 UTC (rev 407) @@ -88,7 +88,7 @@ if(cs2.lt.0.0d0) then !$OMP CRITICAL if (abs(cs2) .gt. 1.0d-4) then - write(warnline,'(a50,6g16.7)') 'rho, dpdrho, press*dpdeps/rho**2, eps, press/rho: ', abs(cs2), rho, dpdrho, press * dpdeps / (rho**2), eps, press/rho + write(warnline,'(a50,6g16.7)') 'abs(cs2), rho, dpdrho, press*dpdeps/rho**2, eps, press/rho: ', abs(cs2), rho, dpdrho, press * dpdeps / (rho**2), eps, press/rho call CCTK_WARN(1,warnline) call CCTK_WARN(1,"cs2 < 0! Check speed of sound calculation!") cs2 = 0.0d0
File [modified]: GRHydro_PPMMReconstruct_drv.F90 Delta lines: +1 -0 =================================================================== --- trunk/src/GRHydro_PPMMReconstruct_drv.F90 2012-08-04 07:18:51 UTC (rev 406) +++ trunk/src/GRHydro_PPMMReconstruct_drv.F90 2012-08-04 07:41:17 UTC (rev 407) @@ -240,6 +240,7 @@ end do !$OMP END PARALLEL DO end if + if(evolve_Y_e.ne.0) then !$OMP PARALLEL DO PRIVATE(j, k) do k = GRHydro_stencil, nz - GRHydro_stencil + 1
File [modified]: GRHydro_UpdateMask.F90 Delta lines: +2 -0 =================================================================== --- trunk/src/GRHydro_UpdateMask.F90 2012-08-04 07:18:51 UTC (rev 406) +++ trunk/src/GRHydro_UpdateMask.F90 2012-08-04 07:41:17 UTC (rev 407) @@ -48,6 +48,7 @@ densrhs(i,j,k) = 0.0d0 srhs(i,j,k,:) = 0.0d0 taurhs(i,j,k) = 0.0d0 + ! Set real-valued mask! This will be sync'ed and right after syncing translated to ! our standard integer based mask (so that atmosphere_mask is still valid!). atmosphere_mask_real(i,j,k) = 1 @@ -69,6 +70,7 @@ densrhs(i,j,k) = 0.0d0 srhs(i,j,k,:) = 0.0d0 taurhs(i,j,k) = 0.0d0 + ! Set real-valued mask! This will be sync'ed and right after syncing translated to ! our standard integer based mask (so that atmosphere_mask is still valid!). atmosphere_mask_real(i,j,k) = 1
commits@lists.einsteintoolkit.org