some polishing

src/main.f90

@@ -405,8 +405,7 @@ program cans
         call boundp(cbcpsi,n,bcpre,nb,is_bound,dl,dzc,kappa)
         call boundp(cbcpsi,n,bcpsi,nb,is_bound,dl,dzc,normx)
         call boundp(cbcpsi,n,bcpsi,nb,is_bound,dl,dzc,normy)
-        !call boundp(cbcpsi,n,bcpsi,nb,is_bound,dl,dzc,normz)
-        call boundp(cbcvel(:,:,2),n,bcvel(:,:,2),nb,is_bound,dl,dzc,normz)
+        call boundp(cbcpsi,n,bcpsi,nb,is_bound,dl,dzc,normz)
       end if
 #if defined(_SCALAR)
       call tm_scal(tm_coeff,n,dli,dzci,dzfi,dt,ssource,rho12,ka12,cp12,psi,u,v,w,s)

src/output.f90

@@ -637,6 +637,7 @@ subroutine cmpt_total_area(n,dli,dzci,dzfi,psi,area)
     real(rp), dimension(8) :: mx,my,mz
     real(rp) :: dpsidx,dpsidy,dpsidz
     integer  :: i,j,k
+    area = 0._rp
     do k=1,n(3)
       do j=1,n(2)
         do i=1,n(1)
@@ -793,7 +794,7 @@ end subroutine cmpt_total_mass
   !
   subroutine cmpt_total_energy(n,dl,dzf,rho12,psi,u,v,w,kin12)
     !
-    ! computes total mass of the system
+    ! computes total kinetic energy of the system
     !
     implicit none
     integer , intent(in ), dimension(3)        :: n
@@ -805,7 +806,7 @@ subroutine cmpt_total_energy(n,dl,dzf,rho12,psi,u,v,w,kin12)
     real(rp) :: vcell,ecell,kin1,kin2
     integer :: i,j,k
     !
-    ! calculate the total mass of the system
+    ! calculate the total kinetic of the system
     !
     kin1 = 0._rp
     kin2 = 0._rp