Initial commit

1 files changed, 45 insertions, 0 deletions

diff --git a/OpacHelium.f b/OpacHelium.f
new file mode 100755
index 0000000..b02fc91
--- /dev/null
+++ b/OpacHelium.f
@@ -0,0 +1,45 @@
+
+c******************************************************************************
+c  The subroutines needed to calculate the He- b-f and f-f opacities are in 
+c  this file.  These are from ATLAS9.
+c******************************************************************************
+
+
+
+
+
+      subroutine opacHeminus
+c******************************************************************************
+c  This routine computes the He- bound-free and free-free opacities.
+c******************************************************************************
+
+      implicit real*8 (a-h,o-z)
+      include 'Atmos.com'
+      include 'Kappa.com'
+      save
+      data freq1 /0./
+         
+      if (freq .ne. freq1) then
+         freq1 = freq
+         a1 =  3.397d-46 + (-5.216d-31+7.039d-15/freq)/freq
+         b1 = -4.116d-42 + ( 1.067d-26+8.135d-11/freq)/freq
+         cc =  5.081d-37 + (-8.724d-23-5.659d-08/freq)/freq
+      endif
+
+      do i=1,ntau
+         aHeminus(i) = (a1*t(i) + b1 +c1/t(i))*ne(i)*
+     .                 numdens(2,1,i)/u(2,1,i)
+      enddo
+
+      return
+      end
+
+
+
+
+
+
+
+
+
+