@@ -242,15 +242,12 @@ contains
242242 real(wp), dimension(-offset_x%beg:m + offset_x%end, &
243243 -offset_y%beg:n + offset_y%end, &
244244 -offset_z%beg:p + offset_z%end, 3) :: liutex_axis
245- integer :: i, j, k, l, kx, ky, kz, kf, Nf, j_glb, k_glb, l_glb
245+ integer :: i, j, k, l, kx, ky, kz, kf, j_glb, k_glb, l_glb
246246 real(wp) :: En_tot
247247 character(20) :: filename
248248 logical :: file_exists
249249 integer :: x_beg, x_end, y_beg, y_end, z_beg, z_end
250250
251-
252- print *, proc_rank, proc_coords(1), proc_coords(2), proc_coords(3)
253-
254251 if (output_partial_domain) then
255252 call s_define_output_region
256253 x_beg = -offset_x%beg + x_output_idx%beg
@@ -904,6 +901,8 @@ contains
904901 Nyloc2 = (n_glb + 1 ) / num_procs_z
905902 Nzloc = p + 1
906903
904+ Nf = max (Nx, Ny, Nz)
905+
907906 @:ALLOCATE(data_in(Nx* Nyloc* Nzloc))
908907 @:ALLOCATE(data_out(Nx* Nyloc* Nzloc))
909908
@@ -912,7 +911,6 @@ contains
912911 @:ALLOCATE(data_cmplx_z(Nxloc, Nyloc2, Nz))
913912
914913 @:ALLOCATE(En_real(Nxloc, Nyloc2, Nz))
915- Nf = max (Nx, Ny, Nz)
916914 @:ALLOCATE(En(Nf))
917915
918916 size_n(1 ) = Nx
0 commit comments