xepsi, yepsi and zepsi are local variables

Author: mathrack

src/genepsi3d.f90

@@ -65,6 +65,8 @@ subroutine geomcomplex(epsi, nxi, nxf, ny, nyi, nyf, nzi, nzf, dx, yp, dz, remp)
     REAL(mytype),DIMENSION(ny) :: yp
     REAL(mytype)               :: remp
 
+    epsi = 0._mytype
+
     IF (itype.EQ.itype_cyl) THEN
 
        CALL geomcomplex_cyl(epsi, nxi, nxf, ny, nyi, nyf, nzi, nzf, dx, yp, remp)
@@ -159,7 +161,6 @@ subroutine gene_epsi_3D(ep1,nx,ny,nz,dx,dy,dz,xlx,yly,zlz ,&
     use var, only : ta2, ta3
     use decomp_2d
     use MPI
-    use complex_geometry, only: xepsi, yepsi, zepsi
     implicit none
     !
     real(mytype),dimension(xsize(1),xsize(2),xsize(3)), intent(inout):: ep1
@@ -184,9 +185,9 @@ subroutine gene_epsi_3D(ep1,nx,ny,nz,dx,dy,dz,xlx,yly,zlz ,&
     integer,     dimension(xsize(2),xsize(3))          :: nobjxraf
     integer,     dimension(ysize(1),ysize(3))          :: nobjyraf
     integer,     dimension(zsize(1),zsize(2))          :: nobjzraf
-    !real(mytype),dimension(nxraf,xsize(2),xsize(3))    :: xepsi
-    !real(mytype),dimension(ysize(1),nyraf,ysize(3))    :: yepsi
-    !real(mytype),dimension(zsize(1),zsize(2),nzraf)    :: zepsi
+    real(mytype),dimension(nxraf,1,1)                  :: xepsi
+    real(mytype),dimension(1,nyraf,1)                  :: yepsi
+    real(mytype),dimension(1,1,nzraf)                  :: zepsi
     real(mytype),dimension(nyraf)                      :: ypraf
     real(mytype)                     :: dxraf,dyraf,dzraf
     integer                          :: i,j,k

src/module_param.f90

@@ -551,7 +551,6 @@ module complex_geometry
   integer     ,allocatable,dimension(:,:)   :: nobjx,nobjy,nobjz
   integer     ,allocatable,dimension(:,:,:) :: nxipif,nxfpif,nyipif,nyfpif,nzipif,nzfpif
   real(mytype),allocatable,dimension(:,:,:) :: xi,xf,yi,yf,zi,zf
-  real(mytype),allocatable,dimension(:,:,:) :: xepsi, yepsi, zepsi  
   integer :: nxraf,nyraf,nzraf,nraf,nobjmax
 end module complex_geometry
 !############################################################################

src/variables.f90

@@ -656,12 +656,6 @@ subroutine init_variables
        zi=zero
        allocate(zf(nobjmax,zsize(1),zsize(2)))
        zf=zero
-       allocate(xepsi(nxraf,1,1)) 
-       xepsi=zero
-       allocate(yepsi(1,nyraf,1))
-       yepsi=zero
-       allocate(zepsi(1,1,nzraf)) 
-       zepsi=zero
 
     endif