added average terms
This commit is contained in:
ignis
parent
2cea7e01bb
commit
251dbfc010
2 changed files with 204 additions and 50 deletions
|
|
@ -16,15 +16,24 @@ divn = ddx(nx) + ddy(ny) + ddz(nz)
|
|||
|
||||
absk = abs(divn)
|
||||
|
||||
fu = (vn + sd) * nx
|
||||
|
||||
lapc = d2dx(c) + d2dy(c) + d2dz(c)
|
||||
|
||||
t1 = (vn + sd) * nx
|
||||
|
||||
t1 = sqr(fu')
|
||||
t2 = sqr(absk')
|
||||
t3 = sqr(divn')
|
||||
t2 = t1' * nx'
|
||||
|
||||
t3 = t1' * divn'
|
||||
|
||||
t4 = t1' * absk'
|
||||
|
||||
t5 = sqr(t1')
|
||||
|
||||
t6 = sqr(nx')
|
||||
|
||||
t7 = sqr(divn')
|
||||
|
||||
t8 = sqr(absk')
|
||||
|
||||
avg { c, u, v, w, fsd, lapc, ddx_c, ddy_c, ddz_c, d2dx_c, d2dy_c, d2dz_c }
|
||||
|
||||
avg fsd { u, v, w, nx, ny, nz, t1, t2, fu, absk, divn, sd }
|
||||
avg fsd { u, v, w, nx, ny, nz, t1, t2, t3, t4, t5, t6, t7, t8, absk, divn, sd }
|
||||
|
|
|
|||
233
code/m_terms.f90
233
code/m_terms.f90
|
|
@ -7,8 +7,9 @@ use m_calculate
|
|||
|
||||
implicit none
|
||||
|
||||
character (len = *), parameter :: output_header="x avg_c avg_u avg_v avg_w avg_fsd avg_lapc avg_ddx_c avg_ddy_c avg_ddz_c avg_d2dx_c avg_d2dy_c avg_d2dz_c fsd_avg_u fsd_avg_v fsd_avg_w fsd_avg_nx fsd_avg_ny fsd_avg_nz fsd_avg_t1 fsd_avg_t2 fsd_avg_fu fsd_avg_absk fsd_avg_divn fsd_avg_sd"
|
||||
character (len = *), parameter :: output_header="x avg_c avg_u avg_v avg_w avg_fsd avg_lapc avg_ddx_c avg_ddy_c avg_ddz_c avg_d2dx_c avg_d2dy_c avg_d2dz_c fsd_avg_u fsd_avg_v fsd_avg_w fsd_avg_nx fsd_avg_ny fsd_avg_nz fsd_avg_t1 fsd_avg_t2 fsd_avg_t3 fsd_avg_t4 fsd_avg_t5 fsd_avg_t6 fsd_avg_t7 fsd_avg_t8 fsd_avg_absk fsd_avg_divn fsd_avg_sd"
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_divn
|
||||
real*8, allocatable, dimension(:,:,:) :: t3_fsd
|
||||
real*8, allocatable, dimension(:,:,:) :: ddx_nx
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_ny
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_nx
|
||||
|
|
@ -20,12 +21,16 @@ real*8, allocatable, dimension(:) :: avg_d2dz_c
|
|||
real*8, allocatable, dimension(:) :: avg_d2dx_c
|
||||
real*8, allocatable, dimension(:,:,:) :: ddy_c
|
||||
real*8, allocatable, dimension(:,:,:) :: ddx_c
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_fu
|
||||
real*8, allocatable, dimension(:,:,:) :: t8_fsd
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_t3
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_t2
|
||||
real*8, allocatable, dimension(:) :: avg_ddx_c
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_t1
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_t7
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_v
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_t5
|
||||
real*8, allocatable, dimension(:,:,:) :: ddy_ny
|
||||
real*8, allocatable, dimension(:,:,:) :: wrate
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_t8
|
||||
real*8, allocatable, dimension(:) :: avg_fsd
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_absk
|
||||
real*8, allocatable, dimension(:,:,:) :: nx
|
||||
|
|
@ -36,24 +41,29 @@ real*8, allocatable, dimension(:) :: fsd_avg_sd
|
|||
real*8, allocatable, dimension(:,:,:) :: lapc
|
||||
real*8, allocatable, dimension(:) :: avg_ddy_c
|
||||
real*8, allocatable, dimension(:,:,:) :: divn
|
||||
real*8, allocatable, dimension(:,:,:) :: fu
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_t1
|
||||
real*8, allocatable, dimension(:,:,:) :: t7_fsd
|
||||
real*8, allocatable, dimension(:) :: avg_ddx_c
|
||||
real*8, allocatable, dimension(:) :: avg_lapc
|
||||
real*8, allocatable, dimension(:,:,:) :: t1_fsd
|
||||
real*8, allocatable, dimension(:,:,:) :: t6_fsd
|
||||
real*8, allocatable, dimension(:,:,:) :: absk
|
||||
real*8, allocatable, dimension(:) :: avg_ddz_c
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_t6
|
||||
real*8, allocatable, dimension(:,:,:) :: c
|
||||
real*8, allocatable, dimension(:) :: avg_c
|
||||
real*8, allocatable, dimension(:,:,:) :: fsd
|
||||
real*8, allocatable, dimension(:,:,:) :: t5_fsd
|
||||
real*8, allocatable, dimension(:,:,:) :: d2dy_c
|
||||
real*8, allocatable, dimension(:) :: avg_d2dy_c
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_t4
|
||||
real*8, allocatable, dimension(:,:,:) :: d2dx_c
|
||||
real*8, allocatable, dimension(:,:,:) :: t1
|
||||
real*8, allocatable, dimension(:) :: avg_u
|
||||
real*8, allocatable, dimension(:) :: avg_v
|
||||
real*8, allocatable, dimension(:) :: avg_w
|
||||
real*8, allocatable, dimension(:,:,:) :: d2dz_c
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_w
|
||||
real*8, allocatable, dimension(:) :: fsd_avg_u
|
||||
real*8, allocatable, dimension(:,:,:) :: t4_fsd
|
||||
real*8, allocatable, dimension(:,:,:) :: sd
|
||||
|
||||
|
||||
|
|
@ -64,6 +74,7 @@ subroutine m_terms_init
|
|||
integer :: ierr
|
||||
|
||||
allocate(fsd_avg_divn(nxp), stat=ierr) ; fsd_avg_divn = 0.
|
||||
allocate(t3_fsd(nxp,nyp,nzp), stat=ierr) ; t3_fsd = 0.
|
||||
allocate(ddx_nx(nxp,nyp,nzp), stat=ierr) ; ddx_nx = 0.
|
||||
allocate(fsd_avg_ny(nxp), stat=ierr) ; fsd_avg_ny = 0.
|
||||
allocate(fsd_avg_nx(nxp), stat=ierr) ; fsd_avg_nx = 0.
|
||||
|
|
@ -75,12 +86,16 @@ allocate(avg_d2dz_c(nxp), stat=ierr) ; avg_d2dz_c = 0.
|
|||
allocate(avg_d2dx_c(nxp), stat=ierr) ; avg_d2dx_c = 0.
|
||||
allocate(ddy_c(nxp,nyp,nzp), stat=ierr) ; ddy_c = 0.
|
||||
allocate(ddx_c(nxp,nyp,nzp), stat=ierr) ; ddx_c = 0.
|
||||
allocate(fsd_avg_fu(nxp), stat=ierr) ; fsd_avg_fu = 0.
|
||||
allocate(t8_fsd(nxp,nyp,nzp), stat=ierr) ; t8_fsd = 0.
|
||||
allocate(fsd_avg_t3(nxp), stat=ierr) ; fsd_avg_t3 = 0.
|
||||
allocate(fsd_avg_t2(nxp), stat=ierr) ; fsd_avg_t2 = 0.
|
||||
allocate(avg_ddx_c(nxp), stat=ierr) ; avg_ddx_c = 0.
|
||||
allocate(fsd_avg_t1(nxp), stat=ierr) ; fsd_avg_t1 = 0.
|
||||
allocate(fsd_avg_t7(nxp), stat=ierr) ; fsd_avg_t7 = 0.
|
||||
allocate(fsd_avg_v(nxp), stat=ierr) ; fsd_avg_v = 0.
|
||||
allocate(fsd_avg_t5(nxp), stat=ierr) ; fsd_avg_t5 = 0.
|
||||
allocate(ddy_ny(nxp,nyp,nzp), stat=ierr) ; ddy_ny = 0.
|
||||
allocate(wrate(nxp,nyp,nzp), stat=ierr) ; wrate = 0.
|
||||
allocate(fsd_avg_t8(nxp), stat=ierr) ; fsd_avg_t8 = 0.
|
||||
allocate(avg_fsd(nxp), stat=ierr) ; avg_fsd = 0.
|
||||
allocate(fsd_avg_absk(nxp), stat=ierr) ; fsd_avg_absk = 0.
|
||||
allocate(nx(nxp,nyp,nzp), stat=ierr) ; nx = 0.
|
||||
|
|
@ -91,24 +106,29 @@ allocate(fsd_avg_sd(nxp), stat=ierr) ; fsd_avg_sd = 0.
|
|||
allocate(lapc(nxp,nyp,nzp), stat=ierr) ; lapc = 0.
|
||||
allocate(avg_ddy_c(nxp), stat=ierr) ; avg_ddy_c = 0.
|
||||
allocate(divn(nxp,nyp,nzp), stat=ierr) ; divn = 0.
|
||||
allocate(fu(nxp,nyp,nzp), stat=ierr) ; fu = 0.
|
||||
allocate(fsd_avg_t1(nxp), stat=ierr) ; fsd_avg_t1 = 0.
|
||||
allocate(t7_fsd(nxp,nyp,nzp), stat=ierr) ; t7_fsd = 0.
|
||||
allocate(avg_ddx_c(nxp), stat=ierr) ; avg_ddx_c = 0.
|
||||
allocate(avg_lapc(nxp), stat=ierr) ; avg_lapc = 0.
|
||||
allocate(t1_fsd(nxp,nyp,nzp), stat=ierr) ; t1_fsd = 0.
|
||||
allocate(t6_fsd(nxp,nyp,nzp), stat=ierr) ; t6_fsd = 0.
|
||||
allocate(absk(nxp,nyp,nzp), stat=ierr) ; absk = 0.
|
||||
allocate(avg_ddz_c(nxp), stat=ierr) ; avg_ddz_c = 0.
|
||||
allocate(fsd_avg_t6(nxp), stat=ierr) ; fsd_avg_t6 = 0.
|
||||
allocate(c(nxp,nyp,nzp), stat=ierr) ; c = 0.
|
||||
allocate(avg_c(nxp), stat=ierr) ; avg_c = 0.
|
||||
allocate(fsd(nxp,nyp,nzp), stat=ierr) ; fsd = 0.
|
||||
allocate(t5_fsd(nxp,nyp,nzp), stat=ierr) ; t5_fsd = 0.
|
||||
allocate(d2dy_c(nxp,nyp,nzp), stat=ierr) ; d2dy_c = 0.
|
||||
allocate(avg_d2dy_c(nxp), stat=ierr) ; avg_d2dy_c = 0.
|
||||
allocate(fsd_avg_t4(nxp), stat=ierr) ; fsd_avg_t4 = 0.
|
||||
allocate(d2dx_c(nxp,nyp,nzp), stat=ierr) ; d2dx_c = 0.
|
||||
allocate(t1(nxp,nyp,nzp), stat=ierr) ; t1 = 0.
|
||||
allocate(avg_u(nxp), stat=ierr) ; avg_u = 0.
|
||||
allocate(avg_v(nxp), stat=ierr) ; avg_v = 0.
|
||||
allocate(avg_w(nxp), stat=ierr) ; avg_w = 0.
|
||||
allocate(d2dz_c(nxp,nyp,nzp), stat=ierr) ; d2dz_c = 0.
|
||||
allocate(fsd_avg_w(nxp), stat=ierr) ; fsd_avg_w = 0.
|
||||
allocate(fsd_avg_u(nxp), stat=ierr) ; fsd_avg_u = 0.
|
||||
allocate(t4_fsd(nxp,nyp,nzp), stat=ierr) ; t4_fsd = 0.
|
||||
allocate(sd(nxp,nyp,nzp), stat=ierr) ; sd = 0.
|
||||
|
||||
|
||||
|
|
@ -118,6 +138,7 @@ end subroutine m_terms_init
|
|||
subroutine m_terms_finalize
|
||||
|
||||
deallocate(fsd_avg_divn)
|
||||
deallocate(t3_fsd)
|
||||
deallocate(ddx_nx)
|
||||
deallocate(fsd_avg_ny)
|
||||
deallocate(fsd_avg_nx)
|
||||
|
|
@ -129,12 +150,16 @@ deallocate(avg_d2dz_c)
|
|||
deallocate(avg_d2dx_c)
|
||||
deallocate(ddy_c)
|
||||
deallocate(ddx_c)
|
||||
deallocate(fsd_avg_fu)
|
||||
deallocate(t8_fsd)
|
||||
deallocate(fsd_avg_t3)
|
||||
deallocate(fsd_avg_t2)
|
||||
deallocate(avg_ddx_c)
|
||||
deallocate(fsd_avg_t1)
|
||||
deallocate(fsd_avg_t7)
|
||||
deallocate(fsd_avg_v)
|
||||
deallocate(fsd_avg_t5)
|
||||
deallocate(ddy_ny)
|
||||
deallocate(wrate)
|
||||
deallocate(fsd_avg_t8)
|
||||
deallocate(avg_fsd)
|
||||
deallocate(fsd_avg_absk)
|
||||
deallocate(nx)
|
||||
|
|
@ -145,24 +170,29 @@ deallocate(fsd_avg_sd)
|
|||
deallocate(lapc)
|
||||
deallocate(avg_ddy_c)
|
||||
deallocate(divn)
|
||||
deallocate(fu)
|
||||
deallocate(fsd_avg_t1)
|
||||
deallocate(t7_fsd)
|
||||
deallocate(avg_ddx_c)
|
||||
deallocate(avg_lapc)
|
||||
deallocate(t1_fsd)
|
||||
deallocate(t6_fsd)
|
||||
deallocate(absk)
|
||||
deallocate(avg_ddz_c)
|
||||
deallocate(fsd_avg_t6)
|
||||
deallocate(c)
|
||||
deallocate(avg_c)
|
||||
deallocate(fsd)
|
||||
deallocate(t5_fsd)
|
||||
deallocate(d2dy_c)
|
||||
deallocate(avg_d2dy_c)
|
||||
deallocate(fsd_avg_t4)
|
||||
deallocate(d2dx_c)
|
||||
deallocate(t1)
|
||||
deallocate(avg_u)
|
||||
deallocate(avg_v)
|
||||
deallocate(avg_w)
|
||||
deallocate(d2dz_c)
|
||||
deallocate(fsd_avg_w)
|
||||
deallocate(fsd_avg_u)
|
||||
deallocate(t4_fsd)
|
||||
deallocate(sd)
|
||||
|
||||
|
||||
|
|
@ -186,15 +216,6 @@ call ddz ( ddz_c, c )
|
|||
call ddy ( ddy_c, c )
|
||||
call ddx ( ddx_c, c )
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
avg_ddx_c(i) = avg_ddx_c(i) + ddx_c(i,j,k)
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
|
|
@ -213,6 +234,15 @@ end do
|
|||
end do
|
||||
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
avg_ddx_c(i) = avg_ddx_c(i) + ddx_c(i,j,k)
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
|
|
@ -474,7 +504,7 @@ end do
|
|||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
fu(i,j,k) = ( ( vn(i,j,k) + sd(i,j,k) ) * nx(i,j,k) )
|
||||
t1(i,j,k) = ( ( vn(i,j,k) + sd(i,j,k) ) * nx(i,j,k) )
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
|
@ -483,7 +513,7 @@ end do
|
|||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
fsd_avg_fu(i) = fsd_avg_fu(i) + fu(i,j,k) * fsd(i,j,k)
|
||||
fsd_avg_t1(i) = fsd_avg_t1(i) + t1(i,j,k) * fsd(i,j,k)
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
|
@ -591,9 +621,9 @@ call MPI_ALLREDUCE(MPI_IN_PLACE, fsd_avg_nz, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_T
|
|||
fsd_avg_nz = fsd_avg_nz / denum / avg_fsd
|
||||
|
||||
|
||||
call MPI_ALLREDUCE(MPI_IN_PLACE, fsd_avg_fu, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err)
|
||||
call MPI_ALLREDUCE(MPI_IN_PLACE, fsd_avg_t1, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err)
|
||||
|
||||
fsd_avg_fu = fsd_avg_fu / denum / avg_fsd
|
||||
fsd_avg_t1 = fsd_avg_t1 / denum / avg_fsd
|
||||
|
||||
|
||||
call MPI_ALLREDUCE(MPI_IN_PLACE, fsd_avg_absk, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err)
|
||||
|
|
@ -635,7 +665,7 @@ call ddx ( ddx_c, c )
|
|||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
fsd(i,j,k) = ( sqrt ( ( ( ((ddx_c(i,j,k))*(ddx_c(i,j,k))) + ((ddy_c(i,j,k))*(ddy_c(i,j,k))) ) + ((ddz_c(i,j,k))*(ddz_c(i,j,k))) ) ) )
|
||||
wrate(i,j,k) = ( rxn_rate ( c(i,j,k) ) )
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
|
@ -644,13 +674,11 @@ end do
|
|||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
wrate(i,j,k) = ( rxn_rate ( c(i,j,k) ) )
|
||||
fsd(i,j,k) = ( sqrt ( ( ( ((ddx_c(i,j,k))*(ddx_c(i,j,k))) + ((ddy_c(i,j,k))*(ddy_c(i,j,k))) ) + ((ddz_c(i,j,k))*(ddz_c(i,j,k))) ) ) )
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
||||
call d2dy ( d2dy_c, c )
|
||||
call d2dx ( d2dx_c, c )
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
|
|
@ -699,6 +727,44 @@ end do
|
|||
end do
|
||||
end do
|
||||
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
t7_fsd(i,j,k) = (((divn(i,j,k) - fsd_avg_divn(i)))*((divn(i,j,k) - fsd_avg_divn(i))))
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
fsd_avg_t7(i) = fsd_avg_t7(i) + t7_fsd(i,j,k) * fsd(i,j,k)
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
t6_fsd(i,j,k) = (((nx(i,j,k) - fsd_avg_nx(i)))*((nx(i,j,k) - fsd_avg_nx(i))))
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
fsd_avg_t6(i) = fsd_avg_t6(i) + t6_fsd(i,j,k) * fsd(i,j,k)
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
||||
call d2dy ( d2dy_c, c )
|
||||
call d2dx ( d2dx_c, c )
|
||||
call d2dz ( d2dz_c, c )
|
||||
|
||||
do k = 1, nzp
|
||||
|
|
@ -713,7 +779,7 @@ end do
|
|||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
t2_fsd(i,j,k) = (((absk(i,j,k) - fsd_avg_absk(i)))*((absk(i,j,k) - fsd_avg_absk(i))))
|
||||
t8_fsd(i,j,k) = (((absk(i,j,k) - fsd_avg_absk(i)))*((absk(i,j,k) - fsd_avg_absk(i))))
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
|
@ -722,7 +788,7 @@ end do
|
|||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
fsd_avg_t2(i) = fsd_avg_t2(i) + t2_fsd(i,j,k) * fsd(i,j,k)
|
||||
fsd_avg_t8(i) = fsd_avg_t8(i) + t8_fsd(i,j,k) * fsd(i,j,k)
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
|
@ -740,7 +806,7 @@ end do
|
|||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
fu(i,j,k) = ( ( vn(i,j,k) + sd(i,j,k) ) * nx(i,j,k) )
|
||||
t1(i,j,k) = ( ( vn(i,j,k) + sd(i,j,k) ) * nx(i,j,k) )
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
|
@ -749,7 +815,7 @@ end do
|
|||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
t1_fsd(i,j,k) = (((fu(i,j,k) - fsd_avg_fu(i)))*((fu(i,j,k) - fsd_avg_fu(i))))
|
||||
t2_fsd(i,j,k) = ( (t1(i,j,k) - fsd_avg_t1(i)) * (nx(i,j,k) - fsd_avg_nx(i)) )
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
|
@ -758,7 +824,61 @@ end do
|
|||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
fsd_avg_t1(i) = fsd_avg_t1(i) + t1_fsd(i,j,k) * fsd(i,j,k)
|
||||
fsd_avg_t2(i) = fsd_avg_t2(i) + t2_fsd(i,j,k) * fsd(i,j,k)
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
t5_fsd(i,j,k) = (((t1(i,j,k) - fsd_avg_t1(i)))*((t1(i,j,k) - fsd_avg_t1(i))))
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
fsd_avg_t5(i) = fsd_avg_t5(i) + t5_fsd(i,j,k) * fsd(i,j,k)
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
t3_fsd(i,j,k) = ( (t1(i,j,k) - fsd_avg_t1(i)) * (divn(i,j,k) - fsd_avg_divn(i)) )
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
fsd_avg_t3(i) = fsd_avg_t3(i) + t3_fsd(i,j,k) * fsd(i,j,k)
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
t4_fsd(i,j,k) = ( (t1(i,j,k) - fsd_avg_t1(i)) * (absk(i,j,k) - fsd_avg_absk(i)) )
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
||||
|
||||
do k = 1, nzp
|
||||
do j = 1, nyp
|
||||
do i = 1, nxp
|
||||
fsd_avg_t4(i) = fsd_avg_t4(i) + t4_fsd(i,j,k) * fsd(i,j,k)
|
||||
end do
|
||||
end do
|
||||
end do
|
||||
|
|
@ -776,16 +896,41 @@ real*8 :: denum
|
|||
denum=real(nfiles*nyp*nzp)
|
||||
|
||||
|
||||
call MPI_ALLREDUCE(MPI_IN_PLACE, fsd_avg_t1, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err)
|
||||
|
||||
fsd_avg_t1 = fsd_avg_t1 / denum / avg_fsd
|
||||
|
||||
|
||||
call MPI_ALLREDUCE(MPI_IN_PLACE, fsd_avg_t2, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err)
|
||||
|
||||
fsd_avg_t2 = fsd_avg_t2 / denum / avg_fsd
|
||||
|
||||
|
||||
call MPI_ALLREDUCE(MPI_IN_PLACE, fsd_avg_t3, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err)
|
||||
|
||||
fsd_avg_t3 = fsd_avg_t3 / denum / avg_fsd
|
||||
|
||||
|
||||
call MPI_ALLREDUCE(MPI_IN_PLACE, fsd_avg_t4, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err)
|
||||
|
||||
fsd_avg_t4 = fsd_avg_t4 / denum / avg_fsd
|
||||
|
||||
|
||||
call MPI_ALLREDUCE(MPI_IN_PLACE, fsd_avg_t5, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err)
|
||||
|
||||
fsd_avg_t5 = fsd_avg_t5 / denum / avg_fsd
|
||||
|
||||
|
||||
call MPI_ALLREDUCE(MPI_IN_PLACE, fsd_avg_t6, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err)
|
||||
|
||||
fsd_avg_t6 = fsd_avg_t6 / denum / avg_fsd
|
||||
|
||||
|
||||
call MPI_ALLREDUCE(MPI_IN_PLACE, fsd_avg_t7, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err)
|
||||
|
||||
fsd_avg_t7 = fsd_avg_t7 / denum / avg_fsd
|
||||
|
||||
|
||||
call MPI_ALLREDUCE(MPI_IN_PLACE, fsd_avg_t8, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err)
|
||||
|
||||
fsd_avg_t8 = fsd_avg_t8 / denum / avg_fsd
|
||||
|
||||
|
||||
|
||||
end subroutine m_terms_average_pass2
|
||||
|
||||
|
|
@ -800,7 +945,7 @@ open (200, file="qEdge_X.dat")
|
|||
write (200,*) output_header
|
||||
|
||||
do i=1,nxp
|
||||
write (200,'(25e20.10)') real(i)*hxp, avg_c(i), avg_u(i), avg_v(i), avg_w(i), avg_fsd(i), avg_lapc(i), avg_ddx_c(i), avg_ddy_c(i), avg_ddz_c(i), avg_d2dx_c(i), avg_d2dy_c(i), avg_d2dz_c(i), fsd_avg_u(i), fsd_avg_v(i), fsd_avg_w(i), fsd_avg_nx(i), fsd_avg_ny(i), fsd_avg_nz(i), fsd_avg_t1(i), fsd_avg_t2(i), fsd_avg_fu(i), fsd_avg_absk(i), fsd_avg_divn(i), fsd_avg_sd(i)
|
||||
write (200,'(30e20.10)') real(i)*hxp, avg_c(i), avg_u(i), avg_v(i), avg_w(i), avg_fsd(i), avg_lapc(i), avg_ddx_c(i), avg_ddy_c(i), avg_ddz_c(i), avg_d2dx_c(i), avg_d2dy_c(i), avg_d2dz_c(i), fsd_avg_u(i), fsd_avg_v(i), fsd_avg_w(i), fsd_avg_nx(i), fsd_avg_ny(i), fsd_avg_nz(i), fsd_avg_t1(i), fsd_avg_t2(i), fsd_avg_t3(i), fsd_avg_t4(i), fsd_avg_t5(i), fsd_avg_t6(i), fsd_avg_t7(i), fsd_avg_t8(i), fsd_avg_absk(i), fsd_avg_divn(i), fsd_avg_sd(i)
|
||||
end do
|
||||
|
||||
close (200)
|
||||
|
|
|
|||
Loading…
Add table
Reference in a new issue