cal_d.f

Go to the documentation of this file.

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!SUBROUTINE CAL_D calculates shear velocity
!

      SUBROUTINE cal_d(V_sh)

      USE param
      USE param1
      USE parallel
      USE scales
      USE constant
      USE physprop
      USE fldvar
      USE visc_s
      USE rxns
      USE toleranc
      USE geometry
      USE indices
      USE is
      USE tau_s
      USE bc
      USE vshear
      USE compar
      USE fun_avg
      USE functions

      DOUBLE PRECISION V_sh,dis
!     DOUBLE PRECISION xdist(IMAX2,JMAX2)
!     DOUBLE PRECISION xdist3(IMAX2,JMAX2,KMAX2),cnter3(IMAX2,JMAX2,KMAX2)

!//SP Note the rational behind using global direction in I direction. This is to ensure
!     correctness in the way distances are calculated in the serial version and also
!     to give the capability to perform calculations over additional ghost layers in
!     to avoid if checks and communication
      DOUBLE PRECISION xdist(imin3:imax3,jstart3:jend3)
      DOUBLE PRECISION xdist3(imin3:imax3,jstart3:jend3,kstart3:kend3)

      INTEGER IJK,I1,J1,K1,I,J,K

      IF (no_k) THEN

!calculate distances
      DO  j1= jstart3, jend3
        xdist(1,j1)=1d0/(odx(1))
        if(imin3.ne.imin2) xdist(imin3,j1)=-1d0/(odx(imin3))
        DO  i1 = 2, imax3
        xdist(i1,j1)=1d0/(odx(i1))+xdist((i1-1),j1)
        END DO
      END DO

      DO  ijk= ijkstart3, ijkend3
          i = i_of(ijk)
          j = j_of(ijk)

        dis=xdist(i,j)

!shear velocity alligned u momentum cells

        vshe(ijk)=v_sh-2d0*(1d0-(dis-(1d0/odx(i))/xlength))*v_sh

!shear velocity alligned with scalar, v, w momentum cells

        vsh(ijk)=v_sh-2d0*(1d0-(dis-(1d0/odx(i))/xlength))*v_sh&
        -2d0*(1d0/(2d0*odx(i)*xlength))*v_sh

      END DO


      ELSE

!calculate distances
      DO k1=kstart3,kend3
        DO j1= jstart3, jend3
          IF (dead_cell_at(1,j1,k1)) cycle  ! skip dead cells
          xdist3(1,j1,k1)=1d0/(odx(1))
          if(imin3.ne.imin2) xdist(imin3,j1)=-1d0/(odx(imin3))
          DO  i1 = 2, imax3
           IF (dead_cell_at(i1,j1,k1)) cycle  ! skip dead cells
           xdist3(i1,j1,k1)=1d0/(odx(i1))+xdist3((i1-1),j1,k1)
          END DO
        END DO
      END DO


      DO  ijk= ijkstart3, ijkend3
          i = i_of(ijk)
          j = j_of(ijk)
          k = k_of(ijk)

          dis=xdist3(i,j,k)

!shear velocity alligned u momentum cells

        vshe(ijk)=v_sh-2d0*(1d0-(dis-(1d0/odx(i))/xlength))*v_sh

!shear velocity alligned with scalar, v, w momentum cells

        vsh(ijk)=v_sh-2d0*(1d0-(dis-(1d0/odx(i))/xlength))*v_sh&
        -2d0*(1d0/(2d0*odx(i)*xlength))*v_sh

      END DO

      END IF
      RETURN
      END SUBROUTINE cal_d

!// Comments on the modifications for DMP version implementation
!// 001 Include header file and common declarations for parallelization
!// 350 Changed do loop limits: 1,ijkmax2-> ijkstart3, ijkend3
!// 350 1206 change do loop limits: 1,kmax2->kstart3,kend3