39 lines
928 B
Fortran
39 lines
928 B
Fortran
SUBROUTINE TDPINI
|
|
C =================
|
|
C
|
|
C initialization of only temperature dependent quantities
|
|
C
|
|
INCLUDE 'IMPLIC.FOR'
|
|
INCLUDE 'BASICS.FOR'
|
|
INCLUDE 'ATOMIC.FOR'
|
|
INCLUDE 'MODELQ.FOR'
|
|
INCLUDE 'ODFPAR.FOR'
|
|
INCLUDE 'ALIPAR.FOR'
|
|
PARAMETER (CFF1=1.3727D-25,CFF2=4.3748D-10,CFF3=2.5993D-7)
|
|
PARAMETER (SGFF0 = 3.694D8)
|
|
C
|
|
C temperature-dependent quantities
|
|
C
|
|
DO ID=1,ND
|
|
T=TEMP(ID)
|
|
T1=UN/T
|
|
HKT1(ID)=HK*T1
|
|
HKT21(ID)=HKT1(ID)*T1
|
|
TK1(ID)=HKT1(ID)/H
|
|
SQT1(ID)=SQRT(T)
|
|
TEMP1(ID)=T1
|
|
CALL GFREE0(ID)
|
|
EMEL1(ID)=UN
|
|
END DO
|
|
C
|
|
C delta m (for evaluation of optical depths)
|
|
C
|
|
DO ID=1,ND-1
|
|
DELDM(ID)=HALF*(DM(ID+1)-DM(ID))
|
|
deldmz(id)=deldm(id)
|
|
if(izscal.eq.1) deldmz(id)=half*(zd(id)-zd(id+1))
|
|
END DO
|
|
DEDM1=DM(1)/DENS(1)
|
|
RETURN
|
|
END
|