41 un_m1, un, pn_m1, pn, phin_m1, phin)
44 REAL(KIND=8),
INTENT(OUT):: time
45 REAL(KIND=8),
INTENT(IN) :: dt
46 INTEGER,
DIMENSION(:),
INTENT(IN) :: list_mode
47 REAL(KIND=8),
DIMENSION(:,:,:),
INTENT(OUT):: un_m1, un
48 REAL(KIND=8),
DIMENSION(:,:,:),
INTENT(OUT):: pn_m1, pn, phin_m1, phin
50 REAL(KIND=8),
DIMENSION(mesh_c%np) :: pn_m2
53 DO i= 1,
SIZE(list_mode)
57 un_m1(:,j,i) =
vv_exact(j,mesh_f%rr,mode,time-dt)
58 un(:,j,i) =
vv_exact(j,mesh_f%rr,mode,time)
62 pn_m2(:) =
pp_exact(j,mesh_c%rr,mode,time-2*dt)
63 pn_m1(:,j,i) =
pp_exact(j,mesh_c%rr,mode,time-dt)
64 pn(:,j,i) =
pp_exact(j,mesh_c%rr,mode,time)
65 phin_m1(:,j,i) = pn_m1(:,j,i) - pn_m2(:)
66 phin(:,j,i) = pn(:,j,i) - pn_m1(:,j,i)
119 INTEGER ,
INTENT(IN) :: TYPE
120 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr
121 INTEGER ,
INTENT(IN) :: mode, i
122 REAL(KIND=8),
INTENT(IN) :: time
123 REAL(KIND=8),
INTENT(IN) :: Re
124 CHARACTER(LEN=2),
INTENT(IN) :: ty
125 REAL(KIND=8),
DIMENSION(:,:,:),
OPTIONAL,
INTENT(IN) :: opt_density
126 REAL(KIND=8),
DIMENSION(:,:,:),
OPTIONAL,
INTENT(IN) :: opt_tempn
127 REAL(KIND=8),
DIMENSION(SIZE(rr,2)) :: vv
130 CHARACTER(LEN=2) :: np
132 IF (
PRESENT(opt_density))
CALL error_petsc(
'density should not be present for test 16')
133 IF (
PRESENT(opt_tempn))
CALL error_petsc(
'temperature should not be present for test 16')
139 n=type; n=
SIZE(rr,1); n=mode; n=i; r=time; r=re; np=ty
172 FUNCTION vv_exact(TYPE,rr,m,t) RESULT(vv)
174 INTEGER ,
INTENT(IN) :: TYPE
175 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr
176 INTEGER,
INTENT(IN) :: m
177 REAL(KIND=8),
INTENT(IN) :: t
178 REAL(KIND=8),
DIMENSION(SIZE(rr,2)) :: vv
179 REAL(KIND=8),
DIMENSION(SIZE(rr,2)) :: r, z
217 FUNCTION pp_exact(TYPE,rr,m,t) RESULT (vv)
219 INTEGER ,
INTENT(IN) :: TYPE
220 REAL(KIND=8),
DIMENSION(:,:),
INTENT(IN) :: rr
221 INTEGER ,
INTENT(IN) :: m
222 REAL(KIND=8),
INTENT(IN) :: t
223 REAL(KIND=8),
DIMENSION(SIZE(rr,2)) :: vv
224 REAL(KIND=8),
DIMENSION(SIZE(rr,2)) :: r, z
236 vv = (0.1d0*r)**2/2.d0
real(kind=8) function, dimension(size(rr, 2)), public vv_exact(TYPE, rr, m, t)
subroutine error_petsc(string)
subroutine, public init_velocity_pressure(mesh_f, mesh_c, time, dt, list_mode, un_m1, un, pn_m1, pn, phin_m1, phin)
real(kind=8) function, dimension(size(rr, 2)), public pp_exact(TYPE, rr, m, t)
real(kind=8) function, dimension(size(rr, 2)), public source_in_ns_momentum(TYPE, rr, mode, i, time, Re, ty, opt_density, opt_tempn)