From 65be0d52bf42c74e7f0f58b852c14d8334a2891c Mon Sep 17 00:00:00 2001 From: ignis Date: Sun, 7 Jul 2019 00:11:10 +0900 Subject: [PATCH] pass2 switch --- code/code_gen/post.py | 5 +- code/m_terms.f90 | 3723 +---------------------------------------- code/post.f90 | 7 + 3 files changed, 50 insertions(+), 3685 deletions(-) diff --git a/code/code_gen/post.py b/code/code_gen/post.py index 6c217f3..549e8f3 100644 --- a/code/code_gen/post.py +++ b/code/code_gen/post.py @@ -699,6 +699,9 @@ class Stage4(): avg1 = "\n".join(allvar[v].code_avg() for v in set1) avg2 = "\n".join(allvar[v].code_avg() for v in set2) + ffmt = 'logical, parameter :: pass2_required={}' + declf = ffmt.format('.true.' if len(set2) > 0 else '.false.') + hfmt = 'character (len = *), parameter :: output_header="{}"' declh = hfmt.format(" ".join(["x"] + set1 + set2)) @@ -724,7 +727,7 @@ close (200) md = {} md["module_name"] = "terms" - md["module_data"] = "\n".join((declh, declavg, declarr)) + md["module_data"] = "\n".join((declf, declh, declavg, declarr)) md["module_init"] = "\n".join((allocavg, allocarr)) md["module_finalize"] = "\n".join((freeavg, freearr)) md["module_pass1"] = calc1 diff --git a/code/m_terms.f90 b/code/m_terms.f90 index e31949b..77dfc5e 100644 --- a/code/m_terms.f90 +++ b/code/m_terms.f90 @@ -7,173 +7,16 @@ use m_calculate implicit none -character (len = *), parameter :: output_header="x avg_c avg_d2dx_c avg_d2dy_c avg_d2dz_c avg_ddx_c avg_ddy_c avg_ddz_c avg_fsd1 avg_fsd2 avg_fsd3 avg_fsd4 avg_fsd5 avg_fsd6 avg_fsd7 avg_fsd8 avg_fsd9 avg_lapc avg_u avg_v avg_w fsd1_avg_absk fsd1_avg_divn fsd1_avg_nx fsd1_avg_ny fsd1_avg_nz fsd1_avg_sd fsd1_avg_t1 fsd2_avg_absk fsd2_avg_divn fsd2_avg_nx fsd2_avg_ny fsd2_avg_nz fsd2_avg_sd fsd2_avg_t1 fsd3_avg_absk fsd3_avg_divn fsd3_avg_nx fsd3_avg_ny fsd3_avg_nz fsd3_avg_sd fsd3_avg_t1 fsd4_avg_absk fsd4_avg_divn fsd4_avg_nx fsd4_avg_ny fsd4_avg_nz fsd4_avg_sd fsd4_avg_t1 fsd5_avg_absk fsd5_avg_divn fsd5_avg_nx fsd5_avg_ny fsd5_avg_nz fsd5_avg_sd fsd5_avg_t1 fsd6_avg_absk fsd6_avg_divn fsd6_avg_nx fsd6_avg_ny fsd6_avg_nz fsd6_avg_sd fsd6_avg_t1 fsd7_avg_absk fsd7_avg_divn fsd7_avg_nx fsd7_avg_ny fsd7_avg_nz fsd7_avg_sd fsd7_avg_t1 fsd8_avg_absk fsd8_avg_divn fsd8_avg_nx fsd8_avg_ny fsd8_avg_nz fsd8_avg_sd fsd8_avg_t1 fsd9_avg_absk fsd9_avg_divn fsd9_avg_nx fsd9_avg_ny fsd9_avg_nz fsd9_avg_sd fsd9_avg_t1 fsd1_avg_t2 fsd1_avg_t3 fsd1_avg_t4 fsd1_avg_t5 fsd1_avg_t6 fsd1_avg_t7 fsd1_avg_t8 fsd2_avg_t2 fsd2_avg_t3 fsd2_avg_t4 fsd2_avg_t5 fsd2_avg_t6 fsd2_avg_t7 fsd2_avg_t8 fsd3_avg_t2 fsd3_avg_t3 fsd3_avg_t4 fsd3_avg_t5 fsd3_avg_t6 fsd3_avg_t7 fsd3_avg_t8 fsd4_avg_t2 fsd4_avg_t3 fsd4_avg_t4 fsd4_avg_t5 fsd4_avg_t6 fsd4_avg_t7 fsd4_avg_t8 fsd5_avg_t2 fsd5_avg_t3 fsd5_avg_t4 fsd5_avg_t5 fsd5_avg_t6 fsd5_avg_t7 fsd5_avg_t8 fsd6_avg_t2 fsd6_avg_t3 fsd6_avg_t4 fsd6_avg_t5 fsd6_avg_t6 fsd6_avg_t7 fsd6_avg_t8 fsd7_avg_t2 fsd7_avg_t3 fsd7_avg_t4 fsd7_avg_t5 fsd7_avg_t6 fsd7_avg_t7 fsd7_avg_t8 fsd8_avg_t2 fsd8_avg_t3 fsd8_avg_t4 fsd8_avg_t5 fsd8_avg_t6 fsd8_avg_t7 fsd8_avg_t8 fsd9_avg_t2 fsd9_avg_t3 fsd9_avg_t4 fsd9_avg_t5 fsd9_avg_t6 fsd9_avg_t7 fsd9_avg_t8" -real*8, allocatable, dimension(:) :: fsd6_avg_nx -real*8, allocatable, dimension(:) :: fsd6_avg_ny -real*8, allocatable, dimension(:) :: fsd6_avg_nz -real*8, allocatable, dimension(:) :: fsd3_avg_t1 -real*8, allocatable, dimension(:) :: fsd3_avg_t3 -real*8, allocatable, dimension(:) :: fsd3_avg_t2 -real*8, allocatable, dimension(:) :: fsd3_avg_t5 -real*8, allocatable, dimension(:) :: fsd3_avg_t4 -real*8, allocatable, dimension(:) :: fsd3_avg_t7 -real*8, allocatable, dimension(:) :: fsd3_avg_t6 -real*8, allocatable, dimension(:) :: fsd3_avg_t8 -real*8, allocatable, dimension(:) :: avg_ddy_c -real*8, allocatable, dimension(:) :: fsd8_avg_nx -real*8, allocatable, dimension(:) :: avg_fsd6 -real*8, allocatable, dimension(:) :: avg_fsd7 -real*8, allocatable, dimension(:) :: avg_fsd4 -real*8, allocatable, dimension(:) :: avg_fsd5 -real*8, allocatable, dimension(:) :: avg_fsd2 -real*8, allocatable, dimension(:) :: fsd8_avg_sd -real*8, allocatable, dimension(:) :: avg_fsd1 -real*8, allocatable, dimension(:) :: avg_fsd8 -real*8, allocatable, dimension(:) :: avg_fsd9 -real*8, allocatable, dimension(:) :: fsd6_avg_t5 -real*8, allocatable, dimension(:) :: fsd8_avg_t1 -real*8, allocatable, dimension(:) :: fsd7_avg_nz -real*8, allocatable, dimension(:) :: fsd7_avg_ny -real*8, allocatable, dimension(:) :: fsd7_avg_nx -real*8, allocatable, dimension(:) :: fsd1_avg_sd -real*8, allocatable, dimension(:) :: fsd8_avg_divn -real*8, allocatable, dimension(:) :: fsd1_avg_t8 -real*8, allocatable, dimension(:) :: fsd1_avg_t7 -real*8, allocatable, dimension(:) :: fsd1_avg_t6 -real*8, allocatable, dimension(:) :: fsd1_avg_t5 -real*8, allocatable, dimension(:) :: fsd1_avg_t4 -real*8, allocatable, dimension(:) :: fsd5_avg_ny -real*8, allocatable, dimension(:) :: fsd5_avg_nx -real*8, allocatable, dimension(:) :: fsd1_avg_t1 -real*8, allocatable, dimension(:) :: fsd5_avg_nz -real*8, allocatable, dimension(:) :: fsd7_avg_t8 -real*8, allocatable, dimension(:) :: fsd6_avg_t8 -real*8, allocatable, dimension(:) :: fsd7_avg_t5 -real*8, allocatable, dimension(:) :: fsd7_avg_t4 -real*8, allocatable, dimension(:) :: fsd7_avg_t7 -real*8, allocatable, dimension(:) :: fsd7_avg_t6 -real*8, allocatable, dimension(:) :: fsd7_avg_t1 -real*8, allocatable, dimension(:) :: fsd6_avg_t3 -real*8, allocatable, dimension(:) :: fsd7_avg_t3 -real*8, allocatable, dimension(:) :: fsd2_avg_absk -real*8, allocatable, dimension(:) :: fsd1_avg_ny -real*8, allocatable, dimension(:) :: fsd1_avg_nx -real*8, allocatable, dimension(:) :: fsd1_avg_nz -real*8, allocatable, dimension(:) :: fsd6_avg_divn -real*8, allocatable, dimension(:) :: avg_ddx_c -real*8, allocatable, dimension(:) :: fsd6_avg_sd -real*8, allocatable, dimension(:) :: fsd3_avg_sd -real*8, allocatable, dimension(:) :: fsd9_avg_divn -real*8, allocatable, dimension(:) :: fsd8_avg_absk -real*8, allocatable, dimension(:) :: fsd4_avg_absk -real*8, allocatable, dimension(:) :: avg_ddz_c -real*8, allocatable, dimension(:) :: fsd5_avg_t1 -real*8, allocatable, dimension(:) :: fsd2_avg_divn -real*8, allocatable, dimension(:) :: fsd8_avg_t2 -real*8, allocatable, dimension(:) :: fsd8_avg_t3 -real*8, allocatable, dimension(:) :: fsd8_avg_t4 -real*8, allocatable, dimension(:) :: fsd8_avg_t5 -real*8, allocatable, dimension(:) :: fsd8_avg_t6 -real*8, allocatable, dimension(:) :: fsd8_avg_t7 -real*8, allocatable, dimension(:) :: fsd2_avg_nx -real*8, allocatable, dimension(:) :: fsd2_avg_ny -real*8, allocatable, dimension(:) :: fsd2_avg_nz -real*8, allocatable, dimension(:) :: fsd4_avg_divn -real*8, allocatable, dimension(:) :: fsd9_avg_absk -real*8, allocatable, dimension(:) :: avg_fsd3 -real*8, allocatable, dimension(:) :: fsd2_avg_t2 -real*8, allocatable, dimension(:) :: fsd2_avg_t3 -real*8, allocatable, dimension(:) :: fsd2_avg_t1 -real*8, allocatable, dimension(:) :: fsd2_avg_t6 -real*8, allocatable, dimension(:) :: fsd2_avg_t7 -real*8, allocatable, dimension(:) :: fsd2_avg_t4 -real*8, allocatable, dimension(:) :: fsd2_avg_t5 -real*8, allocatable, dimension(:) :: fsd8_avg_nz -real*8, allocatable, dimension(:) :: fsd2_avg_t8 -real*8, allocatable, dimension(:) :: fsd8_avg_ny -real*8, allocatable, dimension(:) :: fsd7_avg_sd -real*8, allocatable, dimension(:) :: fsd1_avg_t3 -real*8, allocatable, dimension(:) :: fsd2_avg_sd -real*8, allocatable, dimension(:) :: fsd4_avg_sd -real*8, allocatable, dimension(:) :: fsd7_avg_absk -real*8, allocatable, dimension(:) :: fsd6_avg_t1 -real*8, allocatable, dimension(:) :: fsd3_avg_divn -real*8, allocatable, dimension(:) :: fsd1_avg_t2 -real*8, allocatable, dimension(:) :: fsd6_avg_t6 -real*8, allocatable, dimension(:) :: fsd6_avg_t7 -real*8, allocatable, dimension(:) :: fsd6_avg_t4 -real*8, allocatable, dimension(:) :: fsd4_avg_t4 -real*8, allocatable, dimension(:) :: fsd4_avg_t5 -real*8, allocatable, dimension(:) :: fsd4_avg_t6 -real*8, allocatable, dimension(:) :: fsd4_avg_t7 -real*8, allocatable, dimension(:) :: fsd4_avg_t1 -real*8, allocatable, dimension(:) :: fsd4_avg_t2 -real*8, allocatable, dimension(:) :: fsd4_avg_t3 -real*8, allocatable, dimension(:) :: fsd6_avg_t2 -real*8, allocatable, dimension(:) :: fsd4_avg_t8 -real*8, allocatable, dimension(:) :: avg_lapc -real*8, allocatable, dimension(:) :: fsd5_avg_divn -real*8, allocatable, dimension(:) :: fsd3_avg_absk -real*8, allocatable, dimension(:) :: fsd7_avg_t2 -real*8, allocatable, dimension(:) :: fsd9_avg_t7 -real*8, allocatable, dimension(:) :: fsd9_avg_t6 -real*8, allocatable, dimension(:) :: fsd9_avg_t5 -real*8, allocatable, dimension(:) :: fsd9_avg_t4 -real*8, allocatable, dimension(:) :: fsd9_avg_t3 -real*8, allocatable, dimension(:) :: fsd9_avg_t2 -real*8, allocatable, dimension(:) :: fsd9_avg_t1 -real*8, allocatable, dimension(:) :: fsd9_avg_t8 -real*8, allocatable, dimension(:) :: fsd8_avg_t8 -real*8, allocatable, dimension(:) :: fsd5_avg_t8 -real*8, allocatable, dimension(:) :: fsd9_avg_sd -real*8, allocatable, dimension(:) :: fsd5_avg_t3 -real*8, allocatable, dimension(:) :: fsd5_avg_t2 -real*8, allocatable, dimension(:) :: fsd5_avg_sd -real*8, allocatable, dimension(:) :: fsd5_avg_t7 -real*8, allocatable, dimension(:) :: fsd5_avg_t6 -real*8, allocatable, dimension(:) :: fsd5_avg_t5 -real*8, allocatable, dimension(:) :: fsd5_avg_t4 -real*8, allocatable, dimension(:) :: fsd1_avg_divn -real*8, allocatable, dimension(:) :: fsd5_avg_absk -real*8, allocatable, dimension(:) :: avg_d2dz_c -real*8, allocatable, dimension(:) :: avg_d2dx_c -real*8, allocatable, dimension(:) :: fsd4_avg_nz -real*8, allocatable, dimension(:) :: fsd4_avg_nx -real*8, allocatable, dimension(:) :: fsd4_avg_ny -real*8, allocatable, dimension(:) :: fsd6_avg_absk -real*8, allocatable, dimension(:) :: fsd9_avg_ny -real*8, allocatable, dimension(:) :: fsd9_avg_nx -real*8, allocatable, dimension(:) :: fsd9_avg_nz -real*8, allocatable, dimension(:) :: fsd3_avg_nz -real*8, allocatable, dimension(:) :: fsd3_avg_ny -real*8, allocatable, dimension(:) :: fsd3_avg_nx -real*8, allocatable, dimension(:) :: fsd7_avg_divn -real*8, allocatable, dimension(:) :: fsd1_avg_absk +logical, parameter :: pass2_required=.false. +character (len = *), parameter :: output_header="x avg_c avg_divn avg_fsd" real*8, allocatable, dimension(:) :: avg_c -real*8, allocatable, dimension(:) :: avg_d2dy_c -real*8, allocatable, dimension(:) :: avg_u -real*8, allocatable, dimension(:) :: avg_v -real*8, allocatable, dimension(:) :: avg_w +real*8, allocatable, dimension(:) :: avg_fsd +real*8, allocatable, dimension(:) :: avg_divn real*8, allocatable, dimension(:,:,:) :: xyzbuffer0 real*8, allocatable, dimension(:,:,:) :: xyzbuffer1 real*8, allocatable, dimension(:,:,:) :: xyzbuffer2 real*8, allocatable, dimension(:,:,:) :: xyzbuffer3 real*8, allocatable, dimension(:,:,:) :: xyzbuffer4 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer5 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer6 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer7 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer8 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer9 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer10 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer11 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer12 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer13 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer14 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer15 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer16 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer17 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer18 -real*8, allocatable, dimension(:,:,:) :: xyzbuffer19 contains @@ -181,344 +24,28 @@ subroutine m_terms_init integer :: ierr -allocate(fsd6_avg_nx(nxp), stat=ierr) ; fsd6_avg_nx = 0. -allocate(fsd6_avg_ny(nxp), stat=ierr) ; fsd6_avg_ny = 0. -allocate(fsd6_avg_nz(nxp), stat=ierr) ; fsd6_avg_nz = 0. -allocate(fsd3_avg_t1(nxp), stat=ierr) ; fsd3_avg_t1 = 0. -allocate(fsd3_avg_t3(nxp), stat=ierr) ; fsd3_avg_t3 = 0. -allocate(fsd3_avg_t2(nxp), stat=ierr) ; fsd3_avg_t2 = 0. -allocate(fsd3_avg_t5(nxp), stat=ierr) ; fsd3_avg_t5 = 0. -allocate(fsd3_avg_t4(nxp), stat=ierr) ; fsd3_avg_t4 = 0. -allocate(fsd3_avg_t7(nxp), stat=ierr) ; fsd3_avg_t7 = 0. -allocate(fsd3_avg_t6(nxp), stat=ierr) ; fsd3_avg_t6 = 0. -allocate(fsd3_avg_t8(nxp), stat=ierr) ; fsd3_avg_t8 = 0. -allocate(avg_ddy_c(nxp), stat=ierr) ; avg_ddy_c = 0. -allocate(fsd8_avg_nx(nxp), stat=ierr) ; fsd8_avg_nx = 0. -allocate(avg_fsd6(nxp), stat=ierr) ; avg_fsd6 = 0. -allocate(avg_fsd7(nxp), stat=ierr) ; avg_fsd7 = 0. -allocate(avg_fsd4(nxp), stat=ierr) ; avg_fsd4 = 0. -allocate(avg_fsd5(nxp), stat=ierr) ; avg_fsd5 = 0. -allocate(avg_fsd2(nxp), stat=ierr) ; avg_fsd2 = 0. -allocate(fsd8_avg_sd(nxp), stat=ierr) ; fsd8_avg_sd = 0. -allocate(avg_fsd1(nxp), stat=ierr) ; avg_fsd1 = 0. -allocate(avg_fsd8(nxp), stat=ierr) ; avg_fsd8 = 0. -allocate(avg_fsd9(nxp), stat=ierr) ; avg_fsd9 = 0. -allocate(fsd6_avg_t5(nxp), stat=ierr) ; fsd6_avg_t5 = 0. -allocate(fsd8_avg_t1(nxp), stat=ierr) ; fsd8_avg_t1 = 0. -allocate(fsd7_avg_nz(nxp), stat=ierr) ; fsd7_avg_nz = 0. -allocate(fsd7_avg_ny(nxp), stat=ierr) ; fsd7_avg_ny = 0. -allocate(fsd7_avg_nx(nxp), stat=ierr) ; fsd7_avg_nx = 0. -allocate(fsd1_avg_sd(nxp), stat=ierr) ; fsd1_avg_sd = 0. -allocate(fsd8_avg_divn(nxp), stat=ierr) ; fsd8_avg_divn = 0. -allocate(fsd1_avg_t8(nxp), stat=ierr) ; fsd1_avg_t8 = 0. -allocate(fsd1_avg_t7(nxp), stat=ierr) ; fsd1_avg_t7 = 0. -allocate(fsd1_avg_t6(nxp), stat=ierr) ; fsd1_avg_t6 = 0. -allocate(fsd1_avg_t5(nxp), stat=ierr) ; fsd1_avg_t5 = 0. -allocate(fsd1_avg_t4(nxp), stat=ierr) ; fsd1_avg_t4 = 0. -allocate(fsd5_avg_ny(nxp), stat=ierr) ; fsd5_avg_ny = 0. -allocate(fsd5_avg_nx(nxp), stat=ierr) ; fsd5_avg_nx = 0. -allocate(fsd1_avg_t1(nxp), stat=ierr) ; fsd1_avg_t1 = 0. -allocate(fsd5_avg_nz(nxp), stat=ierr) ; fsd5_avg_nz = 0. -allocate(fsd7_avg_t8(nxp), stat=ierr) ; fsd7_avg_t8 = 0. -allocate(fsd6_avg_t8(nxp), stat=ierr) ; fsd6_avg_t8 = 0. -allocate(fsd7_avg_t5(nxp), stat=ierr) ; fsd7_avg_t5 = 0. -allocate(fsd7_avg_t4(nxp), stat=ierr) ; fsd7_avg_t4 = 0. -allocate(fsd7_avg_t7(nxp), stat=ierr) ; fsd7_avg_t7 = 0. -allocate(fsd7_avg_t6(nxp), stat=ierr) ; fsd7_avg_t6 = 0. -allocate(fsd7_avg_t1(nxp), stat=ierr) ; fsd7_avg_t1 = 0. -allocate(fsd6_avg_t3(nxp), stat=ierr) ; fsd6_avg_t3 = 0. -allocate(fsd7_avg_t3(nxp), stat=ierr) ; fsd7_avg_t3 = 0. -allocate(fsd2_avg_absk(nxp), stat=ierr) ; fsd2_avg_absk = 0. -allocate(fsd1_avg_ny(nxp), stat=ierr) ; fsd1_avg_ny = 0. -allocate(fsd1_avg_nx(nxp), stat=ierr) ; fsd1_avg_nx = 0. -allocate(fsd1_avg_nz(nxp), stat=ierr) ; fsd1_avg_nz = 0. -allocate(fsd6_avg_divn(nxp), stat=ierr) ; fsd6_avg_divn = 0. -allocate(avg_ddx_c(nxp), stat=ierr) ; avg_ddx_c = 0. -allocate(fsd6_avg_sd(nxp), stat=ierr) ; fsd6_avg_sd = 0. -allocate(fsd3_avg_sd(nxp), stat=ierr) ; fsd3_avg_sd = 0. -allocate(fsd9_avg_divn(nxp), stat=ierr) ; fsd9_avg_divn = 0. -allocate(fsd8_avg_absk(nxp), stat=ierr) ; fsd8_avg_absk = 0. -allocate(fsd4_avg_absk(nxp), stat=ierr) ; fsd4_avg_absk = 0. -allocate(avg_ddz_c(nxp), stat=ierr) ; avg_ddz_c = 0. -allocate(fsd5_avg_t1(nxp), stat=ierr) ; fsd5_avg_t1 = 0. -allocate(fsd2_avg_divn(nxp), stat=ierr) ; fsd2_avg_divn = 0. -allocate(fsd8_avg_t2(nxp), stat=ierr) ; fsd8_avg_t2 = 0. -allocate(fsd8_avg_t3(nxp), stat=ierr) ; fsd8_avg_t3 = 0. -allocate(fsd8_avg_t4(nxp), stat=ierr) ; fsd8_avg_t4 = 0. -allocate(fsd8_avg_t5(nxp), stat=ierr) ; fsd8_avg_t5 = 0. -allocate(fsd8_avg_t6(nxp), stat=ierr) ; fsd8_avg_t6 = 0. -allocate(fsd8_avg_t7(nxp), stat=ierr) ; fsd8_avg_t7 = 0. -allocate(fsd2_avg_nx(nxp), stat=ierr) ; fsd2_avg_nx = 0. -allocate(fsd2_avg_ny(nxp), stat=ierr) ; fsd2_avg_ny = 0. -allocate(fsd2_avg_nz(nxp), stat=ierr) ; fsd2_avg_nz = 0. -allocate(fsd4_avg_divn(nxp), stat=ierr) ; fsd4_avg_divn = 0. -allocate(fsd9_avg_absk(nxp), stat=ierr) ; fsd9_avg_absk = 0. -allocate(avg_fsd3(nxp), stat=ierr) ; avg_fsd3 = 0. -allocate(fsd2_avg_t2(nxp), stat=ierr) ; fsd2_avg_t2 = 0. -allocate(fsd2_avg_t3(nxp), stat=ierr) ; fsd2_avg_t3 = 0. -allocate(fsd2_avg_t1(nxp), stat=ierr) ; fsd2_avg_t1 = 0. -allocate(fsd2_avg_t6(nxp), stat=ierr) ; fsd2_avg_t6 = 0. -allocate(fsd2_avg_t7(nxp), stat=ierr) ; fsd2_avg_t7 = 0. -allocate(fsd2_avg_t4(nxp), stat=ierr) ; fsd2_avg_t4 = 0. -allocate(fsd2_avg_t5(nxp), stat=ierr) ; fsd2_avg_t5 = 0. -allocate(fsd8_avg_nz(nxp), stat=ierr) ; fsd8_avg_nz = 0. -allocate(fsd2_avg_t8(nxp), stat=ierr) ; fsd2_avg_t8 = 0. -allocate(fsd8_avg_ny(nxp), stat=ierr) ; fsd8_avg_ny = 0. -allocate(fsd7_avg_sd(nxp), stat=ierr) ; fsd7_avg_sd = 0. -allocate(fsd1_avg_t3(nxp), stat=ierr) ; fsd1_avg_t3 = 0. -allocate(fsd2_avg_sd(nxp), stat=ierr) ; fsd2_avg_sd = 0. -allocate(fsd4_avg_sd(nxp), stat=ierr) ; fsd4_avg_sd = 0. -allocate(fsd7_avg_absk(nxp), stat=ierr) ; fsd7_avg_absk = 0. -allocate(fsd6_avg_t1(nxp), stat=ierr) ; fsd6_avg_t1 = 0. -allocate(fsd3_avg_divn(nxp), stat=ierr) ; fsd3_avg_divn = 0. -allocate(fsd1_avg_t2(nxp), stat=ierr) ; fsd1_avg_t2 = 0. -allocate(fsd6_avg_t6(nxp), stat=ierr) ; fsd6_avg_t6 = 0. -allocate(fsd6_avg_t7(nxp), stat=ierr) ; fsd6_avg_t7 = 0. -allocate(fsd6_avg_t4(nxp), stat=ierr) ; fsd6_avg_t4 = 0. -allocate(fsd4_avg_t4(nxp), stat=ierr) ; fsd4_avg_t4 = 0. -allocate(fsd4_avg_t5(nxp), stat=ierr) ; fsd4_avg_t5 = 0. -allocate(fsd4_avg_t6(nxp), stat=ierr) ; fsd4_avg_t6 = 0. -allocate(fsd4_avg_t7(nxp), stat=ierr) ; fsd4_avg_t7 = 0. -allocate(fsd4_avg_t1(nxp), stat=ierr) ; fsd4_avg_t1 = 0. -allocate(fsd4_avg_t2(nxp), stat=ierr) ; fsd4_avg_t2 = 0. -allocate(fsd4_avg_t3(nxp), stat=ierr) ; fsd4_avg_t3 = 0. -allocate(fsd6_avg_t2(nxp), stat=ierr) ; fsd6_avg_t2 = 0. -allocate(fsd4_avg_t8(nxp), stat=ierr) ; fsd4_avg_t8 = 0. -allocate(avg_lapc(nxp), stat=ierr) ; avg_lapc = 0. -allocate(fsd5_avg_divn(nxp), stat=ierr) ; fsd5_avg_divn = 0. -allocate(fsd3_avg_absk(nxp), stat=ierr) ; fsd3_avg_absk = 0. -allocate(fsd7_avg_t2(nxp), stat=ierr) ; fsd7_avg_t2 = 0. -allocate(fsd9_avg_t7(nxp), stat=ierr) ; fsd9_avg_t7 = 0. -allocate(fsd9_avg_t6(nxp), stat=ierr) ; fsd9_avg_t6 = 0. -allocate(fsd9_avg_t5(nxp), stat=ierr) ; fsd9_avg_t5 = 0. -allocate(fsd9_avg_t4(nxp), stat=ierr) ; fsd9_avg_t4 = 0. -allocate(fsd9_avg_t3(nxp), stat=ierr) ; fsd9_avg_t3 = 0. -allocate(fsd9_avg_t2(nxp), stat=ierr) ; fsd9_avg_t2 = 0. -allocate(fsd9_avg_t1(nxp), stat=ierr) ; fsd9_avg_t1 = 0. -allocate(fsd9_avg_t8(nxp), stat=ierr) ; fsd9_avg_t8 = 0. -allocate(fsd8_avg_t8(nxp), stat=ierr) ; fsd8_avg_t8 = 0. -allocate(fsd5_avg_t8(nxp), stat=ierr) ; fsd5_avg_t8 = 0. -allocate(fsd9_avg_sd(nxp), stat=ierr) ; fsd9_avg_sd = 0. -allocate(fsd5_avg_t3(nxp), stat=ierr) ; fsd5_avg_t3 = 0. -allocate(fsd5_avg_t2(nxp), stat=ierr) ; fsd5_avg_t2 = 0. -allocate(fsd5_avg_sd(nxp), stat=ierr) ; fsd5_avg_sd = 0. -allocate(fsd5_avg_t7(nxp), stat=ierr) ; fsd5_avg_t7 = 0. -allocate(fsd5_avg_t6(nxp), stat=ierr) ; fsd5_avg_t6 = 0. -allocate(fsd5_avg_t5(nxp), stat=ierr) ; fsd5_avg_t5 = 0. -allocate(fsd5_avg_t4(nxp), stat=ierr) ; fsd5_avg_t4 = 0. -allocate(fsd1_avg_divn(nxp), stat=ierr) ; fsd1_avg_divn = 0. -allocate(fsd5_avg_absk(nxp), stat=ierr) ; fsd5_avg_absk = 0. -allocate(avg_d2dz_c(nxp), stat=ierr) ; avg_d2dz_c = 0. -allocate(avg_d2dx_c(nxp), stat=ierr) ; avg_d2dx_c = 0. -allocate(fsd4_avg_nz(nxp), stat=ierr) ; fsd4_avg_nz = 0. -allocate(fsd4_avg_nx(nxp), stat=ierr) ; fsd4_avg_nx = 0. -allocate(fsd4_avg_ny(nxp), stat=ierr) ; fsd4_avg_ny = 0. -allocate(fsd6_avg_absk(nxp), stat=ierr) ; fsd6_avg_absk = 0. -allocate(fsd9_avg_ny(nxp), stat=ierr) ; fsd9_avg_ny = 0. -allocate(fsd9_avg_nx(nxp), stat=ierr) ; fsd9_avg_nx = 0. -allocate(fsd9_avg_nz(nxp), stat=ierr) ; fsd9_avg_nz = 0. -allocate(fsd3_avg_nz(nxp), stat=ierr) ; fsd3_avg_nz = 0. -allocate(fsd3_avg_ny(nxp), stat=ierr) ; fsd3_avg_ny = 0. -allocate(fsd3_avg_nx(nxp), stat=ierr) ; fsd3_avg_nx = 0. -allocate(fsd7_avg_divn(nxp), stat=ierr) ; fsd7_avg_divn = 0. -allocate(fsd1_avg_absk(nxp), stat=ierr) ; fsd1_avg_absk = 0. allocate(avg_c(nxp), stat=ierr) ; avg_c = 0. -allocate(avg_d2dy_c(nxp), stat=ierr) ; avg_d2dy_c = 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(avg_fsd(nxp), stat=ierr) ; avg_fsd = 0. +allocate(avg_divn(nxp), stat=ierr) ; avg_divn = 0. allocate(xyzbuffer0(nxp,nyp,nzp), stat=ierr) ; xyzbuffer0 = 0. allocate(xyzbuffer1(nxp,nyp,nzp), stat=ierr) ; xyzbuffer1 = 0. allocate(xyzbuffer2(nxp,nyp,nzp), stat=ierr) ; xyzbuffer2 = 0. allocate(xyzbuffer3(nxp,nyp,nzp), stat=ierr) ; xyzbuffer3 = 0. allocate(xyzbuffer4(nxp,nyp,nzp), stat=ierr) ; xyzbuffer4 = 0. -allocate(xyzbuffer5(nxp,nyp,nzp), stat=ierr) ; xyzbuffer5 = 0. -allocate(xyzbuffer6(nxp,nyp,nzp), stat=ierr) ; xyzbuffer6 = 0. -allocate(xyzbuffer7(nxp,nyp,nzp), stat=ierr) ; xyzbuffer7 = 0. -allocate(xyzbuffer8(nxp,nyp,nzp), stat=ierr) ; xyzbuffer8 = 0. -allocate(xyzbuffer9(nxp,nyp,nzp), stat=ierr) ; xyzbuffer9 = 0. -allocate(xyzbuffer10(nxp,nyp,nzp), stat=ierr) ; xyzbuffer10 = 0. -allocate(xyzbuffer11(nxp,nyp,nzp), stat=ierr) ; xyzbuffer11 = 0. -allocate(xyzbuffer12(nxp,nyp,nzp), stat=ierr) ; xyzbuffer12 = 0. -allocate(xyzbuffer13(nxp,nyp,nzp), stat=ierr) ; xyzbuffer13 = 0. -allocate(xyzbuffer14(nxp,nyp,nzp), stat=ierr) ; xyzbuffer14 = 0. -allocate(xyzbuffer15(nxp,nyp,nzp), stat=ierr) ; xyzbuffer15 = 0. -allocate(xyzbuffer16(nxp,nyp,nzp), stat=ierr) ; xyzbuffer16 = 0. -allocate(xyzbuffer17(nxp,nyp,nzp), stat=ierr) ; xyzbuffer17 = 0. -allocate(xyzbuffer18(nxp,nyp,nzp), stat=ierr) ; xyzbuffer18 = 0. -allocate(xyzbuffer19(nxp,nyp,nzp), stat=ierr) ; xyzbuffer19 = 0. end subroutine m_terms_init subroutine m_terms_finalize -deallocate(fsd6_avg_nx) -deallocate(fsd6_avg_ny) -deallocate(fsd6_avg_nz) -deallocate(fsd3_avg_t1) -deallocate(fsd3_avg_t3) -deallocate(fsd3_avg_t2) -deallocate(fsd3_avg_t5) -deallocate(fsd3_avg_t4) -deallocate(fsd3_avg_t7) -deallocate(fsd3_avg_t6) -deallocate(fsd3_avg_t8) -deallocate(avg_ddy_c) -deallocate(fsd8_avg_nx) -deallocate(avg_fsd6) -deallocate(avg_fsd7) -deallocate(avg_fsd4) -deallocate(avg_fsd5) -deallocate(avg_fsd2) -deallocate(fsd8_avg_sd) -deallocate(avg_fsd1) -deallocate(avg_fsd8) -deallocate(avg_fsd9) -deallocate(fsd6_avg_t5) -deallocate(fsd8_avg_t1) -deallocate(fsd7_avg_nz) -deallocate(fsd7_avg_ny) -deallocate(fsd7_avg_nx) -deallocate(fsd1_avg_sd) -deallocate(fsd8_avg_divn) -deallocate(fsd1_avg_t8) -deallocate(fsd1_avg_t7) -deallocate(fsd1_avg_t6) -deallocate(fsd1_avg_t5) -deallocate(fsd1_avg_t4) -deallocate(fsd5_avg_ny) -deallocate(fsd5_avg_nx) -deallocate(fsd1_avg_t1) -deallocate(fsd5_avg_nz) -deallocate(fsd7_avg_t8) -deallocate(fsd6_avg_t8) -deallocate(fsd7_avg_t5) -deallocate(fsd7_avg_t4) -deallocate(fsd7_avg_t7) -deallocate(fsd7_avg_t6) -deallocate(fsd7_avg_t1) -deallocate(fsd6_avg_t3) -deallocate(fsd7_avg_t3) -deallocate(fsd2_avg_absk) -deallocate(fsd1_avg_ny) -deallocate(fsd1_avg_nx) -deallocate(fsd1_avg_nz) -deallocate(fsd6_avg_divn) -deallocate(avg_ddx_c) -deallocate(fsd6_avg_sd) -deallocate(fsd3_avg_sd) -deallocate(fsd9_avg_divn) -deallocate(fsd8_avg_absk) -deallocate(fsd4_avg_absk) -deallocate(avg_ddz_c) -deallocate(fsd5_avg_t1) -deallocate(fsd2_avg_divn) -deallocate(fsd8_avg_t2) -deallocate(fsd8_avg_t3) -deallocate(fsd8_avg_t4) -deallocate(fsd8_avg_t5) -deallocate(fsd8_avg_t6) -deallocate(fsd8_avg_t7) -deallocate(fsd2_avg_nx) -deallocate(fsd2_avg_ny) -deallocate(fsd2_avg_nz) -deallocate(fsd4_avg_divn) -deallocate(fsd9_avg_absk) -deallocate(avg_fsd3) -deallocate(fsd2_avg_t2) -deallocate(fsd2_avg_t3) -deallocate(fsd2_avg_t1) -deallocate(fsd2_avg_t6) -deallocate(fsd2_avg_t7) -deallocate(fsd2_avg_t4) -deallocate(fsd2_avg_t5) -deallocate(fsd8_avg_nz) -deallocate(fsd2_avg_t8) -deallocate(fsd8_avg_ny) -deallocate(fsd7_avg_sd) -deallocate(fsd1_avg_t3) -deallocate(fsd2_avg_sd) -deallocate(fsd4_avg_sd) -deallocate(fsd7_avg_absk) -deallocate(fsd6_avg_t1) -deallocate(fsd3_avg_divn) -deallocate(fsd1_avg_t2) -deallocate(fsd6_avg_t6) -deallocate(fsd6_avg_t7) -deallocate(fsd6_avg_t4) -deallocate(fsd4_avg_t4) -deallocate(fsd4_avg_t5) -deallocate(fsd4_avg_t6) -deallocate(fsd4_avg_t7) -deallocate(fsd4_avg_t1) -deallocate(fsd4_avg_t2) -deallocate(fsd4_avg_t3) -deallocate(fsd6_avg_t2) -deallocate(fsd4_avg_t8) -deallocate(avg_lapc) -deallocate(fsd5_avg_divn) -deallocate(fsd3_avg_absk) -deallocate(fsd7_avg_t2) -deallocate(fsd9_avg_t7) -deallocate(fsd9_avg_t6) -deallocate(fsd9_avg_t5) -deallocate(fsd9_avg_t4) -deallocate(fsd9_avg_t3) -deallocate(fsd9_avg_t2) -deallocate(fsd9_avg_t1) -deallocate(fsd9_avg_t8) -deallocate(fsd8_avg_t8) -deallocate(fsd5_avg_t8) -deallocate(fsd9_avg_sd) -deallocate(fsd5_avg_t3) -deallocate(fsd5_avg_t2) -deallocate(fsd5_avg_sd) -deallocate(fsd5_avg_t7) -deallocate(fsd5_avg_t6) -deallocate(fsd5_avg_t5) -deallocate(fsd5_avg_t4) -deallocate(fsd1_avg_divn) -deallocate(fsd5_avg_absk) -deallocate(avg_d2dz_c) -deallocate(avg_d2dx_c) -deallocate(fsd4_avg_nz) -deallocate(fsd4_avg_nx) -deallocate(fsd4_avg_ny) -deallocate(fsd6_avg_absk) -deallocate(fsd9_avg_ny) -deallocate(fsd9_avg_nx) -deallocate(fsd9_avg_nz) -deallocate(fsd3_avg_nz) -deallocate(fsd3_avg_ny) -deallocate(fsd3_avg_nx) -deallocate(fsd7_avg_divn) -deallocate(fsd1_avg_absk) deallocate(avg_c) -deallocate(avg_d2dy_c) -deallocate(avg_u) -deallocate(avg_v) -deallocate(avg_w) +deallocate(avg_fsd) +deallocate(avg_divn) deallocate(xyzbuffer0) deallocate(xyzbuffer1) deallocate(xyzbuffer2) deallocate(xyzbuffer3) deallocate(xyzbuffer4) -deallocate(xyzbuffer5) -deallocate(xyzbuffer6) -deallocate(xyzbuffer7) -deallocate(xyzbuffer8) -deallocate(xyzbuffer9) -deallocate(xyzbuffer10) -deallocate(xyzbuffer11) -deallocate(xyzbuffer12) -deallocate(xyzbuffer13) -deallocate(xyzbuffer14) -deallocate(xyzbuffer15) -deallocate(xyzbuffer16) -deallocate(xyzbuffer17) -deallocate(xyzbuffer18) -deallocate(xyzbuffer19) end subroutine m_terms_finalize @@ -528,38 +55,11 @@ subroutine m_terms_calculate_pass1 integer :: i, j, k -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_u(i) = avg_u(i) + u(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_v(i) = avg_v(i) + v(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_w(i) = avg_w(i) + w(i,j,k) -end do -end do -end do - - ! ( 1.0 - y(i,j,k) ) do k = 1, nzp do j = 1, nyp do i = 1, nxp -xyzbuffer19(i,j,k) = ( 1.0 - y(i,j,k) ) +xyzbuffer3(i,j,k) = ( 1.0 - y(i,j,k) ) end do end do end do @@ -568,1017 +68,81 @@ end do do k = 1, nzp do j = 1, nyp do i = 1, nxp -avg_c(i) = avg_c(i) + xyzbuffer19(i,j,k) -end do -end do -end do - -call d2dx ( xyzbuffer18, xyzbuffer19 ) - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_d2dx_c(i) = avg_d2dx_c(i) + xyzbuffer18(i,j,k) -end do -end do -end do - -call d2dy ( xyzbuffer13, xyzbuffer19 ) - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_d2dy_c(i) = avg_d2dy_c(i) + xyzbuffer13(i,j,k) -end do -end do -end do - -call d2dz ( xyzbuffer12, xyzbuffer19 ) - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_d2dz_c(i) = avg_d2dz_c(i) + xyzbuffer12(i,j,k) -end do -end do -end do - -call ddx ( xyzbuffer11, xyzbuffer19 ) - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_ddx_c(i) = avg_ddx_c(i) + xyzbuffer11(i,j,k) -end do -end do -end do - -call ddy ( xyzbuffer10, xyzbuffer19 ) - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_ddy_c(i) = avg_ddy_c(i) + xyzbuffer10(i,j,k) -end do -end do -end do - -call ddz ( xyzbuffer17, xyzbuffer19 ) - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_ddz_c(i) = avg_ddz_c(i) + xyzbuffer17(i,j,k) +avg_c(i) = avg_c(i) + xyzbuffer3(i,j,k) end do end do end do +call ddx ( xyzbuffer2, xyzbuffer3 ) +call ddy ( xyzbuffer1, xyzbuffer3 ) +call ddz ( xyzbuffer0, xyzbuffer3 ) ! ( 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))) ) ) ) do k = 1, nzp do j = 1, nyp do i = 1, nxp -xyzbuffer16(i,j,k) = ( sqrt ( ( ( ((xyzbuffer11(i,j,k))*(xyzbuffer11(i,j,k))) + ((xyzbuffer10(i,j,k))*(xyzbuffer10(i,j,k))) ) + ((xyzbuffer17(i,j,k))*(xyzbuffer17(i,j,k))) ) ) ) +xyzbuffer3(i,j,k) = ( sqrt ( ( ( ((xyzbuffer2(i,j,k))*(xyzbuffer2(i,j,k))) + ((xyzbuffer1(i,j,k))*(xyzbuffer1(i,j,k))) ) + ((xyzbuffer0(i,j,k))*(xyzbuffer0(i,j,k))) ) ) ) end do end do end do -! ( ( - ddx_c(i,j,k) ) / fsd_orig(i,j,k) ) do k = 1, nzp do j = 1, nyp do i = 1, nxp -xyzbuffer15(i,j,k) = ( ( - xyzbuffer11(i,j,k) ) / xyzbuffer16(i,j,k) ) +avg_fsd(i) = avg_fsd(i) + xyzbuffer3(i,j,k) end do end do end do -! ( ( - ddy_c(i,j,k) ) / fsd_orig(i,j,k) ) +! ( ( - ddx_c(i,j,k) ) / fsd(i,j,k) ) do k = 1, nzp do j = 1, nyp do i = 1, nxp -xyzbuffer14(i,j,k) = ( ( - xyzbuffer10(i,j,k) ) / xyzbuffer16(i,j,k) ) +xyzbuffer4(i,j,k) = ( ( - xyzbuffer2(i,j,k) ) / xyzbuffer3(i,j,k) ) end do end do end do +call ddx ( xyzbuffer2, xyzbuffer4 ) -! ( ( - ddz_c(i,j,k) ) / fsd_orig(i,j,k) ) +! ( ( - ddy_c(i,j,k) ) / fsd(i,j,k) ) do k = 1, nzp do j = 1, nyp do i = 1, nxp -xyzbuffer9(i,j,k) = ( ( - xyzbuffer17(i,j,k) ) / xyzbuffer16(i,j,k) ) +xyzbuffer4(i,j,k) = ( ( - xyzbuffer1(i,j,k) ) / xyzbuffer3(i,j,k) ) end do end do end do -call ddx ( xyzbuffer8, xyzbuffer15 ) -call ddy ( xyzbuffer3, xyzbuffer14 ) -call ddz ( xyzbuffer2, xyzbuffer9 ) + +! ( ( - ddz_c(i,j,k) ) / fsd(i,j,k) ) +do k = 1, nzp +do j = 1, nyp +do i = 1, nxp +xyzbuffer1(i,j,k) = ( ( - xyzbuffer0(i,j,k) ) / xyzbuffer3(i,j,k) ) +end do +end do +end do + +call ddy ( xyzbuffer3, xyzbuffer4 ) +call ddz ( xyzbuffer0, xyzbuffer1 ) ! ( ( ddx_nx(i,j,k) + ddy_ny(i,j,k) ) + ddz_nz(i,j,k) ) do k = 1, nzp do j = 1, nyp do i = 1, nxp -xyzbuffer1(i,j,k) = ( ( xyzbuffer8(i,j,k) + xyzbuffer3(i,j,k) ) + xyzbuffer2(i,j,k) ) +xyzbuffer4(i,j,k) = ( ( xyzbuffer2(i,j,k) + xyzbuffer3(i,j,k) ) + xyzbuffer0(i,j,k) ) end do end do end do -! ( dabs ( divn(i,j,k) ) ) do k = 1, nzp do j = 1, nyp do i = 1, nxp -xyzbuffer0(i,j,k) = ( dabs ( xyzbuffer1(i,j,k) ) ) -end do -end do -end do - - -! ( ( d2dx_c(i,j,k) + d2dy_c(i,j,k) ) + d2dz_c(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer7(i,j,k) = ( ( xyzbuffer18(i,j,k) + xyzbuffer13(i,j,k) ) + xyzbuffer12(i,j,k) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_lapc(i) = avg_lapc(i) + xyzbuffer7(i,j,k) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.05, 0.15 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer6(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.05, 0.15 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask1(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer5(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer6(i,j,k) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_fsd1(i) = avg_fsd1(i) + xyzbuffer5(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd1_avg_absk(i) = fsd1_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer5(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd1_avg_divn(i) = fsd1_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer5(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd1_avg_nx(i) = fsd1_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer5(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd1_avg_ny(i) = fsd1_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer5(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd1_avg_nz(i) = fsd1_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer5(i,j,k) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.15, 0.25 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer4(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.15, 0.25 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask2(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer11(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer4(i,j,k) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_fsd2(i) = avg_fsd2(i) + xyzbuffer11(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd2_avg_absk(i) = fsd2_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer11(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd2_avg_divn(i) = fsd2_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer11(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd2_avg_nx(i) = fsd2_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer11(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd2_avg_ny(i) = fsd2_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer11(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd2_avg_nz(i) = fsd2_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer11(i,j,k) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.25, 0.35 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer10(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.25, 0.35 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask3(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer17(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer10(i,j,k) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_fsd3(i) = avg_fsd3(i) + xyzbuffer17(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd3_avg_absk(i) = fsd3_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer17(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd3_avg_divn(i) = fsd3_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer17(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd3_avg_nx(i) = fsd3_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer17(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd3_avg_ny(i) = fsd3_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer17(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd3_avg_nz(i) = fsd3_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer17(i,j,k) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.35, 0.45 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer8(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.35, 0.45 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask4(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer3(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer8(i,j,k) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_fsd4(i) = avg_fsd4(i) + xyzbuffer3(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd4_avg_absk(i) = fsd4_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer3(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd4_avg_divn(i) = fsd4_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer3(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd4_avg_nx(i) = fsd4_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer3(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd4_avg_ny(i) = fsd4_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer3(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd4_avg_nz(i) = fsd4_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer3(i,j,k) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.45, 0.55 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer2(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.45, 0.55 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask5(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer7(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer2(i,j,k) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_fsd5(i) = avg_fsd5(i) + xyzbuffer7(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd5_avg_absk(i) = fsd5_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer7(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd5_avg_divn(i) = fsd5_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer7(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd5_avg_nx(i) = fsd5_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer7(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd5_avg_ny(i) = fsd5_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer7(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd5_avg_nz(i) = fsd5_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer7(i,j,k) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.55, 0.65 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer6(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.55, 0.65 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask6(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer4(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer6(i,j,k) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_fsd6(i) = avg_fsd6(i) + xyzbuffer4(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd6_avg_absk(i) = fsd6_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer4(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd6_avg_divn(i) = fsd6_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer4(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd6_avg_nx(i) = fsd6_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer4(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd6_avg_ny(i) = fsd6_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer4(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd6_avg_nz(i) = fsd6_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer4(i,j,k) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.65, 0.75 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer10(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.65, 0.75 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask7(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer8(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer10(i,j,k) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_fsd7(i) = avg_fsd7(i) + xyzbuffer8(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd7_avg_absk(i) = fsd7_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer8(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd7_avg_divn(i) = fsd7_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer8(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd7_avg_nx(i) = fsd7_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer8(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd7_avg_ny(i) = fsd7_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer8(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd7_avg_nz(i) = fsd7_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer8(i,j,k) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.75, 0.85 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer2(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.75, 0.85 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask8(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer6(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer2(i,j,k) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_fsd8(i) = avg_fsd8(i) + xyzbuffer6(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd8_avg_absk(i) = fsd8_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer6(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd8_avg_divn(i) = fsd8_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer6(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd8_avg_nx(i) = fsd8_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer6(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd8_avg_ny(i) = fsd8_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer6(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd8_avg_nz(i) = fsd8_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer6(i,j,k) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.85, 0.95 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer10(i,j,k) = ( threshold_min_max ( xyzbuffer19(i,j,k), 0.85, 0.95 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask9(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer2(i,j,k) = ( xyzbuffer16(i,j,k) * xyzbuffer10(i,j,k) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd9_avg_absk(i) = fsd9_avg_absk(i) + xyzbuffer0(i,j,k) * xyzbuffer2(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd9_avg_divn(i) = fsd9_avg_divn(i) + xyzbuffer1(i,j,k) * xyzbuffer2(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -avg_fsd9(i) = avg_fsd9(i) + xyzbuffer2(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd9_avg_nx(i) = fsd9_avg_nx(i) + xyzbuffer15(i,j,k) * xyzbuffer2(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd9_avg_ny(i) = fsd9_avg_ny(i) + xyzbuffer14(i,j,k) * xyzbuffer2(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd9_avg_nz(i) = fsd9_avg_nz(i) + xyzbuffer9(i,j,k) * xyzbuffer2(i,j,k) -end do -end do -end do - - -! ( ( ( u(i,j,k) * nx(i,j,k) ) + ( v(i,j,k) * ny(i,j,k) ) ) + ( w(i,j,k) * nz(i,j,k) ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer1(i,j,k) = ( ( ( u(i,j,k) * xyzbuffer15(i,j,k) ) + ( v(i,j,k) * xyzbuffer14(i,j,k) ) ) + ( w(i,j,k) * xyzbuffer9(i,j,k) ) ) -end do -end do -end do - - -! ( rxn_rate ( c(i,j,k) ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer0(i,j,k) = ( rxn_rate ( xyzbuffer19(i,j,k) ) ) -end do -end do -end do - - -! ( ( ( ( ( d2dx_c(i,j,k) + d2dy_c(i,j,k) ) + d2dz_c(i,j,k) ) * rod ) + wrate(i,j,k) ) / fsd_orig(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer19(i,j,k) = ( ( ( ( ( xyzbuffer18(i,j,k) + xyzbuffer13(i,j,k) ) + xyzbuffer12(i,j,k) ) * rod ) + xyzbuffer0(i,j,k) ) / xyzbuffer16(i,j,k) ) -end do -end do -end do - - -! ( ( vn(i,j,k) + sd(i,j,k) ) * nx(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer18(i,j,k) = ( ( xyzbuffer1(i,j,k) + xyzbuffer19(i,j,k) ) * xyzbuffer15(i,j,k) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd1_avg_sd(i) = fsd1_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer5(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd1_avg_t1(i) = fsd1_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer5(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd2_avg_sd(i) = fsd2_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer11(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd2_avg_t1(i) = fsd2_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer11(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd3_avg_sd(i) = fsd3_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer17(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd3_avg_t1(i) = fsd3_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer17(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd4_avg_sd(i) = fsd4_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer3(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd4_avg_t1(i) = fsd4_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer3(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd5_avg_sd(i) = fsd5_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer7(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd5_avg_t1(i) = fsd5_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer7(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd6_avg_sd(i) = fsd6_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer4(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd6_avg_t1(i) = fsd6_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer4(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd7_avg_sd(i) = fsd7_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer8(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd7_avg_t1(i) = fsd7_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer8(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd8_avg_sd(i) = fsd8_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer6(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd8_avg_t1(i) = fsd8_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer6(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd9_avg_sd(i) = fsd9_avg_sd(i) + xyzbuffer19(i,j,k) * xyzbuffer2(i,j,k) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd9_avg_t1(i) = fsd9_avg_t1(i) + xyzbuffer18(i,j,k) * xyzbuffer2(i,j,k) +avg_divn(i) = avg_divn(i) + xyzbuffer4(i,j,k) end do end do end do @@ -1600,414 +164,14 @@ call MPI_ALLREDUCE(MPI_IN_PLACE, avg_c, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, avg_c = avg_c / denum -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_d2dx_c, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) +call MPI_ALLREDUCE(MPI_IN_PLACE, avg_divn, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) -avg_d2dx_c = avg_d2dx_c / denum +avg_divn = avg_divn / denum -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_d2dy_c, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) +call MPI_ALLREDUCE(MPI_IN_PLACE, avg_fsd, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) -avg_d2dy_c = avg_d2dy_c / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_d2dz_c, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_d2dz_c = avg_d2dz_c / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_ddx_c, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_ddx_c = avg_ddx_c / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_ddy_c, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_ddy_c = avg_ddy_c / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_ddz_c, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_ddz_c = avg_ddz_c / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_fsd1, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_fsd1 = avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_fsd2, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_fsd2 = avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_fsd3, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_fsd3 = avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_fsd4, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_fsd4 = avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_fsd5, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_fsd5 = avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_fsd6, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_fsd6 = avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_fsd7, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_fsd7 = avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_fsd8, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_fsd8 = avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_fsd9, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_fsd9 = avg_fsd9 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_lapc, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_lapc = avg_lapc / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_u, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_u = avg_u / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_v, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_v = avg_v / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, avg_w, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -avg_w = avg_w / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd1_avg_absk, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd1_avg_absk = fsd1_avg_absk / avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd1_avg_divn, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd1_avg_divn = fsd1_avg_divn / avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd1_avg_nx, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd1_avg_nx = fsd1_avg_nx / avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd1_avg_ny, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd1_avg_ny = fsd1_avg_ny / avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd1_avg_nz, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd1_avg_nz = fsd1_avg_nz / avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd1_avg_sd, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd1_avg_sd = fsd1_avg_sd / avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd1_avg_t1, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd1_avg_t1 = fsd1_avg_t1 / avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd2_avg_absk, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd2_avg_absk = fsd2_avg_absk / avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd2_avg_divn, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd2_avg_divn = fsd2_avg_divn / avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd2_avg_nx, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd2_avg_nx = fsd2_avg_nx / avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd2_avg_ny, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd2_avg_ny = fsd2_avg_ny / avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd2_avg_nz, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd2_avg_nz = fsd2_avg_nz / avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd2_avg_sd, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd2_avg_sd = fsd2_avg_sd / avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd2_avg_t1, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd2_avg_t1 = fsd2_avg_t1 / avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd3_avg_absk, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd3_avg_absk = fsd3_avg_absk / avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd3_avg_divn, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd3_avg_divn = fsd3_avg_divn / avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd3_avg_nx, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd3_avg_nx = fsd3_avg_nx / avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd3_avg_ny, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd3_avg_ny = fsd3_avg_ny / avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd3_avg_nz, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd3_avg_nz = fsd3_avg_nz / avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd3_avg_sd, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd3_avg_sd = fsd3_avg_sd / avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd3_avg_t1, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd3_avg_t1 = fsd3_avg_t1 / avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd4_avg_absk, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd4_avg_absk = fsd4_avg_absk / avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd4_avg_divn, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd4_avg_divn = fsd4_avg_divn / avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd4_avg_nx, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd4_avg_nx = fsd4_avg_nx / avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd4_avg_ny, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd4_avg_ny = fsd4_avg_ny / avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd4_avg_nz, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd4_avg_nz = fsd4_avg_nz / avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd4_avg_sd, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd4_avg_sd = fsd4_avg_sd / avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd4_avg_t1, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd4_avg_t1 = fsd4_avg_t1 / avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd5_avg_absk, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd5_avg_absk = fsd5_avg_absk / avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd5_avg_divn, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd5_avg_divn = fsd5_avg_divn / avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd5_avg_nx, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd5_avg_nx = fsd5_avg_nx / avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd5_avg_ny, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd5_avg_ny = fsd5_avg_ny / avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd5_avg_nz, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd5_avg_nz = fsd5_avg_nz / avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd5_avg_sd, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd5_avg_sd = fsd5_avg_sd / avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd5_avg_t1, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd5_avg_t1 = fsd5_avg_t1 / avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd6_avg_absk, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd6_avg_absk = fsd6_avg_absk / avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd6_avg_divn, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd6_avg_divn = fsd6_avg_divn / avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd6_avg_nx, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd6_avg_nx = fsd6_avg_nx / avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd6_avg_ny, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd6_avg_ny = fsd6_avg_ny / avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd6_avg_nz, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd6_avg_nz = fsd6_avg_nz / avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd6_avg_sd, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd6_avg_sd = fsd6_avg_sd / avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd6_avg_t1, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd6_avg_t1 = fsd6_avg_t1 / avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd7_avg_absk, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd7_avg_absk = fsd7_avg_absk / avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd7_avg_divn, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd7_avg_divn = fsd7_avg_divn / avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd7_avg_nx, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd7_avg_nx = fsd7_avg_nx / avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd7_avg_ny, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd7_avg_ny = fsd7_avg_ny / avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd7_avg_nz, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd7_avg_nz = fsd7_avg_nz / avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd7_avg_sd, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd7_avg_sd = fsd7_avg_sd / avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd7_avg_t1, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd7_avg_t1 = fsd7_avg_t1 / avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd8_avg_absk, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd8_avg_absk = fsd8_avg_absk / avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd8_avg_divn, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd8_avg_divn = fsd8_avg_divn / avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd8_avg_nx, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd8_avg_nx = fsd8_avg_nx / avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd8_avg_ny, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd8_avg_ny = fsd8_avg_ny / avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd8_avg_nz, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd8_avg_nz = fsd8_avg_nz / avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd8_avg_sd, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd8_avg_sd = fsd8_avg_sd / avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd8_avg_t1, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd8_avg_t1 = fsd8_avg_t1 / avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd9_avg_absk, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd9_avg_absk = fsd9_avg_absk / avg_fsd9 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd9_avg_divn, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd9_avg_divn = fsd9_avg_divn / avg_fsd9 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd9_avg_nx, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd9_avg_nx = fsd9_avg_nx / avg_fsd9 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd9_avg_ny, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd9_avg_ny = fsd9_avg_ny / avg_fsd9 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd9_avg_nz, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd9_avg_nz = fsd9_avg_nz / avg_fsd9 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd9_avg_sd, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd9_avg_sd = fsd9_avg_sd / avg_fsd9 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd9_avg_t1, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd9_avg_t1 = fsd9_avg_t1 / avg_fsd9 / denum +avg_fsd = avg_fsd / denum end subroutine m_terms_average_pass1 @@ -2018,1501 +182,6 @@ subroutine m_terms_calculate_pass2 integer :: i, j, k -! ( 1.0 - y(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer18(i,j,k) = ( 1.0 - y(i,j,k) ) -end do -end do -end do - -call d2dx ( xyzbuffer13, xyzbuffer18 ) -call d2dy ( xyzbuffer12, xyzbuffer18 ) -call d2dz ( xyzbuffer11, xyzbuffer18 ) -call ddx ( xyzbuffer10, xyzbuffer18 ) -call ddy ( xyzbuffer17, xyzbuffer18 ) -call ddz ( xyzbuffer16, xyzbuffer18 ) - -! ( 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))) ) ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer15(i,j,k) = ( sqrt ( ( ( ((xyzbuffer10(i,j,k))*(xyzbuffer10(i,j,k))) + ((xyzbuffer17(i,j,k))*(xyzbuffer17(i,j,k))) ) + ((xyzbuffer16(i,j,k))*(xyzbuffer16(i,j,k))) ) ) ) -end do -end do -end do - - -! ( ( - ddx_c(i,j,k) ) / fsd_orig(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer14(i,j,k) = ( ( - xyzbuffer10(i,j,k) ) / xyzbuffer15(i,j,k) ) -end do -end do -end do - - -! ( ( - ddy_c(i,j,k) ) / fsd_orig(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer9(i,j,k) = ( ( - xyzbuffer17(i,j,k) ) / xyzbuffer15(i,j,k) ) -end do -end do -end do - - -! ( ( - ddz_c(i,j,k) ) / fsd_orig(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer8(i,j,k) = ( ( - xyzbuffer16(i,j,k) ) / xyzbuffer15(i,j,k) ) -end do -end do -end do - -call ddx ( xyzbuffer3, xyzbuffer14 ) -call ddy ( xyzbuffer2, xyzbuffer9 ) - -! ( ( ( u(i,j,k) * nx(i,j,k) ) + ( v(i,j,k) * ny(i,j,k) ) ) + ( w(i,j,k) * nz(i,j,k) ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer1(i,j,k) = ( ( ( u(i,j,k) * xyzbuffer14(i,j,k) ) + ( v(i,j,k) * xyzbuffer9(i,j,k) ) ) + ( w(i,j,k) * xyzbuffer8(i,j,k) ) ) -end do -end do -end do - -call ddz ( xyzbuffer0, xyzbuffer8 ) - -! ( ( ddx_nx(i,j,k) + ddy_ny(i,j,k) ) + ddz_nz(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer7(i,j,k) = ( ( xyzbuffer3(i,j,k) + xyzbuffer2(i,j,k) ) + xyzbuffer0(i,j,k) ) -end do -end do -end do - - -! ( dabs ( divn(i,j,k) ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer6(i,j,k) = ( dabs ( xyzbuffer7(i,j,k) ) ) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.05, 0.15 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer5(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.05, 0.15 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask1(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer4(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer5(i,j,k) ) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.15, 0.25 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer10(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.15, 0.25 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask2(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer17(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer10(i,j,k) ) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.25, 0.35 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer16(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.25, 0.35 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask3(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer9(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer16(i,j,k) ) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.35, 0.45 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer8(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.35, 0.45 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask4(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer3(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer8(i,j,k) ) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.45, 0.55 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer2(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.45, 0.55 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask5(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer0(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer2(i,j,k) ) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.55, 0.65 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer5(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.55, 0.65 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask6(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer10(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer5(i,j,k) ) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.65, 0.75 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer16(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.65, 0.75 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask7(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer8(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer16(i,j,k) ) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.75, 0.85 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer2(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.75, 0.85 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask8(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer5(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer2(i,j,k) ) -end do -end do -end do - - -! ( threshold_min_max ( c(i,j,k), 0.85, 0.95 ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer16(i,j,k) = ( threshold_min_max ( xyzbuffer18(i,j,k), 0.85, 0.95 ) ) -end do -end do -end do - - -! ( fsd_orig(i,j,k) * mask9(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer2(i,j,k) = ( xyzbuffer15(i,j,k) * xyzbuffer16(i,j,k) ) -end do -end do -end do - - -! ( rxn_rate ( c(i,j,k) ) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer16(i,j,k) = ( rxn_rate ( xyzbuffer18(i,j,k) ) ) -end do -end do -end do - - -! ( ( ( ( ( d2dx_c(i,j,k) + d2dy_c(i,j,k) ) + d2dz_c(i,j,k) ) * rod ) + wrate(i,j,k) ) / fsd_orig(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer18(i,j,k) = ( ( ( ( ( xyzbuffer13(i,j,k) + xyzbuffer12(i,j,k) ) + xyzbuffer11(i,j,k) ) * rod ) + xyzbuffer16(i,j,k) ) / xyzbuffer15(i,j,k) ) -end do -end do -end do - - -! ( ( vn(i,j,k) + sd(i,j,k) ) * nx(i,j,k) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer13(i,j,k) = ( ( xyzbuffer1(i,j,k) + xyzbuffer18(i,j,k) ) * xyzbuffer14(i,j,k) ) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd1_avg_t1(i)) * (nx(i,j,k) - fsd1_avg_nx(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer12(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd1_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd1_avg_nx(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd1_avg_t2(i) = fsd1_avg_t2(i) + xyzbuffer12(i,j,k) * xyzbuffer4(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd2_avg_t1(i)) * (nx(i,j,k) - fsd2_avg_nx(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer11(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd2_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd2_avg_nx(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd2_avg_t2(i) = fsd2_avg_t2(i) + xyzbuffer11(i,j,k) * xyzbuffer17(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd3_avg_t1(i)) * (nx(i,j,k) - fsd3_avg_nx(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer16(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd3_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd3_avg_nx(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd3_avg_t2(i) = fsd3_avg_t2(i) + xyzbuffer16(i,j,k) * xyzbuffer9(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd4_avg_t1(i)) * (nx(i,j,k) - fsd4_avg_nx(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer15(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd4_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd4_avg_nx(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd4_avg_t2(i) = fsd4_avg_t2(i) + xyzbuffer15(i,j,k) * xyzbuffer3(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd5_avg_t1(i)) * (nx(i,j,k) - fsd5_avg_nx(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer1(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd5_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd5_avg_nx(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd5_avg_t2(i) = fsd5_avg_t2(i) + xyzbuffer1(i,j,k) * xyzbuffer0(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd6_avg_t1(i)) * (nx(i,j,k) - fsd6_avg_nx(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer18(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd6_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd6_avg_nx(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd6_avg_t2(i) = fsd6_avg_t2(i) + xyzbuffer18(i,j,k) * xyzbuffer10(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd7_avg_t1(i)) * (nx(i,j,k) - fsd7_avg_nx(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer12(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd7_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd7_avg_nx(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd7_avg_t2(i) = fsd7_avg_t2(i) + xyzbuffer12(i,j,k) * xyzbuffer8(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd8_avg_t1(i)) * (nx(i,j,k) - fsd8_avg_nx(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer11(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd8_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd8_avg_nx(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd8_avg_t2(i) = fsd8_avg_t2(i) + xyzbuffer11(i,j,k) * xyzbuffer5(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd9_avg_t1(i)) * (nx(i,j,k) - fsd9_avg_nx(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer16(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd9_avg_t1(i)) * (xyzbuffer14(i,j,k) - fsd9_avg_nx(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd9_avg_t2(i) = fsd9_avg_t2(i) + xyzbuffer16(i,j,k) * xyzbuffer2(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd1_avg_t1(i)) * (divn(i,j,k) - fsd1_avg_divn(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer15(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd1_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd1_avg_divn(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd1_avg_t3(i) = fsd1_avg_t3(i) + xyzbuffer15(i,j,k) * xyzbuffer4(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd2_avg_t1(i)) * (divn(i,j,k) - fsd2_avg_divn(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer1(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd2_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd2_avg_divn(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd2_avg_t3(i) = fsd2_avg_t3(i) + xyzbuffer1(i,j,k) * xyzbuffer17(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd3_avg_t1(i)) * (divn(i,j,k) - fsd3_avg_divn(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer18(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd3_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd3_avg_divn(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd3_avg_t3(i) = fsd3_avg_t3(i) + xyzbuffer18(i,j,k) * xyzbuffer9(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd4_avg_t1(i)) * (divn(i,j,k) - fsd4_avg_divn(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer12(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd4_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd4_avg_divn(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd4_avg_t3(i) = fsd4_avg_t3(i) + xyzbuffer12(i,j,k) * xyzbuffer3(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd5_avg_t1(i)) * (divn(i,j,k) - fsd5_avg_divn(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer11(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd5_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd5_avg_divn(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd5_avg_t3(i) = fsd5_avg_t3(i) + xyzbuffer11(i,j,k) * xyzbuffer0(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd6_avg_t1(i)) * (divn(i,j,k) - fsd6_avg_divn(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer16(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd6_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd6_avg_divn(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd6_avg_t3(i) = fsd6_avg_t3(i) + xyzbuffer16(i,j,k) * xyzbuffer10(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd7_avg_t1(i)) * (divn(i,j,k) - fsd7_avg_divn(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer15(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd7_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd7_avg_divn(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd7_avg_t3(i) = fsd7_avg_t3(i) + xyzbuffer15(i,j,k) * xyzbuffer8(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd8_avg_t1(i)) * (divn(i,j,k) - fsd8_avg_divn(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer1(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd8_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd8_avg_divn(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd8_avg_t3(i) = fsd8_avg_t3(i) + xyzbuffer1(i,j,k) * xyzbuffer5(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd9_avg_t1(i)) * (divn(i,j,k) - fsd9_avg_divn(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer18(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd9_avg_t1(i)) * (xyzbuffer7(i,j,k) - fsd9_avg_divn(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd9_avg_t3(i) = fsd9_avg_t3(i) + xyzbuffer18(i,j,k) * xyzbuffer2(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd1_avg_t1(i)) * (absk(i,j,k) - fsd1_avg_absk(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer12(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd1_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd1_avg_absk(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd1_avg_t4(i) = fsd1_avg_t4(i) + xyzbuffer12(i,j,k) * xyzbuffer4(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd2_avg_t1(i)) * (absk(i,j,k) - fsd2_avg_absk(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer11(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd2_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd2_avg_absk(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd2_avg_t4(i) = fsd2_avg_t4(i) + xyzbuffer11(i,j,k) * xyzbuffer17(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd3_avg_t1(i)) * (absk(i,j,k) - fsd3_avg_absk(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer16(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd3_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd3_avg_absk(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd3_avg_t4(i) = fsd3_avg_t4(i) + xyzbuffer16(i,j,k) * xyzbuffer9(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd4_avg_t1(i)) * (absk(i,j,k) - fsd4_avg_absk(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer15(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd4_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd4_avg_absk(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd4_avg_t4(i) = fsd4_avg_t4(i) + xyzbuffer15(i,j,k) * xyzbuffer3(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd5_avg_t1(i)) * (absk(i,j,k) - fsd5_avg_absk(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer1(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd5_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd5_avg_absk(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd5_avg_t4(i) = fsd5_avg_t4(i) + xyzbuffer1(i,j,k) * xyzbuffer0(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd6_avg_t1(i)) * (absk(i,j,k) - fsd6_avg_absk(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer18(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd6_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd6_avg_absk(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd6_avg_t4(i) = fsd6_avg_t4(i) + xyzbuffer18(i,j,k) * xyzbuffer10(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd7_avg_t1(i)) * (absk(i,j,k) - fsd7_avg_absk(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer12(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd7_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd7_avg_absk(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd7_avg_t4(i) = fsd7_avg_t4(i) + xyzbuffer12(i,j,k) * xyzbuffer8(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd8_avg_t1(i)) * (absk(i,j,k) - fsd8_avg_absk(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer11(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd8_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd8_avg_absk(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd8_avg_t4(i) = fsd8_avg_t4(i) + xyzbuffer11(i,j,k) * xyzbuffer5(i,j,k) -end do -end do -end do - - -! ( (t1(i,j,k) - fsd9_avg_t1(i)) * (absk(i,j,k) - fsd9_avg_absk(i)) ) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer16(i,j,k) = ( (xyzbuffer13(i,j,k) - fsd9_avg_t1(i)) * (xyzbuffer6(i,j,k) - fsd9_avg_absk(i)) ) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd9_avg_t4(i) = fsd9_avg_t4(i) + xyzbuffer16(i,j,k) * xyzbuffer2(i,j,k) -end do -end do -end do - - -! (((t1(i,j,k) - fsd1_avg_t1(i)))*((t1(i,j,k) - fsd1_avg_t1(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer15(i,j,k) = (((xyzbuffer13(i,j,k) - fsd1_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd1_avg_t1(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd1_avg_t5(i) = fsd1_avg_t5(i) + xyzbuffer15(i,j,k) * xyzbuffer4(i,j,k) -end do -end do -end do - - -! (((t1(i,j,k) - fsd2_avg_t1(i)))*((t1(i,j,k) - fsd2_avg_t1(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer1(i,j,k) = (((xyzbuffer13(i,j,k) - fsd2_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd2_avg_t1(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd2_avg_t5(i) = fsd2_avg_t5(i) + xyzbuffer1(i,j,k) * xyzbuffer17(i,j,k) -end do -end do -end do - - -! (((t1(i,j,k) - fsd3_avg_t1(i)))*((t1(i,j,k) - fsd3_avg_t1(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer18(i,j,k) = (((xyzbuffer13(i,j,k) - fsd3_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd3_avg_t1(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd3_avg_t5(i) = fsd3_avg_t5(i) + xyzbuffer18(i,j,k) * xyzbuffer9(i,j,k) -end do -end do -end do - - -! (((t1(i,j,k) - fsd4_avg_t1(i)))*((t1(i,j,k) - fsd4_avg_t1(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer12(i,j,k) = (((xyzbuffer13(i,j,k) - fsd4_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd4_avg_t1(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd4_avg_t5(i) = fsd4_avg_t5(i) + xyzbuffer12(i,j,k) * xyzbuffer3(i,j,k) -end do -end do -end do - - -! (((t1(i,j,k) - fsd5_avg_t1(i)))*((t1(i,j,k) - fsd5_avg_t1(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer11(i,j,k) = (((xyzbuffer13(i,j,k) - fsd5_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd5_avg_t1(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd5_avg_t5(i) = fsd5_avg_t5(i) + xyzbuffer11(i,j,k) * xyzbuffer0(i,j,k) -end do -end do -end do - - -! (((t1(i,j,k) - fsd6_avg_t1(i)))*((t1(i,j,k) - fsd6_avg_t1(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer16(i,j,k) = (((xyzbuffer13(i,j,k) - fsd6_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd6_avg_t1(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd6_avg_t5(i) = fsd6_avg_t5(i) + xyzbuffer16(i,j,k) * xyzbuffer10(i,j,k) -end do -end do -end do - - -! (((t1(i,j,k) - fsd7_avg_t1(i)))*((t1(i,j,k) - fsd7_avg_t1(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer15(i,j,k) = (((xyzbuffer13(i,j,k) - fsd7_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd7_avg_t1(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd7_avg_t5(i) = fsd7_avg_t5(i) + xyzbuffer15(i,j,k) * xyzbuffer8(i,j,k) -end do -end do -end do - - -! (((t1(i,j,k) - fsd8_avg_t1(i)))*((t1(i,j,k) - fsd8_avg_t1(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer1(i,j,k) = (((xyzbuffer13(i,j,k) - fsd8_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd8_avg_t1(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd8_avg_t5(i) = fsd8_avg_t5(i) + xyzbuffer1(i,j,k) * xyzbuffer5(i,j,k) -end do -end do -end do - - -! (((t1(i,j,k) - fsd9_avg_t1(i)))*((t1(i,j,k) - fsd9_avg_t1(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer18(i,j,k) = (((xyzbuffer13(i,j,k) - fsd9_avg_t1(i)))*((xyzbuffer13(i,j,k) - fsd9_avg_t1(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd9_avg_t5(i) = fsd9_avg_t5(i) + xyzbuffer18(i,j,k) * xyzbuffer2(i,j,k) -end do -end do -end do - - -! (((nx(i,j,k) - fsd1_avg_nx(i)))*((nx(i,j,k) - fsd1_avg_nx(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer13(i,j,k) = (((xyzbuffer14(i,j,k) - fsd1_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd1_avg_nx(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd1_avg_t6(i) = fsd1_avg_t6(i) + xyzbuffer13(i,j,k) * xyzbuffer4(i,j,k) -end do -end do -end do - - -! (((nx(i,j,k) - fsd2_avg_nx(i)))*((nx(i,j,k) - fsd2_avg_nx(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer12(i,j,k) = (((xyzbuffer14(i,j,k) - fsd2_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd2_avg_nx(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd2_avg_t6(i) = fsd2_avg_t6(i) + xyzbuffer12(i,j,k) * xyzbuffer17(i,j,k) -end do -end do -end do - - -! (((nx(i,j,k) - fsd3_avg_nx(i)))*((nx(i,j,k) - fsd3_avg_nx(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer11(i,j,k) = (((xyzbuffer14(i,j,k) - fsd3_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd3_avg_nx(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd3_avg_t6(i) = fsd3_avg_t6(i) + xyzbuffer11(i,j,k) * xyzbuffer9(i,j,k) -end do -end do -end do - - -! (((nx(i,j,k) - fsd4_avg_nx(i)))*((nx(i,j,k) - fsd4_avg_nx(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer16(i,j,k) = (((xyzbuffer14(i,j,k) - fsd4_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd4_avg_nx(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd4_avg_t6(i) = fsd4_avg_t6(i) + xyzbuffer16(i,j,k) * xyzbuffer3(i,j,k) -end do -end do -end do - - -! (((nx(i,j,k) - fsd5_avg_nx(i)))*((nx(i,j,k) - fsd5_avg_nx(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer15(i,j,k) = (((xyzbuffer14(i,j,k) - fsd5_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd5_avg_nx(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd5_avg_t6(i) = fsd5_avg_t6(i) + xyzbuffer15(i,j,k) * xyzbuffer0(i,j,k) -end do -end do -end do - - -! (((nx(i,j,k) - fsd6_avg_nx(i)))*((nx(i,j,k) - fsd6_avg_nx(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer1(i,j,k) = (((xyzbuffer14(i,j,k) - fsd6_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd6_avg_nx(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd6_avg_t6(i) = fsd6_avg_t6(i) + xyzbuffer1(i,j,k) * xyzbuffer10(i,j,k) -end do -end do -end do - - -! (((nx(i,j,k) - fsd7_avg_nx(i)))*((nx(i,j,k) - fsd7_avg_nx(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer18(i,j,k) = (((xyzbuffer14(i,j,k) - fsd7_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd7_avg_nx(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd7_avg_t6(i) = fsd7_avg_t6(i) + xyzbuffer18(i,j,k) * xyzbuffer8(i,j,k) -end do -end do -end do - - -! (((nx(i,j,k) - fsd8_avg_nx(i)))*((nx(i,j,k) - fsd8_avg_nx(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer13(i,j,k) = (((xyzbuffer14(i,j,k) - fsd8_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd8_avg_nx(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd8_avg_t6(i) = fsd8_avg_t6(i) + xyzbuffer13(i,j,k) * xyzbuffer5(i,j,k) -end do -end do -end do - - -! (((nx(i,j,k) - fsd9_avg_nx(i)))*((nx(i,j,k) - fsd9_avg_nx(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer12(i,j,k) = (((xyzbuffer14(i,j,k) - fsd9_avg_nx(i)))*((xyzbuffer14(i,j,k) - fsd9_avg_nx(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd9_avg_t6(i) = fsd9_avg_t6(i) + xyzbuffer12(i,j,k) * xyzbuffer2(i,j,k) -end do -end do -end do - - -! (((divn(i,j,k) - fsd1_avg_divn(i)))*((divn(i,j,k) - fsd1_avg_divn(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer11(i,j,k) = (((xyzbuffer7(i,j,k) - fsd1_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd1_avg_divn(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd1_avg_t7(i) = fsd1_avg_t7(i) + xyzbuffer11(i,j,k) * xyzbuffer4(i,j,k) -end do -end do -end do - - -! (((divn(i,j,k) - fsd2_avg_divn(i)))*((divn(i,j,k) - fsd2_avg_divn(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer16(i,j,k) = (((xyzbuffer7(i,j,k) - fsd2_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd2_avg_divn(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd2_avg_t7(i) = fsd2_avg_t7(i) + xyzbuffer16(i,j,k) * xyzbuffer17(i,j,k) -end do -end do -end do - - -! (((divn(i,j,k) - fsd3_avg_divn(i)))*((divn(i,j,k) - fsd3_avg_divn(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer15(i,j,k) = (((xyzbuffer7(i,j,k) - fsd3_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd3_avg_divn(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd3_avg_t7(i) = fsd3_avg_t7(i) + xyzbuffer15(i,j,k) * xyzbuffer9(i,j,k) -end do -end do -end do - - -! (((divn(i,j,k) - fsd4_avg_divn(i)))*((divn(i,j,k) - fsd4_avg_divn(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer14(i,j,k) = (((xyzbuffer7(i,j,k) - fsd4_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd4_avg_divn(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd4_avg_t7(i) = fsd4_avg_t7(i) + xyzbuffer14(i,j,k) * xyzbuffer3(i,j,k) -end do -end do -end do - - -! (((divn(i,j,k) - fsd5_avg_divn(i)))*((divn(i,j,k) - fsd5_avg_divn(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer1(i,j,k) = (((xyzbuffer7(i,j,k) - fsd5_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd5_avg_divn(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd5_avg_t7(i) = fsd5_avg_t7(i) + xyzbuffer1(i,j,k) * xyzbuffer0(i,j,k) -end do -end do -end do - - -! (((divn(i,j,k) - fsd6_avg_divn(i)))*((divn(i,j,k) - fsd6_avg_divn(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer18(i,j,k) = (((xyzbuffer7(i,j,k) - fsd6_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd6_avg_divn(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd6_avg_t7(i) = fsd6_avg_t7(i) + xyzbuffer18(i,j,k) * xyzbuffer10(i,j,k) -end do -end do -end do - - -! (((divn(i,j,k) - fsd7_avg_divn(i)))*((divn(i,j,k) - fsd7_avg_divn(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer13(i,j,k) = (((xyzbuffer7(i,j,k) - fsd7_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd7_avg_divn(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd7_avg_t7(i) = fsd7_avg_t7(i) + xyzbuffer13(i,j,k) * xyzbuffer8(i,j,k) -end do -end do -end do - - -! (((divn(i,j,k) - fsd8_avg_divn(i)))*((divn(i,j,k) - fsd8_avg_divn(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer12(i,j,k) = (((xyzbuffer7(i,j,k) - fsd8_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd8_avg_divn(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd8_avg_t7(i) = fsd8_avg_t7(i) + xyzbuffer12(i,j,k) * xyzbuffer5(i,j,k) -end do -end do -end do - - -! (((divn(i,j,k) - fsd9_avg_divn(i)))*((divn(i,j,k) - fsd9_avg_divn(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer11(i,j,k) = (((xyzbuffer7(i,j,k) - fsd9_avg_divn(i)))*((xyzbuffer7(i,j,k) - fsd9_avg_divn(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd9_avg_t7(i) = fsd9_avg_t7(i) + xyzbuffer11(i,j,k) * xyzbuffer2(i,j,k) -end do -end do -end do - - -! (((absk(i,j,k) - fsd1_avg_absk(i)))*((absk(i,j,k) - fsd1_avg_absk(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer16(i,j,k) = (((xyzbuffer6(i,j,k) - fsd1_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd1_avg_absk(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd1_avg_t8(i) = fsd1_avg_t8(i) + xyzbuffer16(i,j,k) * xyzbuffer4(i,j,k) -end do -end do -end do - - -! (((absk(i,j,k) - fsd2_avg_absk(i)))*((absk(i,j,k) - fsd2_avg_absk(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer15(i,j,k) = (((xyzbuffer6(i,j,k) - fsd2_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd2_avg_absk(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd2_avg_t8(i) = fsd2_avg_t8(i) + xyzbuffer15(i,j,k) * xyzbuffer17(i,j,k) -end do -end do -end do - - -! (((absk(i,j,k) - fsd3_avg_absk(i)))*((absk(i,j,k) - fsd3_avg_absk(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer14(i,j,k) = (((xyzbuffer6(i,j,k) - fsd3_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd3_avg_absk(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd3_avg_t8(i) = fsd3_avg_t8(i) + xyzbuffer14(i,j,k) * xyzbuffer9(i,j,k) -end do -end do -end do - - -! (((absk(i,j,k) - fsd4_avg_absk(i)))*((absk(i,j,k) - fsd4_avg_absk(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer9(i,j,k) = (((xyzbuffer6(i,j,k) - fsd4_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd4_avg_absk(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd4_avg_t8(i) = fsd4_avg_t8(i) + xyzbuffer9(i,j,k) * xyzbuffer3(i,j,k) -end do -end do -end do - - -! (((absk(i,j,k) - fsd5_avg_absk(i)))*((absk(i,j,k) - fsd5_avg_absk(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer3(i,j,k) = (((xyzbuffer6(i,j,k) - fsd5_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd5_avg_absk(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd5_avg_t8(i) = fsd5_avg_t8(i) + xyzbuffer3(i,j,k) * xyzbuffer0(i,j,k) -end do -end do -end do - - -! (((absk(i,j,k) - fsd6_avg_absk(i)))*((absk(i,j,k) - fsd6_avg_absk(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer1(i,j,k) = (((xyzbuffer6(i,j,k) - fsd6_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd6_avg_absk(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd6_avg_t8(i) = fsd6_avg_t8(i) + xyzbuffer1(i,j,k) * xyzbuffer10(i,j,k) -end do -end do -end do - - -! (((absk(i,j,k) - fsd7_avg_absk(i)))*((absk(i,j,k) - fsd7_avg_absk(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer0(i,j,k) = (((xyzbuffer6(i,j,k) - fsd7_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd7_avg_absk(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd7_avg_t8(i) = fsd7_avg_t8(i) + xyzbuffer0(i,j,k) * xyzbuffer8(i,j,k) -end do -end do -end do - - -! (((absk(i,j,k) - fsd8_avg_absk(i)))*((absk(i,j,k) - fsd8_avg_absk(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer7(i,j,k) = (((xyzbuffer6(i,j,k) - fsd8_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd8_avg_absk(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd8_avg_t8(i) = fsd8_avg_t8(i) + xyzbuffer7(i,j,k) * xyzbuffer5(i,j,k) -end do -end do -end do - - -! (((absk(i,j,k) - fsd9_avg_absk(i)))*((absk(i,j,k) - fsd9_avg_absk(i)))) -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -xyzbuffer5(i,j,k) = (((xyzbuffer6(i,j,k) - fsd9_avg_absk(i)))*((xyzbuffer6(i,j,k) - fsd9_avg_absk(i)))) -end do -end do -end do - - -do k = 1, nzp -do j = 1, nyp -do i = 1, nxp -fsd9_avg_t8(i) = fsd9_avg_t8(i) + xyzbuffer5(i,j,k) * xyzbuffer2(i,j,k) -end do -end do -end do - end subroutine m_terms_calculate_pass2 @@ -3525,320 +194,6 @@ real*8 :: denum denum=real(nfiles*nyp*nzp) -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd1_avg_t2, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd1_avg_t2 = fsd1_avg_t2 / avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd1_avg_t3, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd1_avg_t3 = fsd1_avg_t3 / avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd1_avg_t4, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd1_avg_t4 = fsd1_avg_t4 / avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd1_avg_t5, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd1_avg_t5 = fsd1_avg_t5 / avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd1_avg_t6, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd1_avg_t6 = fsd1_avg_t6 / avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd1_avg_t7, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd1_avg_t7 = fsd1_avg_t7 / avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd1_avg_t8, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd1_avg_t8 = fsd1_avg_t8 / avg_fsd1 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd2_avg_t2, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd2_avg_t2 = fsd2_avg_t2 / avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd2_avg_t3, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd2_avg_t3 = fsd2_avg_t3 / avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd2_avg_t4, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd2_avg_t4 = fsd2_avg_t4 / avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd2_avg_t5, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd2_avg_t5 = fsd2_avg_t5 / avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd2_avg_t6, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd2_avg_t6 = fsd2_avg_t6 / avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd2_avg_t7, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd2_avg_t7 = fsd2_avg_t7 / avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd2_avg_t8, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd2_avg_t8 = fsd2_avg_t8 / avg_fsd2 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd3_avg_t2, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd3_avg_t2 = fsd3_avg_t2 / avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd3_avg_t3, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd3_avg_t3 = fsd3_avg_t3 / avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd3_avg_t4, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd3_avg_t4 = fsd3_avg_t4 / avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd3_avg_t5, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd3_avg_t5 = fsd3_avg_t5 / avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd3_avg_t6, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd3_avg_t6 = fsd3_avg_t6 / avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd3_avg_t7, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd3_avg_t7 = fsd3_avg_t7 / avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd3_avg_t8, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd3_avg_t8 = fsd3_avg_t8 / avg_fsd3 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd4_avg_t2, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd4_avg_t2 = fsd4_avg_t2 / avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd4_avg_t3, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd4_avg_t3 = fsd4_avg_t3 / avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd4_avg_t4, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd4_avg_t4 = fsd4_avg_t4 / avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd4_avg_t5, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd4_avg_t5 = fsd4_avg_t5 / avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd4_avg_t6, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd4_avg_t6 = fsd4_avg_t6 / avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd4_avg_t7, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd4_avg_t7 = fsd4_avg_t7 / avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd4_avg_t8, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd4_avg_t8 = fsd4_avg_t8 / avg_fsd4 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd5_avg_t2, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd5_avg_t2 = fsd5_avg_t2 / avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd5_avg_t3, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd5_avg_t3 = fsd5_avg_t3 / avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd5_avg_t4, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd5_avg_t4 = fsd5_avg_t4 / avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd5_avg_t5, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd5_avg_t5 = fsd5_avg_t5 / avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd5_avg_t6, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd5_avg_t6 = fsd5_avg_t6 / avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd5_avg_t7, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd5_avg_t7 = fsd5_avg_t7 / avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd5_avg_t8, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd5_avg_t8 = fsd5_avg_t8 / avg_fsd5 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd6_avg_t2, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd6_avg_t2 = fsd6_avg_t2 / avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd6_avg_t3, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd6_avg_t3 = fsd6_avg_t3 / avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd6_avg_t4, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd6_avg_t4 = fsd6_avg_t4 / avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd6_avg_t5, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd6_avg_t5 = fsd6_avg_t5 / avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd6_avg_t6, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd6_avg_t6 = fsd6_avg_t6 / avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd6_avg_t7, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd6_avg_t7 = fsd6_avg_t7 / avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd6_avg_t8, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd6_avg_t8 = fsd6_avg_t8 / avg_fsd6 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd7_avg_t2, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd7_avg_t2 = fsd7_avg_t2 / avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd7_avg_t3, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd7_avg_t3 = fsd7_avg_t3 / avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd7_avg_t4, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd7_avg_t4 = fsd7_avg_t4 / avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd7_avg_t5, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd7_avg_t5 = fsd7_avg_t5 / avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd7_avg_t6, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd7_avg_t6 = fsd7_avg_t6 / avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd7_avg_t7, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd7_avg_t7 = fsd7_avg_t7 / avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd7_avg_t8, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd7_avg_t8 = fsd7_avg_t8 / avg_fsd7 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd8_avg_t2, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd8_avg_t2 = fsd8_avg_t2 / avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd8_avg_t3, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd8_avg_t3 = fsd8_avg_t3 / avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd8_avg_t4, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd8_avg_t4 = fsd8_avg_t4 / avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd8_avg_t5, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd8_avg_t5 = fsd8_avg_t5 / avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd8_avg_t6, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd8_avg_t6 = fsd8_avg_t6 / avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd8_avg_t7, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd8_avg_t7 = fsd8_avg_t7 / avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd8_avg_t8, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd8_avg_t8 = fsd8_avg_t8 / avg_fsd8 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd9_avg_t2, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd9_avg_t2 = fsd9_avg_t2 / avg_fsd9 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd9_avg_t3, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd9_avg_t3 = fsd9_avg_t3 / avg_fsd9 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd9_avg_t4, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd9_avg_t4 = fsd9_avg_t4 / avg_fsd9 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd9_avg_t5, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd9_avg_t5 = fsd9_avg_t5 / avg_fsd9 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd9_avg_t6, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd9_avg_t6 = fsd9_avg_t6 / avg_fsd9 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd9_avg_t7, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd9_avg_t7 = fsd9_avg_t7 / avg_fsd9 / denum - - -call MPI_ALLREDUCE(MPI_IN_PLACE, fsd9_avg_t8, nxp, MPI_REAL8, MPI_SUM, MPI_COMM_TASK, mpi_err) - -fsd9_avg_t8 = fsd9_avg_t8 / avg_fsd9 / denum - end subroutine m_terms_average_pass2 @@ -3853,7 +208,7 @@ open (200, file="qEdge_X.dat") write (200,*) output_header do i=1,nxp - write (200,'(147e20.10)') real(i)*hxp, avg_c(i), avg_d2dx_c(i), avg_d2dy_c(i), avg_d2dz_c(i), avg_ddx_c(i), avg_ddy_c(i), avg_ddz_c(i), avg_fsd1(i), avg_fsd2(i), avg_fsd3(i), avg_fsd4(i), avg_fsd5(i), avg_fsd6(i), avg_fsd7(i), avg_fsd8(i), avg_fsd9(i), avg_lapc(i), avg_u(i), avg_v(i), avg_w(i), fsd1_avg_absk(i), fsd1_avg_divn(i), fsd1_avg_nx(i), fsd1_avg_ny(i), fsd1_avg_nz(i), fsd1_avg_sd(i), fsd1_avg_t1(i), fsd2_avg_absk(i), fsd2_avg_divn(i), fsd2_avg_nx(i), fsd2_avg_ny(i), fsd2_avg_nz(i), fsd2_avg_sd(i), fsd2_avg_t1(i), fsd3_avg_absk(i), fsd3_avg_divn(i), fsd3_avg_nx(i), fsd3_avg_ny(i), fsd3_avg_nz(i), fsd3_avg_sd(i), fsd3_avg_t1(i), fsd4_avg_absk(i), fsd4_avg_divn(i), fsd4_avg_nx(i), fsd4_avg_ny(i), fsd4_avg_nz(i), fsd4_avg_sd(i), fsd4_avg_t1(i), fsd5_avg_absk(i), fsd5_avg_divn(i), fsd5_avg_nx(i), fsd5_avg_ny(i), fsd5_avg_nz(i), fsd5_avg_sd(i), fsd5_avg_t1(i), fsd6_avg_absk(i), fsd6_avg_divn(i), fsd6_avg_nx(i), fsd6_avg_ny(i), fsd6_avg_nz(i), fsd6_avg_sd(i), fsd6_avg_t1(i), fsd7_avg_absk(i), fsd7_avg_divn(i), fsd7_avg_nx(i), fsd7_avg_ny(i), fsd7_avg_nz(i), fsd7_avg_sd(i), fsd7_avg_t1(i), fsd8_avg_absk(i), fsd8_avg_divn(i), fsd8_avg_nx(i), fsd8_avg_ny(i), fsd8_avg_nz(i), fsd8_avg_sd(i), fsd8_avg_t1(i), fsd9_avg_absk(i), fsd9_avg_divn(i), fsd9_avg_nx(i), fsd9_avg_ny(i), fsd9_avg_nz(i), fsd9_avg_sd(i), fsd9_avg_t1(i), fsd1_avg_t2(i), fsd1_avg_t3(i), fsd1_avg_t4(i), fsd1_avg_t5(i), fsd1_avg_t6(i), fsd1_avg_t7(i), fsd1_avg_t8(i), fsd2_avg_t2(i), fsd2_avg_t3(i), fsd2_avg_t4(i), fsd2_avg_t5(i), fsd2_avg_t6(i), fsd2_avg_t7(i), fsd2_avg_t8(i), fsd3_avg_t2(i), fsd3_avg_t3(i), fsd3_avg_t4(i), fsd3_avg_t5(i), fsd3_avg_t6(i), fsd3_avg_t7(i), fsd3_avg_t8(i), fsd4_avg_t2(i), fsd4_avg_t3(i), fsd4_avg_t4(i), fsd4_avg_t5(i), fsd4_avg_t6(i), fsd4_avg_t7(i), fsd4_avg_t8(i), fsd5_avg_t2(i), fsd5_avg_t3(i), fsd5_avg_t4(i), fsd5_avg_t5(i), fsd5_avg_t6(i), fsd5_avg_t7(i), fsd5_avg_t8(i), fsd6_avg_t2(i), fsd6_avg_t3(i), fsd6_avg_t4(i), fsd6_avg_t5(i), fsd6_avg_t6(i), fsd6_avg_t7(i), fsd6_avg_t8(i), fsd7_avg_t2(i), fsd7_avg_t3(i), fsd7_avg_t4(i), fsd7_avg_t5(i), fsd7_avg_t6(i), fsd7_avg_t7(i), fsd7_avg_t8(i), fsd8_avg_t2(i), fsd8_avg_t3(i), fsd8_avg_t4(i), fsd8_avg_t5(i), fsd8_avg_t6(i), fsd8_avg_t7(i), fsd8_avg_t8(i), fsd9_avg_t2(i), fsd9_avg_t3(i), fsd9_avg_t4(i), fsd9_avg_t5(i), fsd9_avg_t6(i), fsd9_avg_t7(i), fsd9_avg_t8(i) + write (200,'(4e20.10)') real(i)*hxp, avg_c(i), avg_divn(i), avg_fsd(i) end do close (200) diff --git a/code/post.f90 b/code/post.f90 index 9ee0d92..546cb8f 100644 --- a/code/post.f90 +++ b/code/post.f90 @@ -59,6 +59,7 @@ CALL m_terms_average_pass1(countnum) + optional_pass2: IF ( pass2_required ) THEN secondloop: DO fread=startnum,endnum,skipnum @@ -86,6 +87,8 @@ CALL m_terms_average_pass2(countnum) + ENDIF optional_pass2 + if (iammaster) then CALL m_terms_write_result WRITE(*,*) ' Avergaing RAW data is FINISHED' @@ -174,6 +177,8 @@ ! PASS 2 + optional_pass2: IF ( pass2_required ) THEN + tag = 2 nprogress = 0 @@ -226,6 +231,8 @@ CALL m_terms_average_pass2(countnum) + ENDIF optional_pass2 + if (iammaster) then CALL m_terms_write_result WRITE(*,*) ' Avergaing RAW data is FINISHED'