unremembered modi

Author: fangjian19

src/bc.F90

@@ -692,10 +692,10 @@ subroutine immbody(timerept)
         vel(i,j,k,1)= -1.d0*var_u(1)*pb%dis2ghost/pb%dis2image
         vel(i,j,k,2)= -1.d0*var_u(2)*pb%dis2ghost/pb%dis2image
         vel(i,j,k,3)= -1.d0*var_u(3)*pb%dis2ghost/pb%dis2image
-          ! tmp(i,j,k)=tinf-(var_t-tinf)*pb%dis2ghost/pb%dis2image
+        ! isothernal
+          tmp(i,j,k)=tinf-(var_t-tinf)*pb%dis2ghost/pb%dis2image
         ! adiabatic
-        tmp(i,j,k)=var_t
-        ! tmp(i,j,k)=twall(3)
+        ! tmp(i,j,k)=var_t
         prs(i,j,k)=var_p
         if(nondimen) then 
           rho(i,j,k)=thermal(pressure=prs(i,j,k),temperature=tmp(i,j,k))
@@ -714,8 +714,9 @@ subroutine immbody(timerept)
         vel(i,j,k,2)=0.d0
         vel(i,j,k,3)=0.d0
         ! aidabatic
-        tmp(i,j,k)=var_t
-        ! tmp(i,j,k)=twall(3)
+        ! tmp(i,j,k)=var_t
+        ! isothermal
+        tmp(i,j,k)=twall(3)
         prs(i,j,k)=var_p
         !
         if(nondimen) then 
@@ -767,8 +768,8 @@ subroutine immbody(timerept)
         vel(i,j,k,1)=0.d0
         vel(i,j,k,2)=0.d0
         vel(i,j,k,3)=0.d0
-        ! tmp(i,j,k)  =twall(3)
-        tmp(i,j,k)  =tinf
+        tmp(i,j,k)  =twall(3)
+        ! tmp(i,j,k)  =tinf
         prs(i,j,k)  =pinf
         if(nondimen) then 
           rho(i,j,k)  =thermal(pressure=prs(i,j,k),temperature=tmp(i,j,k))

src/geom.F90

@@ -4674,7 +4674,7 @@ subroutine solid_udf_cavity(xp,inside,bnode)
     step_right_x=xcav_right
     step_upper_y=ycav_upper
     !
-    if(xp(2)<step_upper_y .and. xp(1)<=step_left_x) then
+    if(xp(2)<=step_upper_y .and. xp(1)<=step_left_x) then
       !
       if(present(inside)) inside=.true.
       !
@@ -4703,7 +4703,7 @@ subroutine solid_udf_cavity(xp,inside,bnode)
         !
       endif
       !
-    elseif(xp(2)<step_upper_y .and. xp(1)>=step_right_x) then
+    elseif(xp(2)<=step_upper_y .and. xp(1)>=step_right_x) then
       !
       if(present(inside)) inside=.true.
       !

src/parallel.F90

@@ -1747,6 +1747,8 @@ subroutine updatable_rel2d_a2a_v(datasend,datarecv,                 &
     !
     deallocate(senddispls,recvdispls)
     !
+    call mpi_type_free(newtype,ierr)
+    !
     if(present(timerept) .and. timerept) then
       !
       subtime=subtime+ptime()-time_beg