floortemperature Subroutine

public subroutine floortemperature(qdt, wCT, wnew, x, fl)

Arguments

Type IntentOptional Attributes Name
double precision, intent(in) :: qdt
double precision, intent(in) :: wCT(nw_phys)
double precision, intent(inout) :: wnew(nw_phys)
double precision, intent(in) :: x(1:ndim)
type(rc_fluid), intent(in) :: fl

Calls

proc~~floortemperature~~CallsGraph proc~floortemperature floortemperature proc~get_pthermal get_pthermal proc~floortemperature->proc~get_pthermal proc~get_rfactor get_Rfactor proc~floortemperature->proc~get_rfactor proc~get_rho get_rho proc~floortemperature->proc~get_rho

Called by

proc~~floortemperature~~CalledByGraph proc~floortemperature floortemperature proc~radiative_cooling_add_source radiative_cooling_add_source proc~radiative_cooling_add_source->proc~floortemperature

Source Code

    subroutine floortemperature(qdt,wCT,wnew,x,fl)
      
    !  Force minimum temperature to a fixed temperature
      double precision, intent(in)    :: qdt, wCT(nw_phys)
      double precision, intent(inout) :: wnew(nw_phys)
      double precision, intent(in)    :: x(1:ndim)
      type(rc_fluid), intent(in) :: fl
      double precision :: etherm, rho, emin, Rfactor

      etherm = get_pthermal(wCT,x)  
      rho = get_rho(wCT,x)  
      Rfactor = get_Rfactor()
      emin = rho*fl%tlow*Rfactor
      if(etherm < emin) then
        wnew(iw_e)=wnew(iw_e)+(emin-etherm)*invgam
      end if
    end subroutine floortemperature