9SI20= Cdtau ) THEN$K8 PK!& + TF = ( EXP( -ABS( tau1 - tau3 )*3/2 ) &XK"1 1 - EXP( -ABS( tau1 - tau2 )*3/2 ) )/dtau*(2./3.),K #   ELSE PK\$& + TF = ( EXP( -ABS( tau1 - tau2 )*3/2 ) &LK%' ' + EXP( -ABS( tau1 - tau3 )*3/2 ) )/2.,K&   END IF$K'4K,( END FUNCTION TFhK)A!----------------------------------------------------------------hK*`#ASUBROUTINE get_urf( tau, temp, Fu) ! ؤξ«̩٤η׻$K +XKx,^% 2 REAL,DIMENSION(0:),INTENT(in) :: tau ! ŪTK-^& - REAL,DIMENSION(0:),INTENT(in) :: temp ! \K(.^$ 7 REAL,DIMENSION(:),INTENT(out) :: Fu ! «̩\K/^$ 5 REAL,DIMENSION(:),ALLOCATABLE :: bf ! «̩8K0^  REAL :: dbdtauHK1^ ! REAL :: dtau ! Ū<K@2Q  INTEGER :: i,j,nnl$Kd38K4  nnl = SIZE(tau)-2 $K5<K6  ALLOCATE( bf(0:nnl+1) )@K<7  CALL get_bbf( temp, bf )$K`80K9 Fu = bf(1:)$K:4K;  DO i = 2,nnl4K<   DO j = 2,iDK`=  dbdtau = bf(j-1) - bf(j)HK0 > ! dtau = tau(j-1) - tau(j)hK0t?= D Fu(i) = Fu(i) + dbdtau * TF( tau(i), tau(j-1), tau(j), dtau )0K0@ END DO,K0A  END DO$K0B4K1(C  DEALLOCATE(bf)$K1LD<K1E END SUBROUTINE get_urfhK1FA!----------------------------------------------------------------dK2TG` >SUBROUTINE get_drf(tau,temp,Fd) ! ؤβ«̩٤η׻$K2xHXK2I^% 2 REAL,DIMENSION(0:),INTENT(in) :: tau ! ŪTK3$J^& - REAL,DIMENSION(0:),INTENT(in) :: temp ! \K3K^$ 7 REAL,DIMENSION(:),INTENT(out) :: Fd ! «̩$K3L\K4M^$ 5 REAL,DIMENSION(:),ALLOCATABLE :: bf ! «̩8K48N^  REAL :: dbdtauHK4O^ ! REAL :: dtau ! Ū<K4PQ  INTEGER :: i,j,nnl$K4Q8K5R  nnl = SIZE(tau) - 2<K5TS  ALLOCATE( bf(0:nnl+1) )@K5T  CALL get_bbf( temp, bf )$K5U4K5V  DO i = 1,nnlhK6TW 8C Fd(i) = bf(i)- bf(nnl)*EXP( -( tau(i) - tau(nnl) )*3/2 )8K6X   DO j = i+1, nnlDK6Y  dbdtau = bf(j-1) - bf(j)HK7Z ! dtau = tau(j-1) - tau(j)hK7[= D Fd(i) = Fd(i) - dbdtau * TF( tau(i), tau(j), tau(j-1), dtau )0K7\ END DO,K7]  END DO$K8^4K84_  DEALLOCATE(bf)$K8X`<K8a END SUBROUTINE get_drfhK8bA!----------------------------------------------------------------dK9`c`)?SUBROUTINE get_cvf( Fu, Fd, Fsol, Fcnv ) ! ήǮή®̩٤η׻$K9dXK9e^# 4 REAL,DIMENSION(:),INTENT(in) :: Fu ! «̩XK:4f^# 4 REAL,DIMENSION(:),INTENT(in) :: Fd ! «̩`K:g^( 9 REAL,DIMENSION(:),INTENT(in) :: Fsol ! «̩\K:h^& 7 REAL,DIMENSION(:),INTENT(out) :: Fcnv ! ήǮή«̩$K;i<K;Pj  Fcnv = Fd + Fsol - Fu$K;tk<K;l END SUBROUTINE get_cvfhK<mA!----------------------------------------------------------------`K REAL,DIMENSION(0:),INTENT(out) :: Fnet ! ̣«̩$K=s@K>t  Fnet(1:) = Fu(1:) - Fd(1:)8K>Lu  Fnet(0) = Fnet(1)LK>v' ( Fnet(SIZE(Fnet)-1) = Fnet(SIZE(Fnet)-2)$K>w<K>x END SUBROUTINE get_fnethK?`yA!----------------------------------------------------------------lK?z`5GSUBROUTINE rdeqiv( TP, tau, bf, temp, Ftop, Fu, Fd ) ! ʿղη׻$K@ { (TK@`|Q / INTEGER,INTENT(in) :: TP ! طüֹXK@}^% 2 REAL,DIMENSION(0:),INTENT(in) :: tau ! Ū\KA~^& 7 REAL,DIMENSION(0:),INTENT(inout):: bf ! «̩TKAh^( / REAL,DIMENSION(0:),INTENT(inout):: temp ! \KA^ 5 REAL,INTENT(out) :: Ftop ! 絤ü«̩\KB ^& 7 REAL,DIMENSION(0:),INTENT(inout):: Fu ! «̩\KB|^& 7 REAL,DIMENSION(0:),INTENT(inout):: Fd ! «̩$KBDKB !-絤ü«̩٤η׻$KCTKC\. / Ftop = 2.*sig*temp(TP)**4/( 1.5*tau(TP) + 1. )$KC@KC!-طβ; ʿղ$KCXKD<. 1 bf(TP+1:) = 0.5*Ftop*( 1.5*tau(TP+1:) + 1. )PKD( ) temp(TP+1:) = ( bf(TP+1:)/sig )**(0.25) TKD. / Fu(TP+1:) = 0.5*Ftop*( 1.5*tau(TP+1:) + 2. )LKE,& ' Fd(TP+1:) = 0.5*Ftop*1.5*tau(TP+1:1)$KEP<KE END SUBROUTINE rdeqivhKEA!----------------------------------------------------------------\KFP`"8SUBROUTINE fnetmaxl( Fnet, LTRP ) ! ̣«̩٤η׻$KFtdKF^& = REAL,DIMENSION(0:),INTENT(in) :: Fnet ! ̣«̩TKG,Q 0 INTEGER,INTENT(inout) :: LTRP ! ֹ̹<KGhQ  INTEGER :: nnl ! 8KGQ  INTEGER :: i ! $KG8KG  nnl = SIZE(Fnet)-1$KH 4KHT   DO i = nnl,1,-1LKH! % IF ( Fnet(i) >= Fnet(LTRP) ) THEN4KH  LTRP = i0KI  END IF,KI0  END DO$KIT<KI END SUBROUTINE fnetmaxlhKIA!----------------------------------------------------------------<KP  END MODULE rad_module