Commit a7851d6f authored by m.tavelli's avatar m.tavelli

tested in 3d

parent b6af85e9
Pipeline #168631 failed with stage
......@@ -45,20 +45,16 @@ RECURSIVE SUBROUTINE InitParameters(STRLEN,PARSETUP)
end select
END SUBROUTINE InitParameters
RECURSIVE SUBROUTINE InitialData(xGPin, tGp, Q)
RECURSIVE SUBROUTINE InitialData(xGP, tGp, Q)
USE, INTRINSIC :: ISO_C_BINDING
USE MainVariables, ONLY : nVar, nDim, EQN, ICType
IMPLICIT NONE
! Argument list
REAL, INTENT(IN) :: xGPin(3), tGp !
REAL, INTENT(IN) :: xGP(3), tGp !
REAL, INTENT(OUT) :: Q(nVar) !
REAL :: up(nVar),V0(nVar),r,xGP(3)
REAL :: up(nVar),V0(nVar),r
xGP=xGPin
IF(nDim==2) THEN
xGP(3) = 0.
ENDIF
select case(ICType)
case('CCZ4MinkowskiSrc')
!
......@@ -72,7 +68,7 @@ RECURSIVE SUBROUTINE InitialData(xGPin, tGp, Q)
!
!V0(54) = 0.2*EXP(-0.5*( xGP(1)**2+xGP(2)**2)/1.0**2 )
!
r = SQRT( xGP(1)**2 + xGP(2)**2 + xGP(3)**2 )
r = SQRT( xGP(1)**2 + xGP(2)**2 + xGP(3)**2)
!
V0(60) = 0.001*EXP(-0.5*( (xGP(1)-2.0)**2+xGP(2)**2+xGP(3)**2)/1.0**2 ) + 0.001*EXP(-0.5*( (xGP(1)+2.0)**2+xGP(2)**2+xGP(3)**2)/1.0**2 )
IF( r>5.0 .AND. r<10. ) THEN
......
......@@ -4030,11 +4030,17 @@ RECURSIVE SUBROUTINE pderefinecriteria(refine_flag, max_luh,min_luh,x)
! return
!end if
#ifdef CCZ4EINSTEIN
if(abs(max_luh(60)-min_luh(60))>1.e-4 .or. abs(max_luh(54)-min_luh(54))>1.e-4) then
if(abs(max_luh(60)-min_luh(60))>1.e-4 .or. abs(max_luh(54)-min_luh(54))>1.e-3) then
refine_flag=2
else
refine_flag=0
end if
if(nDim .eq. 3) then
if(sqrt(x(1)**2+x(2)**2+x(3)**2)<15.0) then
refine_flag=2
end if
end if
#endif
END SUBROUTINE pderefinecriteria
......
......@@ -55,7 +55,7 @@
"name": "FOCCZ4Solver",
"order": 4,
"maximum_mesh_size": 10.0,
"maximum_mesh_depth": 2,
"maximum_mesh_depth": 1,
"time_stepping": "global",
"aderdg_kernel": {
"language": "C",
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment