32 lines
795 B
Fortran
32 lines
795 B
Fortran
SUBROUTINE INIBLM
|
|
C =================
|
|
C
|
|
C driving procedure for treating a partial molecular line list for the
|
|
C current wavelength region
|
|
C
|
|
INCLUDE 'PARAMS.FOR'
|
|
INCLUDE 'MODELP.FOR'
|
|
INCLUDE 'SYNTHP.FOR'
|
|
INCLUDE 'LINDAT.FOR'
|
|
C
|
|
PARAMETER (DP0=3.33564E-11, DP1=1.651E8, UN=1.)
|
|
C
|
|
XX=FREQ(1)
|
|
IF(NFREQ.GE.2) XX=0.5*(FREQ(1)+FREQ(2))
|
|
BNU=BN*(XX*1.E-15)**3
|
|
HKF=HK*XX
|
|
DO ID=1,ND
|
|
T=TEMP(ID)
|
|
EXH=EXP(HKF/T)
|
|
EXHK(ID)=UN/EXH
|
|
PLAN(ID)=BNU/(EXH-UN)
|
|
STIM(ID)=UN-EXHK(ID)
|
|
DO IMOL=1,NMOLEC
|
|
IF(AMMOL(IMOL).GT.0.)
|
|
* DOPMOL(IMOL,ID)=UN/(XX*DP0*SQRT(DP1*T/AMMOL(IMOL)+
|
|
* VTURB(ID)))
|
|
END DO
|
|
END DO
|
|
RETURN
|
|
END
|