You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
On the one hand, in the subroutine of "Grid_Mod_Load_Geo(grid, this_proc, domain)", there are:
read(fu) (grid % xc(c), c = -grid % n_bnd_cells, grid % n_cells)
read(fu) (grid % yc(c), c = -grid % n_bnd_cells, grid % n_cells)
read(fu) (grid % zc(c), c = -grid % n_bnd_cells, grid % n_cells)
read(fu) (grid % dx(s), s = 1, grid % n_faces + grid % n_shadows)
read(fu) (grid % dy(s), s = 1, grid % n_faces + grid % n_shadows)
read(fu) (grid % dz(s), s = 1, grid % n_faces + grid % n_shadows)
On the other hand, in the subroutine of " Grid_Mod_Calculate_Face_Geometry(grid)", there are:
!-------------------------------------------------------!
! Calculate the distance between neighbouring cells !
! (For normal faces, including the periodic ones, !
! dx, dy and dz are distanes between cell centers) !
!-------------------------------------------------------!
do s = 1, grid % n_faces
c1 = grid % faces_c(1,s)
c2 = grid % faces_c(2,s)
On the one hand, in the subroutine of "Grid_Mod_Load_Geo(grid, this_proc, domain)", there are:
read(fu) (grid % xc(c), c = -grid % n_bnd_cells, grid % n_cells)
read(fu) (grid % yc(c), c = -grid % n_bnd_cells, grid % n_cells)
read(fu) (grid % zc(c), c = -grid % n_bnd_cells, grid % n_cells)
read(fu) (grid % dx(s), s = 1, grid % n_faces + grid % n_shadows)
read(fu) (grid % dy(s), s = 1, grid % n_faces + grid % n_shadows)
read(fu) (grid % dz(s), s = 1, grid % n_faces + grid % n_shadows)
On the other hand, in the subroutine of " Grid_Mod_Calculate_Face_Geometry(grid)", there are:
!-------------------------------------------------------!
! Calculate the distance between neighbouring cells !
! (For normal faces, including the periodic ones, !
! dx, dy and dz are distanes between cell centers) !
!-------------------------------------------------------!
do s = 1, grid % n_faces
c1 = grid % faces_c(1,s)
c2 = grid % faces_c(2,s)
end do ! faces
Thus,I think these lines mean "grid % dx(s) = 2 * grid % dx(s) ", "grid % dy(s) = 2 * grid % dy(s) " and "grid % dz(s) = 2 * grid % dz(s) " . Is it correct ?
The text was updated successfully, but these errors were encountered: