radiative_cooling_add_source Subroutine

public subroutine radiative_cooling_add_source(qdt, wCT, wCTprim, wnew, x)

Arguments

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

Calls

proc~~radiative_cooling_add_source~~CallsGraph proc~radiative_cooling_add_source radiative_cooling_add_source proc~cool_exact cool_exact proc~radiative_cooling_add_source->proc~cool_exact proc~floortemperature floortemperature proc~radiative_cooling_add_source->proc~floortemperature proc~calc_l_extended calc_l_extended proc~cool_exact->proc~calc_l_extended proc~findt findT proc~cool_exact->proc~findt proc~findy findY proc~cool_exact->proc~findy proc~get_pthermal get_pthermal proc~cool_exact->proc~get_pthermal proc~get_rfactor get_Rfactor proc~cool_exact->proc~get_rfactor proc~get_rho get_rho proc~cool_exact->proc~get_rho proc~floortemperature->proc~get_pthermal proc~floortemperature->proc~get_rfactor proc~floortemperature->proc~get_rho

Source Code

    subroutine radiative_cooling_add_source(qdt,wCT,wCTprim,wnew,x)
      

    ! w[iw]=w[iw]+qdt*S[wCT,x] where S is the source based on wCT within ixO
      double precision, intent(in) :: qdt, wCT(nw_phys), wCTprim(nw_phys)
      double precision, intent(inout) :: wnew(nw_phys)
      double precision, intent(in) :: x(1:ndim)

      call cool_exact(qdt,wCT,wCTprim,wnew,x,rc_fl)
      if( rc_fl%Tfix ) call floortemperature(qdt,wCT,wnew,x,rc_fl)

    end subroutine radiative_cooling_add_source