get_whitelight_Thomson Subroutine

public subroutine get_whitelight_Thomson(Rl, Rin, Ne, A, B, C, D, fluxTB, fluxPB)

Arguments

Type IntentOptional Attributes Name
double precision, intent(in) :: Rl
double precision, intent(in) :: Rin
double precision, intent(in) :: Ne
double precision, intent(in) :: A
double precision, intent(in) :: B
double precision, intent(in) :: C
double precision, intent(in) :: D
double precision, intent(inout) :: fluxTB
double precision, intent(inout) :: fluxPB

Called by

proc~~get_whitelight_thomson~~CalledByGraph proc~get_whitelight_thomson get_whitelight_Thomson proc~integrate_whitelight_spherical integrate_whitelight_spherical proc~integrate_whitelight_spherical->proc~get_whitelight_thomson proc~get_image get_image proc~get_image->proc~integrate_whitelight_spherical proc~get_euv_image get_EUV_image proc~get_euv_image->proc~get_image proc~get_sxr_image get_SXR_image proc~get_sxr_image->proc~get_image proc~get_whitelight_image get_whitelight_image proc~get_whitelight_image->proc~get_image

Source Code

    subroutine get_whitelight_Thomson(Rl,Rin,Ne,A,B,C,D,fluxTB,fluxPB)
      ! use the method in SSW/eltheory
      double precision, intent(in) :: Rl,Rin,Ne,A,B,C,D
      double precision, intent(inout) :: fluxTB,fluxPB

      double precision :: const,u,Bt,Br,PB,TB,sinchi2

      u=0.63d0
      const=1.24878d-25/(1.d0-u/3.d0)
      sinchi2=(Rin/Rl)**2
      Bt=const*(C+u*(D-C))
      PB=const*sinchi2*((A+u*(B-A)))
      Br=Bt-pB
      TB=Bt+Br
      fluxTB=TB*Ne
      fluxPB=PB*Ne

    end subroutine get_whitelight_Thomson