destroy_viscosity Subroutine

public subroutine destroy_viscosity()

Deallocate and nullify viscosity (turbulent/molecular)

Arguments

None

Calls

proc~~destroy_viscosity~~CallsGraph proc~destroy_viscosity destroy_viscosity interface~dealloc dealloc proc~destroy_viscosity->interface~dealloc proc~dealloc_rank2_real dealloc_rank2_real interface~dealloc->proc~dealloc_rank2_real proc~dealloc_rank4_real dealloc_rank4_real interface~dealloc->proc~dealloc_rank4_real proc~dealloc_rank1_integer dealloc_rank1_integer interface~dealloc->proc~dealloc_rank1_integer proc~dealloc_rank2_integer dealloc_rank2_integer interface~dealloc->proc~dealloc_rank2_integer proc~dealloc_rank1_real dealloc_rank1_real interface~dealloc->proc~dealloc_rank1_real proc~dealloc_rank3_integer dealloc_rank3_integer interface~dealloc->proc~dealloc_rank3_integer proc~dealloc_rank3_real dealloc_rank3_real interface~dealloc->proc~dealloc_rank3_real proc~dealloc_rank5_real dealloc_rank5_real interface~dealloc->proc~dealloc_rank5_real proc~dealloc_rank6_real dealloc_rank6_real interface~dealloc->proc~dealloc_rank6_real

Called by

proc~~destroy_viscosity~~CalledByGraph proc~destroy_viscosity destroy_viscosity proc~destroy_solver destroy_solver proc~destroy_solver->proc~destroy_viscosity proc~finish_run finish_run proc~finish_run->proc~destroy_solver proc~abort_run abort_run proc~abort_run->proc~destroy_solver program~main main program~main->proc~finish_run

Contents

Source Code


Source Code

    subroutine destroy_viscosity()
      !< Deallocate and nullify viscosity (turbulent/molecular)

      ! destroy_molecular_viscosity ---!
      if (mu_ref/=0.) then
        call dealloc(mu)
      end if

      !--- destroy_turbulent_viscosity ---!
      if (turbulence/='none') then
        select case (trim(turbulence))

          case ('none')
            !do nothing
            continue

          case('sa', 'saBC')
            nullify(sa_mu)

          case ('sst', 'sst2003')
            nullify(sst_mu)
            !-- blending funciton F1 --!
            call dealloc(sst_F1)
            !--- sst blending funciton destoryed--!

          case ('kkl')
            nullify(kkl_mu)

          case DEFAULT 
            !call turbulence_read_error()
            Fatal_error

        end select
        call dealloc(mu_t)
      end if
      !--- end of turublent viscosity destruction ---!

    end subroutine destroy_viscosity