gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

PROBLEMS: Simplify the Rayleigh-Taylor problem setup.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2017-03-06 12:23:17 -03:00
parent 251ac66fde
commit dbdb5b3b4a
1 changed files with 9 additions and 20 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

29
src/problems.F90
View File

@ -1259,7 +1259,6 @@ module problems
! local variables ! local variables
! !
integer :: i, j, k integer :: i, j, k
real(kind=8) :: yl, yu, dy, dyh
real(kind=8) :: sn, cs real(kind=8) :: sn, cs
! local arrays ! local arrays
@ -1303,11 +1302,6 @@ module problems
! !
y(1:jm) = pdata%meta%ymin + ay(pdata%meta%level,1:jm) y(1:jm) = pdata%meta%ymin + ay(pdata%meta%level,1:jm)
! calculate mesh intervals and areas
!
dy = ady(pdata%meta%level)
dyh = 0.5d+00 * dy
! set the ambient density and pressure ! set the ambient density and pressure
! !
q(idn,:) = dens q(idn,:) = dens
@ -1337,33 +1331,28 @@ module problems
do k = 1, km do k = 1, km
do j = 1, jm do j = 1, jm
! calculate the corner Y coordinates
!
yl = abs(y(j)) - dyh
yu = abs(y(j)) + dyh
! set the primitive variables for two regions ! set the primitive variables for two regions
! !
if (yu <= ycut) then if (y(j) <= ycut) then
q(idn,1:im) = dens q(idn,1:im) = dens
else if (yl >= ycut) then
q(idn,1:im) = dens * drat
else else
q(idn,1:im) = dens * ((yu - ycut) + (ycut - yl) * drat) / dy q(idn,1:im) = dens * drat
end if end if
! set the pressure ! set the pressure
! !
if (ipr > 0) q(ipr,1:im) = pres + q(idn,1:im) * gacc * y(j) if (ipr > 0) q(ipr,1:im) = pres + q(idn,1:im) * gacc * y(j)
! reset the velocity components
!
q(ivx,1:im) = 0.0d+00
q(ivy,1:im) = 0.0d+00
q(ivz,1:im) = 0.0d+00
! add a random seed velocity component ! add a random seed velocity component
! !
do i = 1, im do i = 1, im
q(ivx,i) = q(ivx,i) + vper * randomn()
q(ivy,i) = q(ivy,i) + vper * randomn() q(ivy,i) = q(ivy,i) + vper * randomn()
#if NDIMS == 3
q(ivz,i) = q(ivz,i) + vper * randomn()
#endif /* NDIMS == 3 */
end do end do
! convert the primitive variables to conservative ones ! convert the primitive variables to conservative ones