Demo entry 6749382

sezioni

   

Submitted by anonymous on Jun 12, 2018 at 19:54
Language: Fortran. Code size: 2.8 kB.

subroutine sezionidurto(sez11,sez33,sez34,sez112,sez114A,sez114O,sez116,Y_H,Y_3HE,Y_4HE,Y_12C,Y_14N,Y_16O)
IMPLICIT NONE

REAL(KIND=8) :: T, T9, WLOT9, X_3HE, X_4HE, X_12C, X_14N,X_16O, T9A, WT9A,z, rho!CONDIZIONI INIZIALI
REAL(KIND=8) :: T912,T913, T923, T932, T943, T953, T9A13,T9A56, S5,S32,S7
REAL(KIND=8) :: sez11,sez33, sez34,sez112,sez114A,sez114O,sez116
REAL(KIND=8) :: Y_H,Y_3HE,A_3HE,Y_4HE, A_4HE,Y_12C,A_12C,Y_14N,A_14N,Y_16O,A_16O

T=1.7d7
T9=T*1.D-9
WLOT9=DLOG(T9)


!METALLICITA' INIZIALE
z=0.02d0

!densita iniziale
rho=1.6e2

!abbondanze iniziali in frazione di massa
X_3HE=1.D-6
X_4HE=0.28
X_12C=3.4714*1e-3
X_14N=1.0652*1e-3
X_16O=9.6702*1e-3

!numero di massa
A_3HE=3.d0
A_4HE=4.d0
A_12C=12.d0
A_14N=14.d0
A_16O=16.d0

!abbondanze iniziali
Y_3HE=X_3HE/A_3HE
Y_4HE=X_4HE/A_4HE
Y_12C=X_12C/A_12C
Y_14N=X_14N/A_14N
Y_16O=X_16O/A_16O
Y_H=1.d0-X_3HE-X_4HE-X_12C-X_14N-X_16O-z



!***************************************************
!*******calcolo tutte le sezioni d'urto già ********
!******* moltiplicate per la densità rho ***********
!***************************************************
T912=T9**(1.D0/2.D0)
T913=T9**(1.D0/3.D0)
T923=T9**(2.D0/3.D0)
T932=T9**(3.D0/2.D0)
T943=T9**(4.D0/3.D0)
T953=T9**(5.D0/3.D0)

!SEZIONE D'URTO REAZIONE P+P-->D+P
sez11=(3.82d-15/T923*DEXP(-3.38/T913)*(1.D0+0.123*T913+1.09*T923+0.938*T9))*rho

!SEZIONE D'URTO HE3 + HE3 --> HE4+2P
sez33=(5.96E10/T923*DEXP(-12.276/T913)*(1.+0.034*T913-0.199*T923-.047*T9+.032*T943+0.019*T953))*rho

!SEZIONE D'URTO HE3+HE4-->BE7+G
T9A=T9/(1.D0+4.95D-2*T9)
WT9A=DLOG(T9A)
T9A13=DEXP(1.D0/3.D0*WT9A)
T9A56=DEXP(5.D0/6.D0*WT9A)
sez34=(5.97E6*T9A56/T932*DEXP(-12.826/T9A13))*rho

!SEZIONE D'URTO C12+P-->C13+G
sez112=(2.04D7/T923*DEXP(-13.690/T913-(T9/1.500)**2)&
	*(1.D0+.030*T913+1.19*T923+.254*T9+2.06*T943+1.12*T953)&
	+1.08D5/T932*DEXP(-4.925/T9)+2.15D5/T932*DEXP(-18.179D0/T9))*rho

!SEZIONE D'URTO N14+P-->O15 +G
S5=(5.08D7/T923*DEXP(-15.228/T913-(T9/3.090)**2)*(1.+.027*T913-.778*T923-.149*T9+.261*T943 &
      +.127*T953)+2.28D3/T932*DEXP(-3.011/T9)+1.65D4*T913*DEXP(-12.007/T9))
S32=(9.78D8/T923*DEXP(-15.251/T913-(T9/0.450)**2)*(1.+.027*T913+.219*T923+0.042*T9+6.83*T943+3.32*T953) &
      +1.11D4/T932*DEXP(-3.328/T9)+1.49D4/T932*DEXP(-4.665/T9) &
      +3.80D6/T932*DEXP(-11.048/T9))
S7=(1.08D12/T923*DEXP(-15.251/T913-(T9/0.522)**2) &
      *(1.+.027*T913+2.62*T923+0.501*T9+5.36*T943+2.60*T953) &
      +1.19D8/T932*DEXP(-3.676/T9)+5.41D8/T912*DEXP(-8.926/T9)+0.1d0* &
      4.72D8/T932*DEXP(-7.721/T9)+2.20D9/T932*DEXP(-11.418/T9))

!SEZIONE D'URTO N15+P-->O16+G

sez114A=S5*(1.D0-S7/(S7+S32))*rho

!SEZIONE D'URTO N15+P-->A+C12
sez114O=S5*S7/(S7+S32)*rho

!SEZIONE D'URTO O16+P-->F17+G
sez116=(1.50D8/(T923*(1.+2.13*(1.-DEXP(-.728*T923))))*DEXP(-16.692/T913))*rho


return
end subroutine sezionidurto

This snippet took 0.01 seconds to highlight.

Back to the Entry List or Home.

Delete this entry (admin only).