diff options
Diffstat (limited to 'src/GRHydro_UpdateMask.F90')
-rw-r--r-- | src/GRHydro_UpdateMask.F90 | 7 |
1 files changed, 7 insertions, 0 deletions
diff --git a/src/GRHydro_UpdateMask.F90 b/src/GRHydro_UpdateMask.F90 index b8c3208..b5a84a5 100644 --- a/src/GRHydro_UpdateMask.F90 +++ b/src/GRHydro_UpdateMask.F90 @@ -201,6 +201,7 @@ subroutine GRHydro_InitAtmosMask(CCTK_ARGUMENTS) CCTK_INT :: i,j,k + !$OMP PARALLEL DO PRIVATE(i,j,k) do k = 1, cctk_lsh(3) do j = 1, cctk_lsh(2) do i = 1, cctk_lsh(1) @@ -211,6 +212,7 @@ subroutine GRHydro_InitAtmosMask(CCTK_ARGUMENTS) end do end do end do + !$OMP END PARALLEL DO call CCTK_INFO("Setting up the atmosphere mask: points with rho<rho_min are set to be atmosphere") @@ -445,6 +447,7 @@ subroutine GRHydro_InitialAtmosphereReset(CCTK_ARGUMENTS) rho_min = rho_min * initial_atmosphere_factor endif +!$OMP PARALLEL DO PRIVATE(det,keytemp,i,j,k,anyerr,keyerr) do k = 1, cctk_lsh(3) do j = 1, cctk_lsh(2) do i = 1, cctk_lsh(1) @@ -477,6 +480,7 @@ subroutine GRHydro_InitialAtmosphereReset(CCTK_ARGUMENTS) temperature(i,j,k),y_e(i,j,k)) y_e_con(i,j,k) = dens(i,j,k) * y_e(i,j,k) else + keytemp = 0 call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& rho(i,j,k),eps(i,j,k),xtemp,xye,press(i,j,k),keyerr,anyerr) call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& @@ -518,6 +522,7 @@ subroutine GRHydro_InitialAtmosphereReset(CCTK_ARGUMENTS) temperature_p(i,j,k),y_e_p(i,j,k)) y_e_con_p(i,j,k) = dens_p(i,j,k) * y_e_p(i,j,k) else + keytemp = 0 call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& rho_p(i,j,k),eps_p(i,j,k),xtemp,xye,press_p(i,j,k),keyerr,anyerr) call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& @@ -559,6 +564,7 @@ subroutine GRHydro_InitialAtmosphereReset(CCTK_ARGUMENTS) temperature_p_p(i,j,k),y_e_p_p(i,j,k)) y_e_con_p_p(i,j,k) = dens_p_p(i,j,k) * y_e_p_p(i,j,k) else + keytemp = 0 call EOS_Omni_press(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& rho_p_p(i,j,k),eps_p_p(i,j,k),xtemp,xye,press_p_p(i,j,k),keyerr,anyerr) call EOS_Omni_EpsFromPress(GRHydro_polytrope_handle,keytemp,GRHydro_eos_rf_prec,n,& @@ -576,6 +582,7 @@ subroutine GRHydro_InitialAtmosphereReset(CCTK_ARGUMENTS) end do end do end do +!$OMP END PARALLEL DO end subroutine GRHydro_InitialAtmosphereReset |