1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
|
subroutine calmod
c******************************************************************************
c this program re-calculates BEGN models on a tau5000 scale.
c******************************************************************************
implicit real*8 (a-h,o-z)
include 'Atmos.com'
include 'Pstuff.com'
c*****open the files for standard output and the re-calculated model
nf1out = 20
lscreen = 4
array = 'STANDARD OUTPUT'
nchars = 15
call infile ('output ',nf1out,'formatted ',0,nchars,
. f1out,lscreen)
nf2out = 21
lscreen = lscreen + 2
array = 'RE-CALCULATED MODEL OUTPUT'
nchars = 26
call infile ('output ',nf2out,'formatted ',0,nchars,
. f2out,lscreen)
c*****open and read the model atmosphere
nfmodel = 30
lscreen = lscreen + 2
array = 'THE MODEL ATMOSPHERE'
nchars = 20
call infile ('input ',nfmodel,'formatted ',0,nchars,
. fmodel,lscreen)
call inmodel
c*****do the tau scale conversion
call opacit(2,wavref)
write (nf2out,1001) wavref,(taulam(i),t(i),pgas(i),
. ne(i),kaplam(i),i=1,ntau)
return
c*****format statements
1001 format ('output model for wavelength = ',f5.0/
. (1pe11.4,0pf7.0,1p3e11.4))
end
|
