Initial commit

1 files changed, 49 insertions, 0 deletions

diff --git a/Calmod.f b/Calmod.f
new file mode 100755
index 0000000..3938bef
--- /dev/null
+++ b/Calmod.f
@@ -0,0 +1,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